Freebayes
FreeBayes is a genetic variant detector designed to find small polymorphisms (SNPs, indels, MNPs and complex events).
FreeBayes is haplotype-based, in the sense that it calls variants based on the literal sequences of reads aligned to a particular target, not their precise alignment. This model is a straightforward generalization of previous ones (e.g. PolyBayes, samtools, GATK) which detect or report variants based on alignments. This method avoids one of the core problems with alignment-based variant detection, that identical sequences may have multiple possible alignments.
FreeBayes uses short-read alignments (BAM files) for any number of individuals from a population and a reference genome to determine the most-likely combination of genotypes for the population at each position in the reference. It reports positions which it finds putatively polymorphic in variant call file (VCF) format. It can also use an input set of variants (VCF) as a source of prior information, and a copy number variant map (BED) to define non-uniform ploidy variation across the samples under analysis.
License
Free to use and open source under MIT License.
Available
- Puhti: 1.3.6, 1.3.7
Usage
First load the FreeBayes module.
module load freebayes
After this you can launch Freebayes. For example:
freebayes -f reference.fa input.bam > results.vcf
Note that FreeBayes requires a BAM file that is indexed. A BAM file can be indexed with command:
samtools index input.bam
FreeBayes analysis jobs can be computationally heavy and should be run as batch jobs on Puhti.
On Puhti, you can use freebayes-puhti to automatically submit a Freebayes job to the batch job system.
This tool also speeds up the analysis by running the analysis as several simultaneous tasks in parallel.
To be able to use freebayes-puhti, you first need to define a regions file for your reference fasta file.
This can be done with the command:
fasta_generate_regions.py reference.fa.fai 100000 > regions.txt
For small datasets you may decrease the region size in the command above so that you will get more than 100 regions in the regions file.
Once you have the regions file created, you can launch your analysis task with the command:
freebayes-puhti -regions regions.txt -f reference.fa input.bam -out results.vcf
freebayes-puhti will execute your FreeBayes analysis as an automatically generated array batch job. The results will also be automatically merged and sorted once the batch jobs have finished. By default, freebayes-puhti allows each sub-job to use 16 GB of memory and to run for 24 hours. For massive FreeBayes jobs, this may not be sufficient. In that case, you can try to use options -mem and -time to extend the limits. -mem option
defines the memory reservation in gigabytes while the -time option defines the time reservation in hours. For example, extending the task to 64 GB of memory and 48 hours of running time could be done with the command:
freebayes-puhti -mem 64 -time 48 -regions regions.txt -f reference.fa input.bam -out results.vcf
Once launched, FreeBayes starts monitoring the progress of the job. As the job may take several days, the connection
may break, or you may need to close the connection. This does not harm the actual computing task. Once all sub-jobs have completed, you can use command freebayes-puhti-recover to collect the results. For example:
freebayes-puhti-recover freebayes_jobnum_tmp
Where freebayes_jobnum_tmp is the temporary FreeBayes directory that was created by the freebayes-puhti command in the same directory where the command was launched.