I want to narrow down the blastn search against
nt database using
I already got all taxids of bacteria (
taxid 2) and extacted GIs with
csvtk (Please refer to this).
The next step was to proceed bacterial species identification.
When I run
blastn -query query.fasta -db /path/to/nt -gilist bacteria.taxid.gi.txt -evalue 1e-6 -outfmt 6 -out sequences.txt
An error occured:
BLAST Database error: Specified file is not a valid GI/TI list.
Please refer to the attached file.
bacteria.taxid.gi.txt (Number of taxids: 309,264,110)
What am I doing wrong? Thanks for the help in advance.