Hi, I am mapping DNA short reads to UCSC human genome build Hg19. The read length is 101 bp with single end generated from Illumina machine. I have no problem to run match step, but after many hours of localalign job (~1570000 reads have been processed), it is aborted with the error message:
AlignNTSpace.c:40: AlignNTSpaceUngapped: Assertion `readLength <= referenceLength' failed.
I check the reads and all of them are 101 bp long. Could you please help me on this problem?
Thanks a lot.
Log in to post a comment.