Update strand information in read.modkit.R

hansenlab · Apr 18, 2024 · 39ac79f · 39ac79f
1 parent 94069c4
commit 39ac79f
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/R/read.modkit.R b/R/read.modkit.R
@@ -11,6 +11,7 @@ read.modkit <- function(files,
     for (i in seq_along(files)){
         data <- read.table(files[i], header = FALSE, sep="\t",
                     stringsAsFactors=FALSE, quote="")
+        data$V6[data$V6 == "."] <- "*"
 
         if (length(unique(data$V4)) == 2){
             gr <- GRanges(seqnames = data[data$V4 == "m", ]$V1,