| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 2,159,426 | 0 | G | A | 73.1% | 34.2 / 14.4 | 26 | intergenic (‑2179/‑1302) | USA300HOU_2042/mutS2 | hypothetical membrane protein/possible DNA mismatch repair protein MutS |
| Reads supporting (aligned to +/- strand): ref base G (3/4); new base A (11/8); total (14/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTACCCGTGCTTCAACAAGTGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTGGAAAATCGACCTTAATTAAAA > CP000730/2159289‑2159558 | ttACCCGTGCTTCAACAAGTGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATca > 1:288446/1‑141 (MQ=255) tACCCGTGCTTCAACAAGTGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTATATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGAcatcat > 2:420768/1‑141 (MQ=255) ccGTGCTTCAACAAGTGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTcatgca > 1:288/1‑141 (MQ=255) aacaaGTGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTaaaaaa < 2:58253/141‑1 (MQ=255) tGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAATCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGa < 1:139624/141‑1 (MQ=255) aaTGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAAt > 2:86182/1‑141 (MQ=255) ggCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCgg > 2:13628/1‑141 (MQ=255) aaaaGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTAt > 1:433679/1‑141 (MQ=255) aGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACAACAAGAAAGAAAAAAATTTAGAAATCGGTTATCAAACt < 2:323582/141‑1 (MQ=255) gACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTg > 1:249101/1‑141 (MQ=255) aCGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGc > 1:306798/1‑141 (MQ=255) gTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCaa < 2:381465/141‑1 (MQ=255) ttttAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCa < 2:288446/141‑1 (MQ=255) gAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCt > 1:437421/1‑141 (MQ=255) aaaaaTGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTAt < 1:98134/141‑1 (MQ=255) aaaaaTGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTAt < 1:98131/141‑1 (MQ=255) aaTGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAg > 1:95615/1‑141 (MQ=255) ttGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGtct < 1:72742/141‑1 (MQ=255) aCCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAgg < 2:451507/141‑1 (MQ=255) aaGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCAt > 2:217390/1‑141 (MQ=255) gTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTg > 1:448937/1‑141 (MQ=255) aaTATTCAATTTGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCa > 2:209426/1‑141 (MQ=255) ttGGCTTTGACCGGAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTg < 2:95615/141‑1 (MQ=255) gACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTGGAAAATCg > 1:283811/1‑141 (MQ=255) ggAATACAGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAAACTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTGGAAAATCGACCt < 2:323708/141‑1 (MQ=255) aGGTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTGGAAAATCGACCttaattaa < 2:398398/141‑1 (MQ=255) gTAATGACATCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATATAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTGGAAAATCGACCTTAATTaaaa > 2:216754/1‑141 (MQ=255) | TTACCCGTGCTTCAACAAGTGGAATGGCAAAAAGTAGACGTAAGATTTTAGAAAAAATGGAACGCATTGATAAACCAATGTTAGATGCCAAAAGTGCAAATATTCAATTTGGCTTTGACCGGAATACAGGTAATGACGTCATGCATGTAAAAAATTTAGAAATCGGTTATCAAACTGCAATTACCAAACCTATGAGTATAGAGGTCTCTAAAGGCGATCATATAGCAATCATTGGGCCAAATGGTATTGGAAAATCGACCTTAATTAAAA > CP000730/2159289‑2159558 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |