trapical: trapical rapidly aligns Pulsenet isolate core alleles

installation

Run the following commands to install trapical

git clone https://github.com/ncezid-biome/trapical.git
conda env create -f trapical/environment.yml
conda activate trapical

usage

usage:
    python3 trapical.py [-ifaocnvh]

required arguments:
    -i, --in            [directory] the directory containing the Pulsenet2.0 WGMLST results

optional arguments:
    -f, --fna_dir       [directory] the directory where locus fasta files should be written (default: ./loci_fastas)
    -a, --aln_dir       [directory] the directory where aligned locus files should be written (default: ./loci_alignments)
    -o, --out           [file] the name of the output alignment file (default: ./msa.aln)
    -c, --count         [file] the name of the file for conserved character counts (default: ./conserved_counts.csv)
    -n, --num_threads   [int] the number of processors for parallel processing (default: 1)
    -v, --version       prints the version
    -h, --help          prints this meessage
    --force             if specified, overwrite existing files and directories

workflow

flowchart TB
    %% variables
    input(["Pulsenet2.0 wgMLST results directory"])
    import1["import allele sequences"]
    import2["import core allele hashes"]
    count1["count conserved characters
            in loci with one allele"]
    write1[/"write a fasta for
             each core locus"/]
    align["align fastas with >1
           allele with clustalo"]
    extract["extract variable sites
             from each alignment"]
    concat["concatenate variable
            sites for each isolate"]
    count2["count conserved characters
           in each alignment"]
    write2[/"write concatenated
             alignment to file"/]
    write3[/"write conserved 
            character counts to file"/]
    
    %% workflow
    input --> import1
    input --> import2
    import1 --> write1
    import2 --> write1
    import1 --> count1
    write1 --> align
    align --> extract
    align --> count2
    extract --> concat
    concat --> write2
    count1 --> write3
    count2 --> write3

Loading

Notices

Public Domain Notice

This repository constitutes a work of the United States Government and is not subject to domestic copyright protection under 17 USC § 105. This repository is in the public domain within the United States, and copyright and related rights in the work worldwide are waived through the CC0 1.0 Universal public domain dedication. All contributions to this repository will be released under the CC0 dedication. By submitting a pull request you are agreeing to comply with this waiver of copyright interest.

License Standard Notice

Privacy Notice

This repository contains only non-sensitive, publicly available data and information. All material and community participation is covered by the Disclaimer and Code of Conduct. For more information about CDC's privacy policy, please visit http://www.cdc.gov/other/privacy.html.

Contributing Notice

Anyone is encouraged to contribute to the repository by forking and submitting a pull request. (If you are new to GitHub, you might start with a basic tutorial.) By contributing to this project, you grant a world-wide, royalty-free, perpetual, irrevocable, non-exclusive, transferable license to all users under the terms of the Apache Software License v2 or later.

All comments, messages, pull requests, and other submissions received through CDC including this GitHub page may be subject to applicable federal law, including but not limited to the Federal Records Act, and may be archived. Learn more at http://www.cdc.gov/other/privacy.html.

Records Management Notice

This repository is not a source of government records, but is a copy to increase collaboration and collaborative potential. All government records will be published through the CDC web site.