samtools mpileup and samtools tview give different reads quality results
0
0
Entering edit mode
6 months ago
cocoboxt • 0

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

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

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