forked from p2pu/course-in-a-box
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge branch 'gh-pages' of github.com:quadram-institute-bioscience/eb…
…ame-viromics into gh-pages
- Loading branch information
Showing
6 changed files
with
45 additions
and
103 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,5 @@ | ||
--- | ||
title: Linking to Anvio' | ||
--- | ||
|
||
You can do the Anvi'o data processing yourself on your own computer using the anvio install that you have used earlier in the week. We will provide the necessary files for you to download from [FigShare](https://doi.org/10.6084/m9.figshare.27231678) so you don't have to worry about downloading your files from the VM. |
This file was deleted.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,40 @@ | ||
--- | ||
title: Add Links to the Top Menu | ||
--- | ||
|
||
### Running the Anvi'o metagenomics workflow on the virome | ||
|
||
In the first step you will make the contigs database from the vOTU file and annotate it with information from the default hmms, the COGs database and identify tRNAs. | ||
|
||
```bash | ||
# CHANGE FILENAMES, PATHS AND THREADS AS APPROPRIATE FOR YOUR OWN COMPUTER | ||
anvi-gen-contigs-database -f votus.fna -o anvio/CONTIGS.db -T 4 | ||
anvi-run-hmms -c anvio/CONTIGS.d -T 4 | ||
anvi-run-ncbi-cogs -c anvio/CONTIGS.db -T 4 | ||
anvi-scan-trnas -c anvio/CONTIGS.db | ||
``` | ||
Check the anvio warnings. What do you see for the `anvi-run-hmms` output? | ||
|
||
<details> | ||
<summary>:green_book: Answer</summary> | ||
|
||
There are not that many hmms that return genes. But don't worry, that means that our virus predictions are of good quality, because we don't want the bacterial, archaeal and eukaryotic gene markers. | ||
</details> | ||
|
||
In the second step, you will create the anvio profiles from the read mapping files. We have provided subsampled files to limit download and computational time. | ||
|
||
```bash | ||
# CHANGE FILENAMES, PATHS AND THREADS AS APPROPRIATE FOR YOUR OWN COMPUTER | ||
anvi-profile -i sample_1.bam -c anvio/CONTIGS.db -T 4 | ||
anvi-profile -i sample_2.bam -c anvio/CONTIGS.db -T 4 | ||
anvi-profile -i sample_3.bam -c anvio/CONTIGS.db -T 4 | ||
``` | ||
Merge profiles to generate the profile database. | ||
```bash | ||
anvi-merge */PROFILE.db -o SAMPLES-MERGED -c CONTIGS.db | ||
```` | ||
Now you can look at your samples with `anvi-interactive`. | ||
|
||
## For the imatients | ||
|
||
You can use the `CONTIGS.db`, `PROFILE.db` and `AUXILIARY-DATA.db` as provided in the [FigShare](https://doi.org/10.6084/m9.figshare.27231678) folder! |
This file was deleted.
Oops, something went wrong.
This file was deleted.
Oops, something went wrong.
This file was deleted.
Oops, something went wrong.