Thread: maq map problem
View Single Post
Old 03-12-2010, 09:12 AM   #1
prm36
Member
 
Location: Edinburgh

Join Date: Mar 2010
Posts: 16
Default maq map problem

Hi,

I am attempting to assemble SOLEXA reads from a bacterial genome.
I have been able to assemble the reads from the raw fastq file.
I have filtered and trimmed the fastq file to remove low quality reads to see if this improves my assembly. However, I receive the following message when attempting to run maq map:

paul@bioinf:~/Desktop/TEST/MAQ/filtered$ maq map out.aln.map ref.bfa 1filtered.bfq 2filtered.bfq > out.map.log
-- maq-0.7.1
[ma_load_reads] loading reads...
[ma_load_reads] set length of the first read as 33.
[ma_load_reads] set length of the second read as 33.
maq: read.cc:61: longreads_t* ma_load_reads(void*, int, void*, int): Assertion `strncmp(name, lr->name[j], tl-1) == 0' failed.
Aborted

Can anyone shed some light on what the problem may be?

Cheers

Paul
prm36 is offline   Reply With Quote