Which Bioperl Module Can I Use For Parsing/Mining Blast Output?

Question

[Raghul]

Hi to all, I have blastx results (default output). There are some blast hits starting at frame 1, others starting from frame2,frame3, frame -1 etc. I want to collect all query IDs that blast in Frame=1, Frame=2, etc separately.

I think it is possible using Bioperl. All kind of suggestions are welcome. If somebody has similar bioperl programs, u can post it. So I will get an idea how to start write my own script! I am giving you a portion of blastx output if u want to know more.

THANK YOU raghul

Query= isotig00043|appr-1-p processing domain protein Length=2377

contig0391611122816230- 144 5.6e-34 gi|20178254|sp|Q96XY5|Y2383_SULTO Length=230

Score = 62.4 bits (150), Expect(2) = 2e-016 Identities = 37/123 (30%), Positives = 60/123 (48%), Gaps = 10/123 (8%) Frame = -2

Query 474 KYVFHTVGPMCQKGKLDLRKESKQLKSCVEEVLKLMIKKNCKSICLLAISTGIFNFPLTN 295 KY H VGP+ G KQL+ + K + C+S+ + AIS+GIF +PL
Sbjct 59 KYCIHAVGPIY--GYSGGENAEKQLELTIVNTYKRAEELKCESVSVPAISSGIFGYPLKE 116

Answer 1

[Michael Schubert]

http://www.bioperl.org/wiki/HOWTOs

Answer 2

[Yogesh Pandit]

Link to the SearchIO module from Michael's answer http://www.bioperl.org/wiki/HOWTO:SearchIO

Similar question was asked in the SEQanswers forum: http://seqanswers.com/forums/showthread.php?t=6484

Answer 3

[Scott Cain]

And, should Michael's excellent suggestion not work out: http://lists.open-bio.org/mailman/listinfo/bioperl-l

[Michael Schubert]

should there really be a problem that is not easily resolvable with the basic tutorials, the question is welcome here as well ;-)