foreach run (ERR145651_t ERR145651_t_norm)
/home/sigenae/bin/runSampe -r ${run}_oases/merge/transcripts.fa -1 ${run}_R1.fastq.gz -2 ${run}_R2.fastq.gz -o ${run}_oases/RMBT
/home/sigenae/bin/runExonerate -q ../Data/Danio_rerio.Zv9.pep.chr3.fa -t ${run}_oases/merge/transcripts.fa -o ${run}_oases/exo -p '-m protein2dna'
end
qarray -N blat -sync y -e err_log -o err_log post_asm.sh
cat ERR145651_t_oases/fasta_vs_genome.psl | /home/sigenae/bin/blat_parser.pl -o psl > ERR145651_t_oases_vs_genome.psl
cat ERR145651_t_norm_oases/fasta_vs_genome.psl | /home/sigenae/bin/blat_parser.pl -o psl > ERR145651_t_norm_oases_vs_genome.psl