Skip to content

Commit

Permalink
Add script for sketching
Browse files Browse the repository at this point in the history
  • Loading branch information
rpetit3 committed Feb 7, 2023
1 parent 88d193e commit c7e0522
Show file tree
Hide file tree
Showing 2 changed files with 30 additions and 0 deletions.
13 changes: 13 additions & 0 deletions CITATION.cff
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
cff-version: 1.2.0
message: "If you use Bactopia, please cite it as below."
authors:
- family-names: "Petit III"
given-names: "Robert A. "
orcid: "https://orcid.org/0000-0002-1350-9426"
- family-names: "Read"
given-names: "Timothy D."
orcid: "https://orcid.org/0000-0001-8966-9680"
title: "Bactopia: a Flexible Pipeline for Complete Analysis of Bacterial Genomes. mSystems. 5 (2020)"
doi: 10.1128/mSystems.00190-20
url: "https://github.com/bactopia/bactopia"
version: 1.4.1
17 changes: 17 additions & 0 deletions bin/bactopia-sketcher
Original file line number Diff line number Diff line change
@@ -0,0 +1,17 @@
#! /bin/bash
VERSION=1.0.0
PREFIX=$1
FASTQ=$2

# If no user input, print usage
if [[ $# == 0 ]]; then
echo "bactopia-sketcher - v${VERSION}"
echo ""
echo "bactopia-sketcher <PREFIX> <FASTQ> <OPT1> ... <OPTN>"
echo ""
exit
fi

gzip -cd ${FASTQ} | mash sketch -o ${PREFIX}-k21 -k 21 -s 21 -r -I ${PREFIX} -
gzip -cd ${FASTQ} | mash sketch -o ${PREFIX}-k31 -k 31 -s 31 -r -I ${PREFIX} -
sourmash sketch dna -p k=21,k=31,k=51,abund,scaled=1000 --merge ${PREFIX} -o ${PREFIX}.sig ${FASTQ}

0 comments on commit c7e0522

Please sign in to comment.