BLAST

BLAST = Base L Aligment S T 
blast -query <




(base) $ cd ~/Ostei/cds.fna/Single_line 
### Activate in my case
(base) $ module load Bio/BLAST+/2.15.0 

### Start constructing the own DataBase
(base) $ cd <dir>/genes/prot.faa/
(base) $ makeblastdb -in protein.fasta -db 'prot'
### Now its work
(base) $ for file in [AS]*;
        do blastx -query $file -db ~/<directory>/prot.faa/<genes>.fasta -evalue 1e-15 -outfmt 6 -out <direcoty>/$file,genes.out;
        done

### After that, see the #13 collum (e-value) and separate just the 
(base) $ for file in Genes_of_*; 
         do awk '{print $1}' $file | uniq > gene_$file;
         done
### Check the list and count how much lines there is each file
(base) $ wc -l gene_*