GT-SNPPER(1)

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

πŸ’– GT-SNPPER(1) πŸ’–

okay bestie, welcome to the GenomeTools Manual page for gt-snpper β€” yes it's a manpage, no it does not have to be boring, let's go ✨


🌸 NAME 🌸

gt-snpper β€” annotates SNPs based on how they affect the genome, as described by a genomic annotation. it's giving "what did this tiny letter swap actually do" energy πŸ’…


✨ SYNOPSIS ✨

gt snpper [option ...] GFF3_file [GVF_file]

πŸ¦‹ DESCRIPTION πŸ¦‹

πŸŽ€ OPTIONS πŸŽ€

  • -trans_table [value] β€” the NCBI translation table number, and you literally get to pick from this whole menu:
    • 1: Standard
    • 2: Vertebrate Mitochondrial
    • 3: Yeast Mitochondrial
    • 4: Mold Mitochondrial; Protozoan Mitochondrial; Coelenterate Mitochondrial; Mycoplasma; Spiroplasma
    • 5: Invertebrate Mitochondrial
    • 6: Ciliate Nuclear; Dasycladacean Nuclear; Hexamita Nuclear
    • 9: Echinoderm Mitochondrial; Flatworm Mitochondrial
    • 10: Euplotid Nuclear
    • 11: Bacterial, Archaeal and Plant Plastid
    • 12: Alternative Yeast Nuclear
    • 13: Ascidian Mitochondrial
    • 14: Alternative Flatworm Mitochondrial
    • 15: Blepharisma Macronuclear
    • 16: Chlorophycean Mitochondrial
    • 21: Trematode Mitochondrial
    • 22: Scenedesmus obliquus Mitochondrial
    • 23: Thraustochytrium Mitochondrial
    • 24: Pterobranchia Mitochondrial
    • 25: Candidate Division SR1 and Gracilibacteria (default: 1)
  • -seqfile [filename] β€” tells it which sequence file to pull the sequences from (default: undefined) 🧬
  • -encseq [filename] β€” sets the encoded sequence indexname to pull the sequences from (default: undefined)
  • -seqfiles β€” sets the sequence files to extract 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 (in GFF3), and reports the first match (default: no)
  • -matchdescstart [yes|no] β€” matches the sequence descriptions from the input files against the desired sequence IDs (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 actual sequence entries. if a description has a sequence range in it (like III:1000001..2000000), then 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) πŸ—ΊοΈ
  • -o [filename] β€” redirects output into the file you specify (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 β€” shows help and dips ✌️
  • -version β€” shows version info and dips ✌️

🌈 File format for option -regionmapping 🌈

the file you hand to option -regionmapping defines a β€œmapping”. a mapping takes the sequence-region entries in the GFF3_file and maps them to a sequence file that holds the matching sequence. you can define mappings in one of these two forms, pick your fighter:

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

the first form defines a Lua (<http://www.lua.org>) table named β€œmapping” that maps each sequence region to its corresponding sequence file. the second one defines a Lua function called β€œmapping”, and that function has to return the sequence file name when you call it with the sequence_region as its argument. iykyk πŸ’«


😭 REPORTING BUGS 😭

found a bug? spill it here: <https://github.com/genometools/genometools/issues>.


ngl this page is from GenomeTools 1.6.6, dated 11/09/2025. that's the whole footer bestie, we did it, slay πŸ’…βœ¨πŸ¦‹
πŸ“¦ where this page comes from