# Different result from lefse at genus level when used metaphlan output profile upto genus and species level

**URL:** <https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807>\
**Category:** LEfSe\
**Created:** [August 2, 2020, 3:53pm UTC](https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807 "2020-08-02T15:53:56Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![DEEPCHANDA7](https://avatars.discourse-cdn.com/v4/letter/d/eb9ed0/32.png) [@DEEPCHANDA7](https://forum.biobakery.org/u/DEEPCHANDA7)\
**Post date:** [August 2, 2020, 3:53pm UTC](https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807/1 "2020-08-02T15:53:56Z")

</div>

Hi there!!!  
I have used same metaphlan output as input for lefse in two formats. In first case, extracted out all the OTUs upto species level like:

> type  
> clade\_name  
> k\_\_Bacteria|p\_\_Actinobacteria|c\_\_Actinobacteria|o\_\_Actinomycetales|f\_\_Actinomycetaceae|g\_\_Actinomyces|s\_\_Actinomyces\_odontolyticus

And, in second case, using only upto genus level from the same metaphlan profile output:

> ```
> type
> clade_name
> k __Bacteria|p__ Actinobacteria|c __Actinobacteria|o__ Actinomycetales|f __Actinomycetaceae|g__ Actinomyces
> 
> ```

But, I found both gives different results. I found some genus in the second case that were not found from the first case and vice-versa. Why am I getting this differences in output? Can anyone please help?

Thanks  
dc

---

<div class="post-metadata">

**Author:** ![sma](https://avatars.discourse-cdn.com/v4/letter/s/dec6dc/32.png) [@sma](https://forum.biobakery.org/u/sma)\
**Post date:** [August 7, 2020, 3:33pm UTC](https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807/2 "2020-08-07T15:33:15Z")

</div>

Hi -

Could you check if the result statistics are them same (log LDA score, KW p-values) between the two? If not I’d be worried - it might be miscalculations occurred when LEfSe collapses species level abundance in to genus. If they are the same, then there might be filtering considerations we could tweak to make the plots comparable.

Thanks,  
Siyuan

---

<div class="post-metadata">

**Author:** ![qianguo](https://avatars.discourse-cdn.com/v4/letter/q/839c29/32.png) [@qianguo](https://forum.biobakery.org/u/qianguo)\
**Post date:** [November 10, 2020, 6:27pm UTC](https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807/3 "2020-11-10T18:27:37Z")

</div>

Hi Siyuan,

```
    I met a problem when I using the web version of LEfSe. I got different results when I used two formats of taxonomy information, one is only the species level, and the other is 'genus|species' format. The former got more biomarkers than the latter, and the same biomarker appeared in the two practices got different LDA scores and KW p-values.

   I wonder have you figured out the reason? And which format will get more accurate results?

```

Thanks,  
Qian

---

<div class="post-metadata">

**Author:** ![leishenggit](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/leishenggit/32/814_2.png) [@leishenggit](https://forum.biobakery.org/u/leishenggit)\
**Post date:** [May 24, 2021, 10:21am UTC](https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807/4 "2021-05-24T10:21:46Z")

</div>

I also have this problem.

---

<div class="post-metadata">

**Author:** ![mishort](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/mishort/32/556_2.png) [@mishort](https://forum.biobakery.org/u/mishort)\
**Post date:** [June 4, 2021, 7:37pm UTC](https://forum.biobakery.org/t/different-result-from-lefse-at-genus-level-when-used-metaphlan-output-profile-upto-genus-and-species-level/807/5 "2021-06-04T19:37:26Z")

</div>

Hello,  
Do you happen to have a reproducible example you could send me? That would help me troubleshoot.  
Thanks,  
Meg