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