Aligning Multiple paired end files together

Aligning Multiple paired end files together

1

Hi All,

I have 72 paired end .fastq file for which i need to do Alignment using BWA. Since its a paired end data and my files are named as

  1. sam_001_1.fastq

  2. sam_001_2.fastq

  3. sam_002_1.fastq

  4. sam_002_2.fastq & so on

Since its a paired end data i am not able to generate a single .sam file using a shell script. my script is generating individual .sam file for each .fastq files.
Please let me know whats the best strategy for handling this.

EDIT:1
My script

for i in *.fastq;

do

bwa mem -t 10 /data2/Exome/HG19/BWA/hg19 1.$i 2.$i > /data2/validation_samples/fastq2sam/$i.sam; 

done

Sincerely,

Dave


assembly


alignment


bwa


shell

• 2.1k views

Read more here: Source link