Number of fSNP is 0 but lead QTL within peak
Zepeng-Mu opened this issue · 0 comments
Zepeng-Mu commented
Hi, this is one example I got from my QTL mapping results using ATAC-seq with some of my own annotation
peak chr peak.start peak.end sid position ref alt allele.frq HWE ia.score log10.BH.q chi.sqr Pi Delta Phi overdispersion sid.within.the.region no.fsnp no.rsnp no.iter.null no.iter.alt location.ties Log.lik.null convergence.status sqr.corr.fsnp sqr.corr.rsnp perm.log10.BH.q perm.chi.sqr q.val emp.p.val emp.q.val midPoint dist qtl.type
chr4_87860220-87860720 chr4 87860220 87860720 rs79243012 87860512 T C 0.214286 4.640955 0.999705 -1.873293 9.048585 0.305516 0.01 0.5 11.88503 3.654591 0 20 6 3 87860512 3.824137 0 NA 0.994061 -0.0668922 0.8656526 0.9999999 0.0001357602 0.06012426 87860470 42 within
Basically the number of fSNP is 0 in the output, but the lead SNP found by RASQUAL is in the peak. Does this mean the fSNP is filtered out because there is not sufficient coverage, and then this test is only the linear QTL results just like SNPs that are outside of the peak?
Thanks