VELVETG_127(1)

πŸ‘‘ VELVETG_127 β€” the manual, but make it cute

okay bestie, listen. this is velvetg_127 from the User Commands section of the manual (section 1, iykyk), and it's the one that does de Bruijn graph construction, error removal and repeat resolution, with MAXKMERLENGTH = 127. that's the whole vibe. 127 is a lot of kmer and she is NOT shy about it πŸ’…


🌸 SYNOPSIS 🌸

./velvetg directory [options]

✨ DESCRIPTION ✨

so basically you run it like this: ./velvetg directory [options], and here's what the pieces mean:

  • directory β€” the working directory name. that's it. that's the whole required argument. she needs a home 🏠

πŸŽ€ Standard options (the main-character ones)

  • -cov_cutoff <floating-point|auto> β€” removes low coverage nodes AFTER tour bus, or you can let the system infer it for you. (default: no removal) 🧹
  • -ins_length <integer> β€” the expected distance between two paired end reads. (default: no read pairing)
  • -read_trkg <yes|no> β€” tracking of short read positions in the assembly. (default: no tracking)
  • -min_contig_lgth <integer> β€” minimum contig length that gets exported to the contigs.fa file. (default: hash length * 2)
  • -amos_file <yes|no> β€” export the assembly to an AMOS file. (default: no export)
  • -exp_cov <floating point|auto> β€” expected coverage of unique regions, or let the system infer it. (default: no long or paired-end read resolution)
  • -long_cov_cutoff <floating-point> β€” removal of nodes with low long-read coverage AFTER tour bus. (default: no removal)

πŸ¦‹ Advanced options (for when you're feeling brave)

  • -ins_length* <integer> β€” expected distance between two paired-end reads in the respective short-read dataset. (default: no read pairing)
  • -ins_length_long <integer> β€” expected distance between two long paired-end reads. (default: no read pairing)
  • -ins_length*_sd <integer> β€” est. standard deviation of the respective dataset. (default: 10% of the corresponding length)
[replace '*' by nothing, '2' or '_long' as necessary] β€” fr fr, the * is a placeholder, swap it out properly πŸ™ˆ
  • -scaffolding <yes|no> β€” scaffolding of contigs using paired end information. (default: on)
  • -max_branch_length <integer> β€” maximum length in base pair of a bubble. (default: 100) 🫧
  • -max_divergence <floating-point> β€” maximum divergence rate between two branches in a bubble. (default: 0.2)
  • -max_gap_count <integer> β€” maximum number of gaps allowed in the alignment of the two branches of a bubble. (default: 3)
  • -min_pair_count <integer> β€” minimum number of paired end connections to justify the scaffolding of two long contigs. (default: 5)
  • -max_coverage <floating point> β€” removal of high coverage nodes AFTER tour bus. (default: no removal)
  • -coverage_mask <int> β€” minimum coverage required for confident regions of contigs. (default: 1)
  • -long_mult_cutoff <int> β€” minimum number of long reads required to merge contigs. (default: 2)
  • -unused_reads <yes|no> β€” export unused reads into the UnusedReads.fa file. (default: no)
  • -alignments <yes|no> β€” export a summary of contig alignment to the reference sequences. (default: no)
  • -exportFiltered <yes|no> β€” export the long nodes which were eliminated by the coverage filters. (default: no)
  • -clean <yes|no> β€” remove all the intermediary files which are useless for recalculation. (default : no)
  • -very_clean <yes|no> β€” remove all the intermediary files. (default: no)
⚠️ heads up bestie: with -very_clean, no recalculation is possible. she said goodbye and she meant it πŸ’€
  • -paired_exp_fraction <double> β€” remove all the paired end connections which [have] less than the specified fraction of the expected count. (default: 0.1)
  • -shortMatePaired* <yes|no> β€” for mate-pair libraries, indicate that the library might be contaminated with paired-end reads. (default no)
  • -conserveLong <yes|no> β€” preserve sequences with long reads in them. (default no)

🌈 Output (the lil gifts she leaves behind)

  • directory/contigs.fa β€” fasta file of contigs longer than twice hash length 🧬
  • directory/stats.txt β€” stats file (tab-spaced) that's useful for determining the appropriate coverage cutoff πŸ“Š
  • directory/LastGraph β€” a special formatted file with all the information on the final graph
  • directory/velvet_asm.afg β€” (if requested) the AMOS compatible assembly file

πŸ’– the sign-off πŸ’–

and that's the whole page, bestie. this version is velvetg_127 1.2.10+dfsg1, dated July 2018 β€” a whole icon, aging like fine wine 🍷 no cap, go build your graphs and slay. periodt. 🫢✨