Question: Obtaining Quality score of the mapped sequences - BBMap
1
gravatar for madhu.9124
21 months ago by
madhu.912430
madhu.912430 wrote:

Hi,

I am new to BBMap, so I kindly request to help me to find the **quality score of the mapped sequences. I also want to increase the quality score of the mapped sequence.

Thanks Madhu

alignment • 1.5k views
ADD COMMENTlink modified 4 weeks ago by madzayasodara0 • written 21 months ago by madhu.912430

Are you referring to alignment scores (MAPQ values, column 5)? Those are in the SAM/BAM files. Otherwise this question is missing what it is you are looking for.

ADD REPLYlink modified 21 months ago • written 21 months ago by genomax64k

thanks, I want the MAPQ value. Any idea of the command to get the MAPQ value using SAMtool.

ADD REPLYlink written 21 months ago by madhu.912430
1

What will you do with standalone MAPQ values? You could do something like samtools view test.bam | awk -F "\t" '{print $5}' OR awk -F "\t" '{print $5}' your_file.sam

ADD REPLYlink modified 13 months ago • written 21 months ago by genomax64k

thanks a lot, I just want to retrieve the value and increase the value to test the BBMap application.

ADD REPLYlink written 21 months ago by madhu.912430

missing ' after

awk -F "\t" '{print $5}

ADD REPLYlink written 13 months ago by arup870

Added to post above and here awk -F "\t" '{print $5}' your_file.sam

ADD REPLYlink written 13 months ago by genomax64k

You guys sure there's not a sam/bamtools cmd that straight up does this? I did my alignment using BWA. I could loop through my 23 whole-genome bam files doing this, but It's gonna be real slow...

ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by madzayasodara0
1

the original question is somewhat undefined. What is the purpose of extracting the mapping quality?

If one wanted to filter alignment by a minimal score this can be done with

samtools view -q 10 myfile.bam

where

 -q INT   only include reads with mapping quality >= INT [0]
ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by Istvan Albert ♦♦ 79k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1943 users visited in the last hour