Skip to content

Number of fSNP is 0 but lead QTL within peak #52

@Zepeng-Mu

Description

@Zepeng-Mu

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions