Skip to content

Commit

Permalink
spades: removed references to min contig length, added always succeed…
Browse files Browse the repository at this point in the history
… flag
  • Loading branch information
notestaff committed May 21, 2017
1 parent 00a5186 commit 1c3335f
Show file tree
Hide file tree
Showing 2 changed files with 1 addition and 9 deletions.
7 changes: 0 additions & 7 deletions assembly.py
Expand Up @@ -344,7 +344,6 @@ def parser_assemble_trinity(parser=argparse.ArgumentParser()):
def assemble_spades(
inBam,
outFasta,
min_contig_len=300,
always_succeed=False,
threads=1,
):
Expand Down Expand Up @@ -375,12 +374,6 @@ def assemble_spades(
def parser_assemble_spades(parser=argparse.ArgumentParser()):
parser.add_argument('inBam', help='Input unaligned reads, BAM format.')
parser.add_argument('outFasta', help='Output assembly.')
parser.add_argument(
'--min_contig_len',
default=300,
type=int,
help='Discard contigs shorter than this many basepairs (default %(default)s)'
)
parser.add_argument(
"--always_succeed",
help="""If SPAdes fails (usually because insufficient reads to assemble),
Expand Down
3 changes: 1 addition & 2 deletions pipes/rules/assembly.rules
Expand Up @@ -66,14 +66,13 @@ rule assemble_spades:
cores=int(config.get("number_of_threads", 1))
params: LSF=config.get('LSF_queues', {}).get('short', '-W 4:00'),
UGER=config.get('UGER_queues', {}).get('short', '-q short'),
min_contig_len=str(config["trinity_min_contig_len"]),
logid="{sample}",
numThreads=str(config.get("number_of_threads", 1))
run:
makedirs(expand("{dir}/{subdir}",
dir=[config["data_dir"],config["tmp_dir"]],
subdir=config["subdirs"]["assembly"]))
shell("{config[bin_dir]}/assembly.py assemble_spades {input} {output}")
shell("{config[bin_dir]}/assembly.py assemble_spades {input} {output} --always_succeed")

rule combine_de_novo_contigs:
'''Combine de novo contigs obtained by different assemblers'''
Expand Down

0 comments on commit 1c3335f

Please sign in to comment.