The bioBakery help forum

Metaphlan marker level analysis

Greetings to your nice forum.

I have finished running my analysis of metaphlan 3.0 on my datasets and trying to run all the marker analyses on the bowtie output. i used the code below

metaphlan -t marker_ab_table --bowtie2out metagenome_dataset1_.bowtie2.bz2 --nreads 9533350 --nproc 12 --input_type bowtie2out -o marker_abundance_table_dataset1.txt

but I don’t get any results… (the metagenomic profile has given a lot of taxonomies so I know I have a lot o f things as i have seen the txt output of metaphlan)

thanks

mm the number of reads was not correct that i provided. it is also stated at the end of the bowtie output

but this was not the issue… still doesnt give results

metaphlan -t marker_ab_table --bowtie2out metagenome_dataset1_.bowtie2.bz2 --nreads 9533350 --nproc 12 --input_type bowtie2out -o marker_abundance_table_dataset1.txt

There’s no need to provide nreads for marker_ab_table.

You forgot to provide the input file in the command line you posted here, this is the correct command:

metaphlan -t marker_ab_table --input_type bowtie2out -o marker_abundance_table_dataset1.txt metagenome_dataset1_.bowtie2.bz2

--bowtie2out refers to the output file in which the marker-to-reads BowTie2 mapping file is stored, you have to use this file as input file.

The command you run before may have overwritten the previous MetaPhlAn profile

thank you for your reply. Fortunately it didnt overwrite the original bowtie output. I would like to ask something else also. I got the relative abundances table but the ids of the rel abundance reads to the mapped ones are not giving me the organism name (just the ncbi id). How is it possible to get the actual names of the organisms so as to run in hclust2 ?
or the names are manually written afterwards by the user?

thanks again for the reply

You can add them manually by looking-up the species name using this mapping file attached below mpa_v30_CHOCOPhlAn_201901_markers_per_species.txt (1.9 MB)

thank you for the reply