Commit 7be3d1fe authored by Sander Bollen's avatar Sander Bollen

fix wrong size estimation for gq qualities.

parent 9de1db12
Pipeline #2124 failed with stage
...@@ -50,7 +50,7 @@ def gq_for_gvcf_record(record: cyvcf2.Variant, maxlen: int = 15000) -> List[int] ...@@ -50,7 +50,7 @@ def gq_for_gvcf_record(record: cyvcf2.Variant, maxlen: int = 15000) -> List[int]
""" """
start = record.start start = record.start
end = record.end end = record.end
size = end - end size = end - start
gq = record.gt_quals[0] gq = record.gt_quals[0]
if size < maxlen: if size < maxlen:
return [gq]*(end-start) return [gq]*(end-start)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment