natsuhiko/rasqual

Number of fSNP is 0 but lead QTL within peak

Zepeng-Mu opened this issue · 0 comments

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