[Bioperl-l] Parse BLAST to find mismatches

Ken Graham kjgraham@ucdavis.edu
Fri, 24 May 2002 13:44:29 -0700


I'm stuck on what might be a simple problem. Using a BLAST report I want to 
find the position, on the query sequence, of each mismatch in a HSP. I can 
read in the BLAST report, and get properties on the HSP (such as start, 
percent similarity etc..) but I don't know where to find any info about the 
mismatches. I tried to use the SimpleAlign object but I just get errors.
Am I missing something obvious?
Thanks

Ken Graham
UC Davis School of Medicine
Biological Chemistry