How to present mapping result ?

Question

[BioGeek]

I usually do samtools flagstat to represent the mapping stats. Is there any best way to represent/present the mapping quality/results in details.

[Kevin Blighe]

Did someone ask you for something in particular in relation to alignment?

Answer 1

[michael.ante]

Most aligners have a decent mapping summary log. You could rely on these, or you can try RSeQC's bam_stat.py.

I'm usually not using samtools flagstat, since it is summarising rather alignments not reads.

For comparing different samples, you can use multiQC.

Answer 2

[sophiespo]

I can't tell you specifically but I find it is a good exercise to read current sequencing publications and see how they present their mapping stats. I am always looking for nicer ways to present boring numbers like that.