NCBI BLAST:
Comparing Sequences

You have used BLAST to identify an unknown sequence of nucleotides.

Now let's compare sequences.

Click the Next button below.

Comparing Sequences with BLAST

1 of 7

For this example, we will again look at human mitochondrial DNA, but this time we will compare three different homo sapiens:

  • "Modern" human
  • Neanderthal
  • Denisova

We will use the reference sequences for the mitochondria for these three organisms and compare them using BLAST.

In the Nucleotide database, search for: mitochondrion[ti] AND human[orgn]

Then, because we are looking for the best quality sequences we can find, use the Source databases limit on the left of your screen to limit to RefSeq (curated) records.

Please use the Sort by Default order to follow along with the exercise.

Comparing Sequences with BLAST

2 of 7

The results are three RefSeq records for human mitochondria sequences:

  • One from modern humans (titled simply, "Homo sapiens mitochondrion")
  • One from Neanderthal (titled "Homo sapiens neanderthalensis")
  • One from Denisova (titled "Homo sp. Altai")

You can click on each record to learn more about each. When you are ready, return to the summary display of all three sequences, then use the link under Analyze these sequences in the right panel to Run BLAST.

This link feature from Nucleotide copies the record numbers into the BLAST query box.

Comparing Sequences with BLAST

3 of 7

You should now be on the Standard Nucleotide BLAST page and your three accession numbers should populate the Enter Query Sequence box.

Because we want to align these sequences, click the box: Align two or more sequences.

Then, move the accession numbers "NC_011137.1" and "NC_013993.1" to the new box.

This will compare the sequence in the first box (the modern human - NC_012920.1) with the other two (Neanderthal - NC_011137.1 and Denisova "Homo sp. Altai" - NC_013993.1).

Leave the other options at their default settings and click the BLAST button near the bottom of the screen.

Comparing Sequences with BLAST

4 of 7

Hopefully you won't be surprised that the results page shows two alignments against the modern human sequence: one for "Homo sapiens neanderthalensis mitochondrion, complete genome" and one for "Homo sp. Altai mitochondrion, complete genome."

What is interesting is in the Alignments view, which you can see by clicking on the Alignments tab.

If you scroll down or use the "Next" and "Previous" navigation, you will see that there are two tables:

  • The first comparing modern human with Neanderthal
  • The second comparing modern human to Denisova ("Homo sp. Altai")

Before going further, let’s change the display options so that you can see the differences between these sequences more clearly. Change the Alignment view options at the top to "Pairwise with dots for identities" and select the CDS feature:.

With the "Pairwise with dots for identities" display option, where there is a difference between the two sequences, the different base(s) are shown in the Subject line. Where there is no difference, a dot appears in the Subject line.

Comparing Sequences with BLAST

5 of 7

The CDS display option means that, directly from this screen, you can see the sequence changes that may result in amino acid changes.

Scroll down and find a section with CDS data:

The CDS data appears as two lines:

  • One above, showing the amino acid translation for the Query (modern human) sequence
  • One below, showing the amino acid translation for the Subject (ancient human) sequence

Comparing Sequences with BLAST

6 of 7

Note the section starting around base 3301 in the Query (modern human) sequence as compared to Neanderthal. this section shows CDS data.

Comparing the modern human mitochondrial sequence to the Neanderthal sequence, answer these questions:

A base difference around 3800 in the modern human sequence (an "A" versus a "G") results in a difference in amino acids. In modern humans, it is a T (threonine). In Neanderthal, it was a(n):

There is a base difference around 3308 in the Query (modern human) sequence that has an effect. What seems to be the effect?

Comparing Sequences with BLAST

7 of 7

Scroll down to the second table to view this same section comparing modern human to the Denisova ("Homo sp. Altai") sequence.

Looking at this same region around 3308 in the human (query) sequence, do modern humans differ from Denisova here, like the difference we saw with Neanderthal?

Conclusion

Congratulations! You have reached the end of this tutorial on BLAST.

You might have guessed, from all of the options on the screens, that there's a lot more to learn about BLAST. This exercise gave you a taste of a couple of things BLAST can do. To learn more about BLAST, check out the BLAST tutorials on the NLM YouTube channel.

Want more examples? See Teaching NCBI Resources Through Case Studies "Five Examples for NCBI BLAST."

You can now close the NLM Navigator windows.