Closed Griffan closed 3 months ago
Thanks @Griffan. Trying this out. will report back.
Hi @yfarjoun, I have pushed the fix by bypassing these samples with uninformative GT or PL fields. Let me know if it finish the job successfully on your side. Thanks!
it worked (in that it didn't explode..) I am missing many sites, due to the fact that they are multi-allelic. how did you pre process the 1KG vcfs, as they have lots of multiallelic sites as well?
I would recommend any multi-allelic sites as verifyBamID is using biallelic variants. If the reference data is too big, you may decompose the VCF into bi-allelics first and filter variants by missingness.
This PR aims to resolve segfault reported by multiple users in issue https://github.com/Griffan/VerifyBamID/issues/45 To better trace the error without exchanging testing dataset, I added this simple debugging mode to print backtrace at the crash site.
Debugging Mode
If you encounter abnormal errors, e.g. "Segmentation fault", you can try to recompile the build under debugging mode:
and then rerun your command line to locate the backtrace message, and then post it to issues page, for example: