Skip to content
This repository was archived by the owner on Jul 20, 2021. It is now read-only.

Commit 2c1c492

Browse files
shiffer1Greg Caporaso
authored andcommitted
Update multiple-sequence-alignment.md (#310)
Changed the wording of questions 6 and 7 and added hints to questions 2, 3, and 4 for clarity on the metric.
1 parent be95bf3 commit 2c1c492

File tree

1 file changed

+7
-4
lines changed

1 file changed

+7
-4
lines changed

book/exercises/multiple-sequence-alignment.md

Lines changed: 7 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -150,7 +150,8 @@ Display the guide tree for the sequences in the cell below.
150150
```
151151
152152
```python
153-
>>> # Display the guide tree in this cell.
153+
>>> # Display the guide tree in this cell.
154+
... # hint: Pass metric=kmer_distance
154155
```
155156
156157
## Question 3 <link src="8xf4kF"/>
@@ -165,11 +166,13 @@ What are the differences in the guide tree from *Question 2*, the tree that is g
165166
>>> from skbio.alignment import global_pairwise_align_nucleotide
166167
>>> # add your command for 1 iterations of iterative multiple sequence alignment here
167168
... # hint: pass pairwise_aligner=global_pairwise_align_nucleotide
169+
... # hint: Pass metric=kmer_distance
168170
```
169171
170172
```python
171173
>>> # add your command for 5 iterations of iterative multiple sequence alignment here
172174
... # hint: pass pairwise_aligner=global_pairwise_align_nucleotide
175+
... # hint: Pass metric=kmer_distance
173176
```
174177
175178
## Question 4 <link src="sVsOoD"/>
@@ -183,6 +186,7 @@ Generate and display a tree based on progressive alignment of the sequences from
183186
```python
184187
>>> # Add your command for progressive alignment and tree building here
185188
... # hint: pass pairwise_aligner=global_pairwise_align_nucleotide
189+
... # hint: Pass metric=kmer_distance
186190
```
187191
188192
## Question 5 <link src="6b7fRx"/>
@@ -223,8 +227,7 @@ Discuss your results here.
223227
224228
## Question 6 <link src="UhSjc4"/>
225229
226-
Choose one representative sequence from each of the clusters you defined in question 5. Look these up in `tax_lookup` by their ids to get the taxonomy of each sequence, and include those in the results below. When you see a key that ends with ``__``, that means that there is no known taxonomic assignment for that sequence at that level.
227-
230+
Using `tax_lookup` by their ids, get the taxonomy of each sequence, and include those in the results below. When you see a key that ends with ``__``, that means that there is no known taxonomic assignment for that sequence at that level. Examine the taxonomy for each of your OTUs at 70% and 90%. Does the taxanomic assignment for these OTUs make sense? Additionally, choose one representative sequence from each of the clusters you defined in question 5. Is the taxonomy of the representative sequences consistent with phylogenetic tree you generated in question 4?
228231
```python
229232
>>> print(tax_lookup['4343117'])
230233
```
@@ -241,6 +244,6 @@ Discuss your results here.
241244
242245
## Question 7 <link src="5goFbb"/>
243246
244-
Is the taxonomy of the representative sequences consistent with phylogenetic tree you generated in question 4? For your 90% and 70% OTUs, list three taxa (e.g., at the phylum, class, or species level) that are monophyletic, if any, and three taxa that are not monophyletic, if any. Discuss two specific reasons why some taxa might appear to not be monophyletic based on your tree.
247+
For your 90% and 70% OTUs, list three taxa (e.g., at the phylum, class, or species level) that are monophyletic, if any, and three taxa that are not monophyletic, if any. Discuss two specific reasons why some taxa might appear to not be monophyletic based on your tree.
245248
246249
Discuss your results here.

0 commit comments

Comments
 (0)