how to select a better qulity variant from Spanning alleles

[zhoudreames]

I only want to save one variant possessed the best calling qulity from spanning alleles，how to batch process it？ thanks~

[dancooke]

There's no way to do this with Octopus, but I'm sure you can write a script using pysam that does what you want.

[zhoudreames]

thanks ~