Is there anyway to get some kind of mapping score statistic from GMAP output? A score similar to e-value or mapping quality. It seems like the threshold values used in their oligomer chaining alignment algorithm are hard coded into the software?
Has anyone turned on the --nosplicing flag and used GMAP for aligning genomic sequences instead of mRNA/EST? How does it perform?