Skip to content
This repository was archived by the owner on May 22, 2026. It is now read-only.

Commit f132443

Browse files
committed
Fixed type for variant in case when ref==alt (for pileup mode).
1 parent 03eb543 commit f132443

3 files changed

Lines changed: 96 additions & 92 deletions

File tree

src/main/java/com/astrazeneca/vardict/postprocessmodules/AmpliconPostProcessModule.java

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -164,8 +164,10 @@ public void process(Region rg, List<Map<Integer, Vars>> vars,
164164
if (goodmap.contains(format("%s-%s-%s", amp, vref.refallele, vref.varallele))) {
165165
continue;
166166
}
167+
if (instance().conf.doPileup && vref.refallele.equals(vref.varallele)) {
168+
continue;
169+
}
167170
if (vref.startPosition >= reg.insertStart && vref.endPosition <= reg.insertEnd) {
168-
169171
String regStr = reg.chr + ":" + reg.start + "-" + reg.end;
170172

171173
if (vars.get(amp).containsKey(position) && vars.get(amp).get(position).variants.size() > 0) {

src/main/java/com/astrazeneca/vardict/variations/Variant.java

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -291,7 +291,9 @@ private StringBuilder debugVariantsContent(String n) {
291291
*/
292292
public String varType() {
293293
Matcher mm = ANY_SV.matcher(varallele);
294-
if (refallele.length() == 1 && varallele.length() == 1) {
294+
if (refallele.equals(varallele) && refallele.length() == 1) {
295+
return "";
296+
} else if (refallele.length() == 1 && varallele.length() == 1) {
295297
return "SNV";
296298
} else if (mm.find()) {
297299
return mm.group(1);

0 commit comments

Comments
 (0)