Skip to content
Snippets Groups Projects
Commit 4ccabd55 authored by CAETANO Tomas's avatar CAETANO Tomas
Browse files

fix bug EMOTE_parse_reads pc_valid

parent 1d1cc579
No related branches found
No related tags found
No related merge requests found
......@@ -284,7 +284,7 @@ EMOTE_parse_read = function(fastq_file,
}
close(fq_streamer)
# Write parse report
stat_tb = mutate(stat_tb, pc_valid = is_valid/total_read) %>%
stat_tb = mutate(stat_tb, pc_valid = ifelse(group=="VALID", total_read/sum(stat_tb$total_read),0)) %>%
mutate(demux_filename = paste0(out_dir, "/", fq_basename,"_", tolower(group), ".fastq.gz"))
if(stat_tb[stat_tb$group=="VALID",]$pc_valid == 0){
warning("There is no valid reads, check your EMOTE_read_features to ensure it is correct (an empty file was created anyway)")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment