###########GRAFIC SUBSEQ5 + CFTR +EXONS###########

# dibuixo -> sequencia.gccontent_combo_zoom5.png
mysegments <- function(x,y) {
  polygon(c(x,y,y,x),
          c(0.185,0.185,0.225,0.225),
          col="blue",border="blue",lty="solid",fg="blue",bg="blue");
}
t1 <- read.table("subseq5.gccontent_500.tbl");
png("sequencia.gccontent_combo_zoom5.png",width=4000,height=500);
plot(t1,type="l",col="black",xlim=c(0,350000),ylim=c(0.2,0.8),axes=FALSE);
axis(1,seq(0,350000,25000),col="black",las=1);
axis(2,seq(0.20,0.80,0.05),col="black",las=1);
segments(-50000,0.38,400000,0.38,col="red",lty="dotted");
segments(109628,0.205,296641,0.205,col="blue",lty="solid");
mysegments(109628,109680);
mysegments(133786,133896);
mysegments(138567,138675);
mysegments(160432,160647);
mysegments(163809,163898);
mysegments(164781,164944);
mysegments(166081,166206);
mysegments(169633,169879);
mysegments(171549,171641);
mysegments(178174,178356);
mysegments(188997,189188);
mysegments(217272,217366);
mysegments(219886,219972);
mysegments(221467,222190);
mysegments(224463,224591);
mysegments(232359,232396);
mysegments(233065,233315);
mysegments(236207,236286);
mysegments(240052,240202);
mysegments(241114,241341);
mysegments(244146,244246);
mysegments(257055,257303);
mysegments(271971,272126);
mysegments(282375,282464);
mysegments(294221,294393);
mysegments(294992,295097);
mysegments(296441,296641);
dev.off();