View Single Post
Old 06-09-2010, 09:19 AM   #2
epigen
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 101
Default

The base quality is also a place holder just as the * in the reads column.

I found similar occurrences in my data (created with samtools pileup -vcf), e.g.

chr14 65392038 T A 0 0 0 1 * W

Note that the weird consensus A must result from a floating underflow in the MAQ SNP calling model (http://sourceforge.net/apps/mediawik...?title=SAM_FAQ). Although the FAQ answer says that this only happens in repetitive regions, the reference T is upper case and I just checked again in the UCSC Genome Browser for hg19 that there is no repeat. So I guess the reason is that there is only one read.
epigen is offline   Reply With Quote