SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
aligning the reads kjaja Bioinformatics 2 07-07-2019 06:45 AM
Mosaik Aligning with Solexa Reads JimC Bioinformatics 5 06-06-2012 05:07 PM
how to align 15bp length reads by SHRiMP 54016565 Bioinformatics 2 10-18-2011 05:47 AM
SHRiMP - how to obtain unique mapped reads? idonaldson Bioinformatics 10 07-02-2010 02:48 PM
SHRiMP: Accurate Mapping of Short Color-space Reads doxologist Literature Watch 0 05-27-2009 06:29 AM

Reply
 
Thread Tools
Old 11-04-2009, 09:43 AM   #1
ahabnar
Junior Member
 
Location: Montreal

Join Date: Oct 2009
Posts: 1
Default to aligning all reads with SHRiMP (as automatic)

Hi,

I'm new in SHRiMP. I have 10 library and for each library, i had used the splitreads.py of SHRiMP to 1000 files. Then, i have 1000 files csfasta type for on library. I had test some of them with rmapper-cs... to align and it's working but i need to run the same command as automatically for the 1000 files csfasta first and for all 10 library.
The command line is like:
./rmapper-cs -m 8 -i -30 -g -35 -e -7 -x -25 -h 70% -v 65% -M 35bp,fast -o 100 -U /root/Desktop/SHRiMP_1_3_0/readsT/1000_to_1999.csfasta /root/Desktop/DB.fasta > /root/Desktop/SHRiMP_1_3_0/results/res_1000_to_1999.out

I try with command find like this one:

#!/bin/bash
#

find /root/Desktop/SHRiMP_1_3_0/readsT/ -name "*.csfasta*" -exec /root/Desktop/SHRiMP_1_3_0/bin/rmapper-cs -m 8 -i -30 -g -35 -e -7 -x -25 -h 70% -v 65% -M 35bp,fast -o 100 -U "{}" /root/Desktop/docShrimp191009/DB.fasta >> /root/Desktop/SHRiMP_1_3_0/results/"{}"res_U_lib1.out {} \;

but it said bad file and nothing in result file

Do you know what's the problem

ps: {} is just to have the name of the file in result, it does't working but for the input file it's ok

ravelna
ahabnar is offline   Reply With Quote
Reply

Tags
align, automatic, color code cs, reads, shrimp

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 01:56 AM.


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