Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 90,688 | 0 | T | C | 70.4% | 32.5 / 21.6 | 27 | F346L (TTC→CTC) | murE | UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase |
Reads supporting (aligned to +/- strand): ref base T (3/5); new base C (8/11); total (11/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCG > NZ_CP009273/90562‑90821 | aTGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTgcg < 1:558119/135‑1 (MQ=255) gcgcTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCaa < 1:39042/139‑1 (MQ=255) ttAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGa < 1:242299/138‑1 (MQ=255) ctgctCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTAcg > 1:66430/1‑139 (MQ=255) gCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCgg < 2:66430/139‑1 (MQ=255) tGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACgg > 1:294463/1‑102 (MQ=255) tGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACgg < 2:294463/102‑1 (MQ=255) tGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCACATACGCCGGATGCACTgg > 1:457753/1‑138 (MQ=255) cTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGa > 1:228115/1‑94 (MQ=255) cTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGa < 2:228115/94‑1 (MQ=255) cGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGa < 1:152127/92‑1 (MQ=255) cGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGa > 2:152127/1‑92 (MQ=255) tCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAg > 2:5043/1‑139 (MQ=255) gCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAggcgg > 1:52559/1‑136 (MQ=255) gCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAggcgg < 2:52559/136‑1 (MQ=255) gCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGc < 1:17936/72‑1 (MQ=255) gCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGc > 2:17936/1‑72 (MQ=255) gCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACgg < 2:285067/82‑1 (MQ=255) gCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACgg > 1:285067/1‑82 (MQ=255) gATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTc > 1:155566/1‑139 (MQ=255) gATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTc > 2:19903/1‑139 (MQ=255) ctgAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCa < 1:80947/137‑1 (MQ=255) ctgAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCa < 1:5043/137‑1 (MQ=255) ctgAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCa < 1:330046/137‑1 (MQ=255) aaaaCCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACtgtg < 2:141636/139‑1 (MQ=255) cAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTgtggt < 2:567338/136‑1 (MQ=255) aaGTGCTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATcg < 2:155566/139‑1 (MQ=255) | ATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCG > NZ_CP009273/90562‑90821 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |