Hi,
I am a beginner to Maq. Now, I am mapping reads to transcripts by Maq.
But, there ara some errors. I need some help.
-- CMD: /usr/local/bin/maq fasta2bfa /home/transcripts.fasta ./test/ref.bfa 2> /dev/null
-- CMD: /usr/local/bin/maq fastq2bfq -n 2000000 /home//A_1000.fq ./test/read1
-- finish writing file './test/read1(a)1.bfq'
-- 2500 sequences were loaded.
-- CMD: (cd ./test; /usr/local/bin/maq map -n 2 -e 70 -u unmap1(a)1.txt aln1(a)1.map ref.bfa read1(a)1.bfq 2> aln1(a)1.map.log)
sh: line 1: 29654 error /usr/local/bin/maq map -n 2 -e 70 -u unmap1(a)1.txt aln1(a)1.map ref.bfa read1(a)1.bfq 2> aln1(a)1.map.log
** fail to run command '(cd ./test; /usr/local/bin/maq map -n 2 -e 70 -u unmap1(a)1.txt aln1(a)1.map ref.bfa read1(a)1.bfq 2> aln1(a)1.map.log)' at /usr/local/bin/maq.pl line 842.
Best,
Jingen
gen0224,jingen.li87(a)gmail.com
2013-1-10