BBtool's Fetchnt Error

[Zapages]

Hi Everyone,

I hope to test bbtool's sketch tool our bioinformatics needs. Before we proceed to doing any tests, we are having issues with retrieving the pre made/indexed NT database.

We ran the fetchNT.sh , which was included in bbtools package.


Unfortunately we are receiving the following error:

java -ea -Xmx1g -Xms1g -cp /opt/conda/opt/bbmap-38.86-0/current/ tax.RenameGiToTaxid -Xmx1g in=stdin.fa.gz out=renamed.fa.gz pigz=32 unpigz bgzip zl=8 server ow shrinknames maxbadheaders=5000 badheaders=badHeaders.txt taxpath=/BBSketch/taxonomy
Executing tax.RenameGiToTaxid [-Xmx1g, in=stdin.fa.gz, out=renamed.fa.gz, pigz=32, unpigz, bgzip, zl=8, server, ow, shrinknames, maxbadheaders=5000, badheaders=badHeaders.txt, taxpath=/BBSketch/taxonomy]

 

java.lang.Exception: Maximum bad headers exceeded: 5000
    at shared.KillSwitch.kill(KillSwitch.java:96)
    at tax.RenameGiToTaxid.updateHeadersFromServer_fasta(RenameGiToTaxid.java:629)
    at tax.RenameGiToTaxid.updateHeadersFromServer(RenameGiToTaxid.java:527)
    at tax.RenameGiToTaxid.processInner_server(RenameGiToTaxid.java:737)
    at tax.RenameGiToTaxid.process(RenameGiToTaxid.java:249)
    at tax.RenameGiToTaxid.main(RenameGiToTaxid.java:38)

Is there any way to fix this issue/error?

Thank you for the help.

[GenoMax]

NCBI's FTP site appears to be offline today. So that may be one reason why this is not working.

[Zapages]

Thank you, unfortunately we are still experiencing the issue...

[GenoMax]

I just tried the script and was able to get it to work (did not let it finish). It started making two files. renamed.fa.gz (which has fasta sequences) and badHeaders.txt.

Check to see if you are able to directly access NCBI FTP server from machine you are running this on.

[Zapages]

I just tried the script and was able to get it to work (did not let it finish). It started making two files. renamed.fa.gz (which has fasta sequences) and badHeaders.txt.

Check to see if you are able to directly access NCBI FTP server from machine you are running this on.

It works for me too, but it will eventually error out when the bad headers exceed more than 5000...

[GenoMax]

Create a ticket via BBMap's repo so Brian is aware of the problem.