GAEVAL(1)

✨ gaeval ✨

okay bestie, welcome to gaeval β€” the little tool that calculates coverage and integrity scores for gene models based on transcript alignments. that's the whole vibe. it looks at how well your gene models line up with your transcript alignments and hands you numbers. πŸ’–


🌸 SYNOPSIS 🌸

gaeval [options] alignments.gff3 genes.gff3 [moregenes.gff3 ...]

so you pass it your alignments file first, then one or more gene files, and it does its thing. no cap. ✨


πŸ’… DESCRIPTION πŸ’…

here's everything it can do, section by section, because organization is self-care. πŸŽ€

🫢 Basic options

  • -h, --help β€” prints this help message and then dips. exits right after, she's busy. πŸ’…
  • -v, --version β€” prints the version number and exits. iconic behavior. ✨

⭐ Weights for calculating integrity score

okay so important: these weights must add up to 1.0. like, they literally have to. that's the rule. πŸ™ˆ

  • -a, --alpha: DOUBLE β€” introns confirmed, or % expected CDS length for single-exon genes; default is 0.6. πŸ’–
  • -b, --beta: DOUBLE β€” exon coverage; default is 0.3. πŸ¦‹
  • -g, --gamma: DOUBLE β€” % expected 5' UTR length; default is 0.05. 🌸
  • -e, --epsilon: DOUBLE β€” % expected 3' UTR length; default is 0.05. πŸŽ€

πŸ‘‘ Expected feature lengths for calculating integrity score

these are your "what even counts as normal around here" numbers. πŸ“

  • -c, --exp-cds: INT β€” expected CDS length (in bp); default is 400. ⭐
  • -5, --exp-5putr: INT β€” expected 5' UTR length; default is 200. 🌈
  • -3, --exp-3putr: INT β€” expected 3' UTR length; default is 100. πŸ’œ

heads up bestie: the weights are the one place you cannot freestyle. -a, -b, -g and -e have to sum to exactly 1.0 or it's not gonna be your moment. 🫢

and that's the whole page, fr fr. this is gaeval version 0.15.1, dated January 2016 β€” a lil vintage, still slaying. thanks for reading my TED talk, mwah πŸ’‹