SEQanswers

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 12:58 PM
BWA parameters for mRNA-seq aligning against mRNA refseq kwicher SOLiD 1 09-19-2011 03:45 AM
bwasw scoring parameters when aligning 454 RNA-seq reads on the human genome dgacquer Bioinformatics 0 09-02-2010 05:24 AM
to aligning all reads with SHRiMP (as automatic) ahabnar Bioinformatics 0 11-04-2009 09:43 AM
SHRiMP doubts seq_GA Bioinformatics 1 06-22-2009 08:53 PM

Reply
 
Thread Tools
Old 10-18-2012, 01:07 PM   #1
all_your_base
Member
 
Location: USA

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

Hi,

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/project-db.py --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, 04:21 PM   #2
Giorgio C
Member
 
Location: ITALY

Join Date: Oct 2010
Posts: 89
Default

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
Giorgio C is offline   Reply With Quote
Reply

Tags
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 09:57 AM.


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