diff --git a/src/1_JU28_59vs17_SNP_calling.sh b/src/1_JU28_59vs17_SNP_calling.sh index 61e11d4e0d78d2ad57087c737dffbb6115745cf3..44d2edce575c078bbe26e77d5ec67e07da253c55 100644 --- a/src/1_JU28_59vs17_SNP_calling.sh +++ b/src/1_JU28_59vs17_SNP_calling.sh @@ -29,6 +29,13 @@ src/intersect_SNP.R \ 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 +blastn -query data/RNA5S_belari.fasta -db data/fasta/DBG2OLC_output2.fasta -out results/blastall/RNA5S_2.out +less results/blastall/RNA5S_2.out +makeblastdb -in data/fasta/DBG2OLC_output1.fasta -parse_seqids -dbtype nucl +blastn -query data/RNA5S_belari.fasta -db data/fasta/DBG2OLC_output1.fasta -out results/blastall/RNA5S_1.out +less results/blastall/RNA5S_1.out + +makeblastdb -in data/fasta/nanoport_denovo.fasta -parse_seqids -dbtype nucl +blastn -query data/RNA5S_belari.fasta -db data/fasta/nanoport_denovo.fasta -out results/blastall/RNA5S_nanoport.out +less results/blastall/RNA5S_nanoport.out