Off topic:How to loop over multiple files for BWA MEM in bash for loop?
0
0
Entering edit mode
15 months ago
DNAngel ▴ 120

Hi all,

I've written a large bash script for single-end paired reads which uses bwa mem. It works great. But now I have a new dataset that is paired-end and the only part I need to change is the beginning where bwa mem uses a reference file and sample files to align them. I cannot figure out how to edit my for loop to allow this.

All my sample files have a forward_R1 while and a reverse_R1 file as they are paired files.

My regular bwa mem loop in bash is:

#!/bin/bash
REF=$1
ID=$2

for files in *.fastq
do
     bwa mem -B 2 -M-t 40 $REF $files > ${files%%.fastq}_aln.sam
done

How do I edit this so that it can take two files for the program, but still output the one aligned file?

I tried:

#!/bin/bash
REF=$1
ID=$2

for forfiles, revfiles in *.fastq
do
     bwa mem -B 2 -M -t 40 $REF $forfiles $revfiles > ${forfiles%%.fastq}_paired_aln.sam
done

But this gives an error because the loop doesn't understand two arguments.

bwa bash • 1.5k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 1135 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6