Question: Bulk quality interleaving with bbmap reformat command
0
gravatar for infenit101
3 days ago by
infenit10110
infenit10110 wrote:

I have several fastq.gz files that look like this:

Sample1_R1.fastq.gz Sample1_R2.fastq.gz

Sample2_R1.fastq.gz Sample2_R2.fastq.gz

etc...

The output needs to be: Sample1_interleaved.fastq.gz, Sample2_interleaved.fastq.gz

I am still learning how to loop commands so bare with me. This is what I tried to run:

for i in `ls -1 *_R#.fastq.gz | sed 's/_R#.fastq.gz//‘`
do
reformat.sh in=$i\_R#.fastq.gz out=$i\_interleaved.fastq.gz 
done

However this did not work. Can someone help me? Many Thanks!

bash bbmap bbduk loop reformat • 71 views
ADD COMMENTlink written 3 days ago by infenit10110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1649 users visited in the last hour