Pandaseq in a bash for-loop - Error nofile
I'm merging R1 and R2 files derived from 16S rRNA miseq amplicon-sequencing, and I'd like to perform it on several samples at once. This is my script:
for r1 in *R1*;
do
r2=${r1/R1/R2};
pandaseq -f $r1 -r $r2 -o 10 -w ${r1/R1_001.sic.00.0_0.cor.fastq/pandaseq.fasta};
done
The screen output says this for all my R1 files:
ERR NOFILE Order-89-1-1_S1_L001_R1_001.sic.00.0_0.cor.fastq
Too confused to continue.
Try -h for help.
Is there something missing in that script? Thanks a lot!
1 answer
I found a long solution, but it works. It was a matter of working on that last in-line replacement for the output (w).
#1
for r1 in *R1*;
do
r2=${r1/R1/R2};
pandaseq -f $r1 -r $r2 -o 10 -w ${r1/R1/pandaseq};
done
#2
mkdir pandaseq
#3
mv *pandaseq.fastq pandaseq
#4
for r1 in pandaseq/*fastq;
do
mv $r1 ${r1/fastq/fasta};
done
This generates now a merged file for each sample, kept in a separate folder.
Log in to answer this question.
More posts like this
-
Need help in using cutadapt to trim paired end fastq reads
written by George •
-
EBI-metagenome - Unequal number of reads in introduction and taxonomy
written by agata88
-
News: Workshop: 16S rRNA gene Metabarcoding (Berlin. 3-7 April 2017)
written by Physalia-courses
-
How do i should proceed with 16S rRNA amplicon sequencing data from Illumina MiSeq using QIIME pipe…
written by antonioggsousa
#TIL in-line replace in bash variables (the
{//}syntax). Thank you.Hi Fibar..I have the same problem.. Can you please let me know your email..so i can write in detail
Regards Mukil
So the solution mentioned below did not work for you?
@fibar has not been seen on Biostars for over a year so your best bet may be to create a new post/thread to describe your problem.
@fibar: post example file list. File list you have and code you are using. please don't post screenshots unless it is error.
Hi, I'm sorry for the poor post. @mxm2623 Please read the solution below.
@fibar.. I am confused. Was this thread edited? Looking at my post, it seems it is unnecessary to tag you in my post (9 days ago). If so, I am sorry.
No you did not tag Fibar. You can't tag a user by just using "@"name. You need to paste their biostar profile link URL in a post.
@Fibar probably had "follow via email" set up for this thread and thus came back to respond.
@genomax! My concern was i unnecessarily directed my post @Fibar for some unknown reason. For that I was sorry.