KMERRESISTANCE(1)

πŸ’– KMERRESISTANCE πŸ’–

okay so bestie, welcome to the kmerresistance manual β€” it's a User Commands girlie through and through, and she's here to help you figure out which genes are hanging out in your WGS samples. no cap.


✨ NAME ✨

kmerresistance β€” correlates mapped genes with the predicted species of WGS samples. that's the whole vibe. πŸ¦‹


🌸 SYNOPSIS 🌸

ngl the original page leaves this section absolutely empty β€” not a single line, not even a crumb. so like, there's nothing to show you here, and I'm not gonna make stuff up. moving on bestie πŸ’…


πŸ’œ DESCRIPTION πŸ’œ

KmerResistance correlates mapped genes with the predicted species of WGS samples β€” and the way that this this opens up so many doors?? It lets you identify genes in samples that have been poorly sequenced, and it gives you high accuracy predictions for samples with contamination. literally iconic behaviour.

KmerResistance has exactly one dependency, and it's KMA, which does the mapping for her. KMA is also freely available, so no gatekeeping here 🫢


πŸŽ€ OPTIONS πŸŽ€

  • -i β€” Input/query file name. REQUIRED, no skipping, periodt.
  • -o β€” Output file. REQUIRED too bestie.
  • -t_db β€” Template DB. REQUIRED.
  • -s_db β€” Species DB. REQUIRED.
  • -id β€” ID threshhold. default 70.0
  • -dct β€” Depth correction threshhold. default 0.1
  • -kma β€” alternative KMA. default kma

🌈 KMA options 🌈

these are the KMA settings, all up in here too:

  • -i β€” Input/query file name. default None, REQUIRED
  • -o β€” Output file. default None, REQUIRED
  • -t_db β€” Template DB. default None, REQUIRED
  • -k β€” Kmersize. default 16
  • -e β€” evalue. default 0.05
  • -delta β€” Align in pieces of delta. default 511
  • -mem_mode β€” Use kmers to choose best template, and save memory. default False
  • -ex_mode β€” Searh kmers exhaustively. default False
  • -deCon β€” Remove contamination. default False
  • -dense β€” Do not allow insertions in assembly. default False
  • -ref_fsa β€” Consensus sequnce will have "n" instead of gaps. default False
  • -matrix β€” Print assembly matrix. default False
  • -mp β€” Minimum phred score. default 30
  • -5p β€” Cut a constant number of nucleotides from the 5 prime. default 0
  • -Sparse β€” Run KmerFinder. default False
  • -ID β€” Minimum ID. default 1.0%
  • -ss β€” Sparse sorting (q,c,d). default q
  • -shm β€” Use shared DB made by kma_shm. default 0 (lvl)
  • -swap β€” Swap DB to disk. default 0 (lvl)
  • -1t1 β€” Skip HMM. default False
  • -boot β€” Bootstrap sequence. default False
  • -mrs β€” Minimum alignment score, normalized to alignment length. default 0.5
  • -reward β€” Score for match. default 1
  • -penalty β€” Penalty for mismatch. default -2
  • -gapopen β€” Penalty for gap opening. default -3
  • -gapextend β€” Penalty for gap extension. default -1
  • -h β€” Shows this help message. she's helpful like that πŸ’–

🧚 EXAMPLES 🧚

Mapping reads against resistance genes β€” this is the moment, bestie:

kmerresistance -i sample_1.fastq sample_2.fastq -o out -t_db ResFinder -s_db bacteria

look at her go 😭✨


πŸ‘‘ AUTHOR πŸ‘‘

This manpage was written by Andreas Tille for the Debian distribution, and can be used for any other usage of the program. we love a generous king πŸ’œ


okay so this is kmerresistance version 2.0+git20180205...., dated June 2018, and yeah β€” that's the whole KMERRESISTANCE(1) page. slay, go map some genes, ily bye πŸŽ€πŸ’…πŸ¦‹