.hg_archival.txt

 repo: 092c2fe2278cb7f0b18d81faeb4aab98b89dc096
-node: 9760413b180fc2c68b817c23602541d3a97528af
-branch: default
-tag: 2.7.8
+node: 89633b311684ece67f32a9461aa1567e32dd42f7
+branch: 2.9
+tag: 2.9.20

.hgignore

@@ -6,3 +6,4 @@ tests/
 build/
 dist/
 *.egg-info/
+test/

.hgsubstate

-c168a100f37e23e2c110849a8d91fac8da49f5bd utils/export2graphlan
-69efddea43ae5e37d761cd138fcf083090371d1a utils/hclust2
+266347d7df0960de6e9771fdeb8899fefe9c80a7 utils/export2graphlan
+6d1023617944d9b7c69312bcb364886e7ef0a9f0 utils/hclust2

.hgtags

@@ -19,3 +19,20 @@ cbd7880df400b453b8beb4e62b39e4a23b5523b6 2.7.6
 4d50b9ccd95234a436541db13bcb10741f99138b 2.7.62
 a71e7d6b9d50b89c4b80714af145e10d050be7cf 2.7.63
 3e9c612e0a922b73214c90db51c5e26be17bf8b6 2.7.7
+9760413b180fc2c68b817c23602541d3a97528af 2.7.8
+28653dd67045f037e26828520e3695c0a0696e70 2.9.0
+0000000000000000000000000000000000000000 2.9.0
+f3840f32febc53e5fc2a815a7c6f2f1d81dfe57c 2.9.1
+0000000000000000000000000000000000000000 2.9.1
+a6017724af42d0dc5b27fae26a42a87704cc72c7 2.9.13
+a6017724af42d0dc5b27fae26a42a87704cc72c7 2.9.13
+559c0f70fe06eb0efd758641a26fbb8902debe89 2.9.13
+fae12d288dc0245549ec2f6db824af939717b5eb 2.9.14
+ac64c6e70db20ea5b4080b924872f29728100325 2.9.15
+ac64c6e70db20ea5b4080b924872f29728100325 2.9.15
+c20357c9133be435919db0b2948f9aa164f19a0b 2.9.15
+c20357c9133be435919db0b2948f9aa164f19a0b 2.9.15
+8f88d3957418b667bafe8aa6d48ec2746853c956 2.9.15
+71f40111c849408b60e1bc7aeddfef84c9e55eea 2.9.17
+170242bd646540bfb5521fe7ff2e54bc8a97dd35 2.9.18
+7030c379d395c9eb463ca48943373d52121140d9 2.9.19

__init__.py

-from metaphlan2 import metaphlan2
+from q2_metaphlan2 import metaphlan2
 from ._metaphlan2 import profile_single_fastq
 from ._metaphlan2 import profile_paired_fastq
 

_metaphlan2.py

@@ -23,6 +23,12 @@ def metaphlan2_helper(raw_data, nproc, input_type, output_file, verbose=True):
 
     sb.run(cmd, check=True)
 
+    print('\n\nIf you use MetaPhlAn2 in your work, please cite:\n\nTruong DT, '
+          'Franzosa EA, Tickle TL, Scholz M, Weingart G,\nPasolli E, Tett A, '
+          'Huttenhower C, Segata N.\nMetaPhlAn2 for enhanced metagenomic taxonomic '
+          'profiling\nNature Methods, 2015 Oct 1;12(10):902-3\n\nPMID: 26418763\n'
+          'doi: https://doi.org/10.1038/nmeth.3589', end='\n\n')
+
 
 def profile_single_fastq(raw_data: SingleLanePerSampleSingleEndFastqDirFmt,
                          nproc: int=1) -> biom.Table:

changeset.txt

+== Version 2.9
+- Automatic retrieval and installation of the latest MetaPhlAn2 database
+- New MetaPhlAn2 marker genes extracted with a newer version of ChocoPhlAn based on UniRef
+- Calculation of metagenome size for improved estimation of reads mapped to a given clade
+- Inclusion of NCBI taxonomy ID in the ouput file
+- CAMI (Taxonomic) Profiling Output Format included
 == Version 2.2.0
 - added option "marker_counts" (by Nicola)
 

debian/changelog

+metaphlan2 (2.9.20-1) UNRELEASED; urgency=medium
+
+  [ Andreas Tille ]
+  * New upstream version
+  * debhelper-compat 12
+  * Standards-Version: 4.4.0
+  * Fix watch file
+  * New upstream version
+
+  [ Steve Langasek ]
+  * Patches to port to Python3
+    Closes: #933661
+
+ -- Andreas Tille <tille@debian.org>  Fri, 16 Aug 2019 22:17:00 +0200
+
 metaphlan2 (2.7.8-1) unstable; urgency=medium
 
   * Drop libjs-twitter-bootstrap dependency which is not needed any more

debian/compat

-11

debian/control

@@ -3,24 +3,24 @@ Maintainer: Debian Med Packaging Team <debian-med-packaging@lists.alioth.debian.
 Uploaders: Andreas Tille <tille@debian.org>
 Section: science
 Priority: optional
-Build-Depends: debhelper (>= 11~),
-               python-all,
+Build-Depends: debhelper-compat (= 12),
+               python3-all,
                dh-python,
                pandoc,
                bowtie2
-Standards-Version: 4.2.1
+Standards-Version: 4.4.0
 Vcs-Browser: https://salsa.debian.org/med-team/metaphlan2
 Vcs-Git: https://salsa.debian.org/med-team/metaphlan2.git
 Homepage: https://bitbucket.org/biobakery/metaphlan2
 
 Package: metaphlan2
 Architecture: all
-Depends: ${python:Depends},
+Depends: ${python3:Depends},
          ${misc:Depends},
          metaphlan2-data,
-         python-biom-format,
-         python-msgpack,
-         python-pandas,
+         python3-biom-format,
+         python3-msgpack,
+         python3-pandas,
          bowtie2
 Description: Metagenomic Phylogenetic Analysis
  MetaPhlAn is a computational tool for profiling the composition of

debian/patches/_metaphlan2.py.patch

@@ -9,8 +9,10 @@ Description: Compile with Python2 despite _metaphlan2.py is using Python3 syntax
  support function annotations: These are optional in Python 3, and are
  removed from the function definitions in "_metaphlan2.py" by the patch.
 
---- a/_metaphlan2.py
-+++ b/_metaphlan2.py
+Index: metaphlan2-2.9.19/_metaphlan2.py
+===================================================================
+--- metaphlan2-2.9.19.orig/_metaphlan2.py
++++ metaphlan2-2.9.19/_metaphlan2.py
 @@ -3,7 +3,7 @@
  # This module defines the functions which run MetaPhlAn2 on
  # single and paired fastq data.
@@ -20,8 +22,8 @@ Description: Compile with Python2 despite _metaphlan2.py is using Python3 syntax
  import subprocess as sb
  from q2_types.per_sample_sequences import SingleLanePerSampleSingleEndFastqDirFmt
  from q2_types.per_sample_sequences import SingleLanePerSamplePairedEndFastqDirFmt
-@@ -24,8 +24,7 @@ def metaphlan2_helper(raw_data, nproc, i
-     sb.run(cmd, check=True)
+@@ -30,8 +30,7 @@ def metaphlan2_helper(raw_data, nproc, i
+           'doi: https://doi.org/10.1038/nmeth.3589', end='\n\n')
  
  
 -def profile_single_fastq(raw_data: SingleLanePerSampleSingleEndFastqDirFmt,
@@ -30,7 +32,7 @@ Description: Compile with Python2 despite _metaphlan2.py is using Python3 syntax
      output_biom = None
  
      with tempfile.TemporaryDirectory() as tmp_dir:
-@@ -36,8 +35,7 @@ def profile_single_fastq(raw_data: Singl
+@@ -42,8 +41,7 @@ def profile_single_fastq(raw_data: Singl
      return output_biom
  
  

debian/patches/mpa_dir-is-usr_share_metaphlan2.patch

@@ -7,27 +7,27 @@ Description: Instead of setting mpa_dir bash variable the path to the
 
 --- a/metaphlan2.py
 +++ b/metaphlan2.py
-@@ -417,7 +417,7 @@ def read_params(args):
+@@ -143,7 +143,7 @@ def read_params(args):
  
              "*  You can also provide an externally BowTie2-mapped SAM if you specify this format with \n"
              "   --input_type. Two steps: first apply BowTie2 and then feed MetaPhlAn2 with the obtained sam:\n"
--            "$ bowtie2 --sam-no-hd --sam-no-sq --no-unal --very-sensitive -S metagenome.sam -x ${mpa_dir}/db_v20/mpa_v20_m200 -U metagenome.fastq\n"
-+            "$ bowtie2 --sam-no-hd --sam-no-sq --no-unal --very-sensitive -S metagenome.sam -x /usr/share/metaphlan2/db_v20/mpa_v20_m200 -U metagenome.fastq\n"
-             "$ metaphlan2.py metagenome.sam --input_type sam > profiled_metagenome.txt\n\n"
+-            "$ bowtie2 --sam-no-hd --sam-no-sq --no-unal --very-sensitive -S metagenome.sam -x ${mpa_dir}/metaphlan_databases/mpa_v25_CHOCOPhlAn_201901 -U metagenome.fastq\n"
++            "$ bowtie2 --sam-no-hd --sam-no-sq --no-unal --very-sensitive -S metagenome.sam -x /usr/share/metaphlan2/metaphlan_databases/mpa_v25_CHOCOPhlAn_201901 -U metagenome.fastq\n"
+             "$ metaphlan2.py metagenome.sam --input_type sam -o profiled_metagenome.txt\n\n"
  
-             # "*  Multiple alternative ways to pass the input are also available:\n"
-@@ -1391,7 +1391,7 @@ def metaphlan2():
+             "*  We can also natively handle paired-end metagenomes, and, more generally, metagenomes stored in \n"
+@@ -1159,7 +1159,7 @@ def metaphlan2():
      # check for the mpa_pkl file
      if not os.path.isfile(pars['mpa_pkl']):
          sys.stderr.write("Error: Unable to find the mpa_pkl file at: " + pars['mpa_pkl'] +
 -                         "\nExpecting location ${mpa_dir}/db_v20/map_v20_m200.pkl "
-+                         "\nExpecting location /usr/share/metaphlan2/db_v20/mpa_v20_m200.pkl "
-                          "\nSelect the file location with the option --mpa_pkl.\n"
++                         "\nExpecting location /usr/share/metaphlan2/db_v20/map_v20_m200.pkl "
                           "Exiting...\n\n")
          sys.exit(1)
+ 
 --- a/README.md
 +++ b/README.md
-@@ -86,33 +86,27 @@ In case you moved the `metaphlan2.py` sc
+@@ -124,13 +124,7 @@ $ metaphlan2.py --install --index v29_CH
  
  This section presents some basic usages of MetaPhlAn2, for more advanced usages, please see at [its wiki](https://bitbucket.org/biobakery/biobakery/wiki/metaphlan2).
  
@@ -42,142 +42,7 @@ Description: Instead of setting mpa_dir bash variable the path to the
  
  Here is the basic example to profile a metagenome from raw reads (requires BowTie2 in the system path with execution and read permissions, Perl installed). 
  
- ```
- #!bash
--$ metaphlan2.py metagenome.fastq --input_type fastq > profiled_metagenome.txt
-+$ metaphlan2 metagenome.fastq --input_type fastq > profiled_metagenome.txt
- ```
- 
- It is highly recommended to save the intermediate BowTie2 output for re-running MetaPhlAn extremely quickly (`--bowtie2out`), and use multiple CPUs (`--nproc`) if available:
- 
- ```
- #!bash
--$ metaphlan2.py metagenome.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq > profiled_metagenome.txt
-+$ metaphlan2 metagenome.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq > profiled_metagenome.txt
- ```
- 
- If you already mapped your metagenome against the marker DB (using a previous  MetaPhlAn run), you can obtain the results in few seconds by using the previously saved `--bowtie2out` file and specifying the input (`--input_type bowtie2out`):
- 
- ```
- #!bash
--$ metaphlan2.py metagenome.bowtie2.bz2 --nproc 5 --input_type bowtie2out > profiled_metagenome.txt
-+$ metaphlan2 metagenome.bowtie2.bz2 --nproc 5 --input_type bowtie2out > profiled_metagenome.txt
- ```
- 
- You can also provide an externally BowTie2-mapped SAM if you specify this format with `--input_type`. Two steps here: first map your metagenome with BowTie2 and then feed MetaPhlAn2 with the obtained sam:
-@@ -120,14 +114,14 @@ You can also provide an externally BowTi
- ```
- #!bash
- $ bowtie2 --sam-no-hd --sam-no-sq --no-unal --very-sensitive -S metagenome.sam -x databases/mpa_v20_m200 -U metagenome.fastq
--$ metaphlan2.py metagenome.sam --input_type sam > profiled_metagenome.txt
-+$ metaphlan2 metagenome.sam --input_type sam > profiled_metagenome.txt
- ```
- 
- MetaPhlAn 2 can also natively **handle paired-end metagenomes** (but does not use the paired-end information), and, more generally, metagenomes stored in multiple files (but you need to specify the --bowtie2out parameter):
- 
- ```
- #!bash
--$ metaphlan2.py metagenome_1.fastq,metagenome_2.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq > profiled_metagenome.txt
-+$ metaphlan2 metagenome_1.fastq,metagenome_2.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq > profiled_metagenome.txt
- ```
- 
- For advanced options and other analysis types (such as strain tracking) please refer to the full command-line options.
-@@ -136,7 +130,7 @@ For advanced options and other analysis
- 
- 
- ```
--usage: metaphlan2.py --input_type
-+usage: metaphlan2 --input_type
-                      {fastq,fasta,multifasta,multifastq,bowtie2out,sam}
-                      [--mpa_pkl MPA_PKL] [--bowtie2db METAPHLAN_BOWTIE2_DB]
-                      [--bt2_ps BowTie2 presets] [--bowtie2_exe BOWTIE2_EXE]
-@@ -161,7 +155,7 @@ AUTHORS: Nicola Segata (nicola.segata@un
- 
- COMMON COMMANDS
- 
-- We assume here that metaphlan2.py is in the system path and that mpa_dir bash variable contains the
-+ We assume here that metaphlan2 is in the system path and that mpa_dir bash variable contains the
-  main MetaPhlAn folder. Also BowTie2 should be in the system path with execution and read
-  permissions, and Perl should be installed.
- 
-@@ -172,25 +166,25 @@ strains in particular cases) present in
- relative abundance. This correspond to the default analysis type (--analysis_type rel_ab).
- 
- *  Profiling a metagenome from raw reads:
--$ metaphlan2.py metagenome.fastq --input_type fastq
-+$ metaphlan2 metagenome.fastq --input_type fastq
- 
- *  You can take advantage of multiple CPUs and save the intermediate BowTie2 output for re-running
-    MetaPhlAn extremely quickly:
--$ metaphlan2.py metagenome.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq
-+$ metaphlan2 metagenome.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq
- 
- *  If you already mapped your metagenome against the marker DB (using a previous MetaPhlAn run), you
-    can obtain the results in few seconds by using the previously saved --bowtie2out file and 
-    specifying the input (--input_type bowtie2out):
--$ metaphlan2.py metagenome.bowtie2.bz2 --nproc 5 --input_type bowtie2out
-+$ metaphlan2 metagenome.bowtie2.bz2 --nproc 5 --input_type bowtie2out
- 
- *  You can also provide an externally BowTie2-mapped SAM if you specify this format with 
-    --input_type. Two steps: first apply BowTie2 and then feed MetaPhlAn2 with the obtained sam:
- $ bowtie2 --sam-no-hd --sam-no-sq --no-unal --very-sensitive -S metagenome.sam -x databases/mpa_v20_m200 -U metagenome.fastq
--$ metaphlan2.py metagenome.sam --input_type sam > profiled_metagenome.txt
-+$ metaphlan2 metagenome.sam --input_type sam > profiled_metagenome.txt
- 
- *  We can also natively handle paired-end metagenomes, and, more generally, metagenomes stored in 
-   multiple files (but you need to specify the --bowtie2out parameter):
--$ metaphlan2.py metagenome_1.fastq,metagenome_2.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq
-+$ metaphlan2 metagenome_1.fastq,metagenome_2.fastq --bowtie2out metagenome.bowtie2.bz2 --nproc 5 --input_type fastq
- 
- ------------------------------------------------------------------- 
-  
-@@ -208,23 +202,23 @@ file saved during the execution of the d
- *  The following command will output the abundance of each marker with a RPK (reads per kil-base) 
-    higher 0.0. (we are assuming that metagenome_outfmt.bz2 has been generated before as 
-    shown above).
--$ metaphlan2.py -t marker_ab_table metagenome_outfmt.bz2 --input_type bowtie2out > marker_abundance_table.txt
-+$ metaphlan2 -t marker_ab_table metagenome_outfmt.bz2 --input_type bowtie2out > marker_abundance_table.txt
-    The obtained RPK can be optionally normalized by the total number of reads in the metagenome 
-    to guarantee fair comparisons of abundances across samples. The number of reads in the metagenome
-    needs to be passed with the '--nreads' argument
- 
- *  The list of markers present in the sample can be obtained with '-t marker_pres_table'
--$ metaphlan2.py -t marker_pres_table metagenome_outfmt.bz2 --input_type bowtie2out > marker_abundance_table.txt
-+$ metaphlan2 -t marker_pres_table metagenome_outfmt.bz2 --input_type bowtie2out > marker_abundance_table.txt
-    The --pres_th argument (default 1.0) set the minimum RPK value to consider a marker present
- 
- *  The list '-t clade_profiles' analysis type reports the same information of '-t marker_ab_table'
-    but the markers are reported on a clade-by-clade basis.
--$ metaphlan2.py -t clade_profiles metagenome_outfmt.bz2 --input_type bowtie2out > marker_abundance_table.txt
-+$ metaphlan2 -t clade_profiles metagenome_outfmt.bz2 --input_type bowtie2out > marker_abundance_table.txt
- 
- *  Finally, to obtain all markers present for a specific clade and all its subclades, the 
-    '-t clade_specific_strain_tracker' should be used. For example, the following command
-    is reporting the presence/absence of the markers for the B. fragulis species and its strains
--$ metaphlan2.py -t clade_specific_strain_tracker --clade s__Bacteroides_fragilis metagenome_outfmt.bz2 databases/mpa_v20_m200.pkl --input_type bowtie2out > marker_abundance_table.txt
-+$ metaphlan2 -t clade_specific_strain_tracker --clade s__Bacteroides_fragilis metagenome_outfmt.bz2 databases/mpa_v20_m200.pkl --input_type bowtie2out > marker_abundance_table.txt
-    the optional argument --min_ab specifies the minimum clade abundance for reporting the markers
- 
- ------------------------------------------------------------------- 
-@@ -521,7 +515,7 @@ pickle.dump(db, ofile, pickle.HIGHEST_PR
- ofile.close()
- ```
- 
--* To use the new database, switch to metaphlan2/db_v21 instead of metaphlan2/db\_v20 when running metaphlan2.py with option "--mpa\_pkl".
-+* To use the new database, switch to metaphlan2/db_v21 instead of metaphlan2/db\_v20 when running metaphlan2 with option "--mpa\_pkl".
- 
- 
- ## Metagenomic strain-level population genomics
-@@ -591,7 +585,7 @@ for f in $(ls fastqs/*.bz2)
- do
-     echo "Running metaphlan2 on ${f}"
-     bn=$(basename ${f} | cut -d . -f 1)
--    tar xjfO ${f} | ../metaphlan2.py --bowtie2db ../databases/mpa_v20_m200 --mpa_pkl ../databases/mpa_v20_m200.pkl --input_type multifastq --nproc 10 -s sams/${bn}.sam.bz2 --bowtie2out sams/${bn}.bowtie2_out.bz2 -o sams/${bn}.profile
-+    tar xjfO ${f} | metaphlan2 --bowtie2db /usr/share/metaphlan2/db_v20/mpa_v20_m200 --mpa_pkl /usr/share/metaphlan2/db_v20/mpa_v20_m200.pkl --input_type multifastq --nproc 10 -s sams/${bn}.sam.bz2 --bowtie2out sams/${bn}.bowtie2_out.bz2 -o sams/${bn}.profile
- done
- ```
- 
-@@ -731,4 +725,4 @@ In the output folder, you can find the f
+@@ -752,4 +746,4 @@ In the output folder, you can find the f
  1. clade_name.fasta: the alignment file of all metagenomic strains.
  3. *.marker_pos: this file shows the starting position of each marker in the strains.
  3. *.info: this file shows the general information like the total length of the concatenated markers (full sequence length), number of used markers, etc.

debian/patches/python3.patch

+Description: set interpreter to python3
+Author: Steve Langasek <steve.langasek@ubuntu.com>
+Last-Modified: 2019-08-15
+Bug-Debian: https://bugs.debian.org/933661
+
+--- a/metaphlan2.py
++++ b/metaphlan2.py
+@@ -1,4 +1,4 @@
+-#!/usr/bin/env python
++#!/usr/bin/env python3
+ from __future__ import with_statement
+ __author__ = ('Nicola Segata (nicola.segata@unitn.it), '
+               'Duy Tin Truong, '
+--- a/strainphlan.py
++++ b/strainphlan.py
+@@ -1,4 +1,4 @@
+-#!/usr/bin/env python
++#!/usr/bin/env python3
+ # Author: Duy Tin Truong (duytin.truong@unitn.it)
+ #		at CIBIO, University of Trento, Italy
+ 
+--- a/utils/extract_markers.py
++++ b/utils/extract_markers.py
+@@ -1,4 +1,4 @@
+-#!/usr/bin/env python
++#!/usr/bin/env python3
+ #Author: Duy Tin Truong (duytin.truong@unitn.it)
+ #        at CIBIO, University of Trento, Italy
+ 
+--- a/utils/metaphlan2krona.py
++++ b/utils/metaphlan2krona.py
+@@ -1,4 +1,4 @@
+-#!/usr/bin/env python
++#!/usr/bin/env python3
+ 
+ # ============================================================================== 
+ # Conversion script: from MetaPhlAn output to Krona text input file
+--- a/utils/plot_bug.py
++++ b/utils/plot_bug.py
+@@ -1,4 +1,4 @@
+-#!/usr/bin/env python2
++#!/usr/bin/env python3
+ 
+ import sys
+ import numpy as np
+--- a/utils/read_fastx.py
++++ b/utils/read_fastx.py
+@@ -1,4 +1,4 @@
+-#!/usr/bin/env python
++#!/usr/bin/env python3
+ 
+ 
+ import sys

debian/patches/series

 mpa_dir-is-usr_share_metaphlan2.patch
 spelling.patch
 _metaphlan2.py.patch
+python3.patch

debian/patches/spelling.patch

@@ -2,29 +2,9 @@ Author: Andreas Tille <tille@debian.org>
 Last-Update: Mon, 23 May 2016 16:09:13 +0200
 Description: Spelling
 
---- a/README.md
-+++ b/README.md
-@@ -307,7 +307,7 @@ Post-mapping arguments:
- Additional analysis types and arguments:
-   -t ANALYSIS TYPE      Type of analysis to perform: 
-                          * rel_ab: profiling a metagenomes in terms of relative abundances
--                         * rel_ab_w_read_stats: profiling a metagenomes in terms of relative abundances and estimate the number of reads comming from each clade.
-+                         * rel_ab_w_read_stats: profiling a metagenomes in terms of relative abundances and estimate the number of reads coming from each clade.
-                          * reads_map: mapping from reads to clades (only reads hitting a marker)
-                          * clade_profiles: normalized marker counts for clades with at least a non-null marker
-                          * marker_ab_table: normalized marker counts (only when > 0.0 and normalized by metagenome size if --nreads is specified)
-@@ -713,7 +713,7 @@ python ../strainphlan.py -h
- The default setting can be stringent for some cases where you have very few samples left in the phylogenetic tree. You can relax some parameters to add more samples back:
- 
- 1. *marker\_in\_clade*: In each sample, the clades with the percentage of present markers less than this threshold are removed. Default "0.8". You can set this parameter to "0.5" to add some more samples.
--2. *sample\_in\_marker*: If the percentage of samples that a marker present in is less than this threhold, that marker is removed. Default "0.8". You can set this parameter to "0.5" to add some more samples.
-+2. *sample\_in\_marker*: If the percentage of samples that a marker present in is less than this threshold, that marker is removed. Default "0.8". You can set this parameter to "0.5" to add some more samples.
- 3. *N\_in\_marker*: The consensus markers with the percentage of N nucleotides greater than this threshold are removed. Default "0.2". You can set this parameter to "0.5" to add some more samples.
- 4. *gap\_in\_sample*: The samples with full sequences concatenated from all markers and having the percentage of gaps greater than this threshold will be removed. Default 0.2. You can set this parameter to "0.5" to add some more samples.
- 5. *relaxed\_parameters*: use this option to automatically set the above parameters to add some more samples by accepting some more gaps, Ns, etc. This option is equivalent to set: marker\_in\_clade=0.5, sample\_in\_marker=0.5, N\_in\_marker=0.5, gap\_in\_sample=0.5. Default "False".
 --- a/strainphlan.py
 +++ b/strainphlan.py
-@@ -337,7 +337,7 @@ def read_params():
+@@ -338,7 +338,7 @@ def read_params():
          required=False,
          default=['all'],
          type=str,
@@ -35,7 +15,7 @@ Description: Spelling
                  'file where each clade name is on a line will be read.'
 --- a/metaphlan2.py
 +++ b/metaphlan2.py
-@@ -596,7 +596,7 @@ def read_params(args):
+@@ -316,7 +316,7 @@ def read_params(args):
           default='rel_ab', help =
           "Type of analysis to perform: \n"
           " * rel_ab: profiling a metagenomes in terms of relative abundances\n"
@@ -44,3 +24,23 @@ Description: Spelling
           " * reads_map: mapping from reads to clades (only reads hitting a marker)\n"
           " * clade_profiles: normalized marker counts for clades with at least a non-null marker\n"
           " * marker_ab_table: normalized marker counts (only when > 0.0 and normalized by metagenome size if --nreads is specified)\n"
+--- a/README.md
++++ b/README.md
+@@ -376,7 +376,7 @@ Post-mapping arguments:
+ Additional analysis types and arguments:
+   -t ANALYSIS TYPE      Type of analysis to perform: 
+                          * rel_ab: profiling a metagenomes in terms of relative abundances
+-                         * rel_ab_w_read_stats: profiling a metagenomes in terms of relative abundances and estimate the number of reads comming from each clade.
++                         * rel_ab_w_read_stats: profiling a metagenomes in terms of relative abundances and estimate the number of reads coming from each clade.
+                          * reads_map: mapping from reads to clades (only reads hitting a marker)
+                          * clade_profiles: normalized marker counts for clades with at least a non-null marker
+                          * marker_ab_table: normalized marker counts (only when > 0.0 and normalized by metagenome size if --nreads is specified)
+@@ -734,7 +734,7 @@ python ../strainphlan.py -h
+ The default setting can be stringent for some cases where you have very few samples left in the phylogenetic tree. You can relax some parameters to add more samples back:
+ 
+ 1. *marker\_in\_clade*: In each sample, the clades with the percentage of present markers less than this threshold are removed. Default "0.8". You can set this parameter to "0.5" to add some more samples.
+-2. *sample\_in\_marker*: If the percentage of samples that a marker present in is less than this threhold, that marker is removed. Default "0.8". You can set this parameter to "0.5" to add some more samples.
++2. *sample\_in\_marker*: If the percentage of samples that a marker present in is less than this threshold, that marker is removed. Default "0.8". You can set this parameter to "0.5" to add some more samples.
+ 3. *N\_in\_marker*: The consensus markers with the percentage of N nucleotides greater than this threshold are removed. Default "0.2". You can set this parameter to "0.5" to add some more samples.
+ 4. *gap\_in\_sample*: The samples with full sequences concatenated from all markers and having the percentage of gaps greater than this threshold will be removed. Default 0.2. You can set this parameter to "0.5" to add some more samples.
+ 5. *relaxed\_parameters*: use this option to automatically set the above parameters to add some more samples by accepting some more gaps, Ns, etc. This option is equivalent to set: marker\_in\_clade=0.5, sample\_in\_marker=0.5, N\_in\_marker=0.5, gap\_in\_sample=0.5. Default "False".

debian/rules

@@ -3,7 +3,7 @@
 # DH_VERBOSE := 1
 
 %:
-	dh $@  --with python2
+	dh $@  --with python3
 
 override_dh_auto_build:
 	dh_auto_build

debian/watch

 version=4
 
 #opts="repacksuffix=+ds,dversionmangle=s/\+ds//g,repack,compression=xz" \
-opts="uversionmangle=s/(\d\.\d\.\d)(\d)/$1.$2/" \
+#opts="uversionmangle=s/(\d\.\d\.\d)(\d)/$1.$2/" \
 https://bitbucket.org/biobakery/metaphlan2/downloads?tab=tags .*/(\d[\d.]+)\.tar\.bz2