Question: How to anntate small-RNA seq reads distribution
1
gravatar for i19870503
3.5 years ago by
i1987050310
Australia
i1987050310 wrote:

I use bowtie mapping the smallRNA-seq, but I want to annotate the percent of reads with miRNA, snoRNA, rRNA, tRNA and repeat.

I read the paper is use blast+ and rfam database. Now I don't know which to download in the rfam and how to run the program that make sequence reads classified to major non-coding RNA groups can be viewed as miRNA, snoRNA, rRNA, tRNA and repeat.

rna-seq • 1.6k views
ADD COMMENTlink modified 3.5 years ago by Ashutosh Pandey11k • written 3.5 years ago by i1987050310
0
gravatar for Ashutosh Pandey
3.5 years ago by
Philadelphia
Ashutosh Pandey11k wrote:

If you have already mapped these reads against the reference genome, you can download genomic coordinates for miRNA, tRNA, rRNA, snoRNA etc using UCSC table browser (http://genome.ucsc.edu/cgi-bin/hgTables?command=start) and can use bedtools or HTseq to count or annotate reads aligned to different small RNAs. Though UCSC annotation for some species of small RNAs won't be as comprehensive as some other resources that are exclusively created for that specie of RNA but it should be able to give you close to the optimal results.

ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by Ashutosh Pandey11k

Thanks for reply, My species is pig, so the annotation is not very suffcient. I think blast is more acceptable method.

ADD REPLYlink modified 3.5 years ago • written 3.5 years ago by i1987050310

Hi,

I am planning to do the same requirement as yours. Annotate the percent of reads with miRNA, snoRNA, rRNA, tRNA and repeat for sequenced reads from pig. Would you provide any recommendations?

ADD REPLYlink written 8 days ago by bioinforesearchquestions120
Please log in to add an answer.

Help
Access

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