This is a test version of Biostars. For the public version, visit https://www.biostars.org.
samtools mpileup and samtools tview give different reads quality results

enter code hereI want to use samtools mpileup to achive batch extraction of variant quality information from CRAM files

I start with a CRAM on which i do samtools mpileup -B -f "$reference_genome" -r "$region" -q 0 -Q 0 "$cram_file" And i get :

chr11 | 56094117 | G | 25 | FALSE | ..A,..,A...AAa,,.,..,,,., | ?!??]!]?!???!]!]?]?!]????

But when i use samtools tview --reference -p "$cram_file" Then i get : enter image description here


can clearly see the samtools tview first line :there are 4A 1a (no color and no underline ),it means these 5 reads are Q > 30. But when i caculate the samtools mpileup ASCII code ,there are just 3A 1a (Q > 30). Which one has the truest reads quality results? Has anyone else had this happen to them?

samtools tview mpileup

But when i use samtools tview --reference -p "$cram_file" Then i get :

there is nothing to see here but the coordinates. Furthermore the arg for reference is missing.

Sry!!!I've uploaded a picture to make sure it's visible!!

0 answers

No answers yet.

Log in to answer this question.