Extract sequences from the blast results
After Performing Local Blast against NT database results generated into the *.blast files.
Now i would like to extract the sequences which are aligned and from my original sequences present into input *.fasta file.
I have tried the tool faSomeRecords.py which extracts the results. But the problem is to create the list of sequences
python faSomeRecords.py -f candidate.list -o candidate_sequences.fasta
In other way, if we can create the list of seqeunces from *.blast sequences. Then faSomeRecords.py will extracts the sequences for us.
• 1,800 views
•
link
updated
by
GenoMax
•
written
by
vmevada102 •
•
written
by
vmevada102 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Local blast search against nt db gains less results than online search
written by john
-
Blastn for gRNA characterization in ANopheles
written by Sbrillo •
-
Using 'blastn' to do a remote search on NCBI BLAST against the nt collection database
written by maciwuk •
-
How to perform bidirectional BLAST?
written by Paul
-
Local blast limit query search by GI list?
written by milt0n •
-
"OID not found" error with GI-list and preformatted nt BLAST database
written by ziv84 •
-
downloading nt databases - cannot extract them all?
written by balasink •
-
Why the local blast and online blast produce different results?
written by ngs101.com
-
Local Blast against a specific organism
written by Ritvik •
-
What'S The Easiest Way To Blast 5000 Sequences Against An Exon Fragment?
written by John
Have you tried the solutions suggested here: How to retrieve fasta sequence after local blast? ?
It would be best to have the blast results in columnar output format (like 6) to be able to easily extract the hit ID's. If you have results in html or some other format then parsing them would be more difficult.
Once you get that list using
blastdbcmdis the easiest option to retrieve sequences from a local blast database.