Running BLAST sequence similarity

BLAST is a program that compares a DNA/RNA or protein sequence, which are strings of letters against other reference sequences (i.e., whole non-redundant nucleic acid database, also known as NT).

Pre-requisites

• Have an HPC account on QUT’s lyra. Apply for a new HPC account here.

• Installed conda3 or miniconda3 ( https://docs.conda.io/projects/conda/en/latest/user-guide/install/linux.html ) → Recommend installing miniconda3 in your HPC’s /home directory.

• (optional) Basic unix command line knowledge (example: Learning Resources: the Linux Command Line ; https://swcarpentry.github.io/shell-novice/ )

• (optional) Familiarity with one unix text editors (example Vi/Vim or Nano):

  • VIM ( https://bioinformatics.uconn.edu/vim-guide/ ; https://missing.csail.mit.edu/2020/editors/)

  • Nano (https://engineering.purdue.edu/ECN/Support/KB/Docs/BasictutorialforNanou ; https://www.howtogeek.com/howto/42980/the-beginners-guide-to-nano-the-linux-command-line-text-editor/ )

Installing BLAST

Once miniconda3 has been installed, you will need to log out and login again to enable the conda command.

To install BLAST or other bioinformatics tools go to https://anaconda.org and search for the tool of interest. You can also use the command

to search for tools in the bioconda channel.

If the tool is available click on the tool link, which will open a new window showing the command line needed to install the tool. For example for BLAST the suggested command is:

run the above command to install blast. Conda will evaluate if the tool or necessary dependencies are available and will automatically install all necessary items to run in this case blast.

Note: Follow a similar process as above to install other tools.

Sample Data

Demo sample data to compare the similarity of DNA sequences generated by an RNA-seq approach against a reference Miscanthus sinensis mosaic virus (MsiMV) can be found at:

We want to compare the similarity (from 0 to 100%) of the sequences (also called ‘reads’) inside the query_sample.fa file against the reference MsiMV_genome.fasta sequence. Note: RNA/DNA (and protein) sequences can be stored in a ‘fasta format’. This is a header noted by “>” symbol followed by a sequence identifier on the first row. From the second row onwards the DNA/RNA(protein) sequence is presented.

Running blast on the HPC

We use a PBS Pro submission script to submit jobs to the HPC cluster. Create a file called ‘launch_blastN.pbs’ and fill it with this content, substituting email@host for your email address, and the files used as input to blastn:

Submit the job:

Checking the progression of the submitted job:

How to interpret the result? check this tutorial.