HIV-Isoform-Checker


NameHIV-Isoform-Checker JSON
Version 1.1.9 PyPI version JSON
download
home_pagehttps://github.com/JessicaA2019/HIV_Isoform_Checker
SummaryFilters .gtf file of suspected HIV isoforms and confirms the isoform identities
upload_time2024-11-07 22:47:15
maintainerNone
docs_urlNone
authorJessica Lauren ALbert
requires_python>=3
licenseMIT License Copyright (c) 2023 Jessica Lauren Albert Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions: The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software. THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
keywords hiv isoforms splicing ont sequencing
VCS
bugtrack_url
requirements No requirements were recorded.
Travis-CI No Travis.
coveralls test coverage No coveralls.
            # HIV Isoform Checker
### Overview
This package is a module in the [APHIX pipeline][APHIX]. HIV_Isoform_Checker takes a .gtf file of preliminarily filtered HIV transcripts and filters them to include only correctly assigned transcripts using the following filters. It also calculates the usage counts and percentages for each isoform type, donor site, acceptor site, and pairwise splice site combination.
- FILTER 1: only include class codes =, J, and m
- FILTER 2: only include samples with end values >= min_end_bp and
           start values <= max_start_bp
- FILTER 3: get rid of any samples with read errors/small gaps
- FILTER 4: keep only correct Env samples
- FILTER 5: keep only correct Nef samples(long samples added to possible_misassigned)  
- FILTER 6: keep only correct Rev samples(long samples added to possible_misassigned)  
- FILTER 7: keep only correct Tat samples(long samples added to possible_misassigned)   
- FILTER 8: keep only correct Vif samples
- FILTER 9: keep only correct Vpr samples
- FILTER 10:  check possible_misassigned for fully spliced then partial splice compatibility (tat -> rev -> nef - > vif -> vpr -> unspliced_tat -> env)

> Note: This code currently relies on a very specific setup of the gtf file to work properly. The note must be in the order designated in the [].
transcript entry = ref genome, analysis_pathway, transcript, start, end, ".", "+", ".", [transcript id; gene id; gene_name; xloc; ref_gene_id; contained_in; cmp_ref; class_code; tss_id]
>exon entry = ref genome, analysis_pathway, exon, start, end, ".", "+", ".", [transcript id; gene id; exon number]

******************
## Installation
### Quick Install
```bash
pip install HIV_Isoform_Checker
```
### Manual Install
```bash
git clone https://github.com/JessicaA2019/HIV_Isoform_Checker.git 
cd HIV_Isoform_Checker
python setup.py install
```
### Dependencies
Some dependencies currently do not install with Chimera_Buster. Please pip install the following dependencies:
* mappy
* regex
* argparse

### Testing
To test if the installation was successful run the following command with the test files provided:
```bash
HIV_Isoform_Checker HIC_test.gtf HIV_Isoform_Checker_test NL43.fa
 ```
 ******************
## Usage
```bash
HIV Isoform Checker [options] input_file_name output_file_prefix ref_file
```
### Inputs
To run the pipeline the following input files are required:
| Input | Description |
| ------ | ------ |
|input_file_name  |    Designates .gtf to be filtered. This is required.|
|output_file_prefix  |  Designates output file prefix. This is required.|
| ref_file | The location of the reference genome fasta file.|

The following input files are optional:
| Arguement | Function |
| ------ | ------ |
|-h, --help |  Prints help message to terminal. |
|-g value, --gap value |  Sets gap tolerance. If a read has a gap longer than this value that is not associated with known splice sites, the read will be filtered out as a sequencing error. The default value is 15. |
|-a value, --startBP value |  Sets maximum starting bp in order to insure the presence of the first exon and by proxy a full read. Default is 700. |
|-z value, --endBP value |  Sets minimum ending bp in order to insure the presence of the last exon and by proxy a full read. Default is 9500. |
|-l value, --lengthFS value |  Sets maximum fully spliced transcript length. If a read is longer than this value, then it will be considered partial spliced. This is based off HIV-1 references without any major insertions or deletions. If your sequences have larger deletions (ex: deltaENV or deltaMA) or insertions (ex: multiple markers such as eGFP), this value may need to be adjusted. Default is 2500. |
|-n value, --NCE value | When set to True, the output csv file will have y/n columns for the precence of NCEs. Default is False.|

### Outputs
 The main output files created by the pipeline are:
| Output | Description |
|--------|-------------|
| {output_file_prefix}.csv | A csv containing the sample ID, cluster size, isoform ID, starting and ending base pair, gffcompare class code, exon boundaries, coverage length, size normalized counts, precence of NCEs (optional), and any additional notes for each sample. |
| {output_file_prefix}_altered.txt | A list of all sample IDs with altered exons to fill in small gaps. |
| {output_file_prefix}_pass.txt | A list of all sample IDs that passed all filters. |
| {output_file_prefix}_fail.txt | A list of all sample IDs that failed one or more filters. |
| {output_file_prefix}_isoform_counts.csv | A csv of the usage counts and percentages for each isoform type. |
| {output_file_prefix}_splice_site_usage.csv | A csv of the the usage counts and percentages for each donor site, acceptor site, and pairwise splice site combination. |
| {output_file_prefix}_ref_coordinates.txt | A list of splice site and CDS coordinates for given reference. |
| {output_file_prefix}.log | A log of the arguments given to the script. |

**************************
## Help
For issues or bugs, please report them on the [issues page][issues]. 

## License

MIT - Copyright (c) 2024 Jessica Lauren Albert


[//]: # (These are reference links used in the body of this note and get stripped out when the markdown processor does its job. There is no need to format nicely because it shouldn't be seen. Thanks SO - http://stackoverflow.com/questions/4823468/store-comments-in-markdown-syntax)

   [APHIX]: <https://github.com/JessicaA2019/APHIX>
   [issues]: <https://github.com/JessicaA2019/HIV_Isoform_Checker/issues>

            

Raw data

            {
    "_id": null,
    "home_page": "https://github.com/JessicaA2019/HIV_Isoform_Checker",
    "name": "HIV-Isoform-Checker",
    "maintainer": null,
    "docs_url": null,
    "requires_python": ">=3",
    "maintainer_email": null,
    "keywords": "HIV, Isoforms, splicing, ONT sequencing",
    "author": "Jessica Lauren ALbert",
    "author_email": "Jessica Lauren Albert <jessica.albert2001@gmail.com>",
    "download_url": "https://files.pythonhosted.org/packages/dd/42/98cacbf2d30bf703687a50873de9bc12870b4020421362275ee2a1f3604a/hiv_isoform_checker-1.1.9.tar.gz",
    "platform": null,
    "description": "# HIV Isoform Checker\n### Overview\nThis package is a module in the [APHIX pipeline][APHIX]. HIV_Isoform_Checker takes a .gtf file of preliminarily filtered HIV transcripts and filters them to include only correctly assigned transcripts using the following filters. It also calculates the usage counts and percentages for each isoform type, donor site, acceptor site, and pairwise splice site combination.\n- FILTER 1: only include class codes =, J, and m\n- FILTER 2: only include samples with end values >= min_end_bp and\n           start values <= max_start_bp\n- FILTER 3: get rid of any samples with read errors/small gaps\n- FILTER 4: keep only correct Env samples\n- FILTER 5: keep only correct Nef samples(long samples added to possible_misassigned)  \n- FILTER 6: keep only correct Rev samples(long samples added to possible_misassigned)  \n- FILTER 7: keep only correct Tat samples(long samples added to possible_misassigned)   \n- FILTER 8: keep only correct Vif samples\n- FILTER 9: keep only correct Vpr samples\n- FILTER 10:  check possible_misassigned for fully spliced then partial splice compatibility (tat -> rev -> nef - > vif -> vpr -> unspliced_tat -> env)\n\n> Note: This code currently relies on a very specific setup of the gtf file to work properly. The note must be in the order designated in the [].\ntranscript entry = ref genome, analysis_pathway, transcript, start, end, \".\", \"+\", \".\", [transcript id; gene id; gene_name; xloc; ref_gene_id; contained_in; cmp_ref; class_code; tss_id]\n>exon entry = ref genome, analysis_pathway, exon, start, end, \".\", \"+\", \".\", [transcript id; gene id; exon number]\n\n******************\n## Installation\n### Quick Install\n```bash\npip install HIV_Isoform_Checker\n```\n### Manual Install\n```bash\ngit clone https://github.com/JessicaA2019/HIV_Isoform_Checker.git \ncd HIV_Isoform_Checker\npython setup.py install\n```\n### Dependencies\nSome dependencies currently do not install with Chimera_Buster. Please pip install the following dependencies:\n* mappy\n* regex\n* argparse\n\n### Testing\nTo test if the installation was successful run the following command with the test files provided:\n```bash\nHIV_Isoform_Checker HIC_test.gtf HIV_Isoform_Checker_test NL43.fa\n ```\n ******************\n## Usage\n```bash\nHIV Isoform Checker [options] input_file_name output_file_prefix ref_file\n```\n### Inputs\nTo run the pipeline the following input files are required:\n| Input | Description |\n| ------ | ------ |\n|input_file_name  |    Designates .gtf to be filtered. This is required.|\n|output_file_prefix  |  Designates output file prefix. This is required.|\n| ref_file | The location of the reference genome fasta file.|\n\nThe following input files are optional:\n| Arguement | Function |\n| ------ | ------ |\n|-h, --help |  Prints help message to terminal. |\n|-g value, --gap value |  Sets gap tolerance. If a read has a gap longer than this value that is not associated with known splice sites, the read will be filtered out as a sequencing error. The default value is 15. |\n|-a value, --startBP value |  Sets maximum starting bp in order to insure the presence of the first exon and by proxy a full read. Default is 700. |\n|-z value, --endBP value |  Sets minimum ending bp in order to insure the presence of the last exon and by proxy a full read. Default is 9500. |\n|-l value, --lengthFS value |  Sets maximum fully spliced transcript length. If a read is longer than this value, then it will be considered partial spliced. This is based off HIV-1 references without any major insertions or deletions. If your sequences have larger deletions (ex: deltaENV or deltaMA) or insertions (ex: multiple markers such as eGFP), this value may need to be adjusted. Default is 2500. |\n|-n value, --NCE value | When set to True, the output csv file will have y/n columns for the precence of NCEs. Default is False.|\n\n### Outputs\n The main output files created by the pipeline are:\n| Output | Description |\n|--------|-------------|\n| {output_file_prefix}.csv | A csv containing the sample ID, cluster size, isoform ID, starting and ending base pair, gffcompare class code, exon boundaries, coverage length, size normalized counts, precence of NCEs (optional), and any additional notes for each sample. |\n| {output_file_prefix}_altered.txt | A list of all sample IDs with altered exons to fill in small gaps. |\n| {output_file_prefix}_pass.txt | A list of all sample IDs that passed all filters. |\n| {output_file_prefix}_fail.txt | A list of all sample IDs that failed one or more filters. |\n| {output_file_prefix}_isoform_counts.csv | A csv of the usage counts and percentages for each isoform type. |\n| {output_file_prefix}_splice_site_usage.csv | A csv of the the usage counts and percentages for each donor site, acceptor site, and pairwise splice site combination. |\n| {output_file_prefix}_ref_coordinates.txt | A list of splice site and CDS coordinates for given reference. |\n| {output_file_prefix}.log | A log of the arguments given to the script. |\n\n**************************\n## Help\nFor issues or bugs, please report them on the [issues page][issues]. \n\n## License\n\nMIT - Copyright (c) 2024 Jessica Lauren Albert\n\n\n[//]: # (These are reference links used in the body of this note and get stripped out when the markdown processor does its job. There is no need to format nicely because it shouldn't be seen. Thanks SO - http://stackoverflow.com/questions/4823468/store-comments-in-markdown-syntax)\n\n   [APHIX]: <https://github.com/JessicaA2019/APHIX>\n   [issues]: <https://github.com/JessicaA2019/HIV_Isoform_Checker/issues>\n",
    "bugtrack_url": null,
    "license": "MIT License  Copyright (c) 2023 Jessica Lauren Albert  Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the \"Software\"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:  The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.  THE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. ",
    "summary": "Filters .gtf file of suspected HIV isoforms and confirms the isoform identities",
    "version": "1.1.9",
    "project_urls": {
        "Homepage": "https://github.com/JessicaA2019/HIV_Isoform_Checker"
    },
    "split_keywords": [
        "hiv",
        " isoforms",
        " splicing",
        " ont sequencing"
    ],
    "urls": [
        {
            "comment_text": "",
            "digests": {
                "blake2b_256": "601b518723f3106560c4279dd7218f252f12814eea7592285bfa19f11302269c",
                "md5": "582d8182dc495d06f0f6fe94286fa3a6",
                "sha256": "2a337a25c4d122920275c623dfbec539944dab89b2bc5de1593e270c9a4e4ab8"
            },
            "downloads": -1,
            "filename": "HIV_Isoform_Checker-1.1.9-py3-none-any.whl",
            "has_sig": false,
            "md5_digest": "582d8182dc495d06f0f6fe94286fa3a6",
            "packagetype": "bdist_wheel",
            "python_version": "py3",
            "requires_python": ">=3",
            "size": 17126,
            "upload_time": "2024-11-07T22:47:14",
            "upload_time_iso_8601": "2024-11-07T22:47:14.718004Z",
            "url": "https://files.pythonhosted.org/packages/60/1b/518723f3106560c4279dd7218f252f12814eea7592285bfa19f11302269c/HIV_Isoform_Checker-1.1.9-py3-none-any.whl",
            "yanked": false,
            "yanked_reason": null
        },
        {
            "comment_text": "",
            "digests": {
                "blake2b_256": "dd4298cacbf2d30bf703687a50873de9bc12870b4020421362275ee2a1f3604a",
                "md5": "a0c4f56a365caeb60e0291dfff766822",
                "sha256": "08398d24b697eae023d74df8241b9b521d45e99703696038dff54ec31dba1361"
            },
            "downloads": -1,
            "filename": "hiv_isoform_checker-1.1.9.tar.gz",
            "has_sig": false,
            "md5_digest": "a0c4f56a365caeb60e0291dfff766822",
            "packagetype": "sdist",
            "python_version": "source",
            "requires_python": ">=3",
            "size": 27276,
            "upload_time": "2024-11-07T22:47:15",
            "upload_time_iso_8601": "2024-11-07T22:47:15.740713Z",
            "url": "https://files.pythonhosted.org/packages/dd/42/98cacbf2d30bf703687a50873de9bc12870b4020421362275ee2a1f3604a/hiv_isoform_checker-1.1.9.tar.gz",
            "yanked": false,
            "yanked_reason": null
        }
    ],
    "upload_time": "2024-11-07 22:47:15",
    "github": true,
    "gitlab": false,
    "bitbucket": false,
    "codeberg": false,
    "github_user": "JessicaA2019",
    "github_project": "HIV_Isoform_Checker",
    "travis_ci": false,
    "coveralls": false,
    "github_actions": false,
    "lcname": "hiv-isoform-checker"
}
        
Elapsed time: 0.83806s