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

1.3.2 read group - can use only 1? #311

Open
RichardCorbett opened this issue Jun 9, 2021 · 1 comment
Open

1.3.2 read group - can use only 1? #311

RichardCorbett opened this issue Jun 9, 2021 · 1 comment

Comments

@RichardCorbett
Copy link

RichardCorbett commented Jun 9, 2021

Is your feature request related to a problem? Please describe.
I have samples that are sequenced deeply and the data span multiple nanopore flowcells. When I merge those bams there will be multiple read group (I per flowcell). When I run medaka I get an error saying RuntimeError: The bam F24722_merged.bam contains more than one read group. Please specify --RG to select which read groupto process from {'272423_fastq_runid_9c2b42a90ae458fb58e19950f06b13da5ceede77', '272422_fastq_runid_5b33097e8a403b7ed7d030764f669d6be91d0ef4'}

How can I run medaka on all data in the bam without editing the read group information?

Describe the solution you'd like

Process multiple read groups if they have the same library and/or sample ID. Perhaps this is possible now, but I haven't found a way to do it.

@cjw85
Copy link
Member

cjw85 commented Jun 9, 2021

Apologies, this is not currently possibly. You will have to write a bam without read group information or with a single read group.

I cannot unfortunately give a timeline on implementing an --all_reads option.

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

No branches or pull requests

2 participants