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

Allow using reference fasta with reference cnn in cnvkit batch when processing crams #869

Open
IvantheDugtrio opened this issue Feb 15, 2024 · 1 comment
Labels

Comments

@IvantheDugtrio
Copy link

I get the following error message when trying to process crams. My cram files have @sq tags that point to a reference path that does not exist on my computer, but I have the reference file. Recreating the path to the references is not practical as each cram has different reference paths. I am looking into samtools reheader but ideally I don't modify these crams in any way.

CNVkit 0.9.10 If -r/--reference is given, options to construct a new reference (-f/--fasta) should not be used. (See: cnvkit.py batch -h)

@etal
Copy link
Owner

etal commented Apr 12, 2024

This error is about CNVkit's reference.cnn file, not the FASTA reference genome sequences. If you used -r with an existing file, then you don't need to (and shouldn't) also use -f.

@etal etal added the question label Apr 12, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants