foreach run (ERR145651_s ERR145651_t)
/home/sigenae/bin/runSampe -r ${run}_trinity/Trinity.fasta -1 ${run}_R1.fastq.gz -2 ${run}_R2.fastq.gz -o ${run}_trinity/RMBT
/home/sigenae/bin/runExonerate -q ../Data/Danio_rerio.Zv9.pep.chr3.fa -t ${run}_trinity/Trinity.fasta -o ${run}_trinity/exo -p '-m protein2dna'
end
qarray -N blat -sync y -e err_log -o err_log post_asm.sh
cat ERR145651_s_trinity/fasta_vs_genome.psl | /home/sigenae/bin/blat_parser.pl -o psl > ERR145651_s_trinity_vs_genome.psl
cat ERR145651_t_trinity/fasta_vs_genome.psl | /home/sigenae/bin/blat_parser.pl -o psl > ERR145651_t_trinity_vs_genome.psl