cnvkit icon indicating copy to clipboard operation
cnvkit copied to clipboard

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

Open IvantheDugtrio opened this issue 1 year ago • 1 comments

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)

IvantheDugtrio avatar Feb 15 '24 21:02 IvantheDugtrio

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 avatar Apr 12 '24 19:04 etal