You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: _episodes/1.5.2_Viral_Taxonomy_and_Phylogeny_II.md
+18-19Lines changed: 18 additions & 19 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -208,7 +208,7 @@ We've built two python scripts using this library to shrink the size of tree out
208
208
209
209
> ## Exercise - Post-process your tree
210
210
>
211
-
> Use the below sbatch script make your tree and post-process it. You will have change the contig name in the sbatch script. **Note:** See exercise question #10 below - you might want to plot this 'pet contig' here, or the one from yesterday!
211
+
> Use the below sbatch script make your tree and post-process it. You will have change the contig name in the sbatch script. **Note:** See exercise question #4 below - you might want to plot this 'pet contig' here, or the one from yesterday!
212
212
>
213
213
> Please pause once you have made the trees! We will visualize the tree together with a demo in itol **so that you can make nice graphics to include in your lab books**!
214
214
>
@@ -262,14 +262,13 @@ We've built two python scripts using this library to shrink the size of tree out
262
262
>{: .source}
263
263
{: .solution}
264
264
265
-
### vConTACT3
265
+
## vConTACT3
266
266
267
267
vConTACT3 has an underlying assumption that the fraction of shared genes between two viruses represents their evolutionary relationship. The vConTACT3 gene-sharing network closely correlates with the ICTV taxonomy.
268
268
269
269
>## Exercise - run vConTACT3
270
270
>
271
-
> {:start='6'}
272
-
> 6. Run vcontact3 on your virome dataset
271
+
> Run vcontact3 on your virome dataset
273
272
>
274
273
> You will need the following commands for the conda environment and database
275
274
>
@@ -301,7 +300,6 @@ vConTACT3 has an underlying assumption that the fraction of shared genes between
301
300
>{: .source}
302
301
{: .challenge}
303
302
304
-
305
303
>## sbatch script to run vcontact3
306
304
>```bash
307
305
>#!/bin/bash
@@ -341,25 +339,26 @@ vConTACT3 has an underlying assumption that the fraction of shared genes between
341
339
>
342
340
>### General Questions
343
341
>
344
-
>{:start='7'}
345
-
>7. What percentage of your total viral contigs were classified by vConTACT3?
346
-
>8. How can you judge the accuracy of your taxonomic classifications?
347
-
>9. What is the maximum taxonomic distance between viruses that are still connected in the vConTACT3 graph?
342
+
>{:start='1'}
343
+
>1. What percentage of your total viral contigs were classified by vConTACT3?
344
+
>2. How can you judge the accuracy of your taxonomic classifications?
345
+
>3. What is the maximum taxonomic distance between viruses that are still connected in the vConTACT3 graph?
348
346
>
349
347
>**Pick one pet viral contig that has a terL gene and also a taxonomic classification by vConTACT3**
350
348
>
351
-
>{:start='10'}
352
-
>10. What is the taxonomy classifcation of your *pet contig* from vConTACT3 and the terL tree?
353
-
>11. How many genes are shared with known reference viruses?*(vConTACT3 results)*
354
-
>12. If your classifications are/were different between the terL tree and the vConTACT3 - which classification would you trust more?
349
+
>{:start='4'}
350
+
>4. What is the taxonomy classifcation of your *pet contig* from vConTACT3 and the terL tree?
351
+
>5. How many genes are shared with known reference viruses?*(vConTACT3 results)*
352
+
>6. If your classifications are/were different between the terL tree and the vConTACT3 - which classification would you trust more?
355
353
>
356
-
>### Optional/Group Discussions
354
+
>### Group Discussions
355
+
>
356
+
>{:start='7'}
357
+
>7. Explore the cytoscape network visualization and/or terL tree together - include a graphic of at least 1 tree
358
+
>8. What do long branches on your terL tree mean?
359
+
>9. What does a rooted versus an unrooted tree represent?
360
+
>10. What are the limitations of using the terL tree? and could you overcome these limitations?
357
361
>
358
-
>{:start='12'}
359
-
>13. Explore the cytoscape network visualization and/or terL tree together
360
-
>14. What do long branches on your terL tree mean?
361
-
>15. Where are the contigs assigned to long branches in the network from vConTACT3?
362
-
>16. What are the limitations of using the terL tree? and could you overcome these limitations?
0 commit comments