BioPython : BLAST comparision and parsing output in particular format

Question

[sukeshssk]

I have query sequence, Suppose query: NNNNNNNNNNNNNNNNNN

• Database 1: Homo sapien
• Database 2: Mycobacterium tuberculosis

I compared query sequence with above two Databases using Standalone BLAST individually and I got result as ex. Result1.txt and Result 2.txt. Now, I wanted to write a python programme which will store first 10 results of BLAST result from both files (Result1.txt and Result2.txt) in following format.

Query     Homo sapien                    Mycobacterium tuberculosis
gi ID     gi ID  (1st compared seq.)     gi ID (1st compared seq.)
Length    Length                         Length
          E-value                        E-value
          Score                          Score

Somebody have any suggestions to store this information in above format.

Thank you

[Ram]

Duplicate of: BioPython : comparing same sequence with two different databases

Answer 1

[Bioinformatics_NewComer]

You can have an output in tabular format, that would simplify your work.

[sukeshssk]

Sir, I don't have output in tabular format but i wanted to convert in a tabular format......

[iraun]

You have an option in BLAST to have the output in tabular format. The parameter is -outfmt and you can set it to:

• 6 = tabular,
• 7 = tabular with comment lines,

In this way you would have the output in tabular format.

[Bioinformatics_NewComer]

A: parsing too large blast result with bioperl OR other methods?

This might be of help.