View Single Post
Old 05-19-2012, 02:19 PM   #8
Location: United States

Join Date: Jul 2011
Posts: 10

I am currently trying to deal with exactly this problem: CASAVA alignments that choke mpileup. We use mpileup output directly for analysis with our own SNP finders, and it works very well with output from a BWA-based aligner (stampy). But with CASAVA output, as stated in another response, mpileup thinks everything is a SNP, all reads with phred scores of 0. This is not the case, fortunately, as the same read data run through stampy produces an alignment file that behaves well with mpileup.
frc1230 is offline   Reply With Quote