Creates self-contained html pages for visual variant review with IGV (igv.js).
Project description
igv-reports
Python application to generate self-contained igv.js pages that can be opened within a browser with "file" protocol. The generated html page contains all data neccessary for IGV as uuencoded blobs.
Installation
Prerequisites
igv-reports requires Python 3.6 or greater and pip.
As with all Python projects, use of a virtual environment is recommended.
Instructions for creating a virtual environment using conda
follow.
1. Install Anaconda from https://docs.anaconda.com/anaconda/
2. Create a virtual environment
conda create -n myenv python=3.7.1
conda install -n myenv pip
conda activate myenv
Installing igv-reports
pip install igv-reports
igv-reports requires the package pysam which should be installed automatically. However on OSX this sometimes
fails due to missing dependent libraries. This can be fixed following the procedure below, from the pysam
docs;
"The recommended way to install pysam is through conda/bioconda.
This will install pysam from the bioconda channel and automatically makes sure that dependencies are installed.
Also, compilation flags will be set automatically, which will potentially save a lot of trouble on OS X."
conda config --add channels r
conda config --add channels bioconda
conda install pysam
Creating a report
A report consists of a table of sites or regions and an associated IGV views for each site. Reports are created with
the command line script create_report
. Command line arguments are described below.
Although --tracks is optional, a typical report will include at least an alignment track
(BAM or CRAM) file from which the variants were called.
Arguments:
-
Required
- sites VCF (tabix indexed vcf.gz file), BED, MAF, or generic tab delimited file of genomic variant sites.
- fasta Reference fasta file; must be indexed.
-
Required for generic tab delimited sites file
- --begin INT. Column of start chromosomal position for sites file. Used for generic tab delimited input.
- --end INT. Column of end chromosomal position for sites. Used for generic tab delimited input.
- --sequence INT. Column of sequence (chromosome) name.
-
Optional for generic tab delimited sites file
- --zero-based Specify that the position in the sites file is 0-based (e.g. UCSC files) rather than 1-based. Default is
false
.
- --zero-based Specify that the position in the sites file is 0-based (e.g. UCSC files) rather than 1-based. Default is
-
Optional
- --tracks LIST. Space-delimited list of track files, see below for supported formats. If both tracks and track-config are specified tracks will appear first by default.
- --track-config FILE. File containing array of json configuration objects for igv.js tracks. See the igv.js wiki for more details. This option allows customization of track parameters. When using this option, the track
url
andindexURL
properties should be set to the paths to the respective files. - --ideogram FILE. Ideogram file in UCSC cytoIdeo format.
- --template FILE. HTML template file.
- --output FILE. Output file name; default="igvjs_viewer.html".
- --info-columns LIST. Space delimited list of field names to includ in the variant table. If _sites is a VCF file these are the info field names. If sites is a tab delimited format these are column names.
- --info-columns-prefixes LIST. Space delimited list of prefixes of VCF info field names to include in variant table.
- --sample-columns LIST. Space delimited list of VCF sample/format field names to include in variant table.
- --flanking INT. Genomic region to include either side of variant; default=1000.
- --standalone Embed all JavaScript referenced via
<script>
tags in the page. - --sort Applies to alignment racks only. If specified alignments are initally sorted by the specified option. Supported values include
BASE, STRAND, INSERT_SIZE, MATE_CHR, and NONE
. Default value isBASE
for single nucleotide variants,NONE
(no sorting) otherwise. See the igv.js documentation for more information.
**Tab delimited sites file
Variant sites can be defined from a VCF,
UCSC BED, or a generic tab delimited file.
Track file formats:
Currently supported track file formats are BAM, CRAM, VCF, BED, GFF3, and GTF. FASTA. BAM, CRAM, and VCF files must be indexed. Tabix is supported for other file types and it is recommended that all large files be indexed.
Examples
Data for the examples are available for download.
Creating a variant report from a VCF file:
create_report examples/variants/variants.vcf.gz https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa --ideogram examples/variants/cytoBandIdeo.txt --flanking 1000 --info-columns GENE TISSUE TUMOR COSMIC_ID GENE SOMATIC --tracks examples/variants/variants.vcf.gz examples/variants/recalibrated.bam examples/variants/refGene.sort.bed.gz --output example1.html
Creating a variant report from a "track-config" json file
create_report examples/variants/variants.vcf.gz https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa --ideogram examples/variants/cytoBandIdeo.txt --flanking 1000 --info-columns GENE TISSUE TUMOR COSMIC_ID GENE SOMATIC --track-config examples/variants/trackConfigs.json --output example_config.html
Creating a variant report from a TCGA MAF file
create_report examples/variants/tcga_test.maf https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg19/hg19.fasta --ideogram examples/variants/cytoBandIdeo.txt --flanking 1000 --info-columns Chromosome Start_position End_position Variant_Classification Variant_Type Reference_Allele Tumor_Seq_Allele1 Tumor_Seq_Allele2 dbSNP_RS --tracks examples/variants/refGene.sort.bed.gz --output example_maf.html
Creating a variant report from a generic tab-delimited file
create_report examples/variants/test.maflite.tsv https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg19/hg19.fasta --ideogram examples/variants/cytoBandIdeo.txt --flanking 1000 --sequence 1 --begin 2 --end 3 --info-columns chr start end ref_allele alt_allele --tracks examples/variants/refGene.sort.bed.gz --output example_tab.html
Creating a junction report from a splice-junction bed file
create_report examples/junctions/Introns.38.bed https://s3.dualstack.us-east-1.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa --type junction --ideogram examples/junctions/cytoBandIdeo.txt --output junctions.html --track-config examples/junctions/tracks.json --info-columns TCGA GTEx variant_name --title "Sample A"
Converting genomic files to data URIs for use in igv.js
The script create_datauri
converts the contents of a file to a data uri for use in igv.js. The datauri will be
printed to stdout. NOTE It is not neccessary to run this script explicitly to create a report, it is documented here
for use with stand-alone igv.js.
Convert a gzipped vcf file to a datauri.
create_datauri examples/variants/variants.vcf.gz
Convert a slice of a remote cram file to a datauri.
create_datauri \
--region 8:127,738,322-127,738,508 \
https://s3.amazonaws.com/1000genomes/data/HG00096/alignment/HG00096.alt_bwamem_GRCh38DH.20150718.GBR.low_coverage.cram
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for igv_reports-1.0.3-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 667f11da3fbfafb0e89425ad2d183565f790622fa30857767a01ca47dee0097e |
|
MD5 | 8e32320d14f7b56a7aa1485824dbeca1 |
|
BLAKE2b-256 | 36198b1e270d79a3838be53959e4cbc4e4003dc94df72b92d97060b8f94a227e |