-
Notifications
You must be signed in to change notification settings - Fork 10
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
Readcount fix #273
base: master
Are you sure you want to change the base?
Readcount fix #273
Conversation
if cmd_args['config']: | ||
cmd += ' -c {config}' | ||
p = run_cmd( cmd.format(**cmd_args), stdout=lfile, stderr=subprocess.STDOUT ) | ||
tcmd = "samtools view -F 2048 -b {tmp_bam} > {bamfile}" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
You might want to remove the tmp_bam file to be considerate of space or make that an option
lofreq = bc_cfg['lofreq']['default'] | ||
lofreq_options = bc_cfg['lofreq_options']['default'] or '' | ||
if lofreq: | ||
cmd1 = 'lofreq call -f {reference} {bamfile} -o {vcf} ' + lofreq_options |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Did you install lofreq anywhere? I don't see that you added it to conda requirements or other requirements file?
ngs_mapper/runsample.py
Outdated
if cmd_args['config']: | ||
cmd += ' -c {config}' | ||
p = run_cmd( cmd.format(**cmd_args), stdout=lfile, stderr=subprocess.STDOUT ) | ||
r = p.wait() | ||
if r != 0: | ||
logger.critical( "{0} did not exit sucessfully".format(cmd.format(**cmd_args)) ) | ||
rets.append( r ) | ||
# Variant Calling | ||
from ngs_mapper.config import load_config | ||
bc_cfg = load_config(cmd_args['config'])['base_caller'] |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Only way to use this is to edit config file and use that config file right?
Maybe we can add a flag to run_sample.py as well?
7d1e3ff
to
8d6b3b6
Compare
8d6b3b6
to
2146a5f
Compare
No description provided.