This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract quality scores from aligned BAM files

Hi guys, maybe someone already asked the same but I looked around without success. I have a list of bam files resulting from bwa-mem. I would like to extract the list of scores (quality) of each bam in order to look at the quality of each sample (in order to do some statistics). I also would like to know the min and max of the scores for each bam. I tried to use samtools pileup but it generates a file with a lot of information I do not want and also I don't know how to extract the list of scores. Can anyone help me please?

Thank you in advance

chip-seq
• 4,458 views • link
written by elb elb

You are talking about mapping quality?

• 0 views • link
written by Arup Ghosh Arup Ghosh

Yes MAPQ number. Eg: 60

• 0 views • link
written by elb elb

Those should be in column 5. A combination of samtools view with cutting that column out should give you that info. Find the min and max after sorting the file. Should all be possible with (awk, cut, sort, uniq).

• 0 views • link
written by GenoMax GenoMax

1 answer

From all the quality assessments you can make in ChIP-seq MAPQ is probably the one that is least meaningful. You should take all reads that are above a reasonable threshold like 20 and then call peaks, check fractions of reads per peaks and produce bigwigs to check data on a browser track. This is way more meaningful than any NGS statistics that are related to mapping quality.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this