Skip to content

VarScan #
Find similar titles

Structured data

About
CNV calling from NGS data
SNV calling from NGS data
Date Published
Programming Language
Java
Provider
Washington University

Variant detection in NGS data.

VarScan is a platform-independent mutation caller for targeted, Exome sequencing, WGS data generated on Illumina, SOLiD, Life/PGM, Roche/454, and similar instruments.

It can be used to detect different types of variation:

  • Germline variants (SNPs an dindels) in individual samples or pools of samples.
  • Multi-sample variants (shared or private) in multi-sample datasets (with mpileup).
  • Somatic mutations, LOH events, and germline variants in tumor-normal pairs.
  • Somatic copy number alterations (CNAs) in tumor-normal exome data.

기능 요약 #

Germline variants #

Commands

  • mpileup2snp
  • mpileup2indel
  • mpileup2cns
  • pileup2snp
  • pileup2indel
  • pileup2cns

명령행 예제

$ samtools mpileup -B -f [reference] [bam] | java -jar VarScan.v2.2.jar mpileup2snp

Somatic mutations #

Commands

  • somatic

명령행 예제

$ samtools pileup -f [reference] [normal_bam] > normal.pileup
$ samtools pileup -f [reference] [tumor_bam] > tumor.pileup
$ java -jar VarScan.jar somatic [normal_pileup] [tumor_pileup]

CNV #

Commands

  • copynumber

명령행예제

$ samtools pileup -f [reference] [normal_bam] > normal.pileup
$ samtools pileup -f [reference] [tumor_bam] > tumor.pileup
$ java -jar VarScan.jar copynumber [normal_pileup] [tumor_pileup] [outfile]

Tumor/normal CNA with DNAcopy #

Tumor BAM파일과 정상 BAM 파일을 함께 SAMtools로 mpileup 파일을 만들고, copynumber와 capyCaller 명령을 수행한 후, DNAcopy로 segmentation 한다.

$ samtools mpileup -q 1 -f ref.fa normal.bam tumor.bam | java -jar VarScan.jar copynumber varScan --mpileup 1 --output-file varScan.copynumber
$ java -jar VarScan.jar copyCaller varScan.copynumber --output-file varScan.copynumber.called
$ R
> library(DNAcopy)
> cn <- read.table('varScan.copynumber.called', header=TRUE)
> CNA.object <- CNA(genomdat=cn[,7], chrom=cn[,1], maploc=cn[,2],
                    data.type = 'logratio')
> CNA.smoothed <- smooth.CNA(CNA.object)
> segs <- segment(CNA.smoothed, verbose=0, alpha=0.01, nperm=10000,
                  min.width=2, undo.splits="sdundo", undo.SD=3)
> psegs <- segments.p(segs)
> write.table(psegs, file="varScan.copynumber.called.cbs", 
              row.names=FALSE, col.names=TRUE, quote=FALSE, sep="\t")

VarScan과 함께 제공되는 Perl script 파일 mergeSegments.pl 을 이용하여 머지한다. (스크립트 파일을 일부 수정해야 함. 김형용이 수정한 파일은 mergeSegments-hygkim.pl)

$ perl mergeSegment-hygkim.pl --ref-arm-sizes hg19.arm --output-basename varScam.copynumber.called.cbs varScam.copynumber.called.cbs

위 명령에 사용되는 arm 길이 파일 (hg19.arm)은 다음과 같다.

1    1   125000000   p
1    125000001   249250621   q
2    1   93300000    p
2    93300001    243199373   q
3    1   91000000    p
3    91000001    198022430   q
4    1   50400000    p
4    50400001    191154276   q
5    1   48400000    p
5    48400001    180915260   q
6    1   61000000    p
6    61000001    171115067   q
7    1   59900000    p
7    59900001    159138663   q
8    1   45600000    p
8    45600001    146364022   q
9    1   49000000    p
9    49000001    141213431   q
10   1   40200000    p
10   40200001    135534747   q
11   1   53700000    p
11   53700001    135006516   q
12   1   35800000    p
12   35800001    133851895   q
13   1   17900000    p
13   17900001    115169878   q
14   1   17600000    p
14   17600001    107349540   q
15   1   19000000    p
15   19000001    102531392   q
16   1   36600000    p
16   36600001    90354753    q
17   1   24000000    p
17   24000001    81195210    q
18   1   17200000    p
18   17200001    78077248    q
19   1   26500000    p
19   26500001    59128983    q
20   1   27500000    p
20   27500001    63025520    q
21   1   13200000    p
21   13200001    48129895    q
22   1   14700000    p
22   14700001    51304566    q
X    1   60600000    p
X    60600001    155270560   q
Y    1   12500000    p
Y    12500001    59373566    q

de novo mutations (trios) #

명령행예제

$ samtools mpileup -B -f [reference] [bam] > trio.mpileup
$ java -jar VarScan.v2.2.jar trio trio.mpileup

Incoming Links #

[Medical Scholarly Articles] About #

Related Articles #

Related Education Events #

Related Medical Scholarly Articles #

Suggested Pages #

web biohackers.net
0.0.1_20140628_0