..
This file was automatically generated by docs/create.py.
README
******
.. image:: https://badge.fury.io/py/fuc.svg
:target: https://badge.fury.io/py/fuc
.. image:: https://readthedocs.org/projects/sbslee-fuc/badge/?version=latest
:target: https://sbslee-fuc.readthedocs.io/en/latest/?badge=latest
:alt: Documentation Status
.. image:: https://anaconda.org/bioconda/fuc/badges/version.svg
:target: https://anaconda.org/bioconda/fuc
.. image:: https://anaconda.org/bioconda/fuc/badges/license.svg
:target: https://github.com/sbslee/fuc/blob/main/LICENSE
.. image:: https://anaconda.org/bioconda/fuc/badges/downloads.svg
:target: https://anaconda.org/bioconda/fuc/files
Introduction
============
The main goal of the fuc package (pronounced "eff-you-see") is to wrap some of the most **f**\ requently **u**\ sed **c**\ ommands in the field of bioinformatics into one place.
The package is written in Python, and supports both command line interface (CLI) and application programming interface (API) whose documentations are available at the `Read the Docs <https://sbslee-fuc.readthedocs.io/en/latest/>`_.
Currently, fuc can be used to analyze, summarize, visualize, and manipulate the following file formats:
- Sequence Alignment/Map (SAM)
- Binary Alignment/Map (BAM)
- CRAM
- Variant Call Format (VCF)
- Mutation Annotation Format (MAF)
- Browser Extensible Data (BED)
- FASTQ
- FASTA
- General Feature Format (GFF)
- Gene Transfer Format (GTF)
- delimiter-separated values format (e.g. comma-separated values or CSV format)
Additionally, fuc can be used to parse output data from the following programs:
- `Ensembl Variant Effect Predictor (VEP) <https://asia.ensembl.org/info/docs/tools/vep/index.html>`__
- `SnpEff <http://pcingola.github.io/SnpEff/>`__
- `bcl2fastq and bcl2fastq2 <https://sapac.support.illumina.com/sequencing/sequencing_software/bcl2fastq-conversion-software.html>`__
- `Kallisto <https://pachterlab.github.io/kallisto/>`__
Your contributions (e.g. feature ideas, pull requests) are most welcome.
| Author: Seung-been "Steven" Lee
| Email: sbstevenlee@gmail.com
| License: MIT License
Citation
========
If you use fuc in a published analysis, please report the program version
and cite the following article:
Lee et al., 2022. `ClinPharmSeq: A targeted sequencing panel for clinical pharmacogenetics implementation <https://doi.org/10.1371/journal.pone.0272129>`__. PLOS ONE.
Support fuc
===========
If you find my work useful, please consider becoming a `sponsor <https://github.com/sponsors/sbslee>`__.
Installation
============
The following packages are required to run fuc:
.. parsed-literal::
biopython
lxml
matplotlib
matplotlib-venn
numpy
pandas
pyranges
pysam
scipy
seaborn
statsmodels
There are various ways you can install fuc. The recommended way is via conda (`Anaconda <https://www.anaconda.com/>`__):
.. code-block:: text
$ conda install -c bioconda fuc
Above will automatically download and install all the dependencies as well. Alternatively, you can use pip (`PyPI <https://pypi.org/>`__) to install fuc and all of its dependencies:
.. code-block:: text
$ pip install fuc
Finally, you can clone the GitHub repository and then install fuc locally:
.. code-block:: text
$ git clone https://github.com/sbslee/fuc
$ cd fuc
$ pip install .
The nice thing about this approach is that you will have access to development versions that are not available in Anaconda or PyPI. For example, you can access a development branch with the ``git checkout`` command. When you do this, please make sure your environment already has all the dependencies installed.
Getting help
============
For detailed documentations on the fuc package's CLI and API, please refer to the `Read the Docs <https://sbslee-fuc.readthedocs.io/en/latest/>`_.
For getting help on the fuc CLI:
.. code-block:: text
$ fuc -h
usage: fuc [-h] [-v] COMMAND ...
positional arguments:
COMMAND
bam-aldepth Compute allelic depth from a BAM file.
bam-depth Compute per-base read depth from BAM files.
bam-head Print the header of a BAM file.
bam-index Index a BAM file.
bam-rename Rename the sample in a BAM file.
bam-slice Slice a BAM file.
bed-intxn Find the intersection of BED files.
bed-sum Summarize a BED file.
cov-concat Concatenate depth of coverage files.
cov-rename Rename the samples in a depth of coverage file.
fa-filter Filter sequence records in a FASTA file.
fq-count Count sequence reads in FASTQ files.
fq-sum Summarize a FASTQ file.
fuc-bgzip Write a BGZF compressed file.
fuc-compf Compare the contents of two files.
fuc-demux Parse the Reports directory from bcl2fastq.
fuc-exist Check whether certain files exist.
fuc-find Retrieve absolute paths of files whose name matches a
specified pattern, optionally recursively.
fuc-undetm Compute top unknown barcodes using undertermined FASTQ from
bcl2fastq.
maf-maf2vcf Convert a MAF file to a VCF file.
maf-oncoplt Create an oncoplot with a MAF file.
maf-sumplt Create a summary plot with a MAF file.
maf-vcf2maf Convert a VCF file to a MAF file.
ngs-bam2fq Pipeline for converting BAM files to FASTQ files.
ngs-fq2bam Pipeline for converting FASTQ files to analysis-ready BAM
files.
ngs-hc Pipeline for germline short variant discovery.
ngs-m2 Pipeline for somatic short variant discovery.
ngs-pon Pipeline for constructing a panel of normals (PoN).
ngs-quant Pipeline for running RNAseq quantification from FASTQ files
with Kallisto.
ngs-trim Pipeline for trimming adapters from FASTQ files.
tabix-index Index a GFF/BED/SAM/VCF file with Tabix.
tabix-slice Slice a GFF/BED/SAM/VCF file with Tabix.
tbl-merge Merge two table files.
tbl-sum Summarize a table file.
vcf-call Call SNVs and indels from BAM files.
vcf-filter Filter a VCF file.
vcf-index Index a VCF file.
vcf-merge Merge two or more VCF files.
vcf-rename Rename the samples in a VCF file.
vcf-slice Slice a VCF file for specified regions.
vcf-split Split a VCF file by individual.
vcf-vcf2bed Convert a VCF file to a BED file.
vcf-vep Filter a VCF file by annotations from Ensembl VEP.
optional arguments:
-h, --help Show this help message and exit.
-v, --version Show the version number and exit.
For getting help on a specific command (e.g. vcf-merge):
.. code-block:: text
$ fuc vcf-merge -h
Below is the list of submodules available in the fuc API:
- **common** : The common submodule is used by other fuc submodules such as pyvcf and pybed. It also provides many day-to-day actions used in the field of bioinformatics.
- **pybam** : The pybam submodule is designed for working with sequence alignment files (SAM/BAM/CRAM). It essentially wraps the `pysam <https://pysam.readthedocs.io/en/latest/api.html>`_ package to allow fast computation and easy manipulation. If you are mainly interested in working with depth of coverage data, please check out the pycov submodule which is specifically designed for the task.
- **pybed** : The pybed submodule is designed for working with BED files. It implements ``pybed.BedFrame`` which stores BED data as ``pandas.DataFrame`` via the `pyranges <https://github.com/biocore-ntnu/pyranges>`_ package to allow fast computation and easy manipulation. The submodule strictly adheres to the standard `BED specification <https://genome.ucsc.edu/FAQ/FAQformat.html>`_.
- **pychip** : The pychip submodule is designed for working with annotation or manifest files from the Axiom (Thermo Fisher Scientific) and Infinium (Illumina) array platforms.
- **pycov** : The pycov submodule is designed for working with depth of coverage data from sequence alingment files (SAM/BAM/CRAM). It implements ``pycov.CovFrame`` which stores read depth data as ``pandas.DataFrame`` via the `pysam <https://pysam.readthedocs.io/en/latest/api.html>`_ package to allow fast computation and easy manipulation. The ``pycov.CovFrame`` class also contains many useful plotting methods such as ``CovFrame.plot_region`` and ``CovFrame.plot_uniformity``.
- **pyfq** : The pyfq submodule is designed for working with FASTQ files. It implements ``pyfq.FqFrame`` which stores FASTQ data as ``pandas.DataFrame`` to allow fast computation and easy manipulation.
- **pygff** : The pygff submodule is designed for working with GFF/GTF files. It implements ``pygff.GffFrame`` which stores GFF/GTF data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The submodule strictly adheres to the standard `GFF specification <https://github.com/The-Sequence-Ontology/Specifications/blob/master/gff3.md>`_.
- **pykallisto** : The pykallisto submodule is designed for working with RNAseq quantification data from Kallisto. It implements ``pykallisto.KallistoFrame`` which stores Kallisto's output data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The ``pykallisto.KallistoFrame`` class also contains many useful plotting methods such as ``KallistoFrame.plot_differential_abundance``.
- **pymaf** : The pymaf submodule is designed for working with MAF files. It implements ``pymaf.MafFrame`` which stores MAF data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The ``pymaf.MafFrame`` class also contains many useful plotting methods such as ``MafFrame.plot_oncoplot`` and ``MafFrame.plot_summary``. The submodule strictly adheres to the standard `MAF specification <https://docs.gdc.cancer.gov/Data/File_Formats/MAF_Format/>`_.
- **pysnpeff** : The pysnpeff submodule is designed for parsing VCF annotation data from the `SnpEff <https://pcingola.github.io/SnpEff/>`_ program. It should be used with ``pyvcf.VcfFrame``.
- **pyvcf** : The pyvcf submodule is designed for working with VCF files. It implements ``pyvcf.VcfFrame`` which stores VCF data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The ``pyvcf.VcfFrame`` class also contains many useful plotting methods such as ``VcfFrame.plot_comparison`` and ``VcfFrame.plot_tmb``. The submodule strictly adheres to the standard `VCF specification <https://samtools.github.io/hts-specs/VCFv4.3.pdf>`_.
- **pyvep** : The pyvep submodule is designed for parsing VCF annotation data from the `Ensembl VEP <https://asia.ensembl.org/info/docs/tools/vep/index.html>`_ program. It should be used with ``pyvcf.VcfFrame``.
For getting help on a specific submodule (e.g. pyvcf):
.. code:: python3
>>> from fuc import pyvcf
>>> help(pyvcf)
In Jupyter Notebook and Lab, you can see the documentation for a python
function by hitting ``SHIFT + TAB``. Hit it twice to expand the view.
CLI examples
============
**SAM/BAM/CRAM**
To print the header of a SAM file:
.. code-block:: text
$ fuc bam-head in.sam
To index a CRAM file:
.. code-block:: text
$ fuc bam-index in.cram
To rename the samples in a SAM file:
.. code-block:: text
$ fuc bam-rename in.sam NA12878 > out.sam
To slice a BAM file:
.. code-block:: text
$ fuc bam-slice in.bam chr1:100-200 > out.bam
**BED**
To find intersection between BED files:
.. code-block:: text
$ fuc bed-intxn 1.bed 2.bed 3.bed > intersect.bed
**FASTQ**
To count sequence reads in a FASTQ file:
.. code-block:: text
$ fuc fq-count example.fastq
**FUC**
To check whether a file exists in the operating system:
.. code-block:: text
$ fuc fuc-exist example.txt
To find all VCF files within the current directory recursively:
.. code-block:: text
$ fuc fuc-find .vcf.gz
**TABLE**
To merge two tab-delimited files:
.. code-block:: text
$ fuc tbl-merge left.tsv right.tsv > merged.tsv
**VCF**
To merge VCF files:
.. code-block:: text
$ fuc vcf-merge 1.vcf 2.vcf 3.vcf > merged.vcf
To filter a VCF file annotated by Ensembl VEP:
.. code-block:: text
$ fuc vcf-vep in.vcf 'SYMBOL == "TP53"' > out.vcf
API examples
============
**BAM**
To create read depth profile of a region from a CRAM file:
.. code:: python3
>>> from fuc import pycov
>>> cf = pycov.CovFrame.from_file('HG00525.final.cram', zero=True,
... region='chr12:21161194-21239796', names=['HG00525'])
>>> cf.plot_region('chr12:21161194-21239796')
.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/coverage.png
**VCF**
To filter a VCF file based on a BED file:
.. code:: python3
>>> from fuc import pyvcf
>>> vf = pyvcf.VcfFrame.from_file('original.vcf')
>>> filtered_vf = vf.filter_bed('targets.bed')
>>> filtered_vf.to_file('filtered.vcf')
To remove indels from a VCF file:
.. code:: python3
>>> from fuc import pyvcf
>>> vf = pyvcf.VcfFrame.from_file('with_indels.vcf')
>>> filtered_vf = vf.filter_indel()
>>> filtered_vf.to_file('no_indels.vcf')
To create a Venn diagram showing genotype concordance between groups:
.. code:: python3
>>> from fuc import pyvcf, common
>>> common.load_dataset('pyvcf')
>>> f = '~/fuc-data/pyvcf/plot_comparison.vcf'
>>> vf = pyvcf.VcfFrame.from_file(f)
>>> a = ['Steven_A', 'John_A', 'Sara_A']
>>> b = ['Steven_B', 'John_B', 'Sara_B']
>>> c = ['Steven_C', 'John_C', 'Sara_C']
>>> vf.plot_comparison(a, b, c)
.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/plot_comparison.png
To create various figures for normal-tumor analysis:
.. code:: python3
>>> import matplotlib.pyplot as plt
>>> from fuc import common, pyvcf
>>> common.load_dataset('pyvcf')
>>> vf = pyvcf.VcfFrame.from_file('~/fuc-data/pyvcf/normal-tumor.vcf')
>>> af = pyvcf.AnnFrame.from_file('~/fuc-data/pyvcf/normal-tumor-annot.tsv', sample_col='Sample')
>>> normal = af.df[af.df.Tissue == 'Normal'].index
>>> tumor = af.df[af.df.Tissue == 'Tumor'].index
>>> fig, [[ax1, ax2], [ax3, ax4]] = plt.subplots(2, 2, figsize=(10, 10))
>>> vf.plot_tmb(ax=ax1)
>>> vf.plot_tmb(ax=ax2, af=af, group_col='Tissue')
>>> vf.plot_hist_format('#DP', ax=ax3, af=af, group_col='Tissue')
>>> vf.plot_regplot(normal, tumor, ax=ax4)
>>> plt.tight_layout()
.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/normal-tumor.png
**MAF**
To create an oncoplot with a MAF file:
.. code:: python3
>>> from fuc import common, pymaf
>>> common.load_dataset('tcga-laml')
>>> maf_file = '~/fuc-data/tcga-laml/tcga_laml.maf.gz'
>>> mf = pymaf.MafFrame.from_file(maf_file)
>>> mf.plot_oncoplot()
.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/oncoplot.png
To create a customized oncoplot with a MAF file, see the `Create customized oncoplot <https://sbslee-fuc.readthedocs.io/en/latest/tutorials.html#create-customized-oncoplots>`__ tutorial:
.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/customized_oncoplot.png
To create a summary figure for a MAF file:
.. code:: python3
>>> from fuc import common, pymaf
>>> common.load_dataset('tcga-laml')
>>> maf_file = '~/fuc-data/tcga-laml/tcga_laml.maf.gz'
>>> mf = pymaf.MafFrame.from_file(maf_file)
>>> mf.plot_summary()
.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/maf_summary-2.png
Raw data
{
"_id": null,
"home_page": "https://github.com/sbslee/fuc",
"name": "fuc",
"maintainer": null,
"docs_url": null,
"requires_python": null,
"maintainer_email": null,
"keywords": null,
"author": "Seung-been \"Steven\" Lee",
"author_email": "sbstevenlee@gmail.com",
"download_url": "https://files.pythonhosted.org/packages/73/cd/accbd737f49c9f313d440d17b9562acebece6e5a8214517eab2adfb192bd/fuc-0.38.0.tar.gz",
"platform": null,
"description": "..\n This file was automatically generated by docs/create.py.\n\nREADME\n******\n\n.. image:: https://badge.fury.io/py/fuc.svg\n :target: https://badge.fury.io/py/fuc\n\n.. image:: https://readthedocs.org/projects/sbslee-fuc/badge/?version=latest\n :target: https://sbslee-fuc.readthedocs.io/en/latest/?badge=latest\n :alt: Documentation Status\n\n.. image:: https://anaconda.org/bioconda/fuc/badges/version.svg\n :target: https://anaconda.org/bioconda/fuc\n\n.. image:: https://anaconda.org/bioconda/fuc/badges/license.svg\n :target: https://github.com/sbslee/fuc/blob/main/LICENSE\n\n.. image:: https://anaconda.org/bioconda/fuc/badges/downloads.svg\n :target: https://anaconda.org/bioconda/fuc/files\n\nIntroduction\n============\n\nThe main goal of the fuc package (pronounced \"eff-you-see\") is to wrap some of the most **f**\\ requently **u**\\ sed **c**\\ ommands in the field of bioinformatics into one place.\n\nThe package is written in Python, and supports both command line interface (CLI) and application programming interface (API) whose documentations are available at the `Read the Docs <https://sbslee-fuc.readthedocs.io/en/latest/>`_.\n\nCurrently, fuc can be used to analyze, summarize, visualize, and manipulate the following file formats:\n\n- Sequence Alignment/Map (SAM)\n- Binary Alignment/Map (BAM)\n- CRAM\n- Variant Call Format (VCF)\n- Mutation Annotation Format (MAF)\n- Browser Extensible Data (BED)\n- FASTQ\n- FASTA\n- General Feature Format (GFF)\n- Gene Transfer Format (GTF)\n- delimiter-separated values format (e.g. comma-separated values or CSV format)\n\nAdditionally, fuc can be used to parse output data from the following programs:\n\n- `Ensembl Variant Effect Predictor (VEP) <https://asia.ensembl.org/info/docs/tools/vep/index.html>`__\n- `SnpEff <http://pcingola.github.io/SnpEff/>`__\n- `bcl2fastq and bcl2fastq2 <https://sapac.support.illumina.com/sequencing/sequencing_software/bcl2fastq-conversion-software.html>`__\n- `Kallisto <https://pachterlab.github.io/kallisto/>`__\n\nYour contributions (e.g. feature ideas, pull requests) are most welcome.\n\n| Author: Seung-been \"Steven\" Lee\n| Email: sbstevenlee@gmail.com\n| License: MIT License\n\nCitation\n========\n\nIf you use fuc in a published analysis, please report the program version\nand cite the following article:\n\nLee et al., 2022. `ClinPharmSeq: A targeted sequencing panel for clinical pharmacogenetics implementation <https://doi.org/10.1371/journal.pone.0272129>`__. PLOS ONE.\n\nSupport fuc\n===========\n\nIf you find my work useful, please consider becoming a `sponsor <https://github.com/sponsors/sbslee>`__.\n\nInstallation\n============\n\nThe following packages are required to run fuc:\n\n.. parsed-literal::\n\n biopython\n lxml\n matplotlib\n matplotlib-venn\n numpy\n pandas\n pyranges\n pysam\n scipy\n seaborn\n statsmodels\n\nThere are various ways you can install fuc. The recommended way is via conda (`Anaconda <https://www.anaconda.com/>`__):\n\n.. code-block:: text\n\n $ conda install -c bioconda fuc\n\nAbove will automatically download and install all the dependencies as well. Alternatively, you can use pip (`PyPI <https://pypi.org/>`__) to install fuc and all of its dependencies:\n\n.. code-block:: text\n\n $ pip install fuc\n\nFinally, you can clone the GitHub repository and then install fuc locally:\n\n.. code-block:: text\n\n $ git clone https://github.com/sbslee/fuc\n $ cd fuc\n $ pip install .\n\nThe nice thing about this approach is that you will have access to development versions that are not available in Anaconda or PyPI. For example, you can access a development branch with the ``git checkout`` command. When you do this, please make sure your environment already has all the dependencies installed.\n\nGetting help\n============\n\nFor detailed documentations on the fuc package's CLI and API, please refer to the `Read the Docs <https://sbslee-fuc.readthedocs.io/en/latest/>`_.\n\nFor getting help on the fuc CLI:\n\n.. code-block:: text\n\n $ fuc -h\n usage: fuc [-h] [-v] COMMAND ...\n \n positional arguments:\n COMMAND\n bam-aldepth Compute allelic depth from a BAM file.\n bam-depth Compute per-base read depth from BAM files.\n bam-head Print the header of a BAM file.\n bam-index Index a BAM file.\n bam-rename Rename the sample in a BAM file.\n bam-slice Slice a BAM file.\n bed-intxn Find the intersection of BED files.\n bed-sum Summarize a BED file.\n cov-concat Concatenate depth of coverage files.\n cov-rename Rename the samples in a depth of coverage file.\n fa-filter Filter sequence records in a FASTA file.\n fq-count Count sequence reads in FASTQ files.\n fq-sum Summarize a FASTQ file.\n fuc-bgzip Write a BGZF compressed file.\n fuc-compf Compare the contents of two files.\n fuc-demux Parse the Reports directory from bcl2fastq.\n fuc-exist Check whether certain files exist.\n fuc-find Retrieve absolute paths of files whose name matches a\n specified pattern, optionally recursively.\n fuc-undetm Compute top unknown barcodes using undertermined FASTQ from\n bcl2fastq.\n maf-maf2vcf Convert a MAF file to a VCF file.\n maf-oncoplt Create an oncoplot with a MAF file.\n maf-sumplt Create a summary plot with a MAF file.\n maf-vcf2maf Convert a VCF file to a MAF file.\n ngs-bam2fq Pipeline for converting BAM files to FASTQ files.\n ngs-fq2bam Pipeline for converting FASTQ files to analysis-ready BAM\n files.\n ngs-hc Pipeline for germline short variant discovery.\n ngs-m2 Pipeline for somatic short variant discovery.\n ngs-pon Pipeline for constructing a panel of normals (PoN).\n ngs-quant Pipeline for running RNAseq quantification from FASTQ files\n with Kallisto.\n ngs-trim Pipeline for trimming adapters from FASTQ files.\n tabix-index Index a GFF/BED/SAM/VCF file with Tabix.\n tabix-slice Slice a GFF/BED/SAM/VCF file with Tabix.\n tbl-merge Merge two table files.\n tbl-sum Summarize a table file.\n vcf-call Call SNVs and indels from BAM files.\n vcf-filter Filter a VCF file.\n vcf-index Index a VCF file.\n vcf-merge Merge two or more VCF files.\n vcf-rename Rename the samples in a VCF file.\n vcf-slice Slice a VCF file for specified regions.\n vcf-split Split a VCF file by individual.\n vcf-vcf2bed Convert a VCF file to a BED file.\n vcf-vep Filter a VCF file by annotations from Ensembl VEP.\n \n optional arguments:\n -h, --help Show this help message and exit.\n -v, --version Show the version number and exit.\n\nFor getting help on a specific command (e.g. vcf-merge):\n\n.. code-block:: text\n\n $ fuc vcf-merge -h\n\nBelow is the list of submodules available in the fuc API:\n\n- **common** : The common submodule is used by other fuc submodules such as pyvcf and pybed. It also provides many day-to-day actions used in the field of bioinformatics.\n- **pybam** : The pybam submodule is designed for working with sequence alignment files (SAM/BAM/CRAM). It essentially wraps the `pysam <https://pysam.readthedocs.io/en/latest/api.html>`_ package to allow fast computation and easy manipulation. If you are mainly interested in working with depth of coverage data, please check out the pycov submodule which is specifically designed for the task.\n- **pybed** : The pybed submodule is designed for working with BED files. It implements ``pybed.BedFrame`` which stores BED data as ``pandas.DataFrame`` via the `pyranges <https://github.com/biocore-ntnu/pyranges>`_ package to allow fast computation and easy manipulation. The submodule strictly adheres to the standard `BED specification <https://genome.ucsc.edu/FAQ/FAQformat.html>`_.\n- **pychip** : The pychip submodule is designed for working with annotation or manifest files from the Axiom (Thermo Fisher Scientific) and Infinium (Illumina) array platforms.\n- **pycov** : The pycov submodule is designed for working with depth of coverage data from sequence alingment files (SAM/BAM/CRAM). It implements ``pycov.CovFrame`` which stores read depth data as ``pandas.DataFrame`` via the `pysam <https://pysam.readthedocs.io/en/latest/api.html>`_ package to allow fast computation and easy manipulation. The ``pycov.CovFrame`` class also contains many useful plotting methods such as ``CovFrame.plot_region`` and ``CovFrame.plot_uniformity``.\n- **pyfq** : The pyfq submodule is designed for working with FASTQ files. It implements ``pyfq.FqFrame`` which stores FASTQ data as ``pandas.DataFrame`` to allow fast computation and easy manipulation.\n- **pygff** : The pygff submodule is designed for working with GFF/GTF files. It implements ``pygff.GffFrame`` which stores GFF/GTF data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The submodule strictly adheres to the standard `GFF specification <https://github.com/The-Sequence-Ontology/Specifications/blob/master/gff3.md>`_.\n- **pykallisto** : The pykallisto submodule is designed for working with RNAseq quantification data from Kallisto. It implements ``pykallisto.KallistoFrame`` which stores Kallisto's output data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The ``pykallisto.KallistoFrame`` class also contains many useful plotting methods such as ``KallistoFrame.plot_differential_abundance``.\n- **pymaf** : The pymaf submodule is designed for working with MAF files. It implements ``pymaf.MafFrame`` which stores MAF data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The ``pymaf.MafFrame`` class also contains many useful plotting methods such as ``MafFrame.plot_oncoplot`` and ``MafFrame.plot_summary``. The submodule strictly adheres to the standard `MAF specification <https://docs.gdc.cancer.gov/Data/File_Formats/MAF_Format/>`_.\n- **pysnpeff** : The pysnpeff submodule is designed for parsing VCF annotation data from the `SnpEff <https://pcingola.github.io/SnpEff/>`_ program. It should be used with ``pyvcf.VcfFrame``.\n- **pyvcf** : The pyvcf submodule is designed for working with VCF files. It implements ``pyvcf.VcfFrame`` which stores VCF data as ``pandas.DataFrame`` to allow fast computation and easy manipulation. The ``pyvcf.VcfFrame`` class also contains many useful plotting methods such as ``VcfFrame.plot_comparison`` and ``VcfFrame.plot_tmb``. The submodule strictly adheres to the standard `VCF specification <https://samtools.github.io/hts-specs/VCFv4.3.pdf>`_.\n- **pyvep** : The pyvep submodule is designed for parsing VCF annotation data from the `Ensembl VEP <https://asia.ensembl.org/info/docs/tools/vep/index.html>`_ program. It should be used with ``pyvcf.VcfFrame``.\n\nFor getting help on a specific submodule (e.g. pyvcf):\n\n.. code:: python3\n\n >>> from fuc import pyvcf\n >>> help(pyvcf)\n\nIn Jupyter Notebook and Lab, you can see the documentation for a python\nfunction by hitting ``SHIFT + TAB``. Hit it twice to expand the view.\n\nCLI examples\n============\n\n**SAM/BAM/CRAM**\n\nTo print the header of a SAM file:\n\n.. code-block:: text\n\n $ fuc bam-head in.sam\n\nTo index a CRAM file:\n\n.. code-block:: text\n\n $ fuc bam-index in.cram\n\nTo rename the samples in a SAM file:\n\n.. code-block:: text\n\n $ fuc bam-rename in.sam NA12878 > out.sam\n\nTo slice a BAM file:\n\n.. code-block:: text\n\n $ fuc bam-slice in.bam chr1:100-200 > out.bam\n\n**BED**\n\nTo find intersection between BED files:\n\n.. code-block:: text\n\n $ fuc bed-intxn 1.bed 2.bed 3.bed > intersect.bed\n\n**FASTQ**\n\nTo count sequence reads in a FASTQ file:\n\n.. code-block:: text\n\n $ fuc fq-count example.fastq\n\n**FUC**\n\nTo check whether a file exists in the operating system:\n\n.. code-block:: text\n\n $ fuc fuc-exist example.txt\n\nTo find all VCF files within the current directory recursively:\n\n.. code-block:: text\n\n $ fuc fuc-find .vcf.gz\n\n**TABLE**\n\nTo merge two tab-delimited files:\n\n.. code-block:: text\n\n $ fuc tbl-merge left.tsv right.tsv > merged.tsv\n\n**VCF**\n\nTo merge VCF files:\n\n.. code-block:: text\n\n $ fuc vcf-merge 1.vcf 2.vcf 3.vcf > merged.vcf\n\nTo filter a VCF file annotated by Ensembl VEP:\n\n.. code-block:: text\n\n $ fuc vcf-vep in.vcf 'SYMBOL == \"TP53\"' > out.vcf\n\nAPI examples\n============\n\n**BAM**\n\nTo create read depth profile of a region from a CRAM file:\n\n.. code:: python3\n\n >>> from fuc import pycov\n >>> cf = pycov.CovFrame.from_file('HG00525.final.cram', zero=True,\n ... region='chr12:21161194-21239796', names=['HG00525'])\n >>> cf.plot_region('chr12:21161194-21239796')\n\n.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/coverage.png\n\n**VCF**\n\nTo filter a VCF file based on a BED file:\n\n.. code:: python3\n\n >>> from fuc import pyvcf\n >>> vf = pyvcf.VcfFrame.from_file('original.vcf')\n >>> filtered_vf = vf.filter_bed('targets.bed')\n >>> filtered_vf.to_file('filtered.vcf')\n\nTo remove indels from a VCF file:\n\n.. code:: python3\n\n >>> from fuc import pyvcf\n >>> vf = pyvcf.VcfFrame.from_file('with_indels.vcf')\n >>> filtered_vf = vf.filter_indel()\n >>> filtered_vf.to_file('no_indels.vcf')\n\nTo create a Venn diagram showing genotype concordance between groups:\n\n.. code:: python3\n\n >>> from fuc import pyvcf, common\n >>> common.load_dataset('pyvcf')\n >>> f = '~/fuc-data/pyvcf/plot_comparison.vcf'\n >>> vf = pyvcf.VcfFrame.from_file(f)\n >>> a = ['Steven_A', 'John_A', 'Sara_A']\n >>> b = ['Steven_B', 'John_B', 'Sara_B']\n >>> c = ['Steven_C', 'John_C', 'Sara_C']\n >>> vf.plot_comparison(a, b, c)\n\n.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/plot_comparison.png\n\nTo create various figures for normal-tumor analysis:\n\n.. code:: python3\n\n >>> import matplotlib.pyplot as plt\n >>> from fuc import common, pyvcf\n >>> common.load_dataset('pyvcf')\n >>> vf = pyvcf.VcfFrame.from_file('~/fuc-data/pyvcf/normal-tumor.vcf')\n >>> af = pyvcf.AnnFrame.from_file('~/fuc-data/pyvcf/normal-tumor-annot.tsv', sample_col='Sample')\n >>> normal = af.df[af.df.Tissue == 'Normal'].index\n >>> tumor = af.df[af.df.Tissue == 'Tumor'].index\n >>> fig, [[ax1, ax2], [ax3, ax4]] = plt.subplots(2, 2, figsize=(10, 10))\n >>> vf.plot_tmb(ax=ax1)\n >>> vf.plot_tmb(ax=ax2, af=af, group_col='Tissue')\n >>> vf.plot_hist_format('#DP', ax=ax3, af=af, group_col='Tissue')\n >>> vf.plot_regplot(normal, tumor, ax=ax4)\n >>> plt.tight_layout()\n\n.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/normal-tumor.png\n\n**MAF**\n\nTo create an oncoplot with a MAF file:\n\n.. code:: python3\n\n >>> from fuc import common, pymaf\n >>> common.load_dataset('tcga-laml')\n >>> maf_file = '~/fuc-data/tcga-laml/tcga_laml.maf.gz'\n >>> mf = pymaf.MafFrame.from_file(maf_file)\n >>> mf.plot_oncoplot()\n\n.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/oncoplot.png\n\nTo create a customized oncoplot with a MAF file, see the `Create customized oncoplot <https://sbslee-fuc.readthedocs.io/en/latest/tutorials.html#create-customized-oncoplots>`__ tutorial:\n\n.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/customized_oncoplot.png\n\nTo create a summary figure for a MAF file:\n\n.. code:: python3\n\n >>> from fuc import common, pymaf\n >>> common.load_dataset('tcga-laml')\n >>> maf_file = '~/fuc-data/tcga-laml/tcga_laml.maf.gz'\n >>> mf = pymaf.MafFrame.from_file(maf_file)\n >>> mf.plot_summary()\n\n.. image:: https://raw.githubusercontent.com/sbslee/fuc-data/main/images/maf_summary-2.png\n\n",
"bugtrack_url": null,
"license": "MIT",
"summary": "Frequently used commands in bioinformatics",
"version": "0.38.0",
"project_urls": {
"Homepage": "https://github.com/sbslee/fuc"
},
"split_keywords": [],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "deef39b9f73bc9dac83aa7c21ff824807884d35c7af333f739813ec1a958f544",
"md5": "e52d83bfb1e25a8b94772aa6f362a0e4",
"sha256": "7b73d5efc220510a735f815eea3ebfa4c775aed45559b325aa390fa3be6490dd"
},
"downloads": -1,
"filename": "fuc-0.38.0-py3-none-any.whl",
"has_sig": false,
"md5_digest": "e52d83bfb1e25a8b94772aa6f362a0e4",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": null,
"size": 154743,
"upload_time": "2024-06-15T22:56:01",
"upload_time_iso_8601": "2024-06-15T22:56:01.220035Z",
"url": "https://files.pythonhosted.org/packages/de/ef/39b9f73bc9dac83aa7c21ff824807884d35c7af333f739813ec1a958f544/fuc-0.38.0-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "73cdaccbd737f49c9f313d440d17b9562acebece6e5a8214517eab2adfb192bd",
"md5": "9e6a7e6a41055e0f64d97340c8ce2363",
"sha256": "379ffdfeaa77ba4d9fd38360e1032986ad3e3da8b380f57652cbcf4b9c1c13dc"
},
"downloads": -1,
"filename": "fuc-0.38.0.tar.gz",
"has_sig": false,
"md5_digest": "9e6a7e6a41055e0f64d97340c8ce2363",
"packagetype": "sdist",
"python_version": "source",
"requires_python": null,
"size": 128592,
"upload_time": "2024-06-15T22:56:03",
"upload_time_iso_8601": "2024-06-15T22:56:03.993472Z",
"url": "https://files.pythonhosted.org/packages/73/cd/accbd737f49c9f313d440d17b9562acebece6e5a8214517eab2adfb192bd/fuc-0.38.0.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2024-06-15 22:56:03",
"github": true,
"gitlab": false,
"bitbucket": false,
"codeberg": false,
"github_user": "sbslee",
"github_project": "fuc",
"travis_ci": false,
"coveralls": false,
"github_actions": false,
"lcname": "fuc"
}
JSON
