SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   SHRiMP aligner options for unaligned (unmapped) reads (http://seqanswers.com/forums/showthread.php?t=46471)

bioWizz 09-08-2014 03:17 AM

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


All times are GMT -8. The time now is 02:42 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.