Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Tablet has a "read concordance" mode for viewing strand direction on reads. Both pairs of reads that are consistent with a mapped direction following the reference sequence are coloured green (customisable), and reads that are consistent with the reverse complement of the reference sequence are coloured red.
-
Originally posted by sowmyai View PostLooks like Jim Robinson had it backwards for the color legends. Actually
red=positive,
blue=negative
He has clarified here:
https://groups.google.com/forum/#!to...lp/YiVzwnUuOZM
I used Illumina dUTP-type stranded library prep and fr-firststrand with Tophat2. I find several places that say dUTP libraries should use fr-firststrand with Tophat, for example: http://onetipperday.sterding.com/201...pe-to-use.html
Is anyone able to clear this up?
Thanks
Leave a comment:
-
Looks like Jim Robinson had it backwards for the color legends. Actually
red=positive,
blue=negative
He has clarified here:
Leave a comment:
-
Originally posted by cliffbeall View PostI am fairly sure fr-firststrand means read one (aka forward read or fr) corresponds to the reference genome sequence (firststrand). fr-secondstrand means read one corresponds to the reverse complement (secondstrand) of the reference genome. FR and RF would be a similar concept - but they are just making a shorthand where if the read one or forward read corresponds to the reference strand it maps to the left of the reverse read on a linear representation of the genome (FR) and if it corresponds to the reverse complement of the reference it maps to the right of the reverse (RF).
I hope that makes sense -- you might need to go back to DNA/genome structure and how Illumina works to really understand the above.
First 'fr' does not mean 'forward read'. What it means is the the two reads are arranged in the FR (Foward/Reverse) as opposed to the RF (Reverse/Forward) orientation. FR means that the read pairs are oriented with the 3' end point toward each other when aligned to the reference and RF means they are pointing away from each other.
Code:FR: ------------------------------------------------------------------------- -----------------> <------------------- RF: ------------------------------------------------------------------------- <----------------- ------------------->
Code:mRNA 5'-----------------------------------------------------------AAAAAAAAAAAA R2 ------------> <------------ R1
Leave a comment:
-
Originally posted by ArghavanAri View PostThanks for the information.
But I still cannot understand the difference between fr-firststrand & fr-secondstrand in Tophat2 (or RF & FR in Hisat).
Can anyone help me with that?
I hope that makes sense -- you might need to go back to DNA/genome structure and how Illumina works to really understand the above.
Leave a comment:
-
Thanks for the information.
But I still cannot understand the difference between fr-firststrand & fr-secondstrand in Tophat2 (or RF & FR in Hisat).
Can anyone help me with that?
Leave a comment:
-
As others have mentioned, you need to consider the direction of mapping as well as the read number. For a standard Illumina strand-specific sample preparation, the reverse read is the first read:
- Reverse mapping, First read -- consistent with a gene model reading left to right on IGV
- Reverse mapping, Second read -- consistent with a gene model reading right to left on IGV
- Forward mapping, First read -- consistent with a gene model reading right to left on IGV
- Forward mapping, second read -- consistent with a gene model reading left to right on IGV
Just in case you weren't aware, Tablet (since v1.14.10.20) has a concordant / directional colouring mode that colours sequences appropriately.Last edited by gringer; 09-14-2015, 07:31 PM.
Leave a comment:
-
The two genes are on different strands, you can see by the little chevrons in the annotation track at bottom, so that is expected.
Leave a comment:
-
Hi,
I am trying to align and count the reads of paired-end strand-specific RNA-seq data using Hisat and htseq-count.
This is how my IGV look like. https://www.dropbox.com/s/snpapn051b...ness.jpeg?dl=0
(Color alignment by "first-of-pair strand") (view as pairs).
I am confused now. Why all the read pairs are red for NFE2L2? and clur for HNRNPA3?
I know that my data is strand specific, and that this color is showing the first pair read. So why both pairs have the same color?
Also can anyone tell me what options I should use for Hisat (or Tophat2) and htseq-count?
Thanks,Last edited by ArghavanAri; 09-14-2015, 07:18 AM.
Leave a comment:
-
So is directional RNA-seq the same as strand-specific RNA-seq? any clarification will be greatly appreciated.
Thanx
Leave a comment:
-
Originally posted by Jim Robinson View Postred = neg strand (first read of the pair is no the neg strand)
blue = pos strand
grey = pair information not available, so "first of pair" is undefined
Leave a comment:
-
No, even in a non-strand-specific library, the two reads of a pair should point towards each other (or away from each other if it was a library-prep protocol with circularization). You should check whether it is predominantly the first read that point in the sense direction of the gene it aligns to.
Leave a comment:
-
Thanks for the info Jim!
So if most of the pairs in a contig contain reads that are like this:
Code:=====>----------<===== (red) (blue)
Am I right?
What if both of them are red or blue?
Leave a comment:
Latest Articles
Collapse
-
by seqadmin
The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...-
Channel: Articles
08-27-2024, 04:44 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:02 AM
|
0 responses
10 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:02 AM
|
||
Started by seqadmin, 09-03-2024, 08:30 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
09-03-2024, 08:30 AM
|
||
Started by seqadmin, 08-27-2024, 04:40 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
08-27-2024, 04:40 AM
|
||
New Single-Molecule Sequencing Platform Introduces Advanced Features for High-Throughput Genomics
by seqadmin
Started by seqadmin, 08-22-2024, 05:00 AM
|
0 responses
360 views
0 likes
|
Last Post
by seqadmin
08-22-2024, 05:00 AM
|
Leave a comment: