remove sequences below threshold (not above)

loosolab · Dec 4, 2018 · 96efc46 · 96efc46
1 parent 966761f
commit 96efc46
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/bin/reduce_bed.R b/bin/reduce_bed.R
@@ -54,7 +54,7 @@ reduce_bed <- function(input, kmer = 10, motif = 10, output = "reduced.bed", thr
 
   # remove sequences below minimum length
   total_rows <- nrow(bed_table)
-  bed_table <- bed_table[nchar(sequence) < min_seq_length]
+  bed_table <- bed_table[nchar(sequence) > min_seq_length]
   if (total_rows > nrow(bed_table)) {
     message("Removed ", total_rows - nrow(bed_table), " sequence(s) below minimum length of ", min_seq_length)
   }