Biostar Beta. Not for public use.
searching tool for calculating switch/phase error rate
2
Entering edit mode
11 months ago
Belgium

I'm looking for a tool which, given a truth vcf file and a test vcf file, calculates the phase/switch error rate. I performed phasing of a vcf using WhatsHap and want to compare the outcome to some ground truth phased vcf I have. I can't find a tool, and I don't want to write it myself. It seems like a rather common thing... so I'd expect it to exist, but my google-fu is weak today.

(Cross posted at https://bioinformatics.stackexchange.com/questions/7242/searching-tool-to-calculate-phase-switch-error-rate )

ADD COMMENTlink
0
Entering edit mode
ADD REPLYlink
0
Entering edit mode

Hi genomax, I took a look at the manual and it wasn't obvious to me which metric is homologous to switch error rate.

ADD REPLYlink
1
Entering edit mode
12 months ago
soulsacross • 30
Brazil

Have you tried vcftools --diff functions?

https://vcftools.github.io/man_0112a.html "COMPARISON OPTIONS These options are used to compare the original variant file to another variant file and output the results. All diff functions cannot be written to standard out.

--diff-switch-error

Used in conjuction with the --diff option to calculate phasing errors (specifically "switch errors"). This option generates two output files describing switch errors found between sites, and the average switch error per individual. These two files have the suffixes ".diff.switch" and ".diff.indv.switch" respectively."

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1