-
Notifications
You must be signed in to change notification settings - Fork 18
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
I have modified the Bowtie mapper class because the index for bowtie2… #361
Conversation
… was not being set and was picking up the bowtie index dir
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.
Thanks @Acribbs
Please see my comments before merging
CGATPipelines/PipelineMapping.py
Outdated
@@ -3044,12 +3044,21 @@ def mapper(self, infiles, outfile): | |||
nfiles -= 2 | |||
else: | |||
raise ValueError("unexpected number of files") | |||
index_prefix = "%(bowtie_index_dir)s/%(genome)s_cs" | |||
if exexutable == "bowtie": |
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.
should this be executable
instead?
CGATPipelines/PipelineMapping.py
Outdated
elif self.datatype == "fasta": | ||
data_options.append("-f") | ||
index_prefix = "%(bowtie_index_dir)s/%(genome)s" | ||
if exexutable == "bowtie": |
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.
should this be executable
instead?
CGATPipelines/PipelineMapping.py
Outdated
else: | ||
index_prefix = "%(bowtie_index_dir)s/%(genome)s" | ||
if exexutable == "bowtie": |
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.
should this be executable
instead?
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.
Additional idea.
CGATPipelines/PipelineMapping.py
Outdated
index_prefix = "%(bowtie_index_dir)s/%(genome)s_cs" | ||
if exexutable == "bowtie": | ||
index_prefix = "%(bowtie_index_dir)s/%(genome)s_cs" | ||
else: |
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.
I would even do elif executable = "bowtie2"
and raise an exception in the else
, as a sanity check
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.
Yes I agree with you. I will make the changes.
@sebastian-luna-valero I have made the changes you suggest. Thanks |
bug fix for ValueError
@sebastian-luna-valero Done, many thanks |
Thanks @Acribbs |
… was not being set and was picking up the bowtie index dir