diff --git a/.Rprofile b/.Rprofile
deleted file mode 100644
index 81b960f..0000000
--- a/.Rprofile
+++ /dev/null
@@ -1 +0,0 @@
-source("renv/activate.R")
diff --git a/R/snpeff_fun.R b/R/snpeff_fun.R
index d8f5367..50bf5b9 100644
--- a/R/snpeff_fun.R
+++ b/R/snpeff_fun.R
@@ -115,7 +115,10 @@ snpeffr <- function(vcf_path,
 
     # get just the samples with mutations
     vcf_long <- vcf_long[value == allele | value == ""]
-    vcf_long[, allele := ifelse(value == "", "undefined", allele)]
+    vcf_long[, allele := ifelse(value == "", "", allele)]
+    vcf_long[, HGVS.p := ifelse(value == "", "undetermined", HGVS.p)]
+    vcf_long[, note := ifelse(value == "", "For this run, MycoSNP was unable to determine the presence/absence of a mutation at this position.", "")]
+    vcf_long <- unique(vcf_long, by = colnames(vcf_long))
     posits <- data.table(region = names(posits), POS = posits)
     vcf_long <- posits[vcf_long, on = "POS"]
 
@@ -124,7 +127,7 @@ snpeffr <- function(vcf_path,
                      "mutation"))
 
     return(vcf_long[, c("sample_id", "snpeff_gene_name", "region", "position", "mutation",
-                        "ref_sequence", "sample_sequence")])
+                        "ref_sequence", "sample_sequence", "note")])
   }
 
 
diff --git a/combined.snpeff (1).vcf.gz b/combined.snpeff (1).vcf.gz
new file mode 100644
index 0000000..c4d7bfa
Binary files /dev/null and b/combined.snpeff (1).vcf.gz differ