how to use bash script to make iterative loop through directory with two file types

[jddavis]

Well you've all been wonderful about answering my novice questions...so here is another. I thank you in advance for generously providing an answer.

What I'm trying to do is loop through a directory, using simple bash script, to a set of files and then process the files to a *.sam format using bwa...here's the code I have so far...

#!bin/bash
#basic syntax below
#bwa samse ref.fa aln_sa.sai short_read.fq > aln-se.sam

module load bwa

FILE1="$WORK/trimmed/*se.fq"
FILE2="$WORK/trimmed/*.sai"

for f in FILE1 and FILE2
do
bwa samse whole_genome.fa "$f2".sai "$f1".fq > "$f".sam
done

[mastal]

if your *.sai and *se.fq files have the same prefix (presumably you have one file of each type for each bwa alignment you want to run), you could try:

FILE2="$WORK/trimmed/*.sai"
WORKDIR="$WORK/trimmed/"

for f in $FILE2
do

      prefix=`basename $f .sai`

      bwa samse whole_genome.fa ${WORKDIR}${prefix}.sai     ${WORKDIR}${prefix}se.fq > ${WORKDIR}${prefix}.sam

done

[jddavis]

Hi Mastal:

Thank you for the reply, and I will certainly try this out...they do have similar prefixes...

-Jennifer