{"id":13704217,"url":"https://github.com/ncbi/SKESA","last_synced_at":"2025-05-05T09:33:28.330Z","repository":{"id":53829916,"uuid":"140623405","full_name":"ncbi/SKESA","owner":"ncbi","description":"SKESA assembler","archived":false,"fork":false,"pushed_at":"2024-01-19T20:14:54.000Z","size":281,"stargazers_count":110,"open_issues_count":16,"forks_count":16,"subscribers_count":17,"default_branch":"master","last_synced_at":"2024-08-03T21:04:45.804Z","etag":null,"topics":[],"latest_commit_sha":null,"homepage":"","language":"C++","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":"other","status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/ncbi.png","metadata":{"files":{"readme":"README.md","changelog":null,"contributing":null,"funding":null,"license":"LICENSE","code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null,"roadmap":null,"authors":null,"dei":null}},"created_at":"2018-07-11T20:13:17.000Z","updated_at":"2024-08-03T00:40:12.000Z","dependencies_parsed_at":"2024-02-03T12:47:30.452Z","dependency_job_id":null,"html_url":"https://github.com/ncbi/SKESA","commit_stats":null,"previous_names":[],"tags_count":7,"template":false,"template_full_name":null,"repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ncbi%2FSKESA","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ncbi%2FSKESA/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ncbi%2FSKESA/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/ncbi%2FSKESA/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/ncbi","download_url":"https://codeload.github.com/ncbi/SKESA/tar.gz/refs/heads/master","host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":224439816,"owners_count":17311531,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2022-07-04T15:15:14.044Z","host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":[],"created_at":"2024-08-02T21:01:05.824Z","updated_at":"2025-05-05T09:33:28.299Z","avatar_url":"https://github.com/ncbi.png","language":"C++","funding_links":[],"categories":["Next Generation Sequencing"],"sub_categories":["Assembly"],"readme":"\nFor questions, please contact\n    Alexandre Souvorov (souvorov@ncbi.nlm.nih.gov)\n    Richa Agarwala (agarwala@ncbi.nlm.nih.gov)\n\n## Citation\n\nAlexandre Souvorov, Richa Agarwala and David J. Lipman. \n**SKESA: strategic k-mer extension for scrupulous assemblies.**\n*Genome Biology* 2018 **19**:153.\n[doi.org/10.1186/s13059-018-1540-z](https://doi.org/10.1186/s13059-018-1540-z)\n\nAlexandre Souvorov, Richa Agarwala.\n**SAUTE: sequence assembly using target enrichment.**\nBMC Bioinformatics **22**, 375 (2021). \n[doi.org/10.1186/s12859-021-04174-9](https://doi.org/10.1186/s12859-021-04174-9)\n\n## Introduction\nThis repositiry is the home for two different short read assemblers and helper programs. These programs are based on de Bruijn graphs and share k-mer related code.\n\n[SKESA](#skesa---strategic-k-mer-extension-for-scrupulous-assemblies) is a de-novo sequence read assembler for microbial genomes. It uses conservative heuristics and is designed to create breaks at repeat regions in the genome. This leads to excellent sequence quality without significantly compromising contiguity. If desired, SKESA contigs could be connected into a GFA graph using [GFA connector](#gfa-connector---create-gfa-graph-starting-from-a-set-of-contigs).\n\n[SAUTE](#saute---sequence-assembly-using-target-enrichment) is a de Bruijn graph based target enriched de-novo assembler designed for assembling genomic and RNA-seq reads sequenced using Illumina. The result is reported as a [GFA graph]( http://gfa-spec.github.io/GFA-spec/) and two nucleotide fasta sequence files for assemblies in the graph.\n\n## Contents\n- [Compilation](#compilation)\n- [SKESA - Strategic K-mer Extension for Scrupulous Assemblies](#skesa---strategic-k-mer-extension-for-scrupulous-assemblies)\n- [SAUTE - Sequence assembly using target enrichment](#saute---sequence-assembly-using-target-enrichment)\n- [GFA connector - create GFA graph starting from a set of contigs](#gfa-connector---create-gfa-graph-starting-from-a-set-of-contigs)\n- [Kmercounter - effectively count k-mers in a read set](#kmercounter---effectively-count-k-mers-in-a-read-set)\n\n## Compilation\n\nDownload current source code\n\n       $ git clone https://github.com/ncbi/SKESA\n\nAlternatively, download last stable release from https://github.com/ncbi/SKESA/releases. Releases also include test data and precompiled binaries. Test data is available in example subdirectory that has commands in file run.test for generating assemblies using the test data.\n\nTo build all executables, do following:\n\n       $ cd SKESA\n\n\nIf you would like to build NGS library for accessing reads from SRA, then do\n\n       $ make\n\nOtherwise, if reading inputs only from files, do\n\n       $ make -f Makefile.nongs\n\nBOOST install is expected by makefiles in the SKESA release. If you do not have BOOST on the system path, please specify BOOST_PATH using a command like\n\n       setenv BOOST_PATH /netopt/ncbi_tools64/boost-1.62.0-ncbi1\n\nbefore running make.\n\nThese make files have been tested with BOOST v 1.72.0 and gcc v 7.3.\n\n# SKESA - Strategic K-mer Extension for Scrupulous Assemblies\nVersion 2.4.0\n\n## Synopsis\n    Running\n       skesa \n    or \n       skesa -h\n    or \n       skesa --help\n    gives information about options and produces the following:\n\n    General options:\n      -h [ --help ]                 Produce help message\n      -v [ --version ]              Print version\n      --cores arg (=0)              Number of cores to use (default all) [integer]\n      --memory arg (=32)            Memory available (GB, only for sorted counter) \n                                    [integer]\n      --hash_count                  Use hash counter [flag]\n      --estimated_kmers arg (=100)  Estimated number of distinct kmers for bloom \n                                    filter (M, only for hash counter) [integer]\n      --skip_bloom_filter           Don't do bloom filter; use --estimated_kmers as\n                                    the hash table size (only for hash counter) \n                                    [flag]\n    \n    Input/output options: at least one input providing reads for assembly must be specified:\n      --reads arg                   Input fasta/fastq file(s) for reads (could be \n                                    used multiple times for different runs, could \n                                    be gzipped) [string]\n      --use_paired_ends             Indicates that single (not comma separated) \n                                    fasta/fastq files contain paired reads [flag]\n      --sra_run arg                 Input sra run accession (could be used multiple\n                                    times for different runs) [string]\n      --contigs_out arg             Output file for contigs (stdout if not \n                                    specified) [string]\n    \n    Assembly options:\n      --kmer arg (=21)              Minimal kmer length for assembly [integer]\n      --min_count arg               Minimal count for kmers retained for comparing \n                                    alternate choices [integer]\n      --max_kmer_count arg          Minimum acceptable average count for estimating\n                                    the maximal kmer length in reads [integer]\n      --vector_percent arg (=0.05)  Percentage of reads containing 19-mer for the\n                                    19-mer to be considered a vector (1. disables) [float (0,1]]\n      --insert_size arg             Expected insert size for paired reads (if not \n                                    provided, it will be estimated) [integer]\n      --steps arg (=11)             Number of assembly iterations from minimal to \n                                    maximal kmer length in reads [integer]\n      --fraction arg (=0.1)         Maximum noise to signal ratio acceptable for \n                                    extension [float [0,1)]\n      --max_snp_len arg (=150)      Maximal snp length [integer]\n      --min_contig arg (=200)       Minimal contig length reported in output \n                                    [integer]\n      --allow_snps                  Allow additional step for snp discovery [flag]\n    \n    Debugging options:\n      --force_single_ends           Don't use paired-end information [flag]\n      --seeds arg                   Input file with seeds [string]\n      --all arg                     Output fasta for each iteration [string]\n      --dbg_out arg                 Output kmer file [string]\n      --hist arg                    File for histogram [string]\n      --connected_reads arg         File for connected paired reads [string]\n\n    Note that --sra_run option is not available if SKESA is built using Makefile.nongs\n    \n\n## Short description\nSKESA is a de-novo sequence read assembler for microbial genomes based on DeBruijn graphs. It uses conservative heuristics and is designed to create breaks at repeat \nregions in the genome. This leads to excellent sequence quality. Using k-mers longer than mate length and up to insert size also allows SKESA to attain good contiguity \nas determined by the N50 statistic. It is a multi-threaded application that scales well with the number of processors. For different runs with the same inputs, \nincluding the order of reads, the order and orientation of contigs in the output is deterministic. \n\nSKESA can process read information by accessing reads from SRA (option `--sra_run`) or from files (option `--reads`) format. Any combination of input streams is \nallowed. Files could be in fasta or fastq format and gzipped, which is recognized automatically. \n    \nWhen accessing reads from SRA, SKESA automatically determines if the read set consists of paired-end or single-end reads. For fasta/fastq input of paired reads \nwith separate files for each mate, filenames separated by a comma for first mate followed by the second mate are listed and in this case, the order of reads \nis expected to be same in files for both mates. Alternatively, a single file with both mates could be specified. In this case the reads are expected to be \ninterleaved with first mate followed by the second, and the option `--use_paired_ends` must be used.\n    \nA limitation of the current release is that in case multiple streams of paired reads are provided, it is assumed that all streams have the same insert size. \nUser can explicitly specify expected insert size for the reads (option --insert_size). Otherwise, a sample of input reads is used to estimate the expected \ninsert size. This sampling may lead to very small differences in assembly of the same read set if the order of reads is different and selected sample gives \na difference in expected insert size.\n    \nTwo additional options users may wish to specify depending on the resources available to them are as follows:\n1. the number of cores (option `--cores`) and\n2. total amount of memory in Gb (option `--memory`)\n\nRemaining options are for debugging or modifying algorithm parameters.\n\nOutput of assembly is contigs in fasta format. The definition line for contig has format Contig_N_cnt where N is consecutive integers starting from one \nfor numbering the contigs and cnt is the average count of k-mers in the contig using minimal k-mer length used in the assembly. Contigs are ordered lexicographically.\n\nLimitations:\n1. SKESA is designed for haploid genomes. If it is used for diploid genomes or RNAseq reads, it should create breaks at all heterozygous sites in the genome \nand sites for alternative splicing, respectively. The allow_snps option can be used to make some joins at well separated heterozygous sites.\n2. SKESA is designed for ILLUMINA reads that do not have systematic homopolymer errors. The assembly for reads that do not have properties similar to ILLUMINA \nis likely to be quite fragmented.\n3. Forward-reverse orientation for paired reads is assumed at this time. If this is not true, steps using paired reads are unlikely to change/improve the assembly.\n4. Requesting expected insert size to be estimated using a sample is guaranteed to give the same result, including the order of contigs, for the same order \nof reads but may give very small differences if read order is changed and insert size estimate is different.\n\n\n## SKESA usage examples\n\n     In all the examples below, we are providing 4 cores and have 48 Gb of memory.\n\n     Example of an assembly that directly accesses SRA for an unpaired read set SRR867211 is:\n\n       $ skesa --sra_run SRR867211 --cores 4 --memory 48 \u003e SRR867211.skesa.fa\n\n     Example of an assembly that directly accesses SRA for a paired read set SRR1960353 is:\n\n       $ skesa --sra_run SRR1960353 --cores 4 --memory 48 \u003e SRR1960353.skesa.fa\n\n     Example of an assembly that uses separate fastq files for each mate of SRR1703350 is:\n\n       $ skesa --reads SRR1703350_1.fq,SRR1703350_2.fq --cores 4 --memory 48 \u003e SRR1703350.skesa.fa\n\n     Example of an assembly that uses interleaved mates for SRR1703350 as fastq input is:\n\n       $ skesa --reads SRR1703350.fq --use_paired_ends --cores 4 --memory 48 \u003e SRR1703350.skesa.fa\n\n     Example of an assembly that uses reads from SRA for SRR1695624 and gzipped fasta for SRR1745628 is:\n\n       $ skesa --sra_run SRR1695624 --reads SRR1745628.fa.gz --use_paired_ends --cores 4 --memory 48 \u003e SAMN03218571.skesa.fa\n\n     Example of the same assembly as above done with both runs accessed from SRA is:\n\n       $ skesa --sra_run SRR1695624 --sra_run SRR1745628 --cores 4 --memory 48 \u003e SAMN03218571.skesa.fa\n       \n# SAUTE - Sequence Assembly Using Target Enrichment\nVersion 1.3.0\n## Short description\nSAUTE is a de Bruijn graph based target enriched de-novo assembler designed for assembling genomic and\nRNA-seq reads sequenced using Illumina. It aligns user-provided reference sequences in the target to the de Bruijn graph\nfor reads and finds all paths with good similarity to each reference. Assembled paths are filtered by SAUTE using full length reads and read pairs to remove possibly chimeric connections. The filtering steps could be partially or completely disabled using options `--no_filter_by_reads` and `--no_filter_by_pairs`.\n\nSimilar to [SKESA](#skesa-usage-examples), the user can provide reads using options `--sra_run` or `--reads`. The fasta file for reference sequences are provided using the option `--targets`. The results are reported as a [GFA graph]( http://gfa-spec.github.io/GFA-spec/) (optin `--gfa`) and fasta files (options `--all_variants` and `--selected_variants`). A GFA graph is a convenient way for viewing and analyzing all variants using [Bandage](https://rrwick.github.io/Bandage/).\n\nSAUTE has two flavors - one uses nucleotide reference sequences and the other protein sequences. In both cases, the output is nucleotide assemblies. Depending on the type of reference sequences, do one of the following:\n- Nucleotide references\n  saute --reads mates1.fa,mates2.fa --targets nuc_references.fa --gfa graph.gfa --all_variants assembly.fa\n- Protein references\n  saute_prot --reads mates1.fa,mates2.fa --targets prot_references.fa --gfa graph.gfa --all_variants assembly.fa --genetic_code 1\n  The genetic code is one of the [NCBI genetic codes](https://www.ncbi.nlm.nih.gov/Taxonomy/Utils/wprintgc.cgi). The protein version translates on the fly the segments of de Bruijn graph to the protein space. Any indel changing the frame will corrupt the translated protein and will stop the assembly. There is an option `--allow_frameshifts` which will resume the correct translation after such an indel. It could be useful, for example, for genes with ribosomal slippage.  \n\n## Command line options\n### Options shared by saute and saute_prot\n    General options:\n    -h [ --help ]                    Produce help message\n    -v [ --version ]                 Print version\n    --cores (=0)                     Number of cores to use (default all) [integer]\n    --estimated_kmers (=1000)        Estimated number of distinct kmers for bloom filter (millions) [integer]\n                                     To avoid expensive rehashing of the de Bruijn graph, SAUTE uses a light-weight counting Bloom filter to find \n                                     the number of different kmers satisfying the --min_count criterion. Parameter --estimated_kmers provides the \n                                     number (in millions) of different kmers with any count present in the reads. Overestimating will result in a \n                                     small waste of memory. Underestimating will trigger iterative Bloom filter recalculation until a good estimate\n                                     is found.\n\n    Input/output options : target file, at least one input for reads and output for gfa must be specified:\n    --gfa                            Output file for GFA graph [string]\n    --all_variants                   Output file for sequences with all variant combinations [string]\n    --max_variants (=1000)           Restricts the number of variants reported in --all_variants per graph [integer]\n    --selected_variants              Output file for selected sequences representing all graph elements [string]\n    --targets                        Input file with reference sequences [string]\n    --sra_run                        Input sra run accession (could be used multiple times for different runs) [string]\n    --reads                          Input fasta/fastq file(s) for reads (could be used multiple times for different runs, could be gzipped) [string]\n    --use_paired_ends                Indicates that single (not comma separated) fasta/fastq files contain paired reads [flag]\n\n    Assembly options:\n    --vector_percent arg (=0.05)  Percentage of reads containing 19-mer for the 19-mer to be considered a vector (1. disables) [float (0,1]]\n                                     Before doing any assembly SAUTE removes vectors from the reads. It finds and clips read after all 19-mers found in a relatively high number of reads.\n\n    --min_count (=2)                 Minimal count for kmers retained in graph [integer]\n                                     This option keeps noise kmers out of the de Bruijn graphs used for assembling. By default, a kmer has to be seen in \n                                     at least 2 reads to be included. It is allowed to use --min_count 1, and it could help in assembling low \n                                     coverage spots. The drawback of this setting is that all read errors satisfying --fraction criterion will be \n                                     recognized as valid SNPs. To alleviate this, the program will remove all sections supported by only 1 read \n                                     if they were not necessary for a connection. \n    --fraction (=0.05)               Maximum noise to signal ratio acceptable for extension [float [0,1)]\n                                     All fork branches with relative abundance less than this are ignored.\n    --kmer                           Primary kmer length for assembly (default automatic) [integer]\n    --secondary_kmer                 Shorter kmer length for low coverage spots (default automatic) [integer]\n                                     For assembling SAUTE uses two de Bruijn graphs with different kmer lengths. It uses the longer kmer most of the \n                                     time and switches to the shorter kmer to assemble the low coverage spots. Both kmers must be odd numbers. In case \n                                     of protein references, both kmers must be divisible by 3. By default, the program will use the closest valid values\n                                     not exceeding half and a fifth of the read length.\n    --secondary_kmer_threshold (=1)  Coverage threshold for using shorter kmer [integer]\n                                     This parameter defines what is considered to be low coverage for the secondary kmer use. With default, these are spots with no \n                                     possible extension or extension supported by only 1 read (could happen only with --min-count 1). Specifying \n                                     a higher number may help with detecting some low coverage forks. With more work shifted toward the shorter \n                                     kmer, there is a higher chance that the program will enter a highly repetitive area and will create a very \n                                     complex graph.\n\n    --word                           Word size for seeds [integer \u003c= 16]\n                                     To start assembling, SAUTE scans the main de Bruijn graph and finds seed kmers with high similarity to the reference.\n                                     These kmers are used as starting points for the assembly. To trigger a kmer comparison to the reference, one of its ends\n                                     has to have an exact match of --word length. For a protein reference, the word must be divisible by 3. Defaults are 8 for \n                                     nucleotide references and 12 for protein references. Latter translates to 4 aa match.\n\n    --kmer_complexity (=2000)        Hard mask reference areas with high number of variant seed possibilities (0 disables masking) [integer]\n                                     Repetitive regions could result in very complex output graph and excessive calculation time. All reference areas\n                                     for which SAUTE found in each position more different seeds than this parameter are internally hard masked and\n                                     will not be assembled. \n    --max_fork_density (=0.1)        Maximal fork density averaged over --buf_length before abandoning assembling (0 disables) [float]\n    --buf_length (=200)              Buffer length for fork density [integer]\n                                     If assembling enters an unmasked repetitive area, the number of forks and paths to follow may become very high. \n                                     SAUTE calculates the average fork density for the last --buf_length long portion of the assembly. If the fork \n                                     density is above the --max_fork_density threshold, assembling for the current seed is stopped and the program \n                                     moves to the next seed.\n\n    --target_coverage (=0.5)         Keep a path in output if it has alignment to the reference that is at least target_coverage*(reference length) long [float (0,1]]\n    --min_hit_len                    If a path is shorter than target_coverage*(reference length), use this length threshold for keeping paths (optional) [integer]\n\n    --extend_ends                    Unambiguously extend graph ends using de-novo assembly [flag]\n    --protect_reference_ends         Near the reference ends, don't check if reads support some minimal extension of the fork's branches [flag]\n                                     Each fork is analysed using the algorithm described in the SKESA paper, and less supported branches are deleted. \n                                     There is one check which could be detrimental to assembling RNAseq -- each valid branch has to be extendable by ~100bp.\n                                     If this is close to a transcript 5' or 3' end, extension is expected to fail. This option disables the extension\n                                     check at the reference ends.\n    --keep_subgraphs                 Don't remove redundant subgraphs [flag]\n                                     If the reference file contains similar references, SAUTE may create multiple graphs which are either identical to\n                                     or are sub-graphs of some other graphs. By default, the program will remove redundant\n                                     graphs and sub-graphs. This option disables sub-graph removal.\n    --use_ambiguous_na               Use ambiguous nucleotide codes for SNPs [flag]\n                                     If the option is used, SAUTE will replace SNPs with ambiguous nucleotide codes. This will simplify the graph\n                                     and will reduce the number of sequence variants. The counts for individual SNPs will be lost.\n\n    --gap_open                       Penalty for gap opening [integer] (default of 5 for nucleotide references and 11 for protein references)\n    --gap_extend                     Penalty for gap extension [integer] (default of 3 for nucleotide references and 1 for protein references)\n    --drop_off (=300)                Maximal decrease of score before abandoning assembly path [integer]\n                                     SAUTE uses a Smith-Waterman type algorithm for aligning references to the de Bruijn graph. For nucleotide references,\n                                     the usual match bonus and mismatch penalty are used. They could be changed from the input (see below Options specific \n                                     for nucleotide targets). For protein references BLOSUM62 substitution matrix is used. In both cases, penalties for gaps \n                                     are specified using --gap_open and --gap_extend parameters. The --drop_off parameter determines how far from the \n                                     reference the assembly may go before the extension is stopped, and the assembly is trimmed to the aligned part.\n                                     Larger gap open penalty is used for opening frameshifts (see below).\n\n    Graph cleaning options:\n    --no_filter_by_reads             Don't use full length reads for variants filtering [flag]\n    --no_filter_by_pairs             Don't use mate pairs for variants filtering [flag]\n                                     After assembling, SAUTE will use full length reads and read pairs to remove chimeric connections from the graph. \n                                     Above options will partially or completely disable these steps.\n    --max_path (=1000)               Maximal number of path extensions allowed for a single filtering check [integer]\n                                     For each fork in graph, SAUTE expands all sequences starting at this fork using either the read length or the insert size.\n                                     If the number of expanded sequences exceed this parameter value, the extension length is reduced.\n    --not_aligned_len (=10)          Not aligned read length for break count [integer]\n    --not_aligned_count (=3)         Number of not aligned reads to make a break [integer]\n    --aligned_count (=2)             Number of aligned reads to confirm a connection [integer]\n                                     For each expanded sequence, SAUTE counts how many reads or pairs confirm this path and how many contradict it. Based on\n                                     these counts, the path is either kept or discarded.\n \n    --remove_homopolymer_indels     Remove homopolymer indels [flag]\n    --homopolymer_len (=4)          Minimal length of homopolymer [integer]\n    --homopolymer_ratio (=0.33)     Coverage ratio threshold for removing homopolymer indels [float]\n\n\n### Options specific for nucleotide targets    \n    --match (=2)                    Bonus for match [integer]\n    --mismatch (=3)                 Penalty for mismatch [integer]\n\n### Options specific for protein targets                                         \n    --genetic_code                   Genetic code [integer]\n                                     Accepted genetic codes:\n                                     1 Standard Code\n                                     2 Vertebrate Mitochondrial Code\n                                     3 Yeast Mitochondrial Code\n                                     4 Mold, Protozoan, and Coelenterate Mitochondrial Code and the Mycoplasma/Spiroplasma Code\n                                     5 Invertebrate Mitochondrial Code\n                                     6 Ciliate, Dasycladacean and Hexamita Nuclear Code\n                                     9 Echinoderm and Flatworm Mitochondrial Code\n                                    10 Euplotid Nuclear Code\n                                    11 Bacterial, Archaeal and Plant Plastid Code\n                                    12 Alternative Yeast Nuclear Code\n                                    13 Ascidian Mitochondrial Code\n                                    14 Alternative Flatworm Mitochondrial Code\n                                    16 Chlorophycean Mitochondrial Code\n                                    21 Trematode Mitochondrial Code\n                                    22 Scenedesmus obliquus Mitochondrial Code\n                                    23 Thraustochytrium Mitochondrial Code\n                                    24 Pterobranchia Mitochondrial Code\n                                    25 Candidate Division SR1 and Gracilibacteria Code\n                                    26 Pachysolen tannophilus Nuclear Code\n                                    27 Karyorelict Nuclear Code\n                                    28 Condylostoma Nuclear Code\n                                    29 Mesodinium Nuclear Code\n                                    30 Peritrich Nuclear Code\n                                    31 Blastocrithidia Nuclear Code\n                                    33 Cephalodiscidae Mitochondrial UAA-Tyr Code\n\n    --allow_frameshifts              Allow frameshifts [flag]\n    --fs_open (=25)                  Penalty for opening a gap that results in a frameshift [integer]\n    --na_targets                     References are in-frame nucleotide sequences that will be translated before use [flag]\n\n\n# GFA connector - create GFA graph starting from a set of contigs\nVersion 1.1.0\n## Synopsis\n    Running\n       gfa_connector \n    or \n       gfa_connector -h\n    or \n       gfa_connector --help\n    gives information about options and produces the following:\n\n    General options:\n      -h [ --help ]                        Produce help message\n      -v [ --version ]                     Print version\n      --cores arg (=0)                     Number of cores to use (default all) \n                                           [integer]\n      --estimated_kmers arg (=100)         Estimated number of distinct kmers for \n                                           bloom filter (millions) [integer]\n      --skip_bloom_filter                  Don't do bloom filter; use \n                                           --estimated_kmers as the hash table size\n                                           [flag]\n\n    Input/output options:\n      --sra_run arg                        Input sra run accession (could be used \n                                           multiple times for different runs) \n                                           [string]\n      --reads arg                          Input fasta/fastq file(s) (could be used\n                                           multiple times for different runs, could\n                                           be gzipped) [string]\n      --use_paired_ends                    Indicates that single (not comma \n                                           separated) fasta/fastq files contain \n                                           paired reads [flag]\n      --contigs arg                        Input file with contigs [string]\n      --gfa arg                            GFA graph output (stdout if not \n                                           specified) [string]\n      --dbg arg                            Input de Bruijn graph (optional) \n                                           [string]\n      --csv arg                            CSV file ouput (optional) [string]\n      --contig_color arg (=Purple)         Color for contigs [string]\n      --connector_color arg (=Lime Green)  Color for connectors [string]\n\n    Assembly options:\n      --kmer arg (=41)                     Kmer length for assembly [integer]\n      --min_count arg (=2)                 Minimal count for kmers retained for \n                                           comparing alternate choices [integer]\n      --vector_percent arg (=0.05)         Percentage of reads containing 19-mer for the\n                                           19-mer to be considered a vector (1. disables) [float (0,1]]\n      --fraction arg (=0.1)                Threshold for extension\n      --entropy arg (=0.51)                Minimal entropy for a seed kmer [float]\n      --ext_len arg (=2000)                Maximal length for extension\n\n    Graph cleaning options:\n      --not_aligned_len arg (=10)          Not aligned read length for break count\n      --not_aligned_count arg (=3)         Number of not aligned reads to make a \n                                           break\n      --aligned_count arg (=2)             Number of aligned reads to confirm a \n                                           connection\n      --max_path (=1000)                   Maximal number of path extensions allowed for a single filtering check [integer]\n      --no_filter_by_reads                 Don't use full length reads for variants\n                                           filtering [flag]\n      --no_filter_by_pairs                 Don't use mate pairs for variants \n                                           filtering [flag]\n\n    Note that --sra_run option is not available if gfa_connector is built using Makefile.nongs\n\n## Short description\nGFA_connector is a SKESA companion program that takes contigs as input, finds pairs of contigs connected by reads or read pairs,\nand produces a [GFA graph]( http://gfa-spec.github.io/GFA-spec/). The graph can be viewed and analyzed in [Bandage](https://rrwick.github.io/Bandage/).\n\nGFA_connector constructs a de Bruijn graph (option `--kmer`) and starts extending from each end of the input contigs. It continues each \npossible path up to the maximal length for extension (option `--ext_len`) or until it is connected to another contig. In the \nlatter case, two contigs are connected by an edge. At the end of process, all paths not connected to contigs are discarded. This process creates a \nnumber of graphs each connecting two or more contigs. \n\nBefore the input contigs and connectors are combined into a single GFA graph for output, the connector graphs are analyzed using the reads. First the input \nreads are trimmed to the longest stretches compatible with the de Bruijn graph, and then they are aligned to the connectors. Chimeric connections which do \nnot agree with read sequences or paired ends are removed from the graphs. This step could be partially or completely disabled using \noptions `--no_filter_by_reads` and `--no_filter_by_pairs`.\n\nTo run GFA_connector, one has to provide reads using options `--sra_run` or `--reads` (see SKESA description) and contigs using the option `--contigs`.\n\n## Usage example\nThe following commands will produce a GFA graph representing a Klebsiella pneumoniae KPR0928 genome (NZ_CP008831.1), two plasmids (NZ_CP008832.1 \nand NZ_CP008833.1) and a few contigs which the program was not able to connect\n```\nskesa --sra_run SRR1510963 --hash_count --contigs_out SRR1510963_skesa.fa \u003e\u0026 SRR1510963_skesa.log\ngfa_connector --sra_run SRR1510963 --contigs SRR1510963_skesa.fa --gfa SRR1510963.gfa --csv SRR1510963.csv --kmer 75 \u003e\u0026 SRR1510963_gfa.log\n```\n# Kmercounter - effectively count k-mers in a read set\nVersion 2.1.0\n## Synopsis\n    Running\n       kmercounter \n    or \n       kmercounter -h\n    or \n       kmercounter --help\n    gives information about options and produces the following:\n```\nOptions:\n  -h [ --help ]                Produce help message\n  -v [ --version ]             Print version\n  --sra_run arg                Input sra run accession (could be used multiple \n                               times for different runs) [string]\n  --reads arg                  Input fasta/fastq file(s) (could be used \n                               multiple times for different runs) [string]\n  --kmer arg (=21)             Kmer length [integer]\n  --min_count arg (=2)         Minimal count for kmers retained [integer]\n  --vector_percent arg (=0.05) Percentage of reads containing 19-mer for the\n                               19-mer to be considered a vector (1. disables) [float (0,1]]\n  --estimated_kmers arg (=100) Estimated number of distinct kmers for bloom \n                               filter (millions) for hash count [integer]\n  --skip_bloom_filter          Don't do bloom filter; use --estimated_kmers as \n                               the hash table size for hash count [flag]\n  --dbg_out arg                de Bruijn graph output\n  --text_out arg               Text kmer output\n  --hist arg                   File for histogram [string]\n  --cores arg (=0)             Number of cores to use (default all) [integer]\n  \n     Note that --sra_run option is not available if kmercounter is built using Makefile.nongs\n```\n## Short description\nKmercounter creates and outputs a de Bruijn graph (option `--dbg_out`) for a specified kmer size (option `--kmer`). This file could be used as an optional \ninput for GFA_connector in which case the kmer counting step will be skipped. \n\nAll accepted kmers and their counts could be output in a text form (option `--text_out`). The tab delimited output file contains:\n- Kmer sequence\n- Number of reads in which this kmer was seen\n- Number of reads in which this kmer was seen in positive direction\n\nThe output order is random and could be different for different runs.\n\nWith option `--hist` one can get a histogram for the k-mer counts.\n\n## Usage example\n```\nkmercounter --sra_run SRR1510963 --kmer 75 --dbg_out SRR1510963_k75.hdbg --text_out SRR1510963_k75.kmers --hist SRR1510963_k75.hist \u003e\u0026 SRR1510963_counter.log\n```\n","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fncbi%2FSKESA","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fncbi%2FSKESA","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fncbi%2FSKESA/lists"}