This is a pre-production deployment of Warehouse. Changes made here affect the production instance of PyPI (pypi.python.org).
Help us improve Python packaging - Donate today!

So U Found A Mutation?

Project Description

Found a mutation in one or more samples? Now you want to check if they are in another sample. Unfortunately mutect, varscan or whatever other variant caller is not calling them. Use SUFAM. The super sensitive validation caller that calls everything on a given position. All you need is a vcf with the mutations that you are interested in and the sam/bam file of the sample where you want to find the same inconsipicuous mutation.

Installation

pip install sufam

Run

usage: sufam [-h] [--sample_name SAMPLE_NAME] [--format {matrix,sufam}]
            [--mpileup-parameters MPILEUP_PARAMETERS] [--version]
            reffa vcf bam

So U Found A Mutation? (SUFAM)

Found a mutation in one or more samples? Now you want to check if they are in
another sample. Unfortunately mutect, varscan or whatever other variant caller
is not calling them. Use SUFAM. The super sensitive validation caller that
calls everything on a given position. All you need is a vcf with the mutations
that you are interested in and the sam/bam file of the sample where you want to
find the same inconsipicuous mutation.

Author: inodb

positional arguments:
reffa                 Reference genome (fasta)
vcf                   VCF with mutations to be validated
bam                   BAM to find mutations in

optional arguments:
-h, --help            show this help message and exit
--sample_name SAMPLE_NAME
                        Set name of sample, used in output [name of bam].
--format {matrix,sufam}
                        Set output format [sufam]
--mpileup-parameters MPILEUP_PARAMETERS
                        Set options for mpileup [--ignore-RG --min-MQ 1 --max-
                        depth 250000 --max-idepth 250000]
--version             show program's version number and exit

Example

VCF file like:

#CHROM      POS     ID      REF     ALT
17  7574012 COSM11286,COSM214290    C       G
17  7574012 COSM11286,COSM214290    C       A

Check if given mutations are in a bam file:

sufam human_g1k_v37_chr17.fa mutations.vcf subset1.bam 2> example/sufam.log > example/sufam.tsv

Output:

Developers

Tests

In root dir run:

nosetests

For individual tests:

nosetests -s tests/test_validation.py:TestValidation.test_validate_mutations_indel
Release History

Release History

This version
History Node

0.4.2

History Node

0.4.1

History Node

0.4

History Node

0.3

History Node

0.2

History Node

0.1

Download Files

Download Files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

File Name & Checksum SHA256 Checksum Help Version File Type Upload Date
sufam-0.4.2-py2.py3-none-any.whl (13.4 kB) Copy SHA256 Checksum SHA256 2.7 Wheel Dec 24, 2015
sufam-0.4.2.tar.gz (25.0 kB) Copy SHA256 Checksum SHA256 Source Dec 24, 2015

Supported By

WebFaction WebFaction Technical Writing Elastic Elastic Search Pingdom Pingdom Monitoring Dyn Dyn DNS Sentry Sentry Error Logging CloudAMQP CloudAMQP RabbitMQ Heroku Heroku PaaS Kabu Creative Kabu Creative UX & Design Fastly Fastly CDN DigiCert DigiCert EV Certificate Rackspace Rackspace Cloud Servers DreamHost DreamHost Log Hosting