+ mkdir -p db
+ esearch -db protein -query PRJNA257197
+ efetch -format fasta
+ esearch -db nuccore -query PRJNA257197
+ efetch -format fasta
+ makeblastdb -in db/proteins.fa -dbtype prot -out db/ebov_prots -parse_seqids
+ makeblastdb -in db/genomes.fa -dbtype nucl -out db/ebov_genomes -parse_seqids
+ blastdbcmd -db db/ebov_prots -info
+ echo '>mystery'
+ echo ATGGACTCTCGTCCTCAGAAAGTCTGGATGACGCCGAGTCTCACTGAATCTGACATGGAT
+ echo TACCACAAGATCTTGACAGCAGGTCTGTCCGTTCAACAGGGGGTTGTTCGGCAAAGAGTC
+ echo ATCCCAGTGTATCAAGTAAACAATCTTGAG
+ F='6 qacc sacc pident qlen length frames'
+ blastn -db db/ebov_genomes -query foo.fa -outfmt '6 qacc sacc pident qlen length frames'
+ sort -k 3rn
+ blastx -db db/ebov_prots -query foo.fa -outfmt '6 qacc sacc pident qlen length frames'
+ sort -k 3rn
+ tblastx -db db/ebov_genomes -query foo.fa -outfmt '6 qacc sacc pident qlen length frames'
+ sort -k 3rn