Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Argument parsing, usage statement, and small improvements for countKmers #11

Open
wants to merge 4 commits into
base: master
Choose a base branch
from

Conversation

standage
Copy link

@standage standage commented Feb 15, 2019

I made a few changes to the countKmers script for my own use, and thought they might be useful for others as well.

  • added option parsing so that users can configure cores, paths, etc from the command line without modifying the script
  • added a usage statement so users can see options by running ./countKmers -h
  • added an option for users to specify the FASTQ file extension; also, it autodetects whether the extension includes .gz or not and dynamically determines whether to use zcat or cat, again so users don't have to modify the script

I also updated README.md to reflect these changes, and the fact that Jellyfish 2 is a suitable alternative to the modified version of Jellyfish bundled with HAWK which seems to be giving some users issues.

@@ -1,72 +0,0 @@
HAWK version 0.9.8
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This file is redundant with README.md, and some parts are out of date. Better to just drop it.

README.md Show resolved Hide resolved
CATCMD="bunzip2 -c"
else
CATCMD=cat
fi
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Here is where the compression is autodetected.

jellyfishDir=/home/atif/jellyfish-Hawk/bin #directory where jellyfish is installed
sortDir=/home/atif/coreutils/deps/bin #directory where parallel sort is installed
#!/usr/bin/env bash
set -eo pipefail
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It's helpful to halt execution if something goes wrong. Makes it easier to isolate the point of failure if it doesn't just try to keep going to the next sample.

@standage
Copy link
Author

Rolled back the text about Jellyfish 2.

@atifrahman
Copy link
Owner

Thank you so much for your modifications! Incorporated a few. We are working on a few things. We will take care of the rest once that's done.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants