# HG changeset patch # User geert-vandeweyer # Date 1423752843 18000 # Node ID 6cb012c8497adab82bdf3642c1c35b5dc2108f06 # Parent fd788f9db899aed260340a826be69d431075684e Added BED format check before collapsing regions. diff -r fd788f9db899 -r 6cb012c8497a CoverageReport.pl --- a/CoverageReport.pl Thu Feb 12 08:51:37 2015 -0500 +++ b/CoverageReport.pl Thu Feb 12 09:54:03 2015 -0500 @@ -78,6 +78,13 @@ ## 0. Collapse overlapping target regions. if (defined($opts{'T'})) { + ## check BED format. Must have 6 cols if using this. + my $head = `head -n 1 $opts{'t'}`; + chomp; + my @c = split(/\t/,$head); + if (scalar(@c) < 6) { + die("Targets BED file must be in 6-column format for collapsings. See tool documentation for more info.\n"); + } my $targets = $opts{'t'}; my $tmptargets = "$wd/collapsedtargets.bed"; system("sort -k1,1 -k2,2n $targets > $wd/sorted.targets.bed");