Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Thanks dpryan!
I do have another question if multiple alignments is specified in the bowtie alignment (-a) will infer_experiment be able to work correctly in predicting strandedness?
Leave a comment:
-
Either the paper was wrong or their kit failed completely. My guess is that the paper was wrong and that they used an unstranded kit.
Leave a comment:
-
Thanks guys! this was very helpful! The paper says an Illumina TruSeq Stranded mRNA LT kit so based on what Michael said it is to be expected that the read are the reverse complement of the gene. Another sequencing set from the same paper ( which was also indicated to use the Illumina TruSeq Stranded mRNA kit) gave me the following RSeqQC results:
Fraction of reads explained by "++,--": 0.4775
Fraction of reads explained by "+-,-+": 0.4956
Does this mean that maybe the strandedness part of the kit didn't work as well? Or maybe a non stranded kit was used?
Leave a comment:
-
The "+-,-+" syntax is described here and has got to be about the most cryptic way of conveying things that's possible. What this means is that you have strand-specific data made with the dUTP method. This is the most common type these days and means that if you need to do counting, that you should use the "-s reverse" option in htseq-count or "-s 2" in featureCounts. If you happen to be using tophat2 for alignment, you want the "fr-firststrand" library type.
Edit: I should have refreshed, Michael beat me to it!
Leave a comment:
-
Hi Leanne,
The result shows you, that you have a stranded library. It also means, as you already said, that the reads are the reverse complement of the transcripts.
A couple of library preps produce reads in that flavour:
TruSeq stranded, ScriptSeq, Sense, ....
Just use the strandedness as detected in your analyses (e.g. htseq-count ... -s reverse ..).
Cheers,
Michael
Leave a comment:
-
RSeqQC infer_experiment strandedness
Hi,
I am using RSeqQC infer_experiment.py to look at whether or not a library is stranded and I got the following results:
This is SingleEnd Data
Fraction of reads failed to determine: 0.0242
Fraction of reads explained by "++,--": 0.0398
Fraction of reads explained by "+-,-+": 0.9360
What does this mean? it appears that the reads are mapping to the reverse complement strand relative to the gene, what kind of sequencing protocol would do this? Does it have a name?
Thanks for any help,
LeanneTags: None
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: