Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Prettyprint and SHRiMP 2.2.2 evilution84 Bioinformatics 1 02-20-2012 01:58 PM
BWA parameters for mRNA-seq aligning against mRNA refseq kwicher SOLiD 1 09-19-2011 04:45 AM
bwasw scoring parameters when aligning 454 RNA-seq reads on the human genome dgacquer Bioinformatics 0 09-02-2010 06:24 AM
to aligning all reads with SHRiMP (as automatic) ahabnar Bioinformatics 0 11-04-2009 10:43 AM
SHRiMP doubts seq_GA Bioinformatics 1 06-22-2009 09:53 PM

Thread Tools
Old 10-18-2012, 02:07 PM   #1
Location: USA

Join Date: Mar 2012
Posts: 40
Default SHRiMP - Which parameters for aligning to miRNAs?


I recently read a few reviews that point to SHRiMP as the short read aligner of choice when mapping small RNA-seq reads to a database of miRNAs.

SHRiMP even includes a miRNA parameter, but is very vague as to what the parameters do and how to optimize. I've emailed the dev team but never received a response.

I'll show everyone what I tried for my first attempt mapping these small RNA-seq reads to a miRNA database, and please offer suggestions for optimization or correct any errors in my command.

I first projected the database:
SHRiMP-2.2.3/utils/ --seed 00111111001111111100,00111111110011111100,00111111111100111100,00111111111111001100,00111111111111110000 --h-flag --shrimp-mode ls refs/all_mature_miRNAs.fa

I then mapped to the database:
SHRiMP-2.2.3/bin/gmapper-ls -L all_mature_miRNAs-ls --fastq --qv-offset 33 my_reads.fastq --mode mirna >map.out

According to the SHRiMP documentation, passing the --mode mirna parameter does the equivalent of entering the following:
"loading the 5 seeds mentioned above; plus '-H -n 1 -w 100% -U -a 0 -g -255 -q -255 -Z'"

If you have any experience using SHRiMP for this purpose, I'm eager to hear your advice. Thanks!
all_your_base is offline   Reply With Quote
Old 03-14-2013, 05:21 PM   #2
Giorgio C
Location: ITALY

Join Date: Oct 2010
Posts: 89

Hi, I know it's an old post but I was wondering if you finally get the right way, in terms of parameter, for the analysis of miRNA in Shrimp.

Since I have to start from the beginning and there is not so much documentation out there, I want to be sure of what I'm going to do.

Hope to have some advices from you.

Thanks in advance,
Giorgio C is offline   Reply With Quote

mirna, shrimp, small rna-seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:05 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO