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@1.bfq'
-- 2500 sequences were loaded.
-- CMD: (cd ./test; /usr/local/bin/maq map  -n 2 -e 70 -u unmap1@1.txt aln1@1.map ref.bfa read1@1.bfq 2> aln1@1.map.log)
sh: line 1: 29654 error                 /usr/local/bin/maq map -n 2 -e 70 -u unmap1@1.txt aln1@1.map ref.bfa read1@1.bfq 2> aln1@1.map.log
** fail to run command '(cd ./test; /usr/local/bin/maq map  -n 2 -e 70 -u unmap1@1.txt aln1@1.map ref.bfa read1@1.bfq 2> aln1@1.map.log)' at /usr/local/bin/maq.pl line 842.

 

Best,

 

Jingen

 

gen0224，jingen.li87@gmail.com

2013-1-10