Commit 1be5585f authored by Alessio Milanese's avatar Alessio Milanese
Browse files

Update README.md

parent 1d1e2cad
......@@ -413,3 +413,38 @@ We use [mOTUs2](https://github.com/motu-tool/mOTUs_v2) to create taxonomic profi
The two profiles are really similar. mOTUs is filtering the reads internally based on how good they map to the marker gene sequences; hence trimming and filtering the reads before will not affect much the profiles. For other analysis (like building metagenome-assembled genomes) trimming the reads improve the result.
</p>
6. (bonus) Download this metagenomic sample:
and profile with mOTUs3 and MetaPhlAn3 (already installed and can be run with `metaphlan`, manual: [link](https://github.com/biobakery/MetaPhlAn/wiki/MetaPhlAn-3.0#basic-usage)).
<details><summary>SOLUTION</summary>
<p>
We can profile the unfiltered reads with:
```
motus profile -f raw_reads_1.fastq -r raw_reads_2.fastq -t 8 -n unfiltered -o unfiltered.motus
```
Let's have a look at the species that are profiled by both:
```
motus merge -i unfiltered.motus,test_sample.motus | grep -v "0.0000000000" | grep -v "#" | head
```
Which results in:
```
Escherichia coli [ref_mOTU_v3_00095] 0.0002054152 0.0001797485
Eggerthella lenta [ref_mOTU_v3_00719] 0.0001531478 0.0001633324
Ruminococcus bromii [ref_mOTU_v3_00853] 0.0046285150 0.0049443856
Bacteroides uniformis [ref_mOTU_v3_00855] 0.0179185757 0.0188356677
Anaerostipes hadrus [ref_mOTU_v3_00857] 0.0000810850 0.0000864773
Roseburia faecis [ref_mOTU_v3_00859] 0.0003944823 0.0004207160
Roseburia inulinivorans [ref_mOTU_v3_00860] 0.0001580503 0.0001690509
Roseburia hominis [ref_mOTU_v3_00861] 0.0004800660 0.0005119912
Bifidobacterium longum [ref_mOTU_v3_01099] 0.0023339573 0.0027482230
Megasphaera elsdenii [ref_mOTU_v3_01516] 0.0006391750 0.0006816812
```
The two profiles are really similar. mOTUs is filtering the reads internally based on how good they map to the marker gene sequences; hence trimming and filtering the reads before will not affect much the profiles. For other analysis (like building metagenome-assembled genomes) trimming the reads improve the result.
</p>
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment