View Single Post
Old 09-17-2019, 12:00 PM   #13
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,030
Default

Output fasta reads were interleaved (put in one file) since you had paired-end reads to begin with. To separate them into two files you should use the command below.

Code:
for i in *.fa; do name=$(basename $(i} .fa); reformat.sh in=${name}.fa out1=${name}_R1.fa out=${name}_R2.fa ; done
You could also do the original conversion by providing "out1=" and "out2=" separate file names during BAM conversion to avoid this step.
GenoMax is offline   Reply With Quote