Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 762,147 | 0 | C | A | 56.0% | ‑0.1 / 22.6 | 25 | L216L (CTC→CTA) | mngA | PTS 2‑O‑a‑mannosyl‑D‑glycerate transporter subunit IIABC |
Reads supporting (aligned to +/- strand): ref base C (8/3); new base A (0/14); total (8/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTC > NZ_CP009273/762059‑762217 | aaaaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCa < 1:1168370/90‑1 (MQ=255) aaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAAt < 2:1372537/90‑1 (MQ=255) gCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTAttt > 2:1610629/1‑90 (MQ=255) aaaTTTTTACGTAAAAAAACAAGGCGCTAACGGCATTGAAGGGCTTTTAACGGGGGATCAACTAAATAGTGCAACCGCCTGTATTTTTGc < 2:542401/90‑1 (MQ=255) ttACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGa > 2:797257/1‑90 (MQ=255) gtaaaaaaaCAAGGCGCTAACGGCATTGAGGGGCTTTAAAGGGGGGACCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTc < 1:911797/86‑1 (MQ=255) gTTGAAAACCAGGGCCCAACGGGGATGGAGGGGGGTTAAAGGGCGAATAAACCAAATAGTGAAACCGCCTGTATTTTTGCGGCTGAAGTc < 2:1847482/90‑1 (MQ=255) taaaaaaCAGGGCGCAAACGGCTTTGAGGGGCTTTAAAGGGGGGATAAACAAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc < 1:716653/88‑1 (MQ=255) taaaaaaCAAGGCGCAAACGGCATTGAGGGGCTTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc < 2:618934/88‑1 (MQ=255) tGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc > 2:1512855/1‑90 (MQ=255) tGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGAATAACTCAATAATGCAAACCCCTGTATTTTTTTGGCTGGAggggc > 1:52832/1‑86 (MQ=255) aaaaaaaaAGGGCGCAACCGGCATTGAGGGGCTTTTAAGGGGGGATAACCTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc < 2:312860/88‑1 (MQ=255) aaaaaaaaAGGGCCCTACCGGAATGGAGGGGTTTTAAACGGGGAACAAACAAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc < 1:1610629/88‑1 (MQ=255) aaaaaaaCAGGGCGCAACCGGCATGGAGGGGCTTTAAGCGGGGGACCACCCAAATGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc < 1:1262775/88‑1 (MQ=255) aaaaaaaCAGGCCGCAAACGGAATGGAGGGGCTTTAAAGGGGGAACCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGc < 1:797257/88‑1 (MQ=255) gAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTc < 1:711208/51‑1 (MQ=255) gAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTc > 2:711208/1‑51 (MQ=255) aaaCAGGGCGCAACGGGCATGGAGGGGTTTTAAGGGGGGGACCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATc < 1:210782/90‑1 (MQ=255) aaGGCGCAACGGGAATGGAGGGGCTTTTAAGGGGGGATCAACTAAATGGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAgg < 2:1072230/90‑1 (MQ=255) aaagggCTTTGAGGGGTTTTAAACGGGGGATCAACTAAATAGTCCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGa < 2:1108087/86‑1 (MQ=255) agggCTTTGAGGGGTTTTAAAGGGGGGATAAACTAAATAGTCAAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGc < 2:1658011/88‑1 (MQ=255) agggCATTGAGGGGCTTTAAACGGGGGATCAACTAAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGc < 2:1414423/88‑1 (MQ=255) ttGAAGGGCGTTTAACGGCGGATCAACTCAAAAATGCAAACCCCCGGATTTTTTCGGCTGGAGGCGCCAACCAGGGGAGGGGGCGgtttt > 1:1370382/1‑89 (MQ=255) tGAGGGGTGTTAAACGGCGAAAAAACTAAATGGTCCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTaa < 1:1764244/90‑1 (MQ=255) gCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGGAGGCGCCCCCAAGGGGAGGGGGGGGTTTTATGGCAt > 1:1215568/1‑90 (MQ=255) gTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGGCGCCATCAAAGAGAGGGAGCGGTTTTATGGCAttt > 1:458426/1‑89 (MQ=255) gTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTc > 2:1309212/1‑90 (MQ=255) | AAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTC > NZ_CP009273/762059‑762217 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |