#!/bin/bash -Ceuo pipefail
samtools idxstats 1029_BSZ8-N1-BDNA-1_235YY7LT3_1.recalibrated.bam > idxstats.txt
x_map=$(grep -E "^chrX\s|^X\s" idxstats.txt | cut -f 3)
x_len=$(grep -E "^chrX\s|^X\s" idxstats.txt | cut -f 2)
y_map=$(grep -E "^chrY\s|^Y\s" idxstats.txt | cut -f 3)
y_len=$(grep -E "^chrY\s|^Y\s" idxstats.txt | cut -f 2)
awk -v xm="$x_map" -v xl="$x_len" -v ym="$y_map" -v yl="$y_len" '
BEGIN {
if (ym == 0 || yl == 0) {
print "XX"
} else {
x_cov = xm / xl
y_cov = ym / yl
ratio = x_cov / y_cov
if (ratio > 50.0) print "XX"
else print "XY"
}
}' > 1029_BSZ8-N1-BDNA-1_235YY7LT3_1.sex.txt