Variant Calling with Multiple Individuals in Freebayes

Variant Calling with Multiple Individuals in Freebayes



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!





Read more here: Source link