diff --git a/src/1_JU28_59vs17_SNP_calling.sh b/src/1_JU28_59vs17_SNP_calling.sh
index 42f5adb34582d014ceb772bbf697cdc01d27d994..61e11d4e0d78d2ad57087c737dffbb6115745cf3 100644
--- a/src/1_JU28_59vs17_SNP_calling.sh
+++ b/src/1_JU28_59vs17_SNP_calling.sh
@@ -27,3 +27,8 @@ src/intersect_SNP.R \
   data/list_of_enzymes.csv
 ~/scripts/sms.sh "SNP analysis done"
 
+mkdir -p results/blastall/
+makeblastdb -in data/fasta/DBG2OLC_output2.fasta -parse_seqids -dbtype nucl
+blastn -query data/RNA5S_belari.fasta -db data/fasta/DBG2OLC_output2.fasta -out results/blastall/RNA5S.out
+less results/blastall/RNA5S.out
+