View Single Post
Old 08-27-2012, 10:28 PM   #1
Location: Asia

Join Date: Aug 2012
Posts: 14
Default read length a problem when mapping with MAQ


I have a file containing 36.8 million reads, as suggested in the MAQ manual I divided the file so that i had files containing 2 million reads. the minimum length of any read is 49bp, and the maximum i cannot approximate because whenever i input the files, if one file has reads with length greater than 65bp then maq returns a "no reasonable reads found" error.

I used the dynamic trim script to trim down the sequences according to their quality scores, which shortened the reads which before were 65bp long to 49-55bp long. These reads when inputted were promptly accepted and yet again when i input files containing sequences with length more than 65bp after trimming these files were rejected with the no reasonable reads found error.

This leads me to believe that maq may have an maximum read length cap of 60bp?

It would be helpful if someone can shed some more light on the matter.
koustavpal is offline   Reply With Quote