BP_SEARCH2TRIBE(1)

πŸ“¦ version 1.7.8, in nixos 26.05, nixpkgs unstable

πŸ¦‹ bp_search2tribe β€” turning SearchIO reports into a TRIBE matrix, bestie πŸ’…

okay so this is the manual for bp_search2tribe, a Perl script for the bioperl crew ✨ (iykyk) β€” and the whole vibe is "turn your search output into something MCL can actually chew on".


✨ NAME ✨

bp_search2tribe β€” Turn SearchIO parseable reports(s) into TRIBE matrix 🧚

no cap, that's literally the whole job.


🌸 SYNOPSIS 🌸

Usage:

bp_search2tribe [-o outputfile] [-f reportformat] [-w/--weight] file1
file2 ..

yes bestie, those files go on separate-ish lines in the original and we're keeping it exactly like that 🫢


πŸ’– DESCRIPTION πŸ’–

okay so real talk first, and I'm gonna be honest with you because that's what besties do: this script is probably too slow for most people's uses. 😭 like, not me being dramatic β€” that's straight from the source.

it's better to use something like scripts/searchio/fastam9_to_table, or -m 9 output from BLAST, or the blast2table from the BLAST O'Reilly book, to get a tabular output from these programs β€” and then feed that table into MCL with the mcxdeblast script and the --m9 option. that's the glow-up route πŸ”₯

but if you're here for what this script does: it will turn a protein Search report (BLASTP, FASTP, SSEARCH) into a Markov Matrix for TribeMCL clustering πŸ‘‘

the options are:

  • -o filename β€” the output filename [default STDOUT] πŸŽ€
  • -f format β€” search result format (blast, fasta) (ssearch is fasta format). default is blast.
  • -w or --weight VALUE β€” Change the default weight for E(0.0) hits to VALUE (default=200 (i.e. 1e-200) ) πŸ’…
  • -h β€” this help menu

additionally, specify the filenames you want to process on the command-line. if no files are specified then STDIN input is assumed. you specify this by doing:

bp_search2tribe < file1 file2 file3

main character energy fr fr 🌈


πŸ‘‘ AUTHOR πŸ‘‘

Jason Stajich, jason-at-bioperl-dot-org πŸ¦‹


πŸ’œ the cute sign-off πŸ’œ

this page is living its best life as part of perl v5.42.3, dated 2021-02-03 β€” and yes, that's the footer, we're mentioning it here instead of copying it because ✨ formatting ✨. mwah πŸ’‹ periodt.

πŸ“¦ where this page comes from