View Single Post
Old 01-24-2019, 11:46 AM   #661
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,022

Have you verified that the $i variable is expanding correctly?

If you have multiple files to align you should first create an index by doing

Code: ref=/users/chutfilz/data/chutfilz/Dm3_Index/dm3.fa
this will create a "ref" directory with all index files. Don't worry about what is in the directory (bbmap uses it own organization).

Then when you run use
in command line, instead of "ref=". This will avoid re-indexing the genome each time.

If you are using a job scheduler then you should submit each alignment job separately (not the way you have the loop setup, which I assume is submitted as a single job?)
GenoMax is offline   Reply With Quote