7
$\begingroup$

I was following a tutorial on "Tuxedo Genome Guided Transcriptome Assembly Workshop" and was wondering how to interpret the following:

enter image description here

From what I understand from 'Color Legends', the color blue represents something that is below normal and red represents something that is above normal.

Also, my understanding is that merged.gtf is the reconstructed transcript and genes.bed is the reference. They look to be similar for the most part, but seems to have breaks between. What's the appropriate way of interpreting this result?

$\endgroup$
6
$\begingroup$

You've apparently colored your alignments by read strand. In this case, red indicates "+ (watson) strand" and blue indicates "- (strand) strand". This strand association is determined by the orientation of read #1 in a pair (or just "the read" if you have single-end data). This isn't actually documented in the user manual, but you can find one of the Broad institute software engineers replying here.

Your next question is likely why the color of the reads and the orientation of the genes are opposite of each other. This is because in modern directional protocols, read #2 in a pair matches the orientation of the gene from which it arose.

$\endgroup$
2
  • $\begingroup$ Thank you for your reply. I'm not sure I understand what you're referring to by read #1 and read #2 based on the picture I attached. $\endgroup$ Sep 6 '17 at 2:41
  • $\begingroup$ Ignore that part if you don't understand it. The direct answer to your question is that the colors represent the strand from which the reads arose. $\endgroup$
    – Devon Ryan
    Sep 6 '17 at 6:54

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.