SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
extrat unaligned reads from unmapped.bam (tophat) vivienne_lovely Bioinformatics 1 05-26-2013 10:58 PM
How to filter Shrimp aligner output Giorgio C Bioinformatics 0 03-14-2013 02:36 PM
Error running Shrimp aligner Mahtab Bioinformatics 2 03-03-2013 04:52 PM
unaligned reads of tophat ae_ucla RNA Sequencing 1 04-07-2011 10:06 AM
Tophat options to report unaligned reads and controlling Bowtie options Siva Bioinformatics 0 10-15-2010 07:38 PM

Reply
 
Thread Tools
Old 09-08-2014, 03:17 AM   #1
bioWizz
Junior Member
 
Location: Manipal

Join Date: Jul 2013
Posts: 6
Smile SHRiMP aligner options for unaligned (unmapped) reads

I am trying to align small RNA reads using SHRiMP using the command below

Code:
$ SHRiMP_2_2_3/bin/gmapper-ls -N 2 -o 1 -E input.fasta reference.fasta > output.sam
This gives me the .sam file output of only aligned reads. I also want the file containing unaligned reads. I saw the documentation and found the options

Code:
Output Control
--------------

  [ -E/--sam ]    Select SAM output format.

  [ --sam-unaligned ]    If SAM output format is also selected dump unaligned reads to the SAM output.

  [ -un <filename> ]    Print unaligned reads to target file.

  [ -al <filename> ]    Print aligned reads to target file.
I have tried using these parameters but have ended up obtaining blank files. Can anyone please guide me on the modification of my existing script to incorporate these parameters to get the unaligned output.

Thank you
bioWizz is offline   Reply With Quote
Reply

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 10:58 AM.


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