number=0
makeblastdb -in "$2" -dbtype nucl -out db -logfile /dev/null
for ((i = 0; i < 100; i++))
do 
    makenucseq -amount 1 -length "$1" -auto
    blastn -task blastn -query makeseq.fasta -db db -out out.blast -evalue 0.1 -word_size 4 -outfmt 6
    x=`wc -l out.blast | cut -f1 -d' '`
    number=$(($number+$x))
    rm out.blast
    rm makeseq.fasta
done
number=`echo "scale = 2; $number/100.0" | bc`
echo The average number of alignmenst with E-value less than 0.1: $number
rm db.*