#######GRAFIC G+C#######


t <- read.table("sequencia.gccontent_500.tbl");
str(t);
plot(t,type="h");
segments(-100000,0.40,2000000,0.40,col="red");
t <- read.table("sequencia.gccontent_1000.tbl");
str(t);
plot(t,type="h");
segments(-100000,0.40,2000000,0.40,col="red");
t <- read.table("sequencia.gccontent_2000.tbl");
str(t);
plot(t,type="h");
segments(-100000,0.40,2000000,0.40,col="red");
t <- read.table("sequencia.gccontent_5000.tbl");
str(t);
plot(t,type="h");
segments(-100000,0.40,2000000,0.40,col="red");
t <- read.table("sequencia.gccontent_10000.tbl");
str(t);
plot(t,type="h");
segments(-100000,0.40,2000000,0.40,col="red");

# dibuixo -> sequencia.gccontent_combo.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("sequencia.gccontent_500.tbl");
t2 <- read.table("sequencia.gccontent_1000.tbl");
t3 <- read.table("sequencia.gccontent_2000.tbl");
t4 <- read.table("sequencia.gccontent_5000.tbl");
t5 <- read.table("sequencia.gccontent_10000.tbl");
png("sequencia.gccontent_combo.png",width=4000,height=500);
plot(t1,type="l",col="black",xlim=c(0,2000000),ylim=c(0.2,0.8),axes=FALSE);
axis(1,seq(0,2000000,100000),col="black",las=1);
axis(2,seq(0.20,0.80,0.05),col="black",las=1);
lines(t2,col="blue");
lines(t3,col="red");
lines(t4,col="green");
lines(t5,col="orange");
segments(-100000,0.38,2000000,0.38,col="red",lty="dotted");
segments(1200000,0.205,1550000,0.205,col="purple",lty="solid");
mysegments(1309628,1496641);
dev.off()