TRF(1)

πŸ’– TRF β€” like, the tandem repeat finder?? yeah bestie, it's a User Command πŸ’–

okay so. hi. you're looking at the trf(1) page, and honestly? iconic of you. trf is a little program that goes snooping through DNA and finds tandem repeats, then shows them to you. that's the whole vibe. it's giving "bioinformatics but make it ✨detective✨". let's get into it πŸ§¬πŸ¦‹


✨ NAME ✨

trf β€” locate and display tandem repeats in DNA sequences πŸ’…


🌸 SYNOPSIS 🌸

trf (File|Match|Mismatch|Delta|PM|PI|Minscore|MaxPeriod) [options]

like, that's the shape of it. one file, a bunch of numbers, and then whatever options you're feeling πŸ’–


πŸ‘‘ OPTIONS πŸ‘‘

okay bestie, listen up, because this part is long and you're gonna want to actually read it. here's the where-the-arguments-go situation:

heads up: all weights, penalties, and scores are positive. no negatives, we don't do that here πŸ™ˆ
  • File β€” the sequences input file. this is your data, your baby, your whole reason for being here πŸ₯Ί
  • Match β€” the matching weight
  • Mismatch β€” the mismatching penalty
  • Delta β€” the indel penalty
  • PM β€” match probability (whole number, pls, no decimals)
  • PI β€” indel probability (whole number, same deal)
  • Minscore β€” the minimum alignment score it needs to actually bother reporting to you
  • MaxPeriod β€” the maximum period size it's willing to report

πŸŽ€ [options] = one or more of the following: πŸŽ€

  • -m β€” masked sequence file
  • -f β€” flanking sequence
  • -d β€” data file
  • -h β€” suppress html output. no html for you bestie, keeping it plain πŸ’€
  • -r β€” no redundancy elimination
  • -l <n> β€” maximum TR length expected (in millions) β€” like, -l 3 or -l=3 both mean 3 million. the human genome HG38 would need -l 6. yes, six million. it's a big girl genome, she contains multitudes 🌈
  • -ngs β€” more compact .dat output on multisequence files, and it returns 0 on success 🫢. the output goes to the screen, not a file β€” so don't go looking for a file, it's not there. you may pipe input in with this option by using - as the file name. also, short 50 flanks get appended to the .dat output. tiny flanks. little guys. πŸ’–

πŸ”₯ SEE ALSO πŸ”₯

more info lives on the TRF Unix Help web page, go read it when you're done here:

https://tandem.bu.edu/trf/trf.unix.help.html

⭐ NOTE ⭐

ok this one actually matters: the sequence file should be in FASTA format. no cap. if it's not FASTA, it's not gonna go well 😭

and FASTA looks like this, exactly like this, don't improvise:

>Name of sequence
aggaaacctgccatggcctcctggtgagctgtcctcatccactgctcgctgcctctccag
atactctgacccatggatcccctgggtgcagccaagccacaatggccatggcgccgctgt
actcccacccgccccaccctcctgatcctgctatggacatggcctttccacatccctgtg

see that > up top? that's the name of the sequence. the rest is the sequence. chef's kiss 🀌


πŸ’œ AUTHOR πŸ’œ

Copyright Β© 1999-2020 Gary Benson β€” absolute legend, honestly πŸ‘‘

this manpage was written by Andreas Tille for the Debian distribution and can be used for any other usage of the program. so generous. so kind. we love that for us πŸ’–


okay that's the whole thing!! you're basically a tandem repeat expert now, no big deal πŸ’…βœ¨

trf 4.09.1 Β· June 2020 Β· TRF(1) β€” that's the version, that's the date, that's the manual section, mwah πŸŽ€ bye bestie, go find some repeats πŸ¦‹πŸŒΈ