This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I creating multisample VCF file ?

Hello,

I want to create a multisample VCF file. I have bam files from various alignment.

I ran the command

bcftools mpileup -d 100000 -f  ~reference  *.bam | bcftools call -c  > concate.vcf

Is this step scientifically correct? Or do I need to merge all bam first before running the above command?

Thanks

vcf
• 976 views • link
updated by Ram Ram • written by eennadi • eennadi

Is this step scientifically correct?

yes , althought most people use -m , not -c http://www.htslib.org/workflow/wgs-call.html

   -c, --consensus-caller          The original calling method (conflicts with -m)
   -m, --multiallelic-caller       Alternative model for multiallelic and rare-variant calling (conflicts with -c)
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this