This is a test version of Biostars. For the public version, visit https://www.biostars.org.
alignment summary from tophat

Hi

I would like to know how can I get the summary of alignment form tophat output. I have looked at the alignment_summary.txt generated by tophat, but I want to know how many reads aligned to two or five locations in the genome and how many aligned to six or more location?

Thanks in advance

rna-seq
• 3,723 views • link
updated by Ram Ram • written by hana hana

2 answers

Not sure there is an 'easy way' to do this. You have to parse the accepted_hit.bam (you can view it using samtools) and count the number of repeating mappings (see http://samtools.github.io/hts-specs/SAMv1.pdf to understand the format.)

• 0 views • link
updated by Ram Ram • written by roy.granit roy.granit

Try using QualiMap to analyze the resulting alignment. There are 2 modes possible: RNA-seq QC or BAM QC with GTF file.

• 0 views • link
updated by Ram Ram • written by madk00k madk00k

Log in to answer this question.

More posts like this