https://github.com/MSanderlab/Pancreatic-progenitor-epigenome-maps-prioritize-type-2-diabetes-risk-genes-with-roles-in-development
Raw File
Tip revision: ba79c687523c2696ea0ef30d8476e28a0d860f18 authored by MSanderlab on 16 January 2021, 09:01:57 UTC
adding files
Tip revision: ba79c68
hmm-state_domains.pl
#!/usr/bin/perl

use strict;

MAIN : {

    my $prev_state = 0;
    my $prev_start;
    my $prev_end;
    while (my $line = <>) {
	chomp $line;
	my ($chr, $start, $end, $state) = split(/\t/,$line);
	if ($prev_state == 0) {
	    unless ($state == 3) {
		next;
	    }
	}
	if ($state != $prev_state) {
	    if (($state == 3) && ($prev_state == 0)) {
		print $chr . "\t" . $start . "\t";
	    }
	    if (($state == 3) && ($prev_state == 1)) {
		print $prev_end . "\n" . $chr . "\t" . $start . "\t";
	    }

	}
	$prev_state = $state;
	$prev_start = $start;
	$prev_end = $end;
    }
    if ($prev_state ==1 ) {
	print $prev_end . "\n";
    }

}

back to top