ABACAS(1)

πŸŽ€ ABACAS(1) β€” User Commands πŸŽ€

✨ NAME ✨

abacas πŸ’– okay so this is the one: Algorithm Based Automatic Contiguation of Assembled Sequences 🧚

πŸ’… SYNOPSIS πŸ’…

Okay bestie, there are two ways to call this thing 🫢

abacas -r ref -q qs -p prog  [OPTIONS]

OR

abacas -r ref -q psf -e

Here's what the bits mean 🌸

  • ref β€” a reference sequence, in a single fasta file πŸ“„
  • qs β€” your contigs, in multi-fasta format, plus the MUMmer program you want to use: nucmer or promer πŸ¦‹
  • psf β€” a pseudomolecule/ordered sequence file in fasta format πŸ‘‘

OPTIONS πŸ’œ

  • -h β€” print usage πŸ™ˆ
  • -d β€” use default nucmer/promer parameters ✨
  • -s int β€” minimum length of exact matching word (nucmer default = 12, promer default = 4) πŸ”₯
  • -m β€” print ordered contigs to file in multifasta format πŸ“
  • -b β€” print contigs in bin to file 🌈
  • -N β€” print a pseudomolecule without "N"s πŸ’€
  • -i int β€” minimum percent identity [default 40] πŸ’–
  • -v int β€” minimum contig coverage [default 40] πŸŽ€
  • -V int β€” minimum contig coverage difference [default 1] 🌸
  • -l int β€” minimum contig length [default 1] ⭐
  • -t β€” run tblastx on contigs that are not mapped 🧚
  • -g string (file name) β€” print uncovered regions (gaps) on reference to file name πŸ¦‹
  • -a β€” append contigs in bin to the pseudomolecule πŸ’…
  • -o prefix β€” output files will have this prefix 🫢
  • -P β€” pick primer sets to close gaps πŸ‘‘
  • -f int β€” number of flanking bases on either side of a gap for primer design (default 350) 😭
  • -R int β€” run mummer [default 1, use -R 0 to avoid running mummer] 🌈
  • -e β€” escape contig ordering, i.e. go to primer design πŸ”₯
  • -c β€” reference sequence is circular πŸ“

🌸 DESCRIPTION 🌸

Okay so ABACAS is literally built to rapidly contiguate (align, order, orientate), visualize and design primers to close gaps on shotgun assembled contigs, all based on a reference sequence. It's giving precision πŸ’…βœ¨

ABACAS uses MUMmer to find alignment positions and identify syntenies of assembled contigs against the reference. Then the output gets processed into a pseudomolecule, taking overlapping contigs and gaps into account πŸ¦‹ The way it generates a comparison file you can use to visualize ordered and oriented contigs in ACT is so iconic 🌈 Synteny shows up as red bars, and the colour intensity drops as the percent identity between comparable blocks gets lower πŸ”₯ Info on contigs (orientation, percent identity, coverage, and overlap with other contigs) can also be visualized by loading the outputted feature file into ACT πŸ’œ


πŸŽ€ AUTHOR πŸŽ€

ABACAS IS Copyright (C) 2008-10 The Wellcome Trust Sanger Institute, Cambridge, UK. πŸ‡¬πŸ‡§πŸ’–

This manual page was written by Andreas Tille <tille@debian.org>, for the Debian project (and may be used by others). Huge thanks to Andreas for the page 🫢🧚


Version 1.3.1, dated 2011-02-11. Bye bestie, go slay your assembly πŸ’…πŸ¦‹πŸŒΈ