BCFtools-Teaching: Difference between revisions
(Created page with "Category:TeachingCategory:SoftwareCategory:Bioinformatics === Category === Bioinformatics === Program On === Teaching === Version === 1.3.1 === A...") |
No edit summary |
||
Line 9: | Line 9: | ||
=== Version === | === Version === | ||
1. | 1.8 | ||
=== Author / Distributor === | === Author / Distributor === | ||
Line 16: | Line 16: | ||
=== Description === | === Description === | ||
" | "Samtools is a suite of programs for interacting with high-throughput sequencing data. BCFtools - Reading/writing BCF2/VCF/gVCF files and calling/filtering/summarising SNP and short indel sequence variants" | ||
More details are at [http://www.htslib.org/ BCFtools] | More details are at [http://www.htslib.org/ BCFtools] | ||
=== Running Program === | === Running Program === | ||
The last version of this application is at /usr/local/apps/eb/BCFtools/1. | The last version of this application is at /usr/local/apps/eb/BCFtools/1.8-foss-2018a | ||
To use this version, please loads the module with | To use this version, please loads the module with | ||
<pre class="gscript"> | <pre class="gscript"> | ||
ml BCFtools/1. | ml BCFtools/1.8-foss-2018a | ||
</pre> | </pre> | ||
Line 44: | Line 42: | ||
cd $SLURM_SUBMIT_DIR<br> | cd $SLURM_SUBMIT_DIR<br> | ||
ml BCFtools/1.8-foss-2018a<br> | |||
ml | bcftools <u>[options]</u><br> | ||
</div> | </div> | ||
In the real submission script, at least all the above underlined values need to be reviewed or to be replaced by the proper values | In the real submission script, at least all the above underlined values need to be reviewed or to be replaced by the proper values. | ||
Please refer to [[Running_Jobs_on_the_teaching_cluster]], [[Running_Jobs_on_the_teaching_cluster#Running_an_X-windows_application | Run X window Jobs]] and [[Running_Jobs_on_the_teaching_cluster#How_to_open_an_interactive_session | Run interactive Jobs]] for more details of running jobs at Teaching cluster. | Please refer to [[Running_Jobs_on_the_teaching_cluster]], [[Running_Jobs_on_the_teaching_cluster#Running_an_X-windows_application | Run X window Jobs]] and [[Running_Jobs_on_the_teaching_cluster#How_to_open_an_interactive_session | Run interactive Jobs]] for more details of running jobs at Teaching cluster. | ||
Line 60: | Line 58: | ||
<pre class="gcommand"> | <pre class="gcommand"> | ||
BCFtools/1. | ml BCFtools/1.8-foss-2018a | ||
bcftools --help | bcftools bcftools --help | ||
Program: bcftools (Tools for variant calling and manipulating VCFs and BCFs) | Program: bcftools (Tools for variant calling and manipulating VCFs and BCFs) | ||
License: GNU GPLv3+, due to use of the GNU Scientific Library | License: GNU GPLv3+, due to use of the GNU Scientific Library | ||
Version: 1. | Version: 1.8 (using htslib 1.8) | ||
Usage: bcftools [--version|--version-only] [--help] <command> <argument> | Usage: bcftools [--version|--version-only] [--help] <command> <argument> | ||
Line 84: | Line 82: | ||
query transform VCF/BCF into user-defined formats | query transform VCF/BCF into user-defined formats | ||
reheader modify VCF/BCF header, change sample names | reheader modify VCF/BCF header, change sample names | ||
sort sort VCF/BCF file | |||
view VCF/BCF conversion, view, subset and filter VCF/BCF files | view VCF/BCF conversion, view, subset and filter VCF/BCF files | ||
Line 90: | Line 89: | ||
consensus create consensus sequence by applying VCF variants | consensus create consensus sequence by applying VCF variants | ||
cnv HMM CNV calling | cnv HMM CNV calling | ||
csq call variation consequences | |||
filter filter VCF/BCF files using fixed thresholds | filter filter VCF/BCF files using fixed thresholds | ||
gtcheck check sample concordance, detect sample swaps and contamination | gtcheck check sample concordance, detect sample swaps and contamination | ||
mpileup multi-way pileup producing genotype likelihoods | |||
polysomy detect number of chromosomal copies | polysomy detect number of chromosomal copies | ||
roh identify runs of autozygosity (HMM) | roh identify runs of autozygosity (HMM) |
Revision as of 11:36, 10 August 2018
Category
Bioinformatics
Program On
Teaching
Version
1.8
Author / Distributor
Description
"Samtools is a suite of programs for interacting with high-throughput sequencing data. BCFtools - Reading/writing BCF2/VCF/gVCF files and calling/filtering/summarising SNP and short indel sequence variants" More details are at BCFtools
Running Program
The last version of this application is at /usr/local/apps/eb/BCFtools/1.8-foss-2018a
To use this version, please loads the module with
ml BCFtools/1.8-foss-2018a
Here is an example of a shell script, sub.sh, to run on at the batch queue:
#!/bin/bash
#SBATCH --job-name=j_BCFtools
#SBATCH --partition=batch
#SBATCH --mail-type=ALL
#SBATCH --mail-user=username@uga.edu
#SBATCH --ntasks=1
#SBATCH --mem=10gb
#SBATCH --time=08:00:00
#SBATCH --output=BCFtools.%j.out
cd $SLURM_SUBMIT_DIR
ml BCFtools/1.8-foss-2018a
bcftools [options]
In the real submission script, at least all the above underlined values need to be reviewed or to be replaced by the proper values.
Please refer to Running_Jobs_on_the_teaching_cluster, Run X window Jobs and Run interactive Jobs for more details of running jobs at Teaching cluster.
Here is an example of job submission command:
sbatch ./sub.sh
Documentation
ml BCFtools/1.8-foss-2018a bcftools bcftools --help Program: bcftools (Tools for variant calling and manipulating VCFs and BCFs) License: GNU GPLv3+, due to use of the GNU Scientific Library Version: 1.8 (using htslib 1.8) Usage: bcftools [--version|--version-only] [--help] <command> <argument> Commands: -- Indexing index index VCF/BCF files -- VCF/BCF manipulation annotate annotate and edit VCF/BCF files concat concatenate VCF/BCF files from the same set of samples convert convert VCF/BCF files to different formats and back isec intersections of VCF/BCF files merge merge VCF/BCF files files from non-overlapping sample sets norm left-align and normalize indels plugin user-defined plugins query transform VCF/BCF into user-defined formats reheader modify VCF/BCF header, change sample names sort sort VCF/BCF file view VCF/BCF conversion, view, subset and filter VCF/BCF files -- VCF/BCF analysis call SNP/indel calling consensus create consensus sequence by applying VCF variants cnv HMM CNV calling csq call variation consequences filter filter VCF/BCF files using fixed thresholds gtcheck check sample concordance, detect sample swaps and contamination mpileup multi-way pileup producing genotype likelihoods polysomy detect number of chromosomal copies roh identify runs of autozygosity (HMM) stats produce VCF/BCF stats Most commands accept VCF, bgzipped VCF, and BCF with the file type detected automatically even when streaming from a pipe. Indexed VCF and BCF will work in all situations. Un-indexed VCF and BCF and streams will work in most but not all situations.
Installation
Source code is obtained from BCFtools
System
64-bit Linux