diff --git a/R/sqtl.seeker.R b/R/sqtl.seeker.R
index f017889..7351ff8 100644
--- a/R/sqtl.seeker.R
+++ b/R/sqtl.seeker.R
@@ -131,7 +131,7 @@ sqtl.seeker <- function(tre.df, genotype.f, gene.loc, covariates = NULL,
                 covariates <- covariates[com.samples, , drop = FALSE]
                 multiclass <- apply(covariates, 2, function(x){length(table(x)) > 1})
                 covariates <- covariates[, multiclass, drop = FALSE]
-                if (verbose){
+                if (verbose & any(!multiclass)){
                   message("\t", "Covariates removed due to only one value: ", 
                           paste(names(multiclass)[!multiclass], collapse = ", "))
                 }
diff --git a/data/genes.bed b/data/genes.bed
index 9e81ffa..598d793 100644
--- a/data/genes.bed
+++ b/data/genes.bed
@@ -1,3 +1,4 @@
+chr	start	end	geneId
 19	38826415	38861589	ENSG00000099338.14
 19	52901102	52921657	ENSG00000167555.8
 19	52956829	53016982	ENSG00000258405.3