I just started to use maq to map the paired-end reads of 108bp.
When I use command "maq map -e 200 -a 500 a1.bfq a2.bfq",
I always got this:
-- maq-0.7.1
[ma_load_reads] loading reads...
[ma_load_reads] set length of the first read as 108.
[ma_load_reads] set length of the second read as 108.
[ma_load_reads] 1456472*2 reads loaded.
[ma_longread2read] encoding reads... 2912944 sequences processed.
[ma_match] set the minimum insert size as 109.
[match_core] Total length of the reference: 1586004629
[match_core] round 1/3...
[match_core] making index...
[match_index_sorted] no reasonable reads are available. Exit!
I checked the data using "maq bfq2fastq a1.bfq". It shows the data is as this:
@IL41_3804:3:10:236:221/1
CAGTGATTTCGTCATTTTTCAAGTCGTTAAGTGGATATTCCTCATTTTCCATGATTTTCAGTTTTCTTGCCATATTCCATGTCCCACAGTGGACATTTCTAAATTTTC
+
CCDCDBACBABE<CCDDDCCACA;=DCCCBCACCBCCCA?BB=BCCCB@CACA<CCCC?@B@@CC?CB@<=@@0-3<6=A@@>77;5978A?<;=@?@=A?;:=8;??
I can't find any problem. Anyone has idea about this? Thanks.
When I use command "maq map -e 200 -a 500 a1.bfq a2.bfq",
I always got this:
-- maq-0.7.1
[ma_load_reads] loading reads...
[ma_load_reads] set length of the first read as 108.
[ma_load_reads] set length of the second read as 108.
[ma_load_reads] 1456472*2 reads loaded.
[ma_longread2read] encoding reads... 2912944 sequences processed.
[ma_match] set the minimum insert size as 109.
[match_core] Total length of the reference: 1586004629
[match_core] round 1/3...
[match_core] making index...
[match_index_sorted] no reasonable reads are available. Exit!
I checked the data using "maq bfq2fastq a1.bfq". It shows the data is as this:
@IL41_3804:3:10:236:221/1
CAGTGATTTCGTCATTTTTCAAGTCGTTAAGTGGATATTCCTCATTTTCCATGATTTTCAGTTTTCTTGCCATATTCCATGTCCCACAGTGGACATTTCTAAATTTTC
+
CCDCDBACBABE<CCDDDCCACA;=DCCCBCACCBCCCA?BB=BCCCB@CACA<CCCC?@B@@CC?CB@<=@@0-3<6=A@@>77;5978A?<;=@?@=A?;:=8;??
I can't find any problem. Anyone has idea about this? Thanks.
Comment