I'm a computer scientist working with biologists at a small school that doesn't have dedicated bioinformatics staff. I apologize if I use incorrect terminology since I have limited bioinformatics background.
Our biology staff is studying a soil sample that has been sequenced. The resulting sequences (in the hundreds of millions) were analyzed using mg-rast, but they didn't seem to have a high confidence in the results.
I've been asked to help analyze these sequences. I setup a local instance of NCBI blast, but it was far too slow to analyze that many sequences. We moved on to using BLAT. They've requested we use the nt database from NCBI. I converted the nt FASTA file to 2bit file, splitting the files into small chunks.
gfServer seems to have limitations however. When I run gfServer, I get the following error while gfServer attempts to start.
Exceeding 4 billion bases, sorry gfServer can't handle that.
Is it possible to use the full nt database using BLAT? Is there an alternative to BLAT that can handle that many sequences?
nt
database? If you're analyzing a soil sample, you already know you don't care about any mammalian sequences, for instance. You can make the search much faster if you choose a subset of nt instead of the entire thing. Especially since your biologists are probably only interested in unicellular species. $\endgroup$