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

params.circularfilter ? '' : '-f true -x false' #638

Closed
zhongxing123 opened this issue Dec 19, 2020 · 5 comments · Fixed by #639
Closed

params.circularfilter ? '' : '-f true -x false' #638

zhongxing123 opened this issue Dec 19, 2020 · 5 comments · Fixed by #639
Labels
question Further information is requested user-support For setup and usage questions and advice rather than code issues

Comments

@zhongxing123
Copy link

zhongxing123 commented Dec 19, 2020

Hello,I'm a little confused about circularmapper. I find script def filter = params.circularfilter ? '' : '-f true -x false' in main.nf. Does this mean that -f true-x false is performed when circularfilter == false? Please give me some guidance. This has been bothering me for a long time. Thank you very much

@apeltzer apeltzer added question Further information is requested user-support For setup and usage questions and advice rather than code issues labels Dec 19, 2020
@apeltzer
Copy link
Member

Hi @zhongxing123 ! Yes, this means when the parameter --circularfilter is set to false, that the pipeline will be run with the -f true -x false option(s).

@zhongxing123
Copy link
Author

zhongxing123 commented Dec 20, 2020 via email

@jfy133
Copy link
Member

jfy133 commented Dec 21, 2020

@apeltzer actually isn't this a bug?

def filter = params.circularfilter ? '' : '-f true -x false'

If --circularfilter is TRUE, you would want filtering to turn, correct? Therefore if TRUE you would want to supply -f (to remove the reads that don't match the circular correct? At least that's what I understand from realignsamfile?

 -f,--filterCircularReads <FILTER>                     filter the reads
                                                       that don't map to a
                                                       circular identifier
                                                       (true/false),
                                                       default false

@jfy133 jfy133 mentioned this issue Dec 21, 2020
11 tasks
@jfy133
Copy link
Member

jfy133 commented Dec 23, 2020

Hi @zhongxing123 -

Thank you for looking so closely at the code. Indeed your confusion may have stemmed from the fact this was indeed a bug.

in the fix #639 we have corrected this. So in the next release:

--circularfilter will turn ON filtering of circularmapper generated BAM file. This will therefore produce BAM files with reads that only mapped to the elongated reference genome, and also now with the non-elongated references to be excluded from the BAM header.

@jfy133 jfy133 closed this as completed Dec 23, 2020
@zhongxing123
Copy link
Author

zhongxing123 commented Dec 24, 2020 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
question Further information is requested user-support For setup and usage questions and advice rather than code issues
Projects
None yet
Development

Successfully merging a pull request may close this issue.

3 participants