[Bioperl-l] Blast > parsing result in Exel

Giorgio C casaburi at ceinge.unina.it
Tue Nov 15 04:25:25 EST 2011


Hy everybody,

in this situation froma blast (-m 1) result file :

Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer,
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997),
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs", Nucleic Acids Res. 25:3389-3402.

Query= 132-291
(59 letters)

Database: Scrivania/orchidea/mature_mirBase.fa
21,643 sequences; 470,608 total letters

Searching..................................................done



Score E
Sequences producing significant alignments: (bits) Value

mtr-miR2644b MIMAT0013413 Medicago truncatula miR2644b 28 0.031
mtr-miR2644a MIMAT0013412 Medicago truncatula miR2644a 28 0.031
gga-miR-1704 MIMAT0007596 Gallus gallus miR-1704 22 1.9
gga-miR-1557 MIMAT0007414 Gallus gallus miR-1557 22 1.9
mmu-miR-880-5p MIMAT0017266 Mus musculus miR-880-5p 22 1.9

132_0 8 cagccgctcagattgatggtgcctacagccttgccagcccgctcagattgat 59
12631 5 .............. 18
12630 5 .............. 18
7826 5 ........... 15
7644 19 ........... 9
5394 3 ........... 13
5394 3 ........... 13
BLASTN 2.2.21 [Jun-14-2009]


Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer,
...
....
..........

______________________________________________________________
I need to parse in an exel sheet :

1)ID 2)Name of the hit 3)E-value 4)Score 5)Species


1) 132-291 2)mir2644b 3) 0,031 4)28 5) Medicago truncatula


Is possible from a big blast result file obtain an exel with 5 columns where
every field is the first hit of the blast result. Can anyone halp me to fix
this problem ??? Also with a little script in perl.


Thank you very much
-- 
View this message in context: http://old.nabble.com/Blast-%3E-parsing-result-in-Exel-tp32846407p32846407.html
Sent from the Perl - Bioperl-L mailing list archive at Nabble.com.



More information about the Bioperl-l mailing list