This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Map Smallrna-Seq On Mirna

Hi,

I've smallRNA-seq data and want to count the number of read for each miRNA. I'm studying samples from a non-well annotated species so I've to use mirbases sequences from a related species. My idea is to use bowtie to map the reads on miRNA mature sequences. But is bowtie a good idea to align small reads on small sequences ? Anyone has an another idea ?

Thanks

N.

mirna

1 answer

See How can I accurately identify miRNA sequences from small RNA seq results

Short answer: For tallying exact matches to known mature miRNAs, a simple homemade string-matching script will suffice. But to also capture isomirs / processing variants / non-templated tailed species in your counts, an aligner is needed. Bowtie 1, Bowtie2, and SHRiMP2 have been reported to be fine.

thanks. I didn't see the other post, sorry ;)

Log in to answer this question.