How to make local BLAST sequentially look through a large number of databases and generate a file for each
0
0
Entering edit mode
13 months ago
nipayrhk • 0

I want to run a tblastn of a single protein sequence against a large number of genomes with each as a separate database so that each one generates a new file even if there is no match. This is important as I'm equally interested in genomes with it as I am for genomes without it. If there is a way to make one database create files for each genome this would work. The main point is just that I need a result file for each genome. The amount of genomes will be in the thousands so doing it by hand would be rather tedious. The genomes I have downloaded from from the ncbi database and have converted to individual databases with a script. If the individual file for genome isn't a feature of BLAST I will just write a script to BLAST the genomes one after the other.

BLAST • 588 views
ADD COMMENT
1
Entering edit mode

so that each one generates a new file even if there is no match

and

If the individual file for genome isn't a feature of BLAST I will just write a script to BLAST the genomes one after the other.

You will have to do that yourself, especially since you want an empty file if there is no match. You scripted creating the databases already so this task should follow along the same lines.

ADD REPLY

Login before adding your answer.

Traffic: 3319 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6