This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to call somatic mutations from vcf files

I have created two vcf.gz file with samtools mpileup and bcftools index, and then i compare them with bacftools isec

it outputs 4 .vcf files. But how can i filter out somatic mutations from these files?

Thanks

cancer vcf bcf bam gene

it outputs 4 .vcf files.

call all your samples in the same cmd. bcftools mpileup -f reference.fa s1.bam s2.bam s3.bam s4.bam

I got this error when running bcftools mpileup

[E::main] unrecognized command 'mpileup'

your version of bcftools is just too old.

Don't delete questions please.

0 answers

No answers yet.

Log in to answer this question.