Update Accessory_functions_consensusNanoMod.R

biocorecrg · Jun 25, 2024 · 6ded4a1 · 6ded4a1
1 parent 6db335f
commit 6ded4a1
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/mop_consensus/bin/scripts/Accessory_functions_consensusNanoMod.R b/mop_consensus/bin/scripts/Accessory_functions_consensusNanoMod.R
@@ -29,7 +29,7 @@ epinano_processing <- function(sample_file, ivt_file, initial_position, final_po
   sample <- subset(sample, pos<=final_position)
   sample$reference <- paste(sample$X.Ref, sample$pos, sep='_')
   sample$Difference <- as.numeric(sample$mis)+as.numeric(sample$ins)+as.numeric(sample$del)
-  sample <- sample[,c(1,2,13,12)]
+  sample <- sample[,c(1,2,14,13)]
   colnames(sample) <- c('Reference', 'Position', 'Difference_sample', 'Merge')
 
   ivt <- read_csv_file(ivt_file)
@@ -38,7 +38,7 @@ epinano_processing <- function(sample_file, ivt_file, initial_position, final_po
   ivt <- subset(ivt, pos<=final_position)
   ivt$reference <- paste(ivt$X.Ref, ivt$pos, sep='_')
   ivt$Difference <- as.numeric(ivt$mis)+as.numeric(ivt$ins)+as.numeric(ivt$del)
-  ivt <- ivt[,c(1,2,13,12)]
+  ivt <- ivt[,c(1,2,14,13)]
   colnames(ivt) <- c('Reference', 'Position', 'Difference_IVT', 'Merge')
 
   if (nrow(sample)!=0 && nrow(ivt)!=0) {