MAPSEMBLER2_KISSREADS_GRAPH(1)

πŸ’… MAPSEMBLER2_KISSREADS_GRAPH β€” the manual page, but make it make sense πŸ’…

okay bestie, so this is the manpage for mapsembler2_kissreads_graph from section 1 of the User Commands manual. it's giving graph-mapping realness and honestly? it ate. let's get into it 🫢


✨ NAME ✨

mapsembler2_kissreads_graph β€” maps the provided reads on the graph πŸ¦‹


🌸 SYNOPSIS 🌸

mapsembler2_kissreads_graph input_graph <readsC1.fasta/fastq[.gz]>
[<readsC2.fasta/fastq[.gz]> [<readsC3.fasta/fastq[.gz] ...] -M -t type
[-k value] [-c min_coverage] [-d max_substitutions] [-o name] [-h]

so like. you hand it a graph, you hand it reads, and it does the thing. the reads can be .fasta or .fastq and they can be gzipped (.gz) because we're not savages πŸ’– and yeah you can stack up as many read files as you want, <readsC2> , <readsC3> , and so on forever, iykyk.


πŸ‘‘ DESCRIPTION πŸ‘‘

mapsembler2_kissreads_graph maps the provided reads on the graph. no cap, that's the whole job. πŸ’…

with the option -t coverage it outputs an equivalent graph that removes uncovered edges and then adds:

  • for each node: the coverage per sample and per position ✨
  • for each edge: the number of mapped reads per sample using this edge 🌈

with the option -t modify it outputs the same simplified graph, going through:

  • removing low covered edges and nodes (less that min_coverage) 🧚
  • then phasing simple non branching paths ⭐

and if you pass -M, the input is considered as a Mapsembler output β€” which means it's composed of multiple independent graphes πŸ¦‹ (yes, "graphes", that's what it says, we respect the original spelling here, periodt).


πŸŽ€ OPTIONS πŸŽ€

πŸ’– MANDATORY πŸ’–

like you literally cannot skip these, bestie, they're mandatory for a reason 😭

  • -t <type> β€” the type. it's giving "pick one":
    • "c" "coverage" β€” the coverage vibe 🌸
    • "m" "modify" β€” the modify vibe πŸ”₯

πŸ¦‹ OPTIONAL πŸ¦‹

these ones are the fun little extras, the accessories to the outfit πŸ’…

  • -M β€” the input is considered as a Mapsembler output πŸ‘‘
  • -o file_name β€” write obtained graph. default: standard output (so if you say nothing, it just spills out at you, main character energy)
  • -k size_seed β€” will use seeds of length size_seed. default: 25. 🌱
  • -c min_coverage β€” will consider an edge as coherent if coverage (number of reads from all sets using this edge) is at least min_coverage. default: 2 ✨
  • -d max_substitutions β€” will map a read on the graph with at most max_substitutions substitutions. default: 1 πŸ™ˆ
also -h is hanging out in the synopsis up there, being mysterious. the page doesn't describe it further and neither will i β€” i'm not out here inventing lore, that's not the vibe πŸ’€

🌈 AUTHOR 🌈

this manpage was written by Andreas Tille for the Debian distribution, and can be used for any other usage of the program. thank you Andreas, you absolute icon 🫢


and that's the whole thing!! mapsembler2_kissreads_graph 2.2.4, October 2018 β€” mwah, go map some reads, bestie, you got this πŸ’–βœ¨πŸŽ€πŸ¦‹πŸŒΈπŸ’