Skip to content

Commit e6afc0a

Browse files
committed
bugfix handling requirePE for mid/frag
The logic for restricting alignment records to properly paired when midpoint or fragment are specified was incorrect. This commit fixes the error #70
1 parent 14e046b commit e6afc0a

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

src/scripts/Read_Analysis/PileupScripts/PileupExtract.java

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -105,7 +105,7 @@ else if(param.getTrans() == 2) {
105105
// Apply insert size filters
106106
if(Math.abs(sr.getInferredInsertSize()) < param.getMinInsert() && param.getMinInsert() != -9999) { continue; } //Test for MIN insert size cutoff here
107107
if(Math.abs(sr.getInferredInsertSize()) > param.getMaxInsert() && param.getMaxInsert() != -9999) { continue; } //Test for MAX insert size cutoff here
108-
} else if(param.getAspect() < 1) { continue; } // Make sure that midpoint/fragment pileup must come from properly paired read
108+
} else if(param.getAspect() == 2 || param.getAspect() == 3) { continue; } // Make sure that midpoint and fragment pileup must come from properly paired read
109109

110110
// Shift as needed
111111
if(sr.getReadNegativeStrandFlag()) { mark -= SHIFT; }

0 commit comments

Comments
 (0)