# Understanding Parameters (stat\_q) for Environmental sample **URL:** \ **Category:** MetaPhlAn\ **Created:** [June 4, 2021, 3:44pm UTC](https://forum.biobakery.org/t/understanding-parameters-stat-q-for-environmental-sample/2204 "2021-06-04T15:44:18Z")\ **Posts on this page:** 3\ **Page:** 1
**Author:** ![liamfriar](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/liamfriar/32/763_2.png) [@liamfriar](https://forum.biobakery.org/u/liamfriar)\ **Post date:** [June 4, 2021, 3:44pm UTC](https://forum.biobakery.org/t/understanding-parameters-stat-q-for-environmental-sample/2204/1 "2021-06-04T15:44:18Z")
I used the default parameters for metaphlan3 to analyze trimmed metagenomic reads from a glacial sample. I do not get many hits (wc -l profiled\_metagenome.txt == 31) and they are dominated (97%) by cyanobacteria, which I do not think corresponds to the composition of the sample very well. I have seen mentioned in other threads to lower --stat\_q and --min\_mapq\_val for environmental samples. I understand --min\_mapq\_val , and I think I do not want to lower that much. I do not understand --stat\_q. Can anyone explain stat\_q to me? Does anyone have any thoughts on minimum values for these parameters for which the results would still be trustworthy? Or maybe the compositioin of my sample is just not well represented in the marker gene database? I do not think data quality is the problem because assembly and binning went very well. Finally, can I use my bowtie2 output from the default run as an input when I change parameters? Or do these parameters affect the creation of the bowtie2 output? Thank you! ---
**Author:** ![fbeghini](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/fbeghini/32/81_2.png) [@fbeghini](https://forum.biobakery.org/u/fbeghini)\ **Post date:** [June 29, 2021, 9:59am UTC](https://forum.biobakery.org/t/understanding-parameters-stat-q-for-environmental-sample/2204/3 "2021-06-29T09:59:14Z")
According to the quantile value chosen, the `stat_q` value, the markers taken into consideration for the relative abundance calculation are the one falling between the 20th and 80th percentile.The default value was chosen as 0.2 as a tradeoff in order to reduce the FP, decreasing it you’ll be able to detect more species but you will also increase the risk of having false positives. > [@liamfriar](#): > > Or maybe the compositioin of my sample is just not well represented in the marker gene database? It is a possibility, we have expanded a lot the genome catalog, but still environmental species are very scarce. > [@liamfriar](#): > > Finally, can I use my bowtie2 output from the default run as an input when I change parameters? Or do these parameters affect the creation of the bowtie2 output? Yes! You can run the first time with default parameters and then use the bowtie2out file as input when the parameters are changed. ---
**Author:** ![liamfriar](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/liamfriar/32/763_2.png) [@liamfriar](https://forum.biobakery.org/u/liamfriar)\ **Post date:** [June 29, 2021, 3:47pm UTC](https://forum.biobakery.org/t/understanding-parameters-stat-q-for-environmental-sample/2204/4 "2021-06-29T15:47:17Z")
Thank you @fbeghini !