| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 887,959 | G→A | *179* (TAG→TAA) | rcdA → | transcriptional regulator of csgD and ybiJI, autoregulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 887,959 | 0 | G | A | 98.1% | 168.6 / ‑3.7 | 53 | coding (537/537 nt) | rcdA | transcriptional regulator of csgD and ybiJI, autoregulator |
| Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (26/26); minor base . (0/1); total (26/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01 | |||||||||||
TGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGC > NC_000913/887836‑888093 | tGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGa > 2:1370564/1‑139 (MQ=255) ggAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaa > 2:271014/1‑139 (MQ=255) gAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaa > 1:11986/1‑139 (MQ=255) gAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaa < 1:1218241/139‑1 (MQ=255) gAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGAATTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaa > 1:1297778/1‑139 (MQ=255) cccgcGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt > 2:280850/1‑115 (MQ=255) cccgcGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt < 1:280850/115‑1 (MQ=255) cccgcGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGAt < 2:78566/126‑1 (MQ=255) cccgcGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGAt > 1:78566/1‑126 (MQ=255) ccgcGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCAt < 2:1062700/139‑1 (MQ=255) ccgcGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCAt > 1:1037694/1‑139 (MQ=255) cgcgcgCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATc < 2:1516834/139‑1 (MQ=255) cGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGa < 1:271014/139‑1 (MQ=255) tATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATg > 1:322100/1‑139 (MQ=255) tATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTA‑AAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATgg < 2:1663989/139‑1 (MQ=255) tGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCg < 2:605113/139‑1 (MQ=255) gAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt > 1:1621319/1‑92 (MQ=255) gAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGAt < 2:1621319/92‑1 (MQ=255) gAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCgg < 2:1037694/139‑1 (MQ=255) accgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg < 1:714655/139‑1 (MQ=255) accgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg < 1:1370564/139‑1 (MQ=255) accgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg < 1:829237/139‑1 (MQ=255) ccgaccgTAAGCCGCTGTCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaaaa < 1:264313/87‑1 (MQ=255) ccgaccgTAAGCCGCTGTCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaaaa > 2:264313/1‑87 (MQ=255) accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACa > 2:1647054/1‑139 (MQ=255) ccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGTTGGGCTCCACAAAATGGGGACAt < 2:322100/139‑1 (MQ=255) tAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCaa > 1:805794/1‑139 (MQ=255) gCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGaa > 2:1715969/1‑94 (MQ=255) gCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGaa < 1:1715969/94‑1 (MQ=255) cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCGGGGTAATAGGTAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAg > 1:1106980/1‑139 (MQ=255) cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGAt > 2:55705/1‑64 (MQ=255) cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGAt < 1:55705/64‑1 (MQ=255) cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCca < 2:686281/115‑1 (MQ=255) cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCca > 1:686281/1‑115 (MQ=255) tATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCAt < 2:813969/77‑1 (MQ=255) tATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCAt > 1:813969/1‑77 (MQ=255) tCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAg < 2:11986/139‑1 (MQ=255) cgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTc < 1:345599/105‑1 (MQ=255) cgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTc > 2:345599/1‑105 (MQ=255) tGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt < 2:269080/82‑1 (MQ=255) tGGTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt > 1:269080/1‑82 (MQ=255) gTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaa < 1:1344812/44‑1 (MQ=255) gTTGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACaaaa > 2:1344812/1‑44 (MQ=255) tGAGAGGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTg < 2:805794/139‑1 (MQ=255) agagGGTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCg > 1:489332/1‑139 (MQ=255) gggTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGc < 1:398622/77‑1 (MQ=255) gggTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGc > 2:398622/1‑77 (MQ=255) gTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa < 1:692382/111‑1 (MQ=255) gTTGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa > 2:692382/1‑111 (MQ=255) ttGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcaca < 1:476087/80‑1 (MQ=255) ttGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCcaca > 2:476087/1‑80 (MQ=255) tGCAGGGTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGt > 1:952307/1‑139 (MQ=255) ggTAATAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGc > 1:1515460/1‑139 (MQ=255) | TGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGC > NC_000913/887836‑888093 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |