Greetings, I am wondering if anyone knows how to pipe sequence data into usearch. I am trying to use GNU parallel to break-up up and distribute multiple smaller ublast jobs over our small server. I can do this with regular blast, but i get fatal errors when i try it with usearch-ublast.
cat CAM_SMPL_001754_RMRNA_derep.fasta | parallel --block 100k --recstart '>' --pipe /public2/Tony/usearch -db /public2/Tony/CAMERA/RefSeqMicrobial/microbial.nonredundant.all.udb -top_hits_only -threads 16 -blast6out DATA
Any suggestions would be helpful. thanks.
cat CAM_SMPL_001754_RMRNA_derep.fasta | parallel --block 100k --recstart '>' --pipe /public2/Tony/usearch -db /public2/Tony/CAMERA/RefSeqMicrobial/microbial.nonredundant.all.udb -top_hits_only -threads 16 -blast6out DATA
Any suggestions would be helpful. thanks.
Comment