---
title: "09-unmap-blast"
output: html_document
date: "2024-11-08"
---
```{bash}
cd ../data/blastdb/
for i in {00..22}; do
    wget ftp://ftp.ncbi.nlm.nih.gov/blast/db/nt_prok.${i}.tar.gz
done


```

```{bash}
cd ../data/blastdb/
for i in {00..22}; do
    tar -xvzf nt_prok.${i}.tar.gz
done
```


```{bash}
fasta="../output/08-trinity-unmapped/Trinity.fasta"

blastn \
-query $fasta \
-db ../data/blastdb/nt_prok \
-out ../output/09-unmap-blast/unmap_blastn_ntprok.tab \
-task blastn \
-evalue 1E-20 \
-num_threads 30 \
-max_target_seqs 1 \
-max_hsps 1 \
-outfmt 6 \
> ../output/09-unmap-blast/unmap_blastn_ntprok.log 2>&1

echo "done"
```


