Skip to content

Commit

Permalink
Update 03_blastExample.md
Browse files Browse the repository at this point in the history
change year
  • Loading branch information
rrbrown98 authored Sep 11, 2023
1 parent 947e035 commit 8732983
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions 2023_edition/03_blastExample.md
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ Blast takes a fasta as an input, so we need to convert the fastq into fasta.

```bash
#recall your folder name
export me=$GROUP_SCRATCH/biochem_minicourse_2022/<your_dir>
export me=$GROUP_SCRATCH/biochem_minicourse_2023/<your_dir>

#move to the folder for the dataset
cd $me/data/external
Expand All @@ -24,7 +24,7 @@ head -n 20 SRR13403380_subset.fasta > SRR13403380.head.fasta

We have multipe databases installed on sherlock, which correspond to some of those available on ncbi blast online. Here, because we know what type of samples we are dealing with we have downloaded just a few databases and combined them into one blastdb. The path of this database is
```
/scratch/groups/astraigh/biochem_minicourse_2022/straightlab/blastdb/combined_db
/scratch/groups/astraigh/biochem_minicourse_2023/straightlab/blastdb/combined_db
```

Let's run blast using this database. We will put the blast output in a `blast` folder
Expand All @@ -34,7 +34,7 @@ Let's run blast using this database. We will put the blast output in a `blast` f
mkdir -p blast

#we can set a variable for this blast database so we don't have to type out the full path everytime
myblastdb=/scratch/groups/astraigh/biochem_minicourse_2022/straightlab/blastdb/combined_db
myblastdb=/scratch/groups/astraigh/biochem_minicourse_2023/straightlab/blastdb/combined_db

#run blast
blastn -query SRR13403380.head.fasta -db $myblastdb -num_threads 2 > blast/SRR13403380_rawoutput.first10.txt
Expand Down

0 comments on commit 8732983

Please sign in to comment.