Richard Corbett

Hi @rizkg , Have you had any luck reproducing my error? I am getting some pressures at my center to have this up and running, so please let me know...

Also, do you think it may help if I try running your binary directly (or in a container?)

Thanks. Trying it out now.

Thanks for the speedy reply! Regarding split reads I have the following flagstats suggesting that the guppy minimap alignment includes both supplementary and secondary alignments. Are these the type of...

I'm looking for some differences in the bams and was surprised to find that the guppy minimap bam doesn't contain either the read sequence or the read quality string. The...

Example read pulled from a standalone minimap2 alignment: ``` f9f1c773-379a-400b-8ba6-680e75efde98 272 1 340778 0 3725S88M[truncated for brevity] * 0 0 TTGTTACGTATTGC[truncated for brevity] %$#%&&&&&&*212**.7889; [truncated for brevity]. RG:Z:272420_fastq_runid_1a9f7de9b4571f6bc9d6ff2d9713f9815f79e33NM:i:763 ms:i:26382 AS:i:26385...

Thats too bad. I was away last week and read your first message this morning so I set up a test myself that is waiting on the cluster and then...

Perfect. Many thanks.

One more question - How would you recommend counting the variants uniquely called in my set or in the GIAB set?

Thanks. Any word on permission to share your code?