diff --git a/genect_rnaseq/dea.R b/genect_rnaseq/dea.R
index d2509a97415eee90ca9a9704e09f5fec12c65f6a..4cd741973a805edf47f941a0ac086f7383d91476 100755
--- a/genect_rnaseq/dea.R
+++ b/genect_rnaseq/dea.R
@@ -118,8 +118,11 @@ for (i in 1:a) {
 
 ###################### Run EdgeR ########################
 ###################### Run QuSage ######################
-
-MSIG.geneSets <- read.gmt('geneset.gmt')
+run.qusage <- 0
+if (file.exists('geneset.gmt')) {
+   MSIG.geneSets <- read.gmt('geneset.gmt')
+   run.qusage <- 1
+}
 
 design <- model.matrix(~grps)
 d <- DGEList(counts=countTable,group=grps,lib.size=libSizes)
@@ -166,10 +169,12 @@ for (i in 1:a) {
       heatmap.2(zscores2, col = bluered(100),Rowv = as.dendrogram(STREE), RowSideColors = col.blocks,dendrogram='row', cexCol=textscale,labCol=s$symbol,srtRow=45,srtCol=45,trace="none", margins=c(5, 5))
       legend("topright",legend=grpnames,col=rainbow(length(grpnames)),pch=20,cex=0.5)
       dev.off()
-      gcont <- paste(cond[j],cond[i],sep='-')
-      qs.results = qusage(logcpm, grps,gcont,MSIG.geneSets)
-      save(qs.results,file=paste(cond[i],'_',cond[j],'.qusage.rda',sep=""))
-      }
+      if (run.qusage > 0) {
+      	 gcont <- paste(cond[j],cond[i],sep='-')
+      	 qs.results = qusage(logcpm, grps,gcont,MSIG.geneSets)
+      	 save(qs.results,file=paste(cond[i],'_',cond[j],'.qusage.rda',sep=""))
+	 }
+	 }
       }
   }
 }
diff --git a/genect_rnaseq/statanal.sh b/genect_rnaseq/statanal.sh
index 9cba1e37976c079754bd830ebf1b38acf12dec03..fd36f0f19d745ac77bb35cec936cbdb4a6100664 100644
--- a/genect_rnaseq/statanal.sh
+++ b/genect_rnaseq/statanal.sh
@@ -41,5 +41,8 @@ else
     module load R/3.2.1-intel
     Rscript  $baseDir/dea.R
     Rscript $baseDir/build_ballgown.R *_stringtie
-    perl $baseDir/concat_edgeR.pl *.edgeR.txt
+    if [[ -n `ls *.edgeR.txt` ]]
+    then
+	perl $baseDir/concat_edgeR.pl *.edgeR.txt
+    fi
 fi
diff --git a/variants/union.sh b/variants/union.sh
index 8ecfe87dbeefad6f691385763a1c9315b6ff2854..a02cd2c2dbd7e3d3b0c95094704bd8cec43cb670 100755
--- a/variants/union.sh
+++ b/variants/union.sh
@@ -48,6 +48,8 @@ for i in ${dir}/*.vcf.gz; do
 done 
 
 echo "##fileformat=VCFv4.2" > header.vcf
+echo "##INFO=<ID=CallSet,Number='.',Type=String,Description=\"Caller\">" >> header.vcf
+echo "##INFO=<ID=CallSetInconsistent,Number=1,Type=Integer,Description=\"Caller Inconsistency Manual Inspection\">" >> header.vcf
 zcat ${dir}/*.vcf.gz |grep "##" |grep -v '#fileformat' |sort -u |grep 'ALT\|FILTER\|FORMAT\|INFO' >> header.vcf
 
 perl $baseDir/unionvcf.pl header.vcf $list2