Hi, I am using humann2 v0.11.2 and MetaPhlAn2 version 2.7.7 for metagenomic data, after running, I found that certain species found from humann2 was not in those found by metaphlan2.
For example, Pseudomonas protegens Pf-5 found from humann2 was absent in the results from metaphlan2. I blast the fasta reads against the legacy kegg database first, then the result was provided as the input of humann2. By the way, I did not provide the --taxonomic-profile
parameter while running humann2. Here is an example of my code.
diamond blastx -p 12 -q A57RS1T1.fa -d /public/home/jkyin/biosoft/database/humann-0.99/data/kegg.reduced.dmnd -o /public/home/jkyin/yjk/tomato_micro/meta/result/humann2/A57RS1T1.txt
humann2 --input A57RS1T1.txt --output /public/home/jkyin/yjk/tomato_micro/meta/result/humann2/func \ --id-mapping /public/home/jkyin/biosoft/database/humann-0.99/data/legacy_kegg_idmapping.tsv \ --pathways-database /public/home/jkyin/biosoft/database/humann-0.99/data/keggc --threads 32 \ --remove-temp-output
Is it necessary to add --taxonomic-profile with metaphlan2 output?