Create sequence tracks for misha database — gseq.create

Create tracks based on the presence of specific sequences.

gseq.create_track(
  s,
  track,
  strand = 1,
  res = NULL,
  intervals = gintervals.all()
)

Arguments

s: sequence to look for. can be a regex.
track: name of the track to create
strand: strand of the sequences
intervals: intervals set

Examples

if (FALSE) {
gseq.create_track("C", "seq.C")
gseq.create_track("T", "seq.T")
gseq.create_track("CG", "seq.CG") # genome CpGs
gseq.create_track("[GC]", "seq.G_or_C") # C or G (for GC content calculations)
}