#!/bin/bash

set -e

blastn -db ../SEC7_yeast/blast_db/SEC7_yeast -query input.fasta -max_target_seqs 300 -outfmt "6 evalue bitscore pident sseqid slen sstart send sframe qstart qend qframe qseq sseq staxid" -out 1_blastn.result -evalue 0.01 -num_threads 5 -max_hsps 1 -task blastn -word_size 5
