How to call somatic mutations from vcf files
0
0
Entering edit mode
10 days ago
DS ▴ 20

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 • 294 views
ADD COMMENT
0
Entering edit mode

samtools mpileup is deprecated for calling VCF. use bcftools mpileup. https://samtools.github.io/bcftools/howtos/variant-calling.html

ADD REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

I got this error when running bcftools mpileup

[E::main] unrecognized command 'mpileup'
ADD REPLY
0
Entering edit mode

your version of bcftools is just too old.

Don't delete questions please.

ADD REPLY

Login before adding your answer.

Traffic: 2091 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6