Beefy Boxes and Bandwidth Generously Provided by pair Networks
XP is just a number
 
PerlMonks  

Re: parsing mismatch from blast output

by BioLion (Curate)
on Aug 09, 2010 at 14:56 UTC ( #853841=note: print w/ replies, xml ) Need Help??


in reply to parsing mismatch from blast output

Boulder::Blast is probably one standard way of parsing BLASTN results (and other tag/val reports), and will do a good job of simplifying the data held in the record, from where you could split the match and reference sequences into individual characters, and compare them, outputting mismatched characters (and their positions within the match).

You might also want to check out bioperl: parse BLAST HSPs which gives some example code on handling the records, to get at the info you want.

Hope this helps - have a go, and if you get into further trouble re-post and I am sure many a Monk will rush to your aid - remember to include input/code/outpur/error messages etc...!

Just a something something...


Comment on Re: parsing mismatch from blast output

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://853841]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others lurking in the Monastery: (13)
As of 2014-09-17 17:14 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    How do you remember the number of days in each month?











    Results (92 votes), past polls