diff --git a/variants/annotvcf.sh b/variants/annotvcf.sh index c538256f148a8b31c92b8513dd670caeb27ad33c..d7626aed2426cf65586831426c333bee12eb5161 100755 --- a/variants/annotvcf.sh +++ b/variants/annotvcf.sh @@ -23,8 +23,11 @@ done shift $(($OPTIND -1)) -source /etc/profile.d/modules.sh -module load bedtools/2.26.0 snpeff/4.3q +if [[ -z $isdocker ]] +then + source /etc/profile.d/modules.sh + module load bedtools/2.26.0 snpeff/4.3q +fi if [[ $index_path == '/project/shared/bicf_workflow_ref/human/grch38_cloud/rnaref' ]] then diff --git a/variants/uniform_vcf_gt.pl b/variants/uniform_vcf_gt.pl index 85b8e06858e4731489ea082957bbf2049e4f821b..d73aad444703cdc2bddd1bcf62801b85bf10109d 100755 --- a/variants/uniform_vcf_gt.pl +++ b/variants/uniform_vcf_gt.pl @@ -93,6 +93,7 @@ while (my $line = <VCF>) { $gtdata{DP} = $gtdata{NR}; $gtdata{AO} = $gtdata{NV}; $gtdata{RO} = $gtdata{DP} - $gtdata{AO}; + $gtdata{AD} = join(',',$gtdata{RO},$gtdata{AO}); } elsif (exists $gtdata{AO} && exists $gtdata{RO}) { $gtdata{AD} = join(',',$gtdata{RO},$gtdata{AO}); $gtdata{DP} = $gtdata{RO};