Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 2,308,304 | 0 | C | G | 30.0% | 54.0 / 11.1 | 40 | A230P (GCC→CCC) | fucA | L‑fuculose phosphate aldolase |
Reads supporting (aligned to +/- strand): ref base C (11/17); new base G (11/1); total (22/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.33e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGC > ACB143/2308168‑2308438 | cAGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGACATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTCGGGGTGGAAATCCAGTCGTGGGc > 1:240730/1‑137 (MQ=255) cggCACCTTTGCCGGTGTTGAGGAAAGAGTACAGCTGCGCTCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc < 4:236050/136‑1 (MQ=255) ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCGCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt < 7:252904/137‑1 (MQ=255) ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCttt > 6:261801/1‑137 (MQ=255) ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCACGCAGTCGCCATGGGCACGCAAGGCACGGCGCCCGTAGTAGCTGAAGGCAGCACCGTGAAGTTAGGGGGTGGAAATCCAGTCGAGGGCCTCttt > 7:185870/1‑137 (MQ=255) gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCTGCACCGTGACGTTTGGGGGTGTAAATCCAGTCGTGGGGCTCTTTTGCCAGTGGTTAgg > 8:216093/1‑137 (MQ=255) gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCCAGGCTTCgg > 8:221269/1‑137 (MQ=255) gAGGAAAGAGTACAGCTCCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGGCGTGGGCCTCTTTGCCCAGTGCTTCgg < 6:201442/137‑1 (MQ=255) aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGTGCTTCGGGGATGGGGCgg > 8:2741/1‑137 (MQ=255) cGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCGTTGCCCAGGGCTTCGGGGATGGGGCGGCTGTCAc > 6:109749/1‑137 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGCCCAGTGCGGCGGGGATGGGGGGGGTGTGACCgg > 1:41199/1‑137 (MQ=255) aGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCACTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGgcg > 5:334522/1‑137 (MQ=255) aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCAt < 8:155121/137‑1 (MQ=255) aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGGTTTGCCCCGTGCGTGGGGGTTGGGGGGGGTGTCTCCGGCGCCCAt > 4:235538/1‑137 (MQ=255) cGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagc < 7:168740/137‑1 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGTGTTGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagcag < 8:202735/137‑1 (MQ=255) cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGTCGGTTGGCAGTGCTTCGGGGGTGGGGCGGGTGTTAACGGCCCCCCAGGCcagcaga > 1:165213/1‑136 (MQ=255) cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGGCCAGGGGTTCGGGGATTGGGCGGATGTCACCGGCGCCCCGGGCCAGCAGt > 8:168740/1‑137 (MQ=255) cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCTTTGGCCAGGGCGTGGGGGCTGGGGCGGGTGTGACCGGCGCCCCAGGCccgccct > 5:117149/1‑134 (MQ=255) aaGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTGGCCAGGGCTTAGGGGCTGGGGCGGCTGTCACCGGCGCCCCAGGCCAGCACTTGCAGCTTtgg > 3:270749/1‑136 (MQ=255) ggCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGGGTTGCCCAGTGGGTGGGGGATGGGGGGGCTGTCAACGGCGCCCATGGCCCGCAGTTGGCTCTTtgctg > 2:346898/1‑137 (MQ=255) cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg < 2:227774/137‑1 (MQ=255) ggCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcgcgc < 5:109749/137‑1 (MQ=255) cgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGTCGTTGTCCACTGCTTGGGGGGTGGGGCGGATGTCACCGGCCCCCAAGGGCAGCAGTTGCATCTTTGCTTCGCGCTCgc > 4:226803/1‑136 (MQ=255) gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGaa < 7:221269/137‑1 (MQ=255) gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGaa < 7:216093/137‑1 (MQ=255) agtagCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCAGGGGTTGGGGGGTGGGGCGGGTGTCACCGGCGCCCCTGGGCCGCAGGTGGATCTTTTCTTCGCGGGCGCAcagc > 6:4002/1‑137 (MQ=255) gCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCagag < 4:142540/137‑1 (MQ=255) aGGCAGCTCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTCTTTGCCGAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGGAGTTGCATCTTTGCTGCGCGCTCGAACAGCAgagcgagc < 5:193345/137‑1 (MQ=255) cagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGTGCGCCGCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACg < 5:258964/137‑1 (MQ=255) cGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGGCCAGTGCTTCCGGGTTGGGGCGGCTGTCCCCGGCGCCCCAGGCCCCCACTTTGATCATTGCTGCGCGCTCGAAAACAACAGCCAGCCCGCAGGCt > 8:381019/1‑137 (MQ=255) cGTTTGGGGGTGGAAATCCAGTCGTGGGCCGCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCCGCGCCCAAGGCCCGCAGTTGGATCTTTTCTGCGCGCTCGAACACCAGAGCGAGCACGCAGGCTTCCtc > 2:100595/1‑137 (MQ=255) cGTTTGGGGGTGGAAATCCAGACGTGGGCCTCTTTGCCCAGAGCTTCGGGGATGGGGCGGATGTAACCGGCGCCCATGGCCAGCAGTTGCATCTTTTCTGCGCGCGCGAACACCACAGCGAGCACCCAGGCTTCCtc > 6:61051/1‑137 (MQ=255) ttGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCtcgat < 4:75284/137‑1 (MQ=255) tGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTGGGGGGTGGGGCGGGTGTCACCGGCGCCCCTGGCCAGCAGTTGCATCTTTTCTGCGCGCTCGCACCGCAGAGCGAGCACGCCGGCTTCCGCGATCGAACGc > 8:123564/1‑137 (MQ=255) gggAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGc < 2:165213/136‑1 (MQ=255) cAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt < 7:72008/137‑1 (MQ=255) aGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATc < 1:100595/137‑1 (MQ=255) tCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGGGGATGTCACCGGCGCCCAAGGCCACCACTTGGATCTTTGCTGCGCCCACGCAACGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATcag > 6:183915/1‑137 (MQ=255) cGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATcagc < 3:241346/137‑1 (MQ=255) tGGGGCTGTTTGGCCAGTGGGTCGGGGGTGGGGCGGATGTCACCGGCGCCCATGGCCCGCAGTTGCATCTTTGCTGCGCGCGAGACCAGCAGAGCGCGCACGCAGGCTTGCTCGATCGAACGCCCCGCGCTcagcag > 2:88285/1‑137 (MQ=255) ggCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGc < 4:184297/137‑1 (MQ=255) | CAGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTCTTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGC > ACB143/2308168‑2308438 |
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: ‑ |