Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 49,765 | C→T | intergenic (+134/‑58) | kefC → / → folA | K(+) : H(+) antiporter KefC/dihydrofolate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 49,765 | 0 | C | T | 98.0% | 179.2 / ‑5.3 | 50 | intergenic (+134/‑58) | kefC/folA | K(+) : H(+) antiporter KefC/dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (26/23); total (27/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA > NC_000913/49703‑49826 | ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga > 1:924134/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga < 1:2465180/65‑1 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga > 1:2330130/1‑65 (MQ=255) ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga > 1:3403084/1‑65 (MQ=255) gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac > 1:1107266/1‑65 (MQ=255) gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga > 1:3445662/1‑65 (MQ=255) agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac > 1:2518516/1‑65 (MQ=255) agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac > 1:1826854/1‑65 (MQ=255) agCGTGGACTCGTCAGCAGAATATAAAATCTTCCTCAACATCATCCTCGCACCAGTTGACGACgg < 1:2476109/65‑1 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCGCCAGTTGACGACgg < 1:852665/65‑1 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg > 1:2133899/1‑65 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg < 1:3568109/65‑1 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg < 1:2062600/65‑1 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg < 1:876923/65‑1 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg < 1:1586325/65‑1 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg > 1:905725/1‑65 (MQ=255) agCGTGGACTCGCCAGCAGAATATAAAAGATTCCTCAACATCATCCTCGCACCAGTTGACGACgg < 1:2453522/65‑1 (MQ=255) gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGANGACGGt < 1:2222536/65‑1 (MQ=255) cGTGGACTCGCCAGCAGACTATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGtt > 1:3486118/1‑65 (MQ=255) ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc > 1:924289/1‑65 (MQ=255) gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg < 1:2945424/65‑1 (MQ=255) aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc > 1:513921/1‑65 (MQ=255) aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc > 1:71573/1‑65 (MQ=255) aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc > 1:2262856/1‑65 (MQ=255) aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc > 1:1865004/1‑65 (MQ=255) aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc > 1:3203408/1‑65 (MQ=255) aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc > 1:1451633/1‑65 (MQ=255) cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt < 1:2762054/65‑1 (MQ=255) cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt < 1:136489/65‑1 (MQ=255) cTCGCCAGCAGAATATAAAATTGTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt > 1:2442179/1‑65 (MQ=255) gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa < 1:3715839/65‑1 (MQ=255) ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGTTTTAc < 1:2187515/65‑1 (MQ=255) ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc < 1:2743213/65‑1 (MQ=255) ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc < 1:573473/65‑1 (MQ=255) ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc < 1:306120/65‑1 (MQ=255) ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc > 1:3029424/1‑65 (MQ=255) cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtat < 1:2886375/65‑1 (MQ=255) agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata > 1:3137606/1‑65 (MQ=255) agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata > 1:2318376/1‑65 (MQ=255) agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata > 1:1800530/1‑65 (MQ=255) agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata > 1:2173054/1‑65 (MQ=255) agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata > 1:3647676/1‑65 (MQ=255) taAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg < 1:772811/65‑1 (MQ=255) aaaTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGAc > 1:2941005/1‑65 (MQ=255) aTTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACaa > 1:1963593/1‑65 (MQ=255) aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa > 1:3470065/1‑65 (MQ=255) aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGAGGCGaaaa > 1:211904/1‑65 (MQ=255) aTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACaa > 1:2111297/1‑65 (MQ=255) ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa > 1:2128020/1‑65 (MQ=255) ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa > 1:2361068/1‑65 (MQ=255) cTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt < 1:861640/65‑1 (MQ=255) tCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc < 1:2713087/65‑1 (MQ=255) tCGCACCAGTTGACGACGGGTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc < 1:1175250/65‑1 (MQ=255) gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc < 1:513326/65‑1 (MQ=255) gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc < 1:1548249/65‑1 (MQ=255) aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTTaa > 1:652439/1‑65 (MQ=255) aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCaa < 1:3475203/65‑1 (MQ=255) aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCaa > 1:2995550/1‑65 (MQ=255) aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa > 1:50520/1‑65 (MQ=255) | GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA > NC_000913/49703‑49826 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |