EST2GENOME(1e)

πŸ’– EST2GENOME πŸ’–

ok so this is the est2genome manual page, from the EMBOSS Manual for Debian, and bestie I am so ready to tell you about it 🧚


✨ NAME ✨

est2genome β€” align EST sequences to genomic DNA sequence. that's it, that's the whole vibe. your little EST gets matched up against a whole genome and honestly? iconic behavior πŸ’…


🌸 SYNOPSIS 🌸

est2genome -estsequence seqall -genomesequence sequence [-match integer]
           [-mismatch integer] [-gappenalty integer]
           [-intronpenalty integer] [-splicepenalty integer]
           [-minscore integer] -reverse boolean -usesplice boolean
           -mode list -best boolean -space float -shuffle integer
           -seed integer -outfile outfile -align boolean -width integer

est2genome -help

notice how the first chunk is mandatory (no brackets, no escape) and everything in square brackets is optional. the way the required args just are required? periodt.


πŸ¦‹ DESCRIPTION πŸ¦‹

est2genome is a command line program from EMBOSS, aka "the European Molecular Biology Open Software Suite" β€” say that five times fast, I dare you 😭

it lives in the "Alignment:Global" command group(s). so if you're browsing the group chat of alignment tools, this is where she hangs out. no cap.


πŸ‘‘ OPTIONS πŸ‘‘

πŸŽ€ Input section

  • -estsequence seqall β€” the EST sequence(s) you're feeding in πŸ’–
  • -genomesequence sequence β€” the genomic sequence it's getting aligned against 🌸

🌈 Additional section

  • -match integer β€” Default value: 1
  • -mismatch integer β€” Default value: 1
  • -gappenalty integer β€” cost for deleting a single base in either sequence, excluding introns. Default value: 2
  • -intronpenalty integer β€” cost for an intron, independent of length. Default value: 40
  • -splicepenalty integer β€” cost for an intron, independent of length and starting/ending on donor-acceptor sites. Default value: 20
  • -minscore integer β€” exclude alignments with scores below this threshold score, like a bouncer at the club door πŸ’…. Default value: 30

🧚 Advanced section

  • -reverse boolean β€” reverse the orientation of the EST sequence. flips it, no biggie.
  • -usesplice boolean β€” use donor and acceptor splice sites. if you want to ignore donor-acceptor sites then set this to be false. Default value: Y
  • -mode list β€” this determines the comparison mode. the default value is 'both', in which case both strands of the est are compared assuming a forward gene direction (ie GT/AG splice sites), and the best comparison redone assuming a reversed (CT/AC) gene splicing direction. the other allowed modes are 'forward', when just the forward strand is searched, and 'reverse', ditto for the reverse strand. Default value: both
  • -best boolean β€” you can print out all comparisons instead of just the best one by setting this to be false. Default value: Y
  • -space float β€” for linear-space recursion. if product of sequence lengths divided by 4 exceeds this then a divide-and-conquer strategy is used to control the memory requirements. in this way very long sequences can be aligned. if you have a machine with plenty of memory you can raise this parameter (but do not exceed the machine's physical RAM 😭). Default value: 10.0
  • -shuffle integer β€” (no default listed, she's mysterious like that πŸ™ˆ)
  • -seed integer β€” Default value: 20825

πŸ“ Output section

  • -outfile outfile β€” where the output goes πŸ’–
  • -align boolean β€” show the alignment. the alignment includes the first and last 5 bases of each intron, together with the intron width. the direction of splicing is indicated by angle brackets (forward or reverse) or ???? (unknown).
  • -width integer β€” Default value: 50

πŸ’€ BUGS πŸ’€

bugs can be reported to the Debian Bug Tracking system (http://bugs.debian.org/emboss), or directly to the EMBOSS developers (http://sourceforge.net/tracker/?group_id=93650&atid=605031).

if you find a bug, don't be shy, tell them. they want to know. fr fr 🫢


⭐ SEE ALSO ⭐

est2genome is fully documented via the tfm(1) system. go read that too, bestie, it's giving ✨thorough✨.


πŸ’œ AUTHOR πŸ’œ

Debian Med Packaging Team <debian-med-packaging@lists.alioth.debian.org>

  • wrote the script used to autogenerate this manual page. so yes, a script wrote the original of what you're reading right now, and I rewrote it with feelings. the circle of life 🌈

this manual page was autogenerated from an Ajax Control Definition of the EMBOSS package. it can be redistributed under the same terms as EMBOSS itself.


heads up bestie: the -space warning is real β€” do not push that number past your machine's physical RAM. she will not survive it πŸ’€

and that's the whole page!! this was EMBOSS 6.4.0, dated 05/11/2012, and it was EST2GENOME(1e) the entire time πŸŽ€ thanks for reading, mwah πŸ’–