This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to separate a BLAST ouput file (tabular form) into files containing the hits for each protein searched.

Hi there,

I was wondering if you could help me. So I have done a multi-sequence blastp search which has generated an output.tsv file. I need to separate that tsv file into separate files containing the hits for each protein search. So for protein 1 have all the information of the hits into one file. And then protein 2 into a separate file. I tried to do this by limiting the target sequences to 10 and then splitting it by line number. So 10 hits go in each file but there are some proteins with 3 or 4 hits so then it messes up the separation. And I have to do this on Python !

I am in dire need of some help!

I know you can parse a blast out file but how would I then direct all the hits for each protein into a different file.

Any help would be really appreciated!! Thank you

blast python
• 1,932 views • link
written by Becca • Becca

Providing an outline for a non-python solution.

cut the first column out and then uniq that list to get sequence ID's that have hits. Then use grep in a loop with -w option to extract lines that contain that ID.

And I have to do this on Python !

Is this an assignment?

• 0 views • link
written by GenoMax GenoMax

No it's not an assignment. But I'm working with someone who only uses python. I could do it if I didn't have to use python but using python is confusing me slightly... Thank you though for your input.

• 0 views • link
written by Becca • Becca

Yeah maybe. I will have a look into that. I usually use pandas to make a dataframe to make plots and things.

Or I was thinking something like this : 

from blast import parse

fh = open('blast.tsv')
for blast_record in parse(fh):
    print('query id: {}'.format(blast_record.qid))
    for hit in blast_record.hits:
        for hsp in hit:
            print('****Alignment****')
            print('sequence:', hsp.sid)
            print('length:', hsp.length)
            print('e value:', hsp.evalue)

Output would look like: 

query id: cgl|CAGL0A00187g
****Alignment****
('sequence:', u'ecy|Ecym_8168')
('length:', 90)
('e value:', 0.0001)
****Alignment****
('sequence:', u'ecy|Ecym_8168')
('length:', 44)
('e value:', 0.0007)
****Alignment****
('sequence:', u'ecy|Ecym_4273')
('length:', 84)
('e value:', 0.64)

But I don't know how to all of that direct that into a file... I'll have a think...

• 0 views • link
updated by GenoMax GenoMax • written by Becca • Becca

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

Thank you!

• 0 views • link
written by GenoMax GenoMax

And to add to this, as soon as you have it in a DataFrame you could use something like the following loop (untested):

for query in df["query_id"].unique():
    df.loc[df"query_id"] == query].to_csv("blast_{}.txt".format(query), sep="\t", index=False)
• 0 views • link
written by WouterDeCoster WouterDeCoster

0 answers

No answers yet.

Log in to answer this question.

More posts like this