Illumina / REViewer

A tool for visualizing alignments of reads in regions containing tandem repeats
GNU General Public License v3.0
78 stars 13 forks source link

metrics.tsv AlleleDepth of head and v0.2.7 deviate significant #59

Closed maarten-k closed 1 year ago

maarten-k commented 1 year ago

I compiled version 0.27 and the head of the default branch. The output of the AlleleDepth differs: the total AlleleDepth is the same, but the first Allele always seems lower. A minor deviation was expected, but the variation is quite significant.

For example(variant names are renamed because of privacy concerns):

==> sampleX.metrics.tsv_0.2.7 <==

VariantId       Genotype        AlleleDepth
var1       4/4     16.44/15.44
var2       6/6     24.11/18.72
var3       4/4     28.38/25.75
var4 9/9     22.39/21.39
var5       4/4     15.75/16.95
var6       2/2     15.45/19.91

==> sampleX.metrics.tsv_head <==

VariantId       Genotype        AlleleDepth
var1       4/4     12.56/19.31
var2       6/6     19.39/23.44
var3       4/4     24.94/29.19
var4 9/9     18.22/25.56
var5       4/4     23.60/9.10
var6       2/2     19.82/15.55

The input files are identical. Both binaries were compiled with the same container(available upon request) and the binary of version V0.2.7 gives the same results as the static compiled version downloaded from Git Hub.

Which version should be considered the "gold standard"?

maarten-k commented 1 year ago

I could not rereplicate this problem, so I think there was another problem.