mgnify-pipelines-toolkit 0.1.1__py3-none-any.whl → 0.1.3__py3-none-any.whl

mgnify_pipelines_toolkit/analysis/amplicon/amplicon_utils.py

@@ -46,18 +46,11 @@ def get_read_count(read_path, type='fastq'):
         ]
         zcat_proc = subprocess.Popen(cmd, stdout=subprocess.PIPE, stderr=subprocess.PIPE)
 
-        cmd = [
-            'sed',
-            '-n',
-            '1~4p',
-        ]
-        sed_proc = subprocess.Popen(cmd, stdin=zcat_proc.stdout, stdout=subprocess.PIPE, stderr=subprocess.PIPE)
-
         cmd = [
             'wc',
             '-l'
         ]
-        wc_proc = subprocess.Popen(cmd, stdin=sed_proc.stdout, stdout=subprocess.PIPE, stderr=subprocess.PIPE)
+        wc_proc = subprocess.Popen(cmd, stdin=zcat_proc.stdout, stdout=subprocess.PIPE, stderr=subprocess.PIPE)
         stdout, stderr = wc_proc.communicate()
 
     elif type == 'fasta':
@@ -78,6 +71,9 @@ def get_read_count(read_path, type='fastq'):
 
     read_count = int(read_count)
 
+    if type == 'fastq':
+        read_count /= 4
+
     return read_count
 
 def build_cons_seq(cons_list, read_count, cons_threshold=0.80, do_not_include=None, counter=1, max_line_count=None):

mgnify_pipelines_toolkit/analysis/amplicon/make_asv_count_table.py

@@ -25,13 +25,14 @@ from mgnify_pipelines_toolkit.constants.tax_ranks import _SILVA_TAX_RANKS, _PR2_
 logging.basicConfig(level=logging.DEBUG)
 
 def parse_args():
-
     parser = argparse.ArgumentParser()
 
-    parser.add_argument("-t", "--taxa", required=True, type=str, help="Path to DADA2 taxa file")
+    parser.add_argument("-t", "--taxa", required=True, type=str, help="Path to taxa file")
     parser.add_argument("-f", "--fwd", required=True, type=str, help="Path to DADA2 forward map file")
     parser.add_argument("-r", "--rev", required=False, type=str, help="Path to DADA2 reverse map file")
-    parser.add_argument("-a", "--amp", required=True, type=str, help="Path to extracted amp_region reads from inference subworkflow")
+    parser.add_argument(
+        "-a", "--amp", required=True, type=str, help="Path to extracted amp_region reads from inference subworkflow"
+    )
     parser.add_argument("-hd", "--headers", required=True, type=str, help="Path to fastq headers")
     parser.add_argument("-s", "--sample", required=True, type=str, help="Sample ID")
 
@@ -48,8 +49,7 @@ def parse_args():
 
 
 def order_df(taxa_df):
-
-    if len(taxa_df.columns) == 8:
+    if len(taxa_df.columns) == 9:
         taxa_df = taxa_df.sort_values(_SILVA_TAX_RANKS, ascending=True)
     elif len(taxa_df.columns) == 10:
         taxa_df = taxa_df.sort_values(_PR2_TAX_RANKS, ascending=True)
@@ -60,70 +60,76 @@ def order_df(taxa_df):
     return taxa_df
 
 def make_tax_assignment_dict_silva(taxa_df, asv_dict):
-
     tax_assignment_dict = defaultdict(int)
 
     for i in range(len(taxa_df)):
-        
         sorted_index = taxa_df.index[i]
-        asv_count = asv_dict[sorted_index]
+        asv_num = taxa_df.iloc[i, 0]
+        asv_count = asv_dict[asv_num]
 
         if asv_count == 0:
             continue
 
+        sk = taxa_df.loc[sorted_index, "Superkingdom"]
         k = taxa_df.loc[sorted_index, "Kingdom"]
         p = taxa_df.loc[sorted_index, "Phylum"]
         c = taxa_df.loc[sorted_index, "Class"]
         o = taxa_df.loc[sorted_index, "Order"]
         f = taxa_df.loc[sorted_index, "Family"]
-        g = taxa_df.loc[sorted_index, "Genus"]        
+        g = taxa_df.loc[sorted_index, "Genus"]
         s = taxa_df.loc[sorted_index, "Species"]
 
         tax_assignment = ""
 
         while True:
 
+            if sk != "0":
+                sk = "_".join(sk.split(" "))
+                tax_assignment += sk
+            else:
+                break
+
             if k != "0":
                 k = "_".join(k.split(" "))
-                if k == "Archaea" or k == "Bacteria":
-                    tax_assignment += f"sk__{k}"
-                elif k == "Eukaryota":
-                    tax_assignment += f"sk__Eukaryota"
-                else:
-                    tax_assignment += f"sk__Eukaryota\tk__{k}"
+                tax_assignment += f"\t{k}"
+            elif sk != "0":
+                tax_assignment += f"\tk__"
             else:
                 break
 
             if p != "0":
-                if k == "Archaea" or k == "Bacteria":
-                    tax_assignment += f"\tk__"
                 p = "_".join(p.split(" "))
-                tax_assignment += f"\tp__{p}"
+                tax_assignment += f"\t{p}"
             else:
                 break
+
             if c != "0":
                 c = "_".join(c.split(" "))
-                tax_assignment += f"\tc__{c}"
+                tax_assignment += f"\t{c}"
             else:
                 break
+
             if o != "0":
                 o = "_".join(o.split(" "))
-                tax_assignment += f"\to__{o}"
+                tax_assignment += f"\t{o}"
             else:
                 break
+
             if f != "0":
                 f = "_".join(f.split(" "))
-                tax_assignment += f"\tf__{f}"
+                tax_assignment += f"\t{f}"
             else:
                 break
+
             if g != "0":
                 g = "_".join(g.split(" "))
-                tax_assignment += f"\tg__{g}"
+                tax_assignment += f"\t{g}"
             else:
                 break
+
             if s != "0":
                 s = "_".join(s.split(" "))
-                tax_assignment += f"\ts__{s}"
+                tax_assignment += f"\t{s}"
             break
 
         if tax_assignment == "":
@@ -134,13 +140,12 @@ def make_tax_assignment_dict_silva(taxa_df, asv_dict):
     return tax_assignment_dict
 
 def make_tax_assignment_dict_pr2(taxa_df, asv_dict):
-
     tax_assignment_dict = defaultdict(int)
 
     for i in range(len(taxa_df)):
-        
         sorted_index = taxa_df.index[i]
-        asv_count = asv_dict[sorted_index]
+        asv_num = taxa_df.iloc[i, 0]
+        asv_count = asv_dict[asv_num]
 
         if asv_count == 0:
             continue
@@ -153,53 +158,62 @@ def make_tax_assignment_dict_pr2(taxa_df, asv_dict):
         o = taxa_df.loc[sorted_index, "Order"]
         f = taxa_df.loc[sorted_index, "Family"]
         g = taxa_df.loc[sorted_index, "Genus"]
-        s = taxa_df.loc[sorted_index, "Species"]         
+        s = taxa_df.loc[sorted_index, "Species"]
 
         tax_assignment = ""
 
         while True:
-
             if d != "0":
                 d = "_".join(d.split(" "))
-                tax_assignment += f"d__{d}"
+                tax_assignment += d
             else:
                 break
 
             if sg != "0":
                 sg = "_".join(sg.split(" "))
-                tax_assignment += f"\tsg__{sg}"
+                tax_assignment += f"\t{sg}"
             else:
                 break
+
             if dv != "0":
                 dv = "_".join(dv.split(" "))
-                tax_assignment += f"\tdv__{dv}"
+                tax_assignment += f"\t{dv}"
+            else:
+                break
 
             if sdv != "0":
                 sdv = "_".join(sdv.split(" "))
-                tax_assignment += f"\tsdv__{sdv}"
+                tax_assignment += f"\t{sdv}"
+            else:
+                break
+
             if c != "0":
                 c = "_".join(c.split(" "))
-                tax_assignment += f"\tc__{c}"
+                tax_assignment += f"\t{c}"
             else:
                 break
+
             if o != "0":
                 o = "_".join(o.split(" "))
-                tax_assignment += f"\to__{o}"
+                tax_assignment += f"\t{o}"
             else:
                 break
+
             if f != "0":
                 f = "_".join(f.split(" "))
-                tax_assignment += f"\tf__{f}"
+                tax_assignment += f"\t{f}"
             else:
                 break
+
             if g != "0":
                 g = "_".join(g.split(" "))
-                tax_assignment += f"\tg__{g}"
+                tax_assignment += f"\t{g}"
             else:
                 break
+
             if s != "0":
                 s = "_".join(s.split(" "))
-                tax_assignment += f"\ts__{s}"
+                tax_assignment += f"\t{s}"
             break
 
         if tax_assignment == "":
@@ -210,9 +224,8 @@ def make_tax_assignment_dict_pr2(taxa_df, asv_dict):
     return tax_assignment_dict
 
 def main():
-
     _TAXA, _FWD, _REV, _AMP, _HEADERS, _SAMPLE = parse_args()
-    
+
     fwd_fr = open(_FWD, "r")
     paired_end = True
 
@@ -225,43 +238,43 @@ def main():
     taxa_df = pd.read_csv(_TAXA, sep="\t", dtype=str)
     taxa_df = taxa_df.fillna("0")
     taxa_df = order_df(taxa_df)
-    
-    amp_reads = [ read.strip() for read in list(open(_AMP, "r")) ]
-    headers = [ read.split(" ")[0][1:] for read in list(open(_HEADERS, "r")) ]
+
+    amp_reads = [read.strip() for read in list(open(_AMP, "r"))]
+    headers = [read.split(" ")[0][1:] for read in 
+               list(open(_HEADERS, "r"))]
     amp_region = ".".join(_AMP.split(".")[1:3])
 
     asv_dict = defaultdict(int)
 
     counter = -1
     for line_fwd in fwd_fr:
-
         counter += 1
         line_fwd = line_fwd.strip()
         fwd_asvs = line_fwd.split(",")
 
         if paired_end:
             line_rev = next(rev_fr).strip()
-            rev_asvs = line_rev.split(",")        
+            rev_asvs = line_rev.split(",")
             asv_intersection = list(set(fwd_asvs).intersection(rev_asvs))
-            
+
             if len(asv_intersection) == 0:
                 continue
-            
+
             if len(asv_intersection) == 1 and asv_intersection[0] == "0":
                 continue
         else:
             asv_intersection = fwd_asvs
 
         if headers[counter] in amp_reads:
-            asv_dict[int(asv_intersection[0]) - 1] += 1
-    
+            asv_dict[f"seq_{int(asv_intersection[0]) - 1}"] += 1
+
     fwd_fr.close()
     if paired_end:
         rev_fr.close()
 
     ref_db = ""
 
-    if len(taxa_df.columns) == 8:
+    if len(taxa_df.columns) == 9:
         tax_assignment_dict = make_tax_assignment_dict_silva(taxa_df, asv_dict)
         ref_db = "silva"
     elif len(taxa_df.columns) == 10:
@@ -271,7 +284,7 @@ def main():
     with open(f"./{_SAMPLE}_{amp_region}_{ref_db}_asv_krona_counts.txt", "w") as fw:
         for tax_assignment, count in tax_assignment_dict.items():
             fw.write(f"{count}\t{tax_assignment}\n")
-    
+
 
 if __name__ == "__main__":
-    main()
+    main()

mgnify_pipelines_toolkit/analysis/amplicon/mapseq_to_asv_table.py

@@ -0,0 +1,139 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+# Copyright 2024 EMBL - European Bioinformatics Institute
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+import argparse
+from collections import defaultdict
+import logging
+
+import pandas as pd
+
+logging.basicConfig(level=logging.DEBUG)
+
+def parse_args():
+    parser = argparse.ArgumentParser()
+    parser.add_argument("-i", "--input", required=True, type=str, help="Input from MAPseq output")
+    parser.add_argument(
+        "-l",
+        "--label",
+        choices=["DADA2-SILVA", "DADA2-PR2"],
+        required=True,
+        type=str,
+        help="Database label - either DADA2-SILVA or DADA2-PR2",
+    )
+    parser.add_argument("-s", "--sample", required=True, type=str, help="Sample ID")
+
+    args = parser.parse_args()
+
+    _INPUT = args.input
+    _LABEL = args.label
+    _SAMPLE = args.sample
+
+    return _INPUT, _LABEL, _SAMPLE
+
+def parse_label(label):
+    silva_short_ranks = ["sk__", "k__", "p__", "c__", "o__", "f__", "g__", "s__"]
+    pr2_short_ranks = ["d__", "sg__", "dv__", "sdv__", "c__", "o__", "f__", "g__", "s__"]
+
+    silva_long_ranks = ["Superkingdom", "Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species"]
+    pr2_long_ranks = ["Domain", "Supergroup", "Division", "Subdivision", "Class", "Order", "Family", "Genus", "Species"]
+
+    chosen_short_ranks = ""
+    chosen_long_ranks = ""
+
+    if label == "DADA2-SILVA":
+        chosen_short_ranks = silva_short_ranks
+        chosen_long_ranks = silva_long_ranks
+    elif label == "DADA2-PR2":
+        chosen_short_ranks = pr2_short_ranks
+        chosen_long_ranks = pr2_long_ranks
+    else:
+        logging.error("Incorrect database label - exiting.")
+        exit(1)
+
+    return chosen_short_ranks, chosen_long_ranks
+
+def parse_mapseq(mseq_df, short_ranks, long_ranks):
+    res_dict = defaultdict(list)
+
+    for i in range(len(mseq_df)):
+        asv_id = mseq_df.iloc[i, 0]
+
+        if pd.isna(mseq_df.iloc[i, 1]):
+            tax_ass = [short_ranks[0]]
+        else:
+            tax_ass = mseq_df.iloc[i, 1].split(";")
+
+        res_dict["ASV"].append(asv_id)
+
+        for j in range(len(short_ranks)):
+            curr_rank = long_ranks[j]
+
+            if j >= len(tax_ass):
+                # This would only be true if the assigned taxonomy is shorter than the total reference database taxononmy
+                # so fill each remaining rank with its respective short rank blank
+                curr_tax = short_ranks[j]
+            else:
+                curr_tax = tax_ass[j]
+
+            res_dict[curr_rank].append(curr_tax)
+    res_df = pd.DataFrame.from_dict(res_dict)
+
+    return(res_df)
+
+def process_blank_tax_ends(res_df, ranks):
+    # Necessary function as we want to replace consecutive blank assignments that start at the last rank as NAs
+    # while avoiding making blanks in the middle as NAs
+
+    for i in range(len(res_df)):
+        last_empty_rank = ""
+        currently_empty = False
+        for j in reversed(
+            range(len(ranks))
+        ):  # Parse an assignment backwards, from Species all the way to Superkingdom/Domain
+            curr_rank = res_df.iloc[i, j + 1]
+            if curr_rank in ranks:
+                if last_empty_rank == "":  # Last rank is empty, start window of consecutive blanks
+                    last_empty_rank = j + 1
+                    currently_empty = True
+                elif (
+                    currently_empty
+                ):  # If we're in a window of consecutive blank assignments that started at the beginning
+                    last_empty_rank = j + 1
+                else:
+                    break
+            else:
+                break
+        if last_empty_rank != "":
+            res_df.iloc[i, last_empty_rank:] = "NA"
+        if last_empty_rank == 1:
+            res_df.iloc[i, 1] = ranks[0]
+
+    return res_df
+
+
+def main():    
+    _INPUT, _LABEL, _SAMPLE = parse_args()
+
+    mseq_df = pd.read_csv(_INPUT, header=1, delim_whitespace=True, usecols=[0, 12])
+
+    short_ranks, long_ranks = parse_label(_LABEL)
+    res_df = parse_mapseq(mseq_df, short_ranks, long_ranks)
+    final_res_df = process_blank_tax_ends(res_df, short_ranks)
+
+    final_res_df.to_csv(f"./{_SAMPLE}_{_LABEL}_asv_taxa.tsv", sep="\t", index=False)
+
+if __name__ == "__main__":
+    main()

mgnify_pipelines_toolkit/analysis/amplicon/primer_val_classification.py

@@ -0,0 +1,116 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+# Copyright 2024 EMBL - European Bioinformatics Institute
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+import argparse
+from collections import defaultdict
+import re
+
+from Bio import SeqIO
+import pandas as pd
+
+from mgnify_pipelines_toolkit.constants.var_region_coordinates import REGIONS_16S_BACTERIA, REGIONS_16S_ARCHAEA, REGIONS_18S
+
+STRAND_FWD = "fwd"
+STRAND_REV = "rev"
+
+def parse_args():
+    parser = argparse.ArgumentParser()
+
+    parser.add_argument("-i", "--input", required=True, type=str, help="Path to cmsearch_deoverlap_tblout file")
+    parser.add_argument("-f", "--fasta", required=True, type=str, help="Path to concatenated primers fasta file")
+    parser.add_argument("-s", "--sample", required=True, type=str, help="Sample ID")
+    args = parser.parse_args()
+
+    _INPUT = args.input
+    _FASTA = args.fasta
+    _SAMPLE = args.sample
+
+    return _INPUT, _FASTA, _SAMPLE
+
+
+def get_amp_region(beg, strand, model):
+    prev_region = ""
+
+    for region, region_coords in model.items():
+
+        region_beg = region_coords[0]
+        beg_diff = region_beg - beg
+
+        if strand == STRAND_FWD:
+            if beg_diff > 0:
+                return region
+        else:
+            if beg_diff > 0:
+                return prev_region
+
+        prev_region = region
+
+    return prev_region
+
+
+def main():
+    _INPUT, _FASTA, _SAMPLE = parse_args()    
+    res_dict = defaultdict(list)
+    fasta_dict = SeqIO.to_dict(SeqIO.parse(_FASTA, "fasta"))
+
+    with open(_INPUT, "r") as fr:
+        for line in fr:
+            line = line.strip()
+            line = re.sub("[ \t]+", "\t", line)
+            line_lst = line.split("\t")
+
+            primer_name = line_lst[0]
+            rfam = line_lst[3]
+            beg = float(line_lst[5])
+
+            if rfam == "RF00177":
+                gene = "16S"
+                model = REGIONS_16S_BACTERIA
+            elif rfam == "RF01959":
+                gene = "16S"
+                model = REGIONS_16S_ARCHAEA
+            elif rfam == "RF01960":
+                gene = "18S"
+                model = REGIONS_18S
+            else:
+                continue
+
+            res_dict["Run"].append(_SAMPLE)
+            res_dict["AssertionEvidence"].append("ECO_0000363")
+            res_dict["AssertionMethod"].append("automatic assertion")
+
+            strand = ""
+
+            if "F" in primer_name:
+                strand = STRAND_FWD
+            elif "R" in primer_name:
+                strand = STRAND_REV
+
+            amp_region = get_amp_region(beg, strand, model)
+            primer_seq = str(fasta_dict[primer_name].seq)
+
+            res_dict["Gene"].append(gene)
+            res_dict["VariableRegion"].append(amp_region)
+            res_dict["PrimerName"].append(primer_name)
+            res_dict["PrimerStrand"].append(strand)
+            res_dict["PrimerSeq"].append(primer_seq)
+
+    res_df = pd.DataFrame.from_dict(res_dict)
+    res_df.to_csv(f"./{_SAMPLE}_primer_validation.tsv", sep="\t", index=False)
+
+
+if __name__ == "__main__":
+    main()

mgnify_pipelines_toolkit/constants/tax_ranks.py

@@ -14,5 +14,5 @@
 # See the License for the specific language governing permissions and
 # limitations under the License.
 
-_SILVA_TAX_RANKS = ["Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species"]
+_SILVA_TAX_RANKS = ["Superkingdom", "Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species"]
 _PR2_TAX_RANKS = ["Domain", "Supergroup", "Division", "Subdivision", "Class", "Order", "Family", "Genus", "Species"]

mgnify_pipelines_toolkit/constants/var_region_coordinates.py

@@ -15,36 +15,36 @@
 # limitations under the License.
 
 REGIONS_16S_BACTERIA = {
-    'V1': [69, 92],
-    'V2': [131, 239],
-    'V3': [430, 487],
-    'V4': [566, 672],
-    'V5': [812, 869],
-    'V6': [976, 1033],
-    'V7': [1107, 1164],
-    'V8': [1234, 1285],
-    'V9': [1426, 1456]
+    "V1": [69, 92],
+    "V2": [131, 239],
+    "V3": [430, 487],
+    "V4": [566, 672],
+    "V5": [812, 869],
+    "V6": [976, 1033],
+    "V7": [1107, 1164],
+    "V8": [1234, 1285],
+    "V9": [1426, 1456]
 }
 
 REGIONS_16S_ARCHAEA = {
-    'V1': [61, 79],
-    'V2': [114, 223],
-    'V3': [397, 436],
-    'V4': [516, 623],
-    'V5': [763, 824],
-    'V6': [932, 982],
-    'V7': [1056, 1119],
-    'V8': [1189, 1240],
-    'V9': [1372, 1410]
+    "V1": [61, 79],
+    "V2": [114, 223],
+    "V3": [397, 436],
+    "V4": [516, 623],
+    "V5": [763, 824],
+    "V6": [932, 982],
+    "V7": [1056, 1119],
+    "V8": [1189, 1240],
+    "V9": [1372, 1410]
 }
 
 REGIONS_18S = {
-    'V1': [69, 109],
-    'V2': [136, 298],
-    'V3': [474, 545],
-    'V4': [627, 873],
-    'V5': [1059, 1102],
-    'V7': [1366, 1454],
-    'V8': [1526, 1608],
-    'V9': [1728, 1795]
+    "V1": [69, 109],
+    "V2": [136, 298],
+    "V3": [474, 545],
+    "V4": [627, 873],
+    "V5": [1059, 1102],
+    "V7": [1366, 1454],
+    "V8": [1526, 1608],
+    "V9": [1728, 1795]
 }

{mgnify_pipelines_toolkit-0.1.1.dist-info → mgnify_pipelines_toolkit-0.1.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: mgnify_pipelines_toolkit
-Version: 0.1.1
+Version: 0.1.3
 Summary: Collection of scripts and tools for MGnify pipelines
 Author-email: MGnify team <metagenomics-help@ebi.ac.uk>
 License: Apache Software License 2.0
@@ -55,13 +55,14 @@ You should then be able to run the packages from the command-line. For example t
 ### New script requirements
 
 There are a few requirements for your script:
+
 - It needs to have a named main function of some kind. See `mgnify_pipelines_toolkit/analysis/shared/get_subunits.py` and the `main()` function for an example
 - Because this package is meant to be run from the command-line, make sure your script can easily pass arguments using tools like `argparse` or `click`
 - A small amount of dependencies. This requirement is subjective, but for example if your script only requires a handful of basic packages like `Biopython`, `numpy`, `pandas`, etc., then it's fine. However if the script has a more extensive list of dependencies, a container is probably a better fit.
 
 ### How to add a new script
 
-To add a new Python script, first copy it over to the `mgnify_pipelines_toolkit` directory in this repository, specifically to the subdirectory that makes the most sense. If none of the subdirectories make sense for your script, create a new one. If your script doesn't have a `main()` type function yet, write one. 
+To add a new Python script, first copy it over to the `mgnify_pipelines_toolkit` directory in this repository, specifically to the subdirectory that makes the most sense. If none of the subdirectories make sense for your script, create a new one. If your script doesn't have a `main()` type function yet, write one.
 
 Then, open `pyproject.toml` as you will need to add some bits. First, add any missing dependencies (include the version) to the `dependencies` field.
 
@@ -73,7 +74,7 @@ Then, scroll down to the `[project.scripts]` line. Here, you will create an alia
 
 - `get_subunits` is the alias
 - `mgnify_pipelines_toolkit.analysis.shared.get_subunits` will link the alias to the script with the path `mgnify_pipelines_toolkit/analysis/shared/get_subunits.py`
-- `:main` will specifically call the function named `main()` when the alias is run. 
+- `:main` will specifically call the function named `main()` when the alias is run.
 
 When you have setup this command, executing `get_subunits` on the command-line will be the equivalent of doing:
 
@@ -86,4 +87,5 @@ Finally, you will need to bump up the version in the `version` line.
 At the moment, these should be the only steps required to setup your script in this package (which is subject to change).
 
 ### Building and uploading to PyPi
+
 The building and pushing of the package is automated by GitHub Actions, which will activate only on a new release. Bioconda should then automatically pick up the new PyPi release and push it to their recipes, though it's worth keeping an eye on their automated pull requests just in case [here](https://github.com/bioconda/bioconda-recipes/pulls).

{mgnify_pipelines_toolkit-0.1.1.dist-info → mgnify_pipelines_toolkit-0.1.3.dist-info}/RECORD

@@ -1,12 +1,14 @@
 mgnify_pipelines_toolkit/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mgnify_pipelines_toolkit/analysis/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-mgnify_pipelines_toolkit/analysis/amplicon/amplicon_utils.py,sha256=OkQsiVTh6Y-bYPZMHOnD9xn0dd9bENySFwgibK6x6CU,6549
+mgnify_pipelines_toolkit/analysis/amplicon/amplicon_utils.py,sha256=EvWTFV4gGn0SkrjwC2hzvNGSXFLeyFDmVj2QDa5DmtE,6402
 mgnify_pipelines_toolkit/analysis/amplicon/are_there_primers.py,sha256=fNjzWpXjqpQNRyWG6CoROyHAvA4lZsvPh8sxDpjyMkY,5141
 mgnify_pipelines_toolkit/analysis/amplicon/assess_inflection_point_mcp.py,sha256=pBNpqHFb0zzWgTD1mY3Q5MslQ5nmT99-pSHpyngVEuo,7159
 mgnify_pipelines_toolkit/analysis/amplicon/assess_mcp_proportions.py,sha256=BCzLlfvRqiPC-YwzET901f_d0anYt1zpf5y0iOCQnvs,5191
 mgnify_pipelines_toolkit/analysis/amplicon/classify_var_regions.py,sha256=89LKH4rvqRydAEGvfaWqIClcitJ1Vbu7b5d4FApzGp4,18392
 mgnify_pipelines_toolkit/analysis/amplicon/find_mcp_inflection_points.py,sha256=5emeZjk8891IgyL5ygVcr8wMP-hGEJoEs2rcBbseWj0,3536
-mgnify_pipelines_toolkit/analysis/amplicon/make_asv_count_table.py,sha256=3ODVM7YxYrPiM95UnbYiKmtN93PdlF_JmjZYTEyoCL8,8468
+mgnify_pipelines_toolkit/analysis/amplicon/make_asv_count_table.py,sha256=dZIygTbMZvVbSVBmFFAZz7x24oQEpvdEOTpTcnYAyoM,8444
+mgnify_pipelines_toolkit/analysis/amplicon/mapseq_to_asv_table.py,sha256=lWVIhDxfsTuDzWjjUlMGx3RL7iD_Yy8m9Ppc9wjfCFg,4765
+mgnify_pipelines_toolkit/analysis/amplicon/primer_val_classification.py,sha256=F4ALSuH8N-0hHUqPCFwHgoAnteb2Ft3tUN9j6DaD5h8,3539
 mgnify_pipelines_toolkit/analysis/amplicon/remove_ambiguous_reads.py,sha256=yhDJFGD3t3TMfUlBCJGwzlK4IjFwm7Bporwp-aIM8uU,3139
 mgnify_pipelines_toolkit/analysis/amplicon/rev_comp_se_primers.py,sha256=VoSbou3BAZL4bJes4FsYJvmd45_PjKj8F2sQDIyLDoI,1680
 mgnify_pipelines_toolkit/analysis/amplicon/standard_primer_matching.py,sha256=odGVde7Ct9dS2aqsySWgdgVLCOqfr_ZGeHFcXcuukxs,10846
@@ -15,12 +17,12 @@ mgnify_pipelines_toolkit/analysis/shared/get_subunits.py,sha256=vvhn8O9t1zzD8rIv
 mgnify_pipelines_toolkit/analysis/shared/get_subunits_coords.py,sha256=hFXUTZb-etmJS7Si3mCCVCXV5ZYN0tP6FSbeiVxG1jo,1879
 mgnify_pipelines_toolkit/analysis/shared/mapseq2biom.py,sha256=glvql2Y-BTyA1hTIZz2nlmST3SE6LJbep9sKdMH-vaI,5565
 mgnify_pipelines_toolkit/constants/regex_ambiguous_bases.py,sha256=oVeeCy33bY1WJ-rffOULZ3ogi48Jz0FfTS73MPTur-A,1095
-mgnify_pipelines_toolkit/constants/tax_ranks.py,sha256=QOVkxZjMTYwIxCN7L4CLYFQjWoxsuZ4WKlKiTblD4tM,851
+mgnify_pipelines_toolkit/constants/tax_ranks.py,sha256=fP97JvlUdxJBakJ694VresIY8-N3pcU99m7kZ9buKys,867
 mgnify_pipelines_toolkit/constants/thresholds.py,sha256=7J3caCikkEcLdKF4zSR0z8qMQw4-h9aSkSbFbS0LNg4,873
-mgnify_pipelines_toolkit/constants/var_region_coordinates.py,sha256=GAW1DdQ9IM3dCBMVMS6ZNtxVFEGwEGdJaDNorVAzR1g,1303
-mgnify_pipelines_toolkit-0.1.1.dist-info/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
-mgnify_pipelines_toolkit-0.1.1.dist-info/METADATA,sha256=D1u8X2ougLo4WNGCmxy_0NZDOyfN5UF5_N65zF2dhR8,4950
-mgnify_pipelines_toolkit-0.1.1.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-mgnify_pipelines_toolkit-0.1.1.dist-info/entry_points.txt,sha256=GH2Vn239CpRUjapHA555c4YQ7aKsxBlaFu-XvloJb4I,1114
-mgnify_pipelines_toolkit-0.1.1.dist-info/top_level.txt,sha256=xA_wC7C01V3VwuDnqwRM2QYeJJ45WtvF6LVav4tYxuE,25
-mgnify_pipelines_toolkit-0.1.1.dist-info/RECORD,,
+mgnify_pipelines_toolkit/constants/var_region_coordinates.py,sha256=E8Cd3A1Hj9M95zw9Ut-2x8sE6_PlH6RJJEoikyZUMaQ,1303
+mgnify_pipelines_toolkit-0.1.3.dist-info/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
+mgnify_pipelines_toolkit-0.1.3.dist-info/METADATA,sha256=D8bYOR2kQZzJPdqtFkHj_Xd4axEHjzJPJXKAHtFj8L0,4950
+mgnify_pipelines_toolkit-0.1.3.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+mgnify_pipelines_toolkit-0.1.3.dist-info/entry_points.txt,sha256=K8nqnyAQG9jqHGgIfMIaCIe20u5a0FFCCqJWi4DoD2U,1306
+mgnify_pipelines_toolkit-0.1.3.dist-info/top_level.txt,sha256=xA_wC7C01V3VwuDnqwRM2QYeJJ45WtvF6LVav4tYxuE,25
+mgnify_pipelines_toolkit-0.1.3.dist-info/RECORD,,

{mgnify_pipelines_toolkit-0.1.1.dist-info → mgnify_pipelines_toolkit-0.1.3.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.42.0)
+Generator: bdist_wheel (0.43.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

{mgnify_pipelines_toolkit-0.1.1.dist-info → mgnify_pipelines_toolkit-0.1.3.dist-info}/entry_points.txt

@@ -8,6 +8,8 @@ get_subunits = mgnify_pipelines_toolkit.analysis.shared.get_subunits:main
 get_subunits_coords = mgnify_pipelines_toolkit.analysis.shared.get_subunits_coords:main
 make_asv_count_table = mgnify_pipelines_toolkit.analysis.amplicon.make_asv_count_table:main
 mapseq2biom = mgnify_pipelines_toolkit.analysis.shared.mapseq2biom:main
+mapseq_to_asv_table = mgnify_pipelines_toolkit.analysis.amplicon.mapseq_to_asv_table:main
+primer_val_classification = mgnify_pipelines_toolkit.analysis.amplicon.primer_val_classification:main
 remove_ambiguous_reads = mgnify_pipelines_toolkit.analysis.amplicon.remove_ambiguous_reads:main
 rev_comp_se_primers = mgnify_pipelines_toolkit.analysis.amplicon.rev_comp_se_primers:main
 standard_primer_matching = mgnify_pipelines_toolkit.analysis.amplicon.standard_primer_matching:main