kssdtree 2.0.1__tar.gz → 2.0.3__tar.gz

{kssdtree-2.0.1/kssdtree.egg-info → kssdtree-2.0.3}/PKG-INFO

@@ -1,8 +1,11 @@
-Metadata-Version: 2.1
+Metadata-Version: 1.1
 Name: kssdtree
-Version: 2.0.1
+Version: 2.0.3
 Summary: Kssdtree is a versatile Python package for phylogenetic analysis. It also provides one-stop tree construction and visualization. It can handle DNA sequences of both fasta or fastq format, whether gzipped or not. 
 Home-page: https://github.com/yhlink/kssdtree
-Download-URL: https://pypi.org/project/kssdtree
 Author: Hang Yang
 Author-email: yhlink1207@gmail.com
+License: UNKNOWN
+Download-URL: https://pypi.org/project/kssdtree
+Description: UNKNOWN
+Platform: UNKNOWN

kssdtree-2.0.3/README.md

@@ -0,0 +1,8 @@
+Kssdtree is a versatile Python package for phylogenetic analysis, offering three distinct pipelines: the Routine Pipeline, the Reference Subtraction Pipeline, and the GTDB-based Phylogenetic Placement Pipeline.
+
+Routine Pipeline: A general-purpose tool for phylogenetic analysis of user genomic data.
+Reference Subtraction Pipeline: Designed for intra-species phylogenomic analysis.
+GTDB-based Phylogenetic Placement Pipeline: Facilitates the search for similar genomes in the Genome Taxonomy Database (GTDB), conducting phylogenetic analysis alongside these genomes and positioning the input genomes within the entire prokaryotic tree of life.
+Kssdtree also provides one-stop tree construction and visualization. It can handle DNA sequences in both fasta and fastq formats, whether gzipped or not. Additionally, Kssdtree is compatible with multiple platforms (Linux, MacOS, and Windows) and can be run using Jupyter notebooks.
+
+More usages about Kssdtree, please see Kssdtree documentation (https://kssdtree.readthedocs.io/en/latest).

{kssdtree-2.0.1 → kssdtree-2.0.3}/command_set.c

@@ -25,7 +25,6 @@
 #include <errno.h>
 #include <math.h>
 #include <unistd.h>
-
 const char skch_prefix[]="combco";
 const char idx_prefix[]="combco.index";
 const char pan_prefix[]="pan";

{kssdtree-2.0.1 → kssdtree-2.0.3/kssdtree.egg-info}/PKG-INFO

@@ -1,8 +1,11 @@
-Metadata-Version: 2.1
+Metadata-Version: 1.1
 Name: kssdtree
-Version: 2.0.1
+Version: 2.0.3
 Summary: Kssdtree is a versatile Python package for phylogenetic analysis. It also provides one-stop tree construction and visualization. It can handle DNA sequences of both fasta or fastq format, whether gzipped or not. 
 Home-page: https://github.com/yhlink/kssdtree
-Download-URL: https://pypi.org/project/kssdtree
 Author: Hang Yang
 Author-email: yhlink1207@gmail.com
+License: UNKNOWN
+Download-URL: https://pypi.org/project/kssdtree
+Description: UNKNOWN
+Platform: UNKNOWN

{kssdtree-2.0.1 → kssdtree-2.0.3}/kssdtree.egg-info/requires.txt

@@ -1,3 +1,4 @@
 pyqt5
 ete3
 requests
+pandas

{kssdtree-2.0.1 → kssdtree-2.0.3}/kssdtree.py

@@ -28,23 +28,18 @@ def sketch(shuf_file=None, genome_files=None, output=None, set_opt=None):
         if not os.path.exists(shuf_file):
             if shuf_file in ['L3K9.shuf', './L3K9.shuf', 'L3K10.shuf', './L3K10.shuf']:
                 print('Downloading...', shuf_file)
-                import http.client
-                http.client.HTTPConnection._http_vsn = 10
-                http.client.HTTPConnection._http_vsn_str = 'HTTP/1.0'
                 if shuf_file == 'L3K9.shuf' or shuf_file == './L3K9.shuf':
                     url = 'https://zenodo.org/records/12699159/files/L3K9.shuf?download=1'
                 else:
                     url = 'https://zenodo.org/records/12699159/files/L3K10.shuf?download=1'
-                start_time = time.time()
-                response = requests.get(url, stream=True)
-                with open(shuf_file, 'wb') as file:
-                    for chunk in response.iter_content(chunk_size=1024):
-                        if chunk:
-                            file.write(chunk)
+                headers = {'Accept-Encoding': 'gzip, deflate'}
+                response = requests.get(url, headers=headers, stream=True)
+                with open(shuf_file, 'wb') as f:
+                    for chunk in response.iter_content(chunk_size=8192):
+                    f.write(chunk)
                 end_time = time.time()
-                if end_time - start_time > 120:
-                    print(
-                        "Network timeout, please manually download from github (https://github.com/yhlink/kssdtree/tree/master/shuffle_file)")
+                if end_time - start_time > 200:
+                    print("Network timeout, please manually download from https://zenodo.org/records/12699159")
                     return False
                 print('Download finished: ', shuf_file)
             elif shuf_file in ['L2K8.shuf', 'L2K9.shuf', 'L3K11.shuf', './L2K8.shuf', './L2K9.shuf', './L3K11.shuf']:
@@ -64,9 +59,9 @@ def sketch(shuf_file=None, genome_files=None, output=None, set_opt=None):
         print('Sketching...')
         start = time.time()
         if set_opt:
-            kssd.dist_dispatch(shuf_file, genome_files, output, 1, 0, 0)
+            kssd.dist_dispatch(shuf_file, genome_files, output, 1, 0, 0, '')
         else:
-            kssd.dist_dispatch(shuf_file, genome_files, output, 0, 0, 0)
+            kssd.dist_dispatch(shuf_file, genome_files, output, 0, 0, 0, '')
         end = time.time()
         print('Sketch spend time：%.2fs' % (end - start))
         print('Sketch finished!')
@@ -76,7 +71,7 @@ def sketch(shuf_file=None, genome_files=None, output=None, set_opt=None):
         return False
 
 
-def dist(genome_sketch=None, output=None, flag=None):
+def dist(genome_sketch=None, output=None, metric = None, flag=None):
     if genome_sketch is not None and output is not None:
         if not os.path.exists(genome_sketch):
             print('No such file or directory: ', genome_sketch)
@@ -86,6 +81,9 @@ def dist(genome_sketch=None, output=None, flag=None):
         #     return False
         if flag is None:
             flag = 0
+        if metric is None:
+            metric = 'mash'
+        
         print('Disting...')
         start = time.time()
         if '/' in output:
@@ -97,11 +95,15 @@ def dist(genome_sketch=None, output=None, flag=None):
         else:
             output_name = output
         if output_name.endswith(".phy") or output_name.endswith(".phylip"):
-            kssd.dist_dispatch(genome_sketch, output, genome_sketch, 2, 0, flag)
-            end = time.time()
-            print('Dist spend time：%.2fs' % (end - start))
-            print('Dist finished!')
-            return True
+            if metric not in ['mash', 'aaf']:
+                print('Metric type error, only supports mash or aaf distance')
+                return False
+            else:
+                kssd.dist_dispatch(genome_sketch, output, genome_sketch, 2, 0, flag, metric)
+                end = time.time()
+                print('Dist spend time：%.2fs' % (end - start))
+                print('Dist finished!')
+                return True
         else:
             print('Output type error, only supports .phylip (.phy) format:', output_name)
             return False
@@ -207,7 +209,12 @@ def visualize(newick=None, taxonomy=None, mode=None):
             return False
         if mode is None:
             mode = 'r'
-        toolutils.view_tree(newick, taxonomy, mode=mode)
+        if taxonomy is not None and mode == 'c':
+            print('Warning: this pipeline only support 'r' (rectangle) mode !!!')
+            mode = 'r'
+            toolutils.view_tree(newick, taxonomy, mode=mode)
+        else:
+            toolutils.view_tree(newick, taxonomy, mode=mode)
     else:
         print('Args error!!!')
         return False
@@ -280,6 +287,9 @@ def subtract(ref_sketch=None, genome_sketch=None, output=None, flag=None):
 def quick(shuf_file=None, genome_files=None, output=None, reference=None, database=None, method='nj', mode='r', N=0):
     if reference is None and database is None:
         if shuf_file is not None and genome_files is not None and output is not None:
+            if toolutils.is_positive_integer(N) or toolutils.is_negative_integer(N):
+                print("N must = 0 !!!")
+                return False
             timeStamp = int(time.mktime(time.localtime(time.time())))
             temp_sketch = toolutils.rs() + '_sketch_' + str(timeStamp)
             temp_phy = toolutils.rs() + '_temp.phy'
@@ -325,7 +335,7 @@ def quick(shuf_file=None, genome_files=None, output=None, reference=None, databa
     elif reference is None and database == 'gtdbr214':
         if shuf_file is not None and genome_files is not None and output is not None:
             if not toolutils.is_positive_integer(N):
-                print("N must >0 !!!")
+                print("N must > 0 !!!")
                 return False
             if shuf_file != 'L3K9.shuf':
                 print("shuf_file must be set to 'L3K9.shuf'")
@@ -354,83 +364,86 @@ def quick(shuf_file=None, genome_files=None, output=None, reference=None, databa
             return False
     elif reference is not None and database is None:
         if shuf_file is not None and genome_files is not None and output is not None and method in ['nj', 'dnj']:
-            if shuf_file is not None and genome_files is not None and output is not None and method in ['nj', 'dnj']:
-                if not toolutils.allowed_file(genome_files):
-                    num = toolutils.get_file_num(genome_files)
-                    if num == 1:
-                        print('genome_files is a folder containing at least two . fasta or .fastq files!!!')
-                        return False
-                else:
-                    print('genome_files is a folder containing at least two . fasta or .fastq files, not a file!!!')
+            if toolutils.is_positive_integer(N) or toolutils.is_negative_integer(N):
+                print("N must = 0 !!!")
+                return False
+            if not toolutils.allowed_file(genome_files):
+                num = toolutils.get_file_num(genome_files)
+                if num == 1:
+                    print('genome_files is a folder containing at least two . fasta or .fastq files!!!')
                     return False
-                timeStamp = int(time.mktime(time.localtime(time.time())))
-                temp_reference_sketch = toolutils.rs() + '_ref_sketch_' + str(timeStamp)
-                temp_genomes_sketch = toolutils.rs() + '_sketch_' + str(timeStamp)
-                if not toolutils.allowed_file(reference):
-                    # cur_path = os.getcwd()
-                    # ref_path = os.path.join(cur_path, reference)
-                    num = toolutils.get_file_num(reference)
-                    if num == 1:
-                        temp_union_sketch = temp_reference_sketch
-                    else:
-                        temp_union_sketch = toolutils.rs() + '_ref_union_sketch_' + str(timeStamp)
-                else:
+            else:
+                print('genome_files is a folder containing at least two . fasta or .fastq files, not a file!!!')
+                return False
+            timeStamp = int(time.mktime(time.localtime(time.time())))
+            temp_reference_sketch = toolutils.rs() + '_ref_sketch_' + str(timeStamp)
+            temp_genomes_sketch = toolutils.rs() + '_sketch_' + str(timeStamp)
+            if not toolutils.allowed_file(reference):
+                # cur_path = os.getcwd()
+                # ref_path = os.path.join(cur_path, reference)
+                num = toolutils.get_file_num(reference)
+                if num == 1:
                     temp_union_sketch = temp_reference_sketch
-                temp_subtract_sketch = toolutils.rs() + '_subtract_sketch_' + str(timeStamp)
-                temp_phy = toolutils.rs() + '_temp.phy'
-                print('Step1...')
-                s1 = sketch(shuf_file=shuf_file, genome_files=reference, output=temp_reference_sketch, set_opt=True)
-                if not s1:
-                    return False
-                s2 = sketch(shuf_file=shuf_file, genome_files=genome_files, output=temp_genomes_sketch, set_opt=True)
-                if not s2:
-                    return False
-                print('Step2...')
-                s3 = union(ref_sketch=temp_reference_sketch, output=temp_union_sketch)
-                if not s3:
-                    return False
-                s4 = subtract(ref_sketch=temp_union_sketch, genome_sketch=temp_genomes_sketch,
-                              output=temp_subtract_sketch, flag=1)
-                if not s4:
-                    return False
-                print('Step3...')
-                if method == 'nj':
-                    s5 = dist(genome_sketch=temp_subtract_sketch, output=temp_phy,
-                              flag=0)
                 else:
-                    s5 = dist(genome_sketch=temp_subtract_sketch, output=temp_phy,
-                              flag=1)
-                if not s5:
-                    return False
-                print('Step4...')
-                s6 = build(phylip=temp_phy, output=output, method=method)
-                if not s6:
-                    return False
-                print('Step5...')
-                print('Tree visualization finished!')
-                visualize(newick=output, mode=mode)
-                if platform.system() == 'Linux':
-                    current_directory = os.getcwd()
-                    temp_dir1 = os.path.join(current_directory, temp_reference_sketch)
-                    temp_dir2 = os.path.join(current_directory, temp_genomes_sketch)
-                    temp_dir3 = os.path.join(current_directory, temp_union_sketch)
-                    temp_dir4 = os.path.join(current_directory, temp_subtract_sketch)
-                    temp_dir5 = os.path.join(current_directory, 'distout')
-                    if os.path.exists(temp_dir1):
-                        shutil.rmtree(temp_dir1)
-                    if os.path.exists(temp_dir2):
-                        shutil.rmtree(temp_dir2)
-                    if os.path.exists(temp_dir3):
-                        shutil.rmtree(temp_dir3)
-                    if os.path.exists(temp_dir4):
-                        shutil.rmtree(temp_dir4)
-                    if os.path.exists(temp_dir5):
-                        shutil.rmtree(temp_dir5)
-                    if os.path.exists(temp_phy):
-                        os.remove(temp_phy)
+                    temp_union_sketch = toolutils.rs() + '_ref_union_sketch_' + str(timeStamp)
             else:
-                print('Args error, please see https://kssdtree.readthedocs.io/en/latest!!!')
+                temp_union_sketch = temp_reference_sketch
+            temp_subtract_sketch = toolutils.rs() + '_subtract_sketch_' + str(timeStamp)
+            temp_phy = toolutils.rs() + '_temp.phy'
+            print('Step1...')
+            s1 = sketch(shuf_file=shuf_file, genome_files=reference, output=temp_reference_sketch, set_opt=True)
+            if not s1:
+                return False
+            s2 = sketch(shuf_file=shuf_file, genome_files=genome_files, output=temp_genomes_sketch, set_opt=True)
+            if not s2:
+                return False
+            print('Step2...')
+            s3 = union(ref_sketch=temp_reference_sketch, output=temp_union_sketch)
+            if not s3:
+                return False
+            s4 = subtract(ref_sketch=temp_union_sketch, genome_sketch=temp_genomes_sketch,
+                          output=temp_subtract_sketch, flag=1)
+            if not s4:
+                return False
+            print('Step3...')
+            if method == 'nj':
+                s5 = dist(genome_sketch=temp_subtract_sketch, output=temp_phy,
+                          flag=0)
+            else:
+                s5 = dist(genome_sketch=temp_subtract_sketch, output=temp_phy,
+                          flag=1)
+            if not s5:
+                return False
+            print('Step4...')
+            s6 = build(phylip=temp_phy, output=output, method=method)
+            if not s6:
                 return False
+            print('Step5...')
+            print('Tree visualization finished!')
+            visualize(newick=output, mode=mode)
+            if platform.system() == 'Linux':
+                current_directory = os.getcwd()
+                temp_dir1 = os.path.join(current_directory, temp_reference_sketch)
+                temp_dir2 = os.path.join(current_directory, temp_genomes_sketch)
+                temp_dir3 = os.path.join(current_directory, temp_union_sketch)
+                temp_dir4 = os.path.join(current_directory, temp_subtract_sketch)
+                temp_dir5 = os.path.join(current_directory, 'distout')
+                if os.path.exists(temp_dir1):
+                    shutil.rmtree(temp_dir1)
+                if os.path.exists(temp_dir2):
+                    shutil.rmtree(temp_dir2)
+                if os.path.exists(temp_dir3):
+                    shutil.rmtree(temp_dir3)
+                if os.path.exists(temp_dir4):
+                    shutil.rmtree(temp_dir4)
+                if os.path.exists(temp_dir5):
+                    shutil.rmtree(temp_dir5)
+                if os.path.exists(temp_phy):
+                    os.remove(temp_phy)
+        else:
+            print('Args error, please see https://kssdtree.readthedocs.io/en/latest!!!')
+            return False
+
     else:
         print('Pipeline error, please see https://kssdtree.readthedocs.io/en/latest!!!')
         return False

{kssdtree-2.0.1 → kssdtree-2.0.3}/mman.c

@@ -1,4 +1,3 @@
-
 #ifdef _WIN32
 #include <windows.h>
 #include <io.h>

{kssdtree-2.0.1 → kssdtree-2.0.3}/pykssd.c

@@ -121,12 +121,18 @@ int create_matrix(char *input_name, char *output_name, int flag) {
             for (int j = 0; j <= i; j++)
                 distances[i][j] = 0.0;
         }
-    } else {
+    } else if (flag == 1) {
         for (int i = 0; i < num_seqs; i++) {
             distances[i] = malloc(i * sizeof(double));
             for (int j = 0; j < i; j++)
                 distances[i][j] = 0.0;
         }
+    } else {
+        for (int i = 0; i < num_seqs; i++) {
+            distances[i] = malloc(num_seqs * sizeof(double));
+            for (int j = 0; j < num_seqs; j++)
+                distances[i][j] = 0.0;
+        }
     }
     rewind(fp);
     fgets(line, PATHLEN, fp);
@@ -144,7 +150,7 @@ int create_matrix(char *input_name, char *output_name, int flag) {
                 j += 1;
             }
         }
-    } else {
+    } else if (flag == 1) {
         while (fgets(line, PATHLEN, fp) != NULL) {
             sscanf(line, "%*s %*s %*s %*s %lf", &distance);
             if (j < i) {
@@ -156,6 +162,16 @@ int create_matrix(char *input_name, char *output_name, int flag) {
                 j += 1;
             }
         }
+    } else {
+        while (fgets(line, PATHLEN, fp) != NULL) {
+            sscanf(line, "%*s %*s %*s %*s %lf", &distance);
+            distances[i][j] = distance;
+            i += 1;
+            if (i == num_seqs && j < num_seqs) {
+                i = 0;
+                j += 1;
+            }
+        }
     }
 
     for (int i = 0; i < num_seqs; i++) {
@@ -187,7 +203,7 @@ int create_matrix(char *input_name, char *output_name, int flag) {
             }
             fprintf(fo, "\n");
         }
-    } else {
+    } else if (flag == 1) {
         for (int i = 0; i < num_seqs; i++) {
             fprintf(fo, "%s\t", seq_names[i]);
             for (int j = 0; j < num_seqs; j++) {
@@ -197,6 +213,14 @@ int create_matrix(char *input_name, char *output_name, int flag) {
             }
             fprintf(fo, "\n");
         }
+    } else {
+        for (int i = 0; i < num_seqs; i++) {
+            fprintf(fo, "%s\t", seq_names[i]);
+            for (int j = 0; j < num_seqs; j++) {
+                fprintf(fo, "%.6f\t", distances[i][j]);
+            }
+            fprintf(fo, "\n");
+        }
     }
     for (int i = 0; i < num_seqs; i++) {
         free(distances[i]);
@@ -235,8 +259,9 @@ static PyObject *py_dist_dispatch(PyObject *self, PyObject *args) {
     char *str3;
     int flag1;
     int flag2;
+    char *str4;
     int N;
-    if (!PyArg_ParseTuple(args, "sssiii", &str1, &str2, &str3, &flag1, &N, &flag2)) {
+    if (!PyArg_ParseTuple(args, "sssiiis", &str1, &str2, &str3, &flag1, &N, &flag2, &str4)) {
         return NULL;
     }
     if (flag1 == 0) {
@@ -394,6 +419,11 @@ static PyObject *py_dist_dispatch(PyObject *self, PyObject *args) {
         dist_opt_val3.num_remaining_args = 1;
         dist_opt_val3.remaining_args = &str3;
         dist_opt_val3.num_neigb = N;
+	if (strcmp(str4, "mash") == 0) {
+		dist_opt_val3.metric = 0;
+	} else {
+		dist_opt_val3.metric = 1;
+	}
 #ifdef _OPENMP
         if(dist_opt_val3.p == 0)
             dist_opt_val3.p = omp_get_num_procs();
@@ -628,4 +658,4 @@ static struct PyModuleDef kssdmodule = {
 
 PyMODINIT_FUNC PyInit_kssd(void) {
     return PyModule_Create(&kssdmodule);
-}
+}

{kssdtree-2.0.1 → kssdtree-2.0.3}/setup.py

@@ -2,14 +2,7 @@ import sys
 from setuptools import setup, Extension, find_packages
 from os import environ
 import os
-import platform
 
-
-# def get_gcc_version():
-#     gcc_version = subprocess.check_output(['gcc', '--version']).decode('utf-8')
-#     version_line = gcc_version.split('\n', 1)[0]
-#     version_str = version_line.split()[-1]
-#     return version_str
 extra_compile_args = []
 extra_link_args = []
 if 'darwin' in sys.platform:
@@ -26,12 +19,12 @@ if 'darwin' in sys.platform:
             gcc_path = "/opt/homebrew/bin/gcc-11"
         elif 'gcc-12' == gcc_version:
             gcc_path = "/opt/homebrew/bin/gcc-12"
+        elif 'gcc-13' == gcc_version:
+            gcc_path = "/opt/homebrew/bin/gcc-13"
         elif 'gcc-14' == gcc_version:
             gcc_path = "/opt/homebrew/bin/gcc-14"
-        elif 'gcc-15' == gcc_version:
-            gcc_path = "/opt/homebrew/bin/gcc-15"
         else:
-            gcc_path = "/opt/homebrew/bin/gcc-13"
+            gcc_path = ""
         extra_compile_args = ['-fopenmp']
         extra_link_args = ['-fopenmp']
         os.environ["CC"] = gcc_path
@@ -85,25 +78,16 @@ include_dirs1 = ['kssdheaders']
 include_dirs2 = ['njheaders']
 include_dirs3 = ['dnjheaders']
 
-if 'darwin' in sys.platform:
-    if platform.machine() == 'arm64':
-        require_pakages = []
-    else:
-        require_pakages = [
-            'pyqt5',
-            'ete3',
-            'requests'
-        ]
-else:
-    require_pakages = [
-        'pyqt5',
-        'ete3',
-        'requests'
-    ]
+require_pakages = [
+    'pyqt5',
+    'ete3',
+    'requests',
+    'pandas'
+]
 
 setup(
     name='kssdtree',
-    version='2.0.1',
+    version='2.0.3',
     author='Hang Yang',
     author_email='yhlink1207@gmail.com',
     description="Kssdtree is a versatile Python package for phylogenetic analysis. It also provides one-stop tree construction and visualization. It can handle DNA sequences of both fasta or fastq format, whether gzipped or not. ",
@@ -126,4 +110,3 @@ setup(
     include_package_data=True
 )
 
-

{kssdtree-2.0.1 → kssdtree-2.0.3}/toolutils.py

@@ -8,7 +8,7 @@ import pandas as pd
 import string
 
 def allowed_file(filename):
-    allowed_extensions = ['.fa', '.fa.gz', '.fasta', '.fasta.gz', '.fna', '.fna.gz', '.fastq', '.fastq.gz']
+    allowed_extensions = ['.fa', '.fa.gz', '.fasta', '.fasta.gz', '.fna', '.fna.gz', '.fastq', '.fastq.gz', '.fq', '.fq.gz']
     return any(filename.endswith(ext) for ext in allowed_extensions)
 
 def rs():
@@ -24,6 +24,11 @@ def is_positive_integer(num):
     else:
         return False
 
+def is_negative_integer(num):
+    if isinstance(num, int) and num < 0:
+        return True
+    else:
+        return False
 
 def randomcolor():
     colorArr = ['1', '2', '3', '4', '5', '6', '7', '8', '9', 'A', 'B', 'C', 'D', 'E', 'F']

kssdtree-2.0.1/README.md

@@ -1,12 +0,0 @@
-# Kssdree: an interactive Python package for phylogenetic analysis based on sketching technique
-
-Kssdtree is a versatile Python package for phylogenetic analysis, including three different pipelines: routine pipeline, reference subtraction pipeline and phylogenetic placement pipeline. The routine pipeline serves as a versatile tool for general-purpose phylogenetic analysis of users' genomic data.  The reference subtraction pipeline designs for population-level phylogenetic analysis. The phylogenetic placement pipeline facilitates the search for similar genomes in the Genome Taxonomy Database (GTDB). It conducts phylogenetic analysis alongside these similar genomes and positions the input genomes within the entire prokaryotic tree of life. 
-
-It also provides one-stop tree construction and visualization. It can handle DNA sequences of both fasta or fastq format, whether gzipped or not. Kssdtree can run on multiple platforms (Linux, Windows, and MacOS) with Jupyter notebooks.
-
-More usages about Kssdtree, please see Kssdtree documentation (https://kssdtree.readthedocs.io/en/latest).
- 
-
-
-
-