@@ -25,7 +25,7 @@ You can run the entire pipeline with **a single command**, or use the scripts **
2525
2626### Total abundance output
2727
28- ` counts_phylum.csv ` parsed from 7 kraken2 reports of metagenomic samples using ` KrakenParser ` :
28+ ` counts_phylum.csv ` parsed from 9 kraken2 reports of metagenomic samples using ` KrakenParser ` :
2929
3030```
3131Sample_id,Calditrichota,Caldisericota,Thermosulfidibacterota,Elusimicrobiota,Candidatus Fervidibacterota,Lentisphaerota,Kiritimatiellota,Vulcanimicrobiota,Thermodesulfobiota,Atribacterota,Dictyoglomota,Nitrospinota,Chrysiogenota,Coprothermobacterota,Aquificota,Thermotogota,Bdellovibrionota,Nitrospirota,Deferribacterota,Synergistota,Myxococcota,Acidobacteriota,Candidatus Bipolaricaulota,Candidatus Saccharibacteria,Candidatus Absconditabacteria,Fusobacteriota,Spirochaetota,Candidatus Omnitrophota,Chlamydiota,Verrucomicrobiota,Planctomycetota,Thermodesulfobacteriota,Campylobacterota,Candidatus Cloacimonadota,Fibrobacterota,Gemmatimonadota,Balneolota,Rhodothermota,Ignavibacteriota,Chlorobiota,Bacteroidota,Deinococcota,Thermomicrobiota,Armatimonadota,Chloroflexota,Cyanobacteriota,Mycoplasmatota,Actinomycetota,Bacillota,Pseudomonadota,Heterolobosea,Parabasalia,Fornicata,Evosea,Bacillariophyta,Cercozoa,Euglenozoa,Apicomplexa,Microsporidia,Basidiomycota,Ascomycota,Nanoarchaeota,Candidatus Micrarchaeota,Candidatus Thermoplasmatota,Candidatus Lokiarchaeota,Nitrososphaerota,Euryarchaeota,Thermoproteota,Hofneiviricota,Artverviricota,Nucleocytoviricota,Cossaviricota,Kitrinoviricota,Negarnaviricota,Lenarviricota,Pisuviricota,Peploviricota,Uroviricota
@@ -39,7 +39,7 @@ X9,0,3,2,16,7,1,23,12,10,9,1,2,134,40,390,289,29,372,27,81,150,90,9,88,32,287,88
3939
4040### Relative abundance output
4141
42- ` ra_phylum.csv ` calculated from 7 kraken2 reports of metagenomic samples using ` KrakenParser ` :
42+ ` ra_phylum.csv ` calculated from 9 kraken2 reports of metagenomic samples using ` KrakenParser ` :
4343
4444```
4545Sample_id,taxon,rel_abund_perc
@@ -61,7 +61,7 @@ X9,Other (<4.0%),1.8979510606688494
6161
6262### α-diversity output
6363
64- ` alpha_div.csv ` calculated from 7 kraken2 reports of metagenomic samples using ` KrakenParser ` :
64+ ` alpha_div.csv ` calculated from 9 kraken2 reports of metagenomic samples using ` KrakenParser ` :
6565
6666```
6767Sample,Shannon,Pielou,Chao1
@@ -74,7 +74,7 @@ X9,4.033664950188261,0.5050385978575492,3492.16
7474
7575### β-diversity output
7676
77- ` beta_div_bray.csv ` calculated from 7 kraken2 reports of metagenomic samples using ` KrakenParser ` :
77+ ` beta_div_bray.csv ` calculated from 9 kraken2 reports of metagenomic samples using ` KrakenParser ` :
7878
7979```
8080,X1,X2,...,X8,X9
@@ -85,7 +85,7 @@ X8,0.61,0.723,...,0.0,0.665
8585X9,0.353,0.388,...,0.665,0.0
8686```
8787
88- ` beta_div_jaccard.csv ` calculated from 7 kraken2 reports of metagenomic samples using ` KrakenParser ` :
88+ ` beta_div_jaccard.csv ` calculated from 9 kraken2 reports of metagenomic samples using ` KrakenParser ` :
8989
9090```
9191,X1,X2,...,X8,X9
@@ -278,8 +278,7 @@ KrakenParser --diversity -i data/counts/counts_species.csv -o data/diversity -s
278278## Example Output Structure
279279After running the full pipeline, the output directory will look like this:
280280```
281- data/
282- ├─ kreports/ # Input Kraken2 reports
281+ results/
283282├─ counts/ # Total abundance CSV output
284283│ ├─ counts_species.csv
285284│ ├─ counts_genus.csv
@@ -296,7 +295,7 @@ data/
296295│ └─ beta_div_jaccard.csv
297296└─ intermediate/ # Intermediate files
298297 ├─ mpa/ # Converted MPA files
299- │ ├─ sample.MPA.TXT
298+ │ ├─ { sample}.txt
300299 │ ├─ ...
301300 ├─ COMBINED.txt # Merged MPA table
302301 └─ txt/ # Extracted taxonomic levels in TXT
0 commit comments