PyamilySeq 0.6.0__py3-none-any.whl → 0.7.1__py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- PyamilySeq/Constants.py +1 -1
- PyamilySeq/PyamilySeq.py +81 -39
- PyamilySeq/PyamilySeq_Genus.py +85 -102
- PyamilySeq/PyamilySeq_Species.py +101 -94
- PyamilySeq/Seq_Combiner.py +26 -7
- PyamilySeq/clusterings.py +111 -73
- PyamilySeq/utils.py +117 -7
- PyamilySeq-0.7.1.dist-info/METADATA +250 -0
- PyamilySeq-0.7.1.dist-info/RECORD +14 -0
- {PyamilySeq-0.6.0.dist-info → PyamilySeq-0.7.1.dist-info}/WHEEL +1 -1
- PyamilySeq/CD-Hit_StORF-Reporter_Cross-Genera_Builder.py +0 -600
- PyamilySeq-0.6.0.dist-info/METADATA +0 -147
- PyamilySeq-0.6.0.dist-info/RECORD +0 -15
- {PyamilySeq-0.6.0.dist-info → PyamilySeq-0.7.1.dist-info}/LICENSE +0 -0
- {PyamilySeq-0.6.0.dist-info → PyamilySeq-0.7.1.dist-info}/entry_points.txt +0 -0
- {PyamilySeq-0.6.0.dist-info → PyamilySeq-0.7.1.dist-info}/top_level.txt +0 -0
|
@@ -1,147 +0,0 @@
|
|
|
1
|
-
Metadata-Version: 2.1
|
|
2
|
-
Name: PyamilySeq
|
|
3
|
-
Version: 0.6.0
|
|
4
|
-
Summary: PyamilySeq - A a tool to look for sequence-based gene families identified by clustering methods such as CD-HIT, DIAMOND, BLAST or MMseqs2.
|
|
5
|
-
Home-page: https://github.com/NickJD/PyamilySeq
|
|
6
|
-
Author: Nicholas Dimonaco
|
|
7
|
-
Author-email: nicholas@dimonaco.co.uk
|
|
8
|
-
Project-URL: Bug Tracker, https://github.com/NickJD/PyamilySeq/issues
|
|
9
|
-
Classifier: Programming Language :: Python :: 3
|
|
10
|
-
Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
|
|
11
|
-
Classifier: Operating System :: OS Independent
|
|
12
|
-
Requires-Python: >=3.6
|
|
13
|
-
Description-Content-Type: text/markdown
|
|
14
|
-
License-File: LICENSE
|
|
15
|
-
|
|
16
|
-
# PyamilySeq - !BETA!
|
|
17
|
-
**PyamilySeq** (Family Seek) is a Python tool for clustering gene sequences into families based on sequence similarity identified by tools such as CD-HIT, BLAST, DIAMOND or MMseqs2.
|
|
18
|
-
This work is an extension of the gene family / pangenome tool developed for the StORF-Reporter publication in NAR (https://doi.org/10.1093/nar/gkad814).
|
|
19
|
-
|
|
20
|
-
## Features
|
|
21
|
-
- **End-to-End**: PyamilySeq can take a directory of GFF+FASTA files, run CD-HIT for clustering and process the results.
|
|
22
|
-
- **Clustering**: Supports input from CD-HIT formatted files as well as CSV and TSV edge lists (-outfmt 6 from BLAST/DIAMOND).
|
|
23
|
-
- **Reclustering**: Allows for the addition of new sequences post-initial clustering.
|
|
24
|
-
- **Output**: Generates a gene 'Roary/Panaroo' formatted presence-absence CSV formatted file for downstream analysis.
|
|
25
|
-
- Align representative sequences using MAFFT.
|
|
26
|
-
- Output concatenated aligned sequences for downstream analysis.
|
|
27
|
-
- Optionally output sequences of identified families.
|
|
28
|
-
|
|
29
|
-
|
|
30
|
-
### Installation
|
|
31
|
-
PyamilySeq requires Python 3.6 or higher. Install using pip:
|
|
32
|
-
|
|
33
|
-
```bash
|
|
34
|
-
pip install PyamilySeq
|
|
35
|
-
```
|
|
36
|
-
|
|
37
|
-
### Examples: Below are two examples of running PyamilySeq in its two main modes.
|
|
38
|
-
#### 'Full Mode': Will conduct clustering of sequences as part of PyamilySeq run
|
|
39
|
-
```bash
|
|
40
|
-
PyamilySeq -run_mode Full -group_mode Species -output_dir ../../test_data/testing -input_type combined -input_dir .../test_data/genomes -name_split _combined.gff3 -pid 0.99 -len_diff 0.99 -clust_tool CD-HIT -gpa True -con True -w 99 -verbose True
|
|
41
|
-
```
|
|
42
|
-
#### 'Partial Mode': Will take the output of a sequence clustering
|
|
43
|
-
```bash
|
|
44
|
-
PyamilySeq -run_mode Partial -group_mode Species -output_dir .../test_data/testing -cluster_file .../test_data/CD-HIT/combined_Ensmbl_pep_CD_90_60.clstr -clust_tool CD-HIT -original_fasta .../test_data/combined_Ensmbl_cds.fasta -gpa True -con True -w 99 -verbose True
|
|
45
|
-
```
|
|
46
|
-
|
|
47
|
-
```bash
|
|
48
|
-
Calculating Groups
|
|
49
|
-
Gene Groups:
|
|
50
|
-
first_core_99: 3103
|
|
51
|
-
first_core_95: 0
|
|
52
|
-
first_core_15: 3217
|
|
53
|
-
first_core_0: 4808
|
|
54
|
-
Total Number of Gene Groups (Including Singletons): 11128
|
|
55
|
-
```
|
|
56
|
-
|
|
57
|
-
|
|
58
|
-
## Usage - Menu
|
|
59
|
-
```
|
|
60
|
-
usage: PyamilySeq.py [-h] -run_mode {Full,Partial} -group_mode {Species,Genus} -clust_tool {CD-HIT} -output_dir OUTPUT_DIR [-input_type {separate,combined}] [-input_dir INPUT_DIR] [-name_split NAME_SPLIT]
|
|
61
|
-
[-pid PIDENT] [-len_diff LEN_DIFF] [-mem CLUSTERING_MEMORY] [-t CLUSTERING_THREADS] [-cluster_file CLUSTER_FILE] [-reclustered RECLUSTERED] [-seq_tag SEQUENCE_TAG]
|
|
62
|
-
[-core_groups CORE_GROUPS] [-genus_groups GENUS_GROUPS] [-w WRITE_FAMILIES] [-con CON_CORE] [-original_fasta ORIGINAL_FASTA] [-gpa GENE_PRESENCE_ABSENCE_OUT] [-verbose {True,False}] [-v]
|
|
63
|
-
|
|
64
|
-
PyamilySeq v0.6.0: PyamilySeq Run Parameters.
|
|
65
|
-
|
|
66
|
-
options:
|
|
67
|
-
-h, --help show this help message and exit
|
|
68
|
-
|
|
69
|
-
Required Arguments:
|
|
70
|
-
-run_mode {Full,Partial}
|
|
71
|
-
Run Mode: Should PyamilySeq be run in "Full" or "Partial" mode?
|
|
72
|
-
-group_mode {Species,Genus}
|
|
73
|
-
Group Mode: Should PyamilySeq be run in "Species" or "Genus" mode?
|
|
74
|
-
-clust_tool {CD-HIT} Clustering tool to use: CD-HIT, DIAMOND, BLAST or MMseqs2.
|
|
75
|
-
-output_dir OUTPUT_DIR
|
|
76
|
-
Directory for all output files.
|
|
77
|
-
|
|
78
|
-
Full-Mode Arguments - Required when "-run_mode Full" is used:
|
|
79
|
-
-input_type {separate,combined}
|
|
80
|
-
Type of input files: 'separate' for separate FASTA and GFF files, 'combined' for GFF files with embedded FASTA sequences.
|
|
81
|
-
-input_dir INPUT_DIR Directory containing GFF/FASTA files.
|
|
82
|
-
-name_split NAME_SPLIT
|
|
83
|
-
substring used to split the filename and extract the genome name ('_combined.gff3' or '.gff').
|
|
84
|
-
-pid PIDENT Default 0.95: Pident threshold for clustering.
|
|
85
|
-
-len_diff LEN_DIFF Default 0.80: Minimum length difference between clustered sequences - (-s) threshold for CD-HIT clustering.
|
|
86
|
-
|
|
87
|
-
Clustering Runtime Arguments - Optional when "-run_mode Full" is used:
|
|
88
|
-
-mem CLUSTERING_MEMORY
|
|
89
|
-
Default 4000: Memory to be allocated for clustering (in MBs).
|
|
90
|
-
-t CLUSTERING_THREADS
|
|
91
|
-
Default 4: Threads to be allocated for clustering.
|
|
92
|
-
|
|
93
|
-
Partial-Mode Arguments - Required when "-run_mode Partial" is used:
|
|
94
|
-
-cluster_file CLUSTER_FILE
|
|
95
|
-
Clustering output file containing CD-HIT, TSV or CSV Edge List
|
|
96
|
-
|
|
97
|
-
Grouping Arguments - Use to fine-tune grouping of genes after clustering:
|
|
98
|
-
-reclustered RECLUSTERED
|
|
99
|
-
Currently only works on Partial Mode: Clustering output file from secondary round of clustering.
|
|
100
|
-
-seq_tag SEQUENCE_TAG
|
|
101
|
-
Default - "StORF": Unique identifier to be used to distinguish the second of two rounds of clustered sequences
|
|
102
|
-
-core_groups CORE_GROUPS
|
|
103
|
-
Default - ('99,95,15'): Gene family groups to use for "Species" mode
|
|
104
|
-
-genus_groups GENUS_GROUPS
|
|
105
|
-
Default - ('1,2,3,4,5,6'): Gene family groups to use for "Genus" mode
|
|
106
|
-
|
|
107
|
-
Output Parameters:
|
|
108
|
-
-w WRITE_FAMILIES Default - No output: Output sequences of identified families (provide levels at which to output "-w 99,95" - Must provide FASTA file with -fasta
|
|
109
|
-
-con CON_CORE Default - No output: Output aligned and concatinated sequences of identified families - used for MSA (provide levels at which to output "-w 99,95" - Must provide FASTA file with -fasta
|
|
110
|
-
-original_fasta ORIGINAL_FASTA
|
|
111
|
-
FASTA file to use in conjunction with "-w" or "-con" when running in Partial Mode.
|
|
112
|
-
-gpa GENE_PRESENCE_ABSENCE_OUT
|
|
113
|
-
Default - False: If selected, a Roary formatted gene_presence_absence.csv will be created - Required for Coinfinder and other downstream tools
|
|
114
|
-
|
|
115
|
-
Misc:
|
|
116
|
-
-verbose {True,False}
|
|
117
|
-
Default - False: Print out runtime messages
|
|
118
|
-
-v Default - False: Print out version number and exit
|
|
119
|
-
|
|
120
|
-
```
|
|
121
|
-
|
|
122
|
-
|
|
123
|
-
## Seq-Combiner: This tool is provided to enable the pre-processing of multiple GFF/FASTA files together ready to be clustered by the user
|
|
124
|
-
### Example:
|
|
125
|
-
```bash
|
|
126
|
-
Seq-Combiner -input_dir .../test_data/genomes -name_split _combined.gff3 -output_dir.../test_data -output_name combine_fasta_seqs.fa -input_type combined
|
|
127
|
-
```
|
|
128
|
-
## Seq-Combiner Menu:
|
|
129
|
-
```bash
|
|
130
|
-
usage: Seq_Combiner.py [-h] -input_dir INPUT_DIR -input_type {separate,combined} -name_split NAME_SPLIT -output_dir OUTPUT_DIR -output_name OUTPUT_FILE
|
|
131
|
-
|
|
132
|
-
Seq-Combiner v0.6.0: Seq-Combiner Run Parameters.
|
|
133
|
-
|
|
134
|
-
options:
|
|
135
|
-
-h, --help show this help message and exit
|
|
136
|
-
|
|
137
|
-
Required Arguments:
|
|
138
|
-
-input_dir INPUT_DIR Directory location where the files are located.
|
|
139
|
-
-input_type {separate,combined}
|
|
140
|
-
Type of input files: 'separate' for separate FASTA and GFF files, 'combined' for GFF files with embedded FASTA sequences.
|
|
141
|
-
-name_split NAME_SPLIT
|
|
142
|
-
substring used to split the filename and extract the genome name ('_combined.gff3' or '.gff').
|
|
143
|
-
-output_dir OUTPUT_DIR
|
|
144
|
-
Directory for all output files.
|
|
145
|
-
-output_name OUTPUT_FILE
|
|
146
|
-
Output file name.
|
|
147
|
-
```
|
|
@@ -1,15 +0,0 @@
|
|
|
1
|
-
PyamilySeq/CD-Hit_StORF-Reporter_Cross-Genera_Builder.py,sha256=UzQ5iOKCNfurxmj1pnkowF11YfWBO5vnBCKxQK6goB8,26538
|
|
2
|
-
PyamilySeq/Constants.py,sha256=LNbn_6kugNL3dl4HFDkq8ZC4cizcXy6MWRz3hP4QJGI,31
|
|
3
|
-
PyamilySeq/PyamilySeq.py,sha256=_Lwk_NYUmOMzoJsCRpiMIx5CFuKYxnYBjZ6NhTdBiCE,13026
|
|
4
|
-
PyamilySeq/PyamilySeq_Genus.py,sha256=dKSM4ZGrOzZPIMkDveyiWD4WqBeYjwAP3g1diR7-0rE,12107
|
|
5
|
-
PyamilySeq/PyamilySeq_Species.py,sha256=urDUDkTu4-5L5rPlaWN2kicXOalqDDVK9lLOrOMdDUU,13219
|
|
6
|
-
PyamilySeq/Seq_Combiner.py,sha256=bHEIR-MZODSGm8n69ZIN-XzEoct5WZ1kH5Xa6uKCu4Y,1972
|
|
7
|
-
PyamilySeq/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
|
|
8
|
-
PyamilySeq/clusterings.py,sha256=Xa0YTvfmI0A8sPdaNHw3j2PVVvGR0JlrebUb_vP2kt8,16183
|
|
9
|
-
PyamilySeq/utils.py,sha256=bC1fpJ8SS14oB4EHvsgZbR0ttn83BiBttlMYeD6IS3g,9863
|
|
10
|
-
PyamilySeq-0.6.0.dist-info/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
|
|
11
|
-
PyamilySeq-0.6.0.dist-info/METADATA,sha256=4O545LBKwlxMTj4q8kCth-IEYJW8Vz4YhaaJ5R3P9EE,7706
|
|
12
|
-
PyamilySeq-0.6.0.dist-info/WHEEL,sha256=Wyh-_nZ0DJYolHNn1_hMa4lM7uDedD_RGVwbmTjyItk,91
|
|
13
|
-
PyamilySeq-0.6.0.dist-info/entry_points.txt,sha256=QtXD1tmnLvRAkIpGWZgXm1lfLH8GGeCwxmgoHZaTp98,102
|
|
14
|
-
PyamilySeq-0.6.0.dist-info/top_level.txt,sha256=J6JhugUQTq4rq96yibAlQu3o4KCM9WuYfqr3w1r119M,11
|
|
15
|
-
PyamilySeq-0.6.0.dist-info/RECORD,,
|
|
File without changes
|
|
File without changes
|
|
File without changes
|