WGLab / NanoCaller

Variant calling tool for long-read sequencing data
MIT License
90 stars 8 forks source link

AD (allelic depth) field #38

Closed jamesdalg closed 11 months ago

jamesdalg commented 11 months ago

To the community: I really think this is a fantastic tool, with incredible speed and great options for phasing. Is there a way to get the AD field to populate with NanoCaller, much like Deepvariant and GATK's HaplotypeCaller? There are multiple downstream tools that require this information and I'm wondering if I can somehow coax NanoCaller into giving this information, it would be a much faster alternative than deepvariant (which has some stages of the code that don't parallelize well and take days).

umahsn commented 11 months ago

Hi, thank you for bringing this to attention. I will add this functionality into NanoCaller.

umahsn commented 11 months ago

I have added information for allelic depth, for overall as well as forward/reverse strands, in the format field for SNPs. I will extend the functionality to indels soon.

jamesdalg commented 11 months ago

If you can, sometime, please update the docker image sometime if it's not too much trouble. I've never been able to get it to work on the cluster without using the image.

umahsn commented 11 months ago

I have pushed a new image version 3.4.0 that has allelic depths.

jamesdalg commented 11 months ago

Thanks so very much!! This is a huge help.