Skip to content
Snippets Groups Projects
Commit 177d794f authored by Venkat Malladi's avatar Venkat Malladi
Browse files

If negative fragmentlength default to 200.

parent cf41e35b
Branches
Tags
1 merge request!62Resolve "callPeaksMACS exits when data is bad"
Pipeline #2842 failed with stages
......@@ -108,6 +108,8 @@ def call_peaks_macs(experiment, xcor, control, prefix, genome_size, chrom_sizes)
frag_lengths = firstline.split()[2] # third column
fragment_length = frag_lengths.split(',')[0] # grab first value
logger.info("Fraglen %s", fragment_length)
if int(fragment_length) < 1:
fragment_length = "200"
# Generate narrow peaks and preliminary signal tracks
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment