This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Variant Calling with Multiple Individuals in Freebayes

Hello,

I am currently a new student in Bioinformatics and was trying to do variant calling with Freebayes software. I have 10 individuals and also have a reference genome. I have already sorted and indexed all my .bam files ready for variant calling. I am able to successfully produce a .vcf output file for each individual with target information.

I want to run the freebayes script in order to perform variant calling on all 10 individuals simultaneously and produce one .vcf output file containing information from all 10 individuals. At this point I am quite lost and was hoping for some direction. Below is the code I use for Freebayes piped with a filter. More specifically, I am hoping to compare number of SNPs and allele frequency between the 10 individuals.

freebayes -f ref.fa aln.bam | vcffilter -f "QUAL > 20" >results.vcf

Any help would be much appreciated!

freebayes variant-calling
• 4,094 views • link
updated by size_t • size_t • written by Laurence • Laurence

try this:

freebayes -f ref.fa -L bamlist --ploidy 2 --vcf result.vcf

-L : A file containing a list of BAM files to be analyzed.

or

freebayes -f ref.fa aln1.bam aln2.bam aln3.bam ...
• 0 views • link
written by size_t • size_t

1 answer

It looks like you can provide multiple bam files and output a gVCF file https://github.com/freebayes/freebayes#usage

• 0 views • link
written by Ming Tommy Tang Ming Tommy Tang

Log in to answer this question.

More posts like this