Skip to content

Added bowtie2 and SNAP aligners (resolves #102)#105

Open
fnothaft wants to merge 1 commit intoBD2KGenomics:masterfrom
fnothaft:issues/102-snap-bowtie2
Open

Added bowtie2 and SNAP aligners (resolves #102)#105
fnothaft wants to merge 1 commit intoBD2KGenomics:masterfrom
fnothaft:issues/102-snap-bowtie2

Conversation

@fnothaft
Copy link
Contributor

No description provided.

def run_bowtie2(job,
read1,
name1, name2, name3, name4, rev1, rev2, ref,
read2=None):
Copy link
Collaborator

Choose a reason for hiding this comment

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

Is this being called directly from the tool that generates the index files? It might be a bit cleaner to pass a dict of IDs, but if you'd prefer it this way that's fine.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

I am OK either way. Let me know if you have strong preferences.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Actually, I've been thinking about this more (i.e., actually wiring up the tools) and you're correct. I'll be refactoring this shortly.

'-o', '-bam', '/data/sample.bam',
'-t', str(job.cores)]
else:

Copy link
Collaborator

Choose a reason for hiding this comment

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

You think whitespace is free or something?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

I do.

if not sort:
return job.fileStore.writeGlobalFile(os.path.join(work_dir, 'sample.bam'))
else:
return (job.fileStore.writeGlobalFile(os.path.join(work_dir, 'sample.bam')),
Copy link
Collaborator

Choose a reason for hiding this comment

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

Minor: outermost parentheses unnecessary as more than one item is always returned as a tuple.

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