Skip to content

Commit 8aa6a08

Browse files
authored
Merge pull request #164 from PolinaBevad/f_fix_tcov_zero
Added fix for total coverage in case of candidate deletion on 3'.
2 parents 4005734 + a1da460 commit 8aa6a08

File tree

6 files changed

+18
-2
lines changed

6 files changed

+18
-2
lines changed

src/main/java/com/astrazeneca/vardict/modules/StructuralVariantsProcessor.java

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1011,6 +1011,13 @@ && isHasAndEquals(complement(ref.get(bp + 1)), ref, p5 - 1)) {
10111011
sv.splits += cnt3;
10121012
sv.clusters += pairs != 0 ? 1 : 0;
10131013

1014+
if (!cov.containsKey(p3)) {
1015+
cov.put(p3, pairs + sc3v.cnt);
1016+
}
1017+
if (cov.containsKey(bp) && cov.get(bp) < cov.get(p3)) {
1018+
cov.put(bp, cov.get(p3));
1019+
}
1020+
10141021
adjCnt(vref, sc3v, conf);
10151022
Variation tmp = new Variation();
10161023
tmp.cnt = pairs;
Binary file not shown.
Binary file not shown.

testdata/fastas/hg38.fa.csv

Lines changed: 3 additions & 1 deletion
Large diffs are not rendered by default.
Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,4 @@
11
Simple,GRCh37.fa,NA24385_tcov_empty.bam,chr8,122660146,122660347,-f 0.001
22
NA24385_tcov_empty testbed 8 122660339 122660339 T A 37 4 21 12 0 4 T/A 0.1081 2;0 56.5 1 40.8 1 60.0 8.000 0.1143 0 0 2.000 1 2.0 4 35 GGTAGGTATTCAAAAGTATT AAGGTCAATAATTGTGATCT 8:122660147-122660347 SNV 0 0
3-
NA24385_tcov_empty testbed 8 122660248 122660248 G A 44 15 9 20 4 11 G/A 0.3409 2;2 42.1 1 39.3 1 60.0 30.000 0.3659 0 0 1.000 1 1.7 15 41 ACTCTTTATAATCATCTCTC ATAGTCCTTTCCACAGATGA 8:122660147-122660347 SNV 0 0
3+
NA24385_tcov_empty testbed 8 122660248 122660248 G A 44 15 9 20 4 11 G/A 0.3409 2;2 42.1 1 39.3 1 60.0 30.000 0.3659 0 0 1.000 1 1.7 15 41 ACTCTTTATAATCATCTCTC ATAGTCCTTTCCACAGATGA 8:122660147-122660347 SNV 0 0
4+
NA24385_tcov_empty testbed 8 122660114 122661194 T <DEL> 3 3 0 0 1 2 -1080/-1080 1.0000 0;2 50.3 1 40.4 1 28.0 6.000 1.0000 1.0000 0 0 0 0.3 3 3 CTCAGGTGATCTACCCACCT GGCCTCCCAAAGTGCTGGGA 8:122660147-122660347 DEL 0 2-1-1
Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
Simple,hg38.fa,Tcov_zero_SV_DEL_candidate3.bam,chr8,90534611,90620070,-f 0.1
2+
Tcov_zero_SV_DEL_candidate3 testbed chr8 90615919 90615919 C T 2 2 0 0 1 1 T/T 1.0000 0;2 55.5 1 37.0 0 60.0 4.000 1.0000 0 0 4.000 1 1.0 2 2 GTGATATGGACAATAAGGTC GGGGTGAGGCGGTGTCAGAT chr8:90534612-90620070 SNV 0 0
3+
Tcov_zero_SV_DEL_candidate3 testbed chr8 90617256 90617256 C T 7 7 0 0 4 3 T/T 1.0000 0;2 42.3 1 37.0 0 60.0 14.000 1.0000 0 1 1.000 1 1.1 7 7 GTTAAGACTTTGAGGGACTG TGGGAAGGCATGATTGGTTT chr8:90534612-90620070 SNV 0 0
4+
Tcov_zero_SV_DEL_candidate3 testbed chr8 90605762 90645105 T <DEL> 6 6 0 0 4 2 -39343/-39343 1.0000 0;2 56.8 1 36.7 1 60.0 12.000 1.0000 1.0000 0 0 0 0.5 6 6 TTAATATGAGAGAACTACAG CTTACCGAAAACACTGCATT chr8:90534612-90620070 DEL 0 2-4-1
5+
Tcov_zero_SV_DEL_candidate3 testbed chr8 90617487 90617487 C T 3 3 0 0 1 2 T/T 1.0000 0;2 29.3 1 24.3 1 60.0 2.000 1.0000 0 1 1.000 1 1.3 2 2 GATGGTTTTATAAAAGGGAG TCCCCTGCACACACTCTCTT chr8:90534612-90620070 SNV 0 0
6+
Tcov_zero_SV_DEL_candidate3 testbed chr8 90615448 90615448 G A 2 2 0 0 1 1 A/A 1.0000 0;2 45.5 1 37.0 0 60.0 4.000 1.0000 0 0 3.000 1 1.0 2 2 TGTCGTGGGAGGGACCTCGC GGAGGTAACTGAATCATGGG chr8:90534612-90620070 SNV 0 0

0 commit comments

Comments
 (0)