CLIPSeqTools
view release on metacpan or search on metacpan
lib/CLIPSeqTools/PlotApp/nucleotide_composition.pm view on Meta::CPAN
ylab="Percent of nucleotides (%)", main="Nucleotide composition per
read position", cex=0.5)});
$R->run(q{lines(idata$position, idata$C_percent, type="o",
col="orange", cex=0.5)});
$R->run(q{lines(idata$position, idata$G_percent, type="o",
col="lightblue", cex=0.5)});
$R->run(q{lines(idata$position, idata$T_percent, type="o",
col="darkblue", cex=0.5)});
$R->run(q{lines(idata$position, idata$N_percent, type="o", col="black",
cex=0.5)});
$R->run(q{legend( "topright" , c("A", "C", "G", "T", "N"),
fill=c("darkred", "orange", "lightblue", "darkblue", "black"))});
$R->run(q{graphics.off()});
# Close R
$R->stop();
}
1;
( run in 1.313 second using v1.01-cache-2.11-cpan-49f99fa48dc )