Skip to content

Commit 20b50a4

Browse files
committed
SexKaryotyper: Fix absolute coverage thresholds for female detection
1 parent 5e941bd commit 20b50a4

File tree

1 file changed

+5
-4
lines changed

1 file changed

+5
-4
lines changed

src/main/groovy/gngs/SexKaryotyper.groovy

Lines changed: 5 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -112,10 +112,6 @@ class SexKaryotyper implements Runnable {
112112

113113
log.info "Means for $sample are: " + [chrX: xCoverage, chrY: yCoverage, autosomes: autosomeCoverage]
114114

115-
if(yCoverage < 5 && xCoverage > 30) {
116-
sex = Sex.FEMALE
117-
}
118-
else
119115
if(autosomeCoverage < 1.0) {
120116
log.warning("Autosomal coverage too low to infer sex: $autosomeCoverage")
121117
sex = Sex.OTHER
@@ -125,6 +121,11 @@ class SexKaryotyper implements Runnable {
125121
log.info "chrX / auto = " + (xCoverage / autosomeCoverage)
126122
sex = Sex.MALE
127123
}
124+
else
125+
if(yCoverage / autosomeCoverage < 0.125) {
126+
log.info "chrY / auto = " + (yCoverage / autosomeCoverage)
127+
sex = Sex.FEMALE
128+
}
128129
else
129130
sex = Sex.OTHER
130131
}

0 commit comments

Comments
 (0)