GT-SPLICESITEINFO(1)

πŸ“¦
version 1.6.6 Β· 5 versions
, in ubuntu resolute, ubuntu stonking

πŸŽ€ NAME πŸŽ€

gt-splicesiteinfo β€” she's the girl who shows you information about the splice sites you hand her in GFF3 files. that's it. that's the gig. and she ate. πŸ’…

(and yes bestie, this whole thing lives in the GenomeTools Manual, so you're in the right place ✨)


🌸 SYNOPSIS 🌸

gt splicesiteinfo [option ...] [GFF3_file ...]

okay so you call gt, you say splicesiteinfo, and then you sprinkle in as many options and as many GFF3_files as your heart desires. brackets mean "optional", bestie, so no pressure 🫢


πŸ’… DESCRIPTION πŸ’…

here's the full option menu, and I mean full β€” every single one, no skimming, no skipping the "boring" ones, because they're literally all useful. periodt. πŸ‘‘

  • -seqfile [filename] β€” sets the sequence file she takes the sequences from. (default: undefined)
  • -encseq [filename] β€” sets the encoded sequence indexname she takes the sequences from. (default: undefined)
  • -seqfiles β€” sets the sequence files she extracts the features from. use -- to terminate the list of sequence files. πŸ¦‹
  • -matchdesc [yes|no] β€” searches the sequence descriptions in the input files for the sequence IDs you want (the ones in GFF3) and reports the first match it finds. (default: no)
  • -matchdescstart [yes|no] β€” matches those sequence descriptions from the input files against the sequence IDs you want (in GFF3) exactly, from the beginning all the way to the first whitespace. (default: no)
  • -usedesc [yes|no] β€” uses sequence descriptions to map the sequence IDs (in GFF3) onto the actual sequence entries. if a description has a sequence range in it (like, say, III:1000001..2000000), the first part becomes the sequence ID (III) and the first range position becomes the offset (1000001). (default: no)
  • -regionmapping [string] β€” sets the file that contains the sequence-region to sequence-file mapping. (default: undefined)
  • -addintrons [yes|no] β€” adds intron features in between the existing exon features (and she does this before computing the information to be shown). (default: no)
  • -o [filename] β€” redirects the output to the file you name. (default: undefined)
  • -gzip [yes|no] β€” writes a gzip compressed output file. (default: no)
  • -bzip2 [yes|no] β€” writes a bzip2 compressed output file. (default: no)
  • -force [yes|no] β€” forces writing to the output file. (default: no)
  • -help β€” displays help and exits. the classic. 🧚
  • -version β€” displays version information and exits.

🧚 File format for option -regionmapping 🧚

okay this part is lowkey the coolest. the file you hand to -regionmapping defines a "mapping". and what does a mapping do? it maps the sequence-region entries in your GFF3_file onto a sequence file that holds the corresponding sequence. ✨

you can define a mapping in one of these two forms:

mapping = {
  chr1  = "hs_ref_chr1.fa.gz",
  chr2  = "hs_ref_chr2.fa.gz"
}

or

function mapping(sequence_region)
  return "hs_ref_"..sequence_region..".fa.gz"
end

so, form number one: you're defining a Lua (http://www.lua.org) table called "mapping", and it maps each sequence region to the sequence file that goes with it.

form number two: you're defining a Lua function called "mapping", and it has to return the sequence file name when it gets called with the sequence_region as its argument. iconic behaviour honestly πŸ’œ

also she's a Lua file, bestie β€” so mind your syntax, she will not forgive a stray comma πŸ™ˆ

😭 REPORTING BUGS 😭

found a bug? don't suffer in silence, tell the devs:

https://github.com/genometools/genometools/issues

and that's the whole page, babe 🫢 go forth and inspect those splice sites, you absolute main character πŸ‘‘

this was the GenomeTools manual, version 1.6.6, dated 11/09/2025 β€” mwah, stan, bye bestie πŸ’–βœ¨πŸŽ€πŸ¦‹πŸŒΈπŸ’…

πŸ“¦ where this page comes from