Genome Variant Analysis

read_NVC.py
Function: This module is used to check the nucleotide composition bias. Due to random priming, certain patterns are over represented at the beginning (5’end) of reads. This bias could be easily examined by NVC (Nucleotide versus cycle) plot. NVC plot is generated by overlaying all reads together, then calculating nucleotide composition for each position of read (or each sequencing cycle). In ideal condition (genome is random and RNA-seq reads is randomly sampled from genome), we expect A%=C%=G%=T%=25% at each position of reads.
Usage: read_NVC.py -i Pairend_nonStrandSpecific_36mer_Human_hg19.bam -o output
java -jar GenomeAnalysisTK.jar
Function: Outputs a list of intervals that are covered to or above a given threshold
Usage: java -jar GenomeAnalysisTK.jar -T FindCoveredIntervals -R reference.fasta -I my_file.bam [-cov 10 \] [-uncovered \] -o output.list
GEMINI interactions
Function: Find genes among variants that are interacting partners.
Usage: gemini interactions -g CTBP2 -r 3 example.db
java -jar GenomeAnalysisTK.jar
Function: Create plots to visualize base recalibration results
Usage: java -jar GenomeAnalysisTK.jar -T AnalyzeCovariates -R myrefernce.fasta -BQSR myrecal.table -plots BQSR.pdf
java -jar GenomeAnalysisTK.jar
Function: Calculates the GC content of the reference sequence for each interval
Usage: java -jar GenomeAnalysisTK.jar -T GCContentByInterval -R reference.fasta -o output.txt -L input.intervals
java -jar GenomeAnalysisTK.jar
Function: Write out sequence read data (for filtering, merging, subsetting etc)
Usage: java -jar GenomeAnalysisTK.jar -T PrintReads -R reference.fasta -I input1.bam -I input2.bam -o output.bam --read_filter MappingQualityZero // Prints the first 2000 reads in the BAM file java -jar GenomeAnalysisTK.jar -T PrintReads -R reference.fasta -I input.bam -o output.bam -n 2000 // Downsamples BAM file to 25% java -jar GenomeAnalysisTK.jar -T PrintReads -R reference.fasta -I input.bam -o output.bam -dfrac 0.25
java -jar GenomeAnalysisTK.jar
Function: Left-align indels in a variant callset
Usage: java -jar GenomeAnalysisTK.jar -T LeftAlignAndTrimVariants -R reference.fasta --variant input.vcf -o output.vcf --splitMultiallelics --dontTrimAlleles --keepOriginalAC
java -jar GenomeAnalysisTK.jar
Function: Convert variants from other file formats to VCF format
Usage: java -jar GenomeAnalysisTK.jar -T VariantsToVCF -R reference.fasta -o output.vcf --variant:RawHapMap input.hapmap
java -jar GenomeAnalysisTK.jar
Function: Haplotype-based resolution of variants in separate callsets.
Usage: java -jar GenomeAnalysisTK.jar -T HaplotypeResolver -R reference.fasta -V:v1 input1.vcf -V:v2 input2.vcf -o output.vcf
java -jar GenomeAnalysisTK.jar
Function: Randomly select variant records according to specified options
Usage: java -jar GenomeAnalysisTK.jar -T ValidationSiteSelectorWalker -R reference.fasta -V input1.vcf -V input2.vcf -sn NA12878 -o output.vcf --numValidationSites 200 -sampleMode POLY_BASED_ON_GT -freqMode KEEP_AF_SPECTRUM
GEMINI autosomal_dominant
Function: Find variants meeting an autosomal dominant model.
Usage: gemini autosomal_dominant test.auto_dom.db --columns "chrom,start,end,gene"
VarScan
Function: Filter somatic mutation calls to remove clusters of false positives and SNV calls near indels. Note: this is a basic filter. More advanced filtering strategies consider mapping quality, read mismatches, soft-trimming, and other factors when deciding whether or not to filter a variant.
Usage: java -jar VarScan.jar somaticFilter [mutations file] OPTIONS
java -jar GenomeAnalysisTK.jar
Function: Convert VCF to binary pedigree file
Usage: java -jar GenomeAnalysisTK.jar -T VariantsToBinaryPed -R reference.fasta -V variants.vcf -m metadata.fam -bed output.bed -bim output.bim -fam output.fam
java -jar GenomeAnalysisTK.jar
Function: Create plots to visualize base recalibration results
Usage: java -jar GenomeAnalysisTK.jar -T AnalyzeCovariates -R myrefernce.fasta -before recal2.table -after recal3.table -plots recalQC.pdf
VarScan
Function: Filter variants in a file by coverage, supporting reads, variant frequency, or average base quality. It is for use with output from pileup2snp or pileup2indel.
Usage: java -jar VarScan.jar filter [variants file] OPTIONS