| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_007779 | 1,164,944 | 0 | T | G | 42.4% | 43.2 / 11.7 | 34 | S36R (AGT→AGG) | ycfN | thiamin kinase |
| Reads supporting (aligned to +/- strand): ref base T (2/17); new base G (14/0); total (17/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.47e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCAC > NC_007779/1164850‑1165044 | gcaATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTtttc > 1:1132066/1‑101 (MQ=255) ccATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACAt < 3:564492/101‑1 (MQ=255) ccATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACAt < 3:104629/101‑1 (MQ=255) gacgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTCTCATTTAACATCAGGGCCAGCg > 3:122237/1‑101 (MQ=255) acgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGt < 3:1531858/101‑1 (MQ=255) gaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTCTCATTGAACATCAGGGCCCGCGgtt > 1:415932/1‑101 (MQ=255) tGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTtgtg < 4:58084/101‑1 (MQ=255) tCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCg < 1:435047/101‑1 (MQ=255) cgcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGGGTTTTTATTTAAGATCAAGAGCACCGCTttgtgttggtgtt > 3:1398922/1‑92 (MQ=255) gcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTCTCCTTTAACATCAGGGCCAGCGGTTTGTTGTGCGTCAg > 1:352618/1‑101 (MQ=255) tttttCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGcc < 3:237353/101‑1 (MQ=255) tttttCCCGCAGTATCATCCCGTCACGACGTCTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGcc < 3:1588786/101‑1 (MQ=255) tCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTGTTCTCATTTGACATCATGACCAACGTATTGTTGTTGGTTAGCCGCAc > 2:74089/1‑101 (MQ=255) tCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGGTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCAc > 2:128660/1‑101 (MQ=255) gCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGATTTCTCCTTTAAACAAAGAGCCAGCGCTTTGTTTTGCGTCTGCCGCAcagcc > 1:1389478/1‑97 (MQ=255) aGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTCATTTAAAATCAAGGCCCGCGGTTTTTTTTGGGTCAGCCGCCCGCTccc > 3:34076/1‑100 (MQ=255) atcCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTCTCATTTTACATCAGGGCCCGCGTTTTGTTGTGCGTCAGGCGCCCCATCCTCGTgcg > 3:662142/1‑101 (MQ=255) ccGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTCTTTTAACAACAGGGCCAGCGTTTTGTTTTGCGTCAGCCGCACGATCCTCATGcgccc > 3:1352042/1‑100 (MQ=255) ccGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTCTCCTTTAACATCAGGGCCAGCGGTTTTTTGTGCGTCAGCCGCACGATCCTGATGcgccc > 3:328351/1‑100 (MQ=255) cGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTCATTGAACATCAGGGCCAGCGGTTTGTTGTGGGTCAGCCGCACCAGCCCGATGcgccgc > 1:37205/1‑101 (MQ=255) tCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGTTTTCTCATTTAAAAACAAGGCCAAAGTTTTGTTTTGCGTCAGCCGCCCGAACCTCCTGCGCCGCAg > 1:1407680/1‑101 (MQ=255) gacgTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTCATTGAAAATCAGGGGCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGGCcg > 3:671215/1‑101 (MQ=255) cgTTTAATAGTGGGCTTAGTGGCGGGAGGTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCcgcg > 1:256893/1‑101 (MQ=255) ttAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTc < 1:179251/101‑1 (MQ=255) tAATAGTGGGCTTAGTGGCGGGGGGTTTTTTATTTAAAAATAGGGCCCGCGCTTTTTTTTTCGGCAGCCGCCCCACCCCGATGCGCCGCAGGcccccgtct > 3:810149/1‑94 (MQ=255) tAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCt < 3:451250/101‑1 (MQ=255) aaTAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCtt < 1:121023/101‑1 (MQ=255) gggCTTAGTGGCGGGATTTTTCTCATTGGACATCAGGGGCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCCAATGCGCCGCACTCCGCGTTCTTGCGCCa < 1:331102/101‑1 (MQ=255) tAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCCGTCCGCGTTCTTGCGCCAGTATc < 3:1596544/101‑1 (MQ=255) gCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCt < 1:1396723/101‑1 (MQ=255) gCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCt < 1:1142114/101‑1 (MQ=255) gCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCt < 1:490758/101‑1 (MQ=255) cGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCtt < 1:805343/101‑1 (MQ=255) ggAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTa < 1:597243/101‑1 (MQ=255) ttttCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATcac < 3:1077265/101‑1 (MQ=255) ttttCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATcac > 3:938952/1‑101 (MQ=255) ttttCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATcac > 3:1041185/1‑101 (MQ=255) | GCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGAGTTTTCTCATTGAACATCAGGGCCAGCGTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCAC > NC_007779/1164850‑1165044 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |