Skip to main content

Project vcf level coordinates from chromosomes to transcript

Project description

VCF projection from chromosome to transcriptome

The program uses ncls, a fast interval tree search program that considerably speeds up the ''stabbing'' problem of finding out the variant that overlaps the transcriptome interval.

GTF_FILE = 'Mus_musculus.GRCm38.91.gtf'
VCF_FILE = 'mgp_v4_indels_NOD_PWK.vcf'

from vcfproj import vcfproj
vcf_gtf = vcfproj.projection(GTF_FILE, VCF_FILE)

vcf_gtf is a dataframe, that would contain the following columns


'chrom_x', - chromosome of variation 
'start', - start of transcript
'end', - end of transcript
'gene', - gene name
'txome', - transcript name
'vcf_index', - ignore
'relative_pos', - relative position of variation on transcript
'transcript_length', - length of transcript
'chrom_y', - same as chrom_x
'pos', - position of variation
'id', - ID
'ref', - REF
'alt', - ALT
'qual', - QUAL
'filter', - FILTER
'info', - INFO
'format', - FORMAT
'samples' - SAMPLES 


Project details

Download files

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

Files for vcfproj, version 0.0.3
Filename, size File type Python version Upload date Hashes
Filename, size vcfproj-0.0.3-py3-none-any.whl (16.9 kB) File type Wheel Python version py3 Upload date Hashes View
Filename, size vcfproj-0.0.3.tar.gz (4.3 kB) File type Source Python version None Upload date Hashes View

Supported by

Pingdom Pingdom Monitoring Google Google Object Storage and Download Analytics Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page