Samtools mpileup vcf file
0
0
Entering edit mode
6.7 years ago
AW ▴ 350

Hi all,

I would like to know how samtools pileup format calls SNPs.

For my pipeline I

  1. Map Illumina reads onto my reference with Bowtie to get a sam file
  2. Have a custom scripts which filter out non uniquely mapping reads and reads > 5 mismatches
  3. Use samtools pileup to produce a vcf file

My custom script removes reads on the basis of NM and AS. Therefore, the sam flags are no longer meaningful. E.g. if both reads map originally then the flag for read 1 and read 2 might be 83 and 163. However, if I remove read 2 because it has > mismatches then the flag for read 1 will remain 83 (i.e. read paired, read mapped in proper pair) when it is now actually an unpaired read. This only has an impact if subsequent programs use the flag information.

My question is, does samtools pileup read/use the flags in the sam file to produce the vcf file?

Thanks!

bam vcf SNPs mpileup samtools • 1.9k views
ADD COMMENT
0
Entering edit mode

ask developers and in practice use NextGene and Ugene

ADD REPLY
0
Entering edit mode

how do you extract uniquely mapping reads with NextGene or Ugene?

ADD REPLY

Login before adding your answer.

Traffic: 811 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