Question: (Closed) How to extract reads belonging to sense intron/exon and antisense intron/exon while working with small RNA sequence
0
gravatar for toralmanvar
11 months ago by
toralmanvar770
toralmanvar770 wrote:

I have Illumina data of 1 x50bp from small RNA library. I am interested in identifying known and novel miRNAs in my plant, whose reference genome is not available. However, first I want to count the reads aligning Exon sense, Exon antisense, Intron sense, Intron antisense. Basically, I want to replicate the information shown in table below. I have gone through the complete publication and they have mentioned the use of Rfam database in "Experimental procedures" section for identifying mRNA. So my question is how can I classify reads as Exon sense, Exon antisense, Intron sense, Intron antisense using Rfam or any other databases? enter image description here

ADD COMMENTlink written 11 months ago by toralmanvar770
1

Hello toralmanvar!

You already has a very similar question posted (Extracting reads belonging to different RNA families from Rfam ). You can edit that question (make no changes if not needed) and it will be bumped to main page.

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

ADD REPLYlink modified 11 months ago • written 11 months ago by genomax68k

According to me, both the questions are totally different. Previous question was about different non-coding RNA classification, whereas this question focuses on type of mRNA classification. However I agree that the reference publication table which I have provided is same for both.

ADD REPLYlink written 11 months ago by toralmanvar770

hello genomax, can you tell me the reason why you think that it is the similar post to my previous post?

ADD REPLYlink modified 11 months ago • written 11 months ago by toralmanvar770

Here is the text of your two posts:

I have Illumina data of 1 x50bp from small RNA library. I am interested in identifying known and novel miRNAs in my plant, whose reference genome is not available. Here first I want to count the reads aligning to non-coding rRNAs like tRNA, rRNA, snoRNA, snRNA, siRNA etc except miRNAs. After removing those reads, I would like to align remaining reads with miRBase using blastn or bowtie to identify known miRNAs. Now problem is that, I want to use Rfam for detecting reads of different RNA class, thus I concatenated the fasta files of 2686 families (2,487,655 sequences) but Rfam v13 sequences does not have proper information of the class they represent in their header. So my concern is how can I count and extract the reads mapping to different classes of RNAs using Rfam.

and

I have Illumina data of 1 x50bp from small RNA library. I am interested in identifying known and novel miRNAs in my plant, whose reference genome is not available. However, first I want to count the reads aligning Exon sense, Exon antisense, Intron sense, Intron antisense. Basically, I want to replicate the information shown in table below. I have gone through the complete publication and they have mentioned the use of Rfam database in "Experimental procedures" section for identifying mRNA. So my question is how can I classify reads as Exon sense, Exon antisense, Intron sense, Intron antisense using Rfam or any other databases?

It looks like you are referring to the same experimental dataset in both posts. You have added the request of how can I classify reads as Exon sense, Exon antisense, Intron sense, Intron antisense.

You can always edit the first post and add this a second requirement. This will automatically bring your post back up to main page. It is my understanding that if you did a smallRNA prep then it is not going to work for normal mRNA discovery.

We make an effort to keep threads focused on topics (with answers in one place) so the information is easy to find for others who search for it in future. That is the only reason this post has been closed with reference to your previous one.

ADD REPLYlink modified 11 months ago • written 11 months ago by genomax68k
Please log in to add an answer.
The thread is closed. No new answers may be added.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1684 users visited in the last hour