-
Notifications
You must be signed in to change notification settings - Fork 4
/
README
38 lines (23 loc) · 2.37 KB
/
README
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
# Copyright 2014 Haley Abel
# This file is part of Copycat2.
# Copycat2 is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
# Copycat2 is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
# You should have received a copy of the GNU General Public License
# along with Copycat2. If not, see <http://www.gnu.org/licenses/>.
copycat2--Normalizes targeted capture sequence coverage against a pooled normal control for CNV detection and plotting. Uses variant allele frequecies to aid in normalization and provide support for CNV calls.
Required input files:
test_cov=case coverage, output from BEDTools coverage bed with target-level (exon or bait) read counts (should be in 7 column format obtained from using 3-col bed)
test_vcf=vcf for case (optional but strongly recommended)
ctl_info=list of control coverage files 1 file per line, with absolute file paths. control coverage files should be in the same format as for cases.
target.bedfile=bedfile (3 column) used to create coverage files
all other files are as described in R documentation.
Example:
library(CC2)
copycat2(test_cov="/gscuser/habel/programs/data/cov/tumor/AML11.RMG.tumor.bam.cov.3.txt", test_vcf="/gscuser/eduncava/my_space/CNV/RMG/vcfs/AML11.germline.RMG.VCF", ctl_info="/gscuser/habel/programs/data/cov/normal/new.info.txt", outdir="/gscuser/habel/programs/data/workdir/test0909",
outprefix="AML11_1", target.bedfile="/gscuser/habel/programs/data/workdir/baits.3col.bed", fa="/gscuser/eduncava/MedSeq/CNV/build37.fa", vcftype="VARSCAN", use_af=TRUE, min_coverage=50, plot.only=FALSE, annot.bedfile="/gscuser/eduncava/MedSeq/CNV/hg19.annot.bed", snp.wingspan=0, plot.bedfile=NULL, plot.starts=1, plot.stops=22, p.af.cutoff=1e-4, p.cov.cutoff=1e-4, af.threshold=0.58, png=FALSE, coverage.min.ratio=0.75, coverage.max.ratio=1.25, minwid=5, segalpha=0.01, sdundo=2, restrict="0.125,6", last.autosome=22, frac.normal=0.5, thin.rate=NULL, plot.id="", min.var.cov=20, vafs_normalize=TRUE, min.normal.corr=0.8, min.num.normals=3, clean.up=TRUE, read_len=100)