diff --git a/variants/cnvkit.sh b/variants/cnvkit.sh index a914b342ad63d2ee687202ba664cb095c0810342..4b2513cf04523fe7ce03cc963cc727d1ae7232d1 100755 --- a/variants/cnvkit.sh +++ b/variants/cnvkit.sh @@ -61,5 +61,5 @@ cnvkit.py coverage ${sbam} ${targets}antitargets.bed -o ${pair_id}.antitargetcov cnvkit.py fix ${pair_id}.targetcoverage.cnn ${pair_id}.antitargetcoverage.cnn ${normals} -o ${pair_id}.cnr cnvkit.py segment ${pair_id}.cnr -o ${pair_id}.cns cnvkit.py call ${pair_id}.cns -o ${pair_id}.call.cns -cnvkit.py scatter ${pair_id}.cnr -s ${pair_id}.call.cns -o ${pair_id}.cnv.scatter.pdf +cnvkit.py scatter ${pair_id}.cnr -s ${pair_id}.call.cns -t --segment-color "blue" -o ${pair_id}.cnv.scatter.pdf perl $baseDir/filter_cnvkit.pl *.call.cns diff --git a/variants/unionvcf.pl b/variants/unionvcf.pl index 6ef347583ab7e9685ed94c06efb4861e18699928..ee84efb14a55e0d4d773c5e6d48fa4131b07d179 100755 --- a/variants/unionvcf.pl +++ b/variants/unionvcf.pl @@ -97,7 +97,7 @@ foreach $vcf (@vcffiles) { } my $mafs = '.'; my @maf = (); - if ($gtdata{DP} && $gtdata{DP} ne '.' && $gtdata{AO} && $gtdata{AO} ne '.') { + if ($gtdata{DP} && $gtdata{DP} ne '.' && exists $gtdata{AO} && $gtdata{AO} ne '.') { foreach $areadct (split(/,/,$gtdata{AO})) { push @maf, sprintf("%.2f",$areadct/$gtdata{DP}); }