Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 2,308,336 | 0 | C | G | 28.9% | 75.4 / 10.7 | 46 | A225P (GCA→CCA) | fucA | L‑fuculose phosphate aldolase |
Reads supporting (aligned to +/- strand): ref base C (10/22); new base G (13/0); total (24/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAG > KT2440/2308212‑2308468 | gagTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATggg < 7:338474/137‑1 (MQ=255) cagTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGGACCGTGTCGTTTGGGGGTGGAACTCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATggg < 4:55179/136‑1 (MQ=255) gTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTTCCCAGTGCTTCGGGGCTtgggg > 4:60352/1‑136 (MQ=255) aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCGTTGGCCAGTGGTTAGGGGATGGGGCgg > 4:278987/1‑137 (MQ=255) aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCAGTGAGTACGGGCTGGGGCgg > 8:184304/1‑137 (MQ=255) cGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCAc < 5:391277/137‑1 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGTGGTTGGGGGGTGGGGGGGGTGTGACCgg > 3:134554/1‑137 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGGGTTGGCCAGGGGGTCGGGGCTGGGGGGGGTGTGCCCgg > 3:62484/1‑137 (MQ=255) aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCTTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCAt < 5:107109/137‑1 (MQ=255) ggCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATg < 3:250839/137‑1 (MQ=255) aGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGTCCCGTGGTTCGGGGGTGGGGCGGCGGGGACCGGCGccccggggc > 3:221866/1‑131 (MQ=255) gTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGCCCCGGGGTTGGGGGATGGGGGGGATGTGACCGGCGCCCCTGGccc > 7:303760/1‑136 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGGGTCGGGGATGGGGCGGCGGTCCCCCGCGCCCCAGGCcaccag > 8:304868/1‑137 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGGCCAGGGGTTCGGGGTTGGGGCGGCTGTCGCCGGCGCCCCAGGcccggcg > 3:384728/1‑132 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCAGTGCTTCGGGGATGGGGGGGGTGTCACCGGCGCCCATGGCcagcag > 8:60320/1‑137 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCGTTTCCCAGTGGTGGGGGGGTGGGGGGGTTGTGACCGGCGCCCACTGCcagcag > 4:269046/1‑137 (MQ=255) cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCGGTGCGCACTGCGTGGGGGGTGGGGCGGCTGTCACCCGCCCCCCAGGCccgcacg > 2:107409/1‑135 (MQ=255) cAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTtg < 3:171026/137‑1 (MQ=255) cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGAGGTTTGGGGGTGGACATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGTGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg < 4:122056/137‑1 (MQ=255) cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg < 5:152488/137‑1 (MQ=255) cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg < 4:353966/137‑1 (MQ=255) cGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGTGGTTGGGGGATGGGGGGGATGTCACCGGCGCCCCAGGCCAGCAGTTGCATAGTTTCTgcgcg > 6:61938/1‑137 (MQ=255) ggCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc < 8:62248/137‑1 (MQ=255) ggCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc < 1:107409/137‑1 (MQ=255) ggCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGACGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc < 4:103991/137‑1 (MQ=255) gcgcgcGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCt < 6:340000/137‑1 (MQ=255) gcgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGCTGTCACCGGCGCCCATGGGCAGCAGTTGGATCTTTTCTGCGCGCTCg > 5:156047/1‑137 (MQ=255) cgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGGTGGCCAGGGCTTGGGGGTTGGGGCGGCTGTCACCGGCGCCCAAGGCCACCAGGTGCATCTTTTATgcgcgggcgc > 3:32606/1‑132 (MQ=255) gcagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGGCGTTGGGCAGTGGGTGGGGGGTGGGGGGGGTGTCACCGGCGCCCCTGGGCAGCACGTGCATCTTTGCTGGGCGCGCGAACAGCCGAGCGCGCAc > 4:23700/1‑137 (MQ=255) gcagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTTGCCAGGGCTTGGGGGATGGGGCGGGTGTCACCGGCGCCCCAGGCCAGCAGGTGCATCATTGCTGAGCGGTCGAACAGCCgcgcgcgccc > 6:89411/1‑135 (MQ=255) agcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGGCCAGTGGTTCGGGGATGGGGCGGGTGTGACCGGCGCCCCAGGCCCGCAGTTGGATCTTTGCTGCGCGCGCGCACAGCAGAGCGCGCACGc > 5:341587/1‑137 (MQ=255) agcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCGTTGCCCACTGGGTCGGGGGTGGGGCGGATGTGACCGGCGCCCCTGGCCAGGCGTTGGATCTTTGCTTCGCGGGCGCACACCACCGCGAGCCCGc > 8:116332/1‑137 (MQ=255) gcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGGGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCa < 1:114403/137‑1 (MQ=255) cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGGCCAGTGCTTCGGGGATTGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTTCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCCGGCt > 5:180908/1‑137 (MQ=255) cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGgcg > 2:67009/1‑77 (MQ=255) cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGgcg < 1:67009/77‑1 (MQ=255) cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTTGCCCGTGGTTCGGGGCTGGGGCGGATGTCACCGGCGCCCATGGCCCGCAGTTGCATCTTTTCTGCGCGCTCGAACCGCAGAGCGAGCACGCCGGcg > 5:134112/1‑136 (MQ=255) gTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCtt < 3:60352/137‑1 (MQ=255) aaaTCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCa > 7:245748/1‑137 (MQ=255) cAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt < 2:126357/137‑1 (MQ=255) gTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATca < 8:245748/137‑1 (MQ=255) gTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATcagca < 7:304868/137‑1 (MQ=255) ctTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCatga < 6:370044/137‑1 (MQ=255) ctTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCatga < 4:350484/137‑1 (MQ=255) tGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATgg > 3:164747/1‑137 (MQ=255) tgTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAg < 4:2372/136‑1 (MQ=255) | GAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAG > KT2440/2308212‑2308468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |