Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,154,068 | 0 | T | G | 52.0% | 12.1 / 21.3 | 26 | G248G (GGT→GGG) | ptsG | PTS glucose transporter subunit IIBC |
Reads supporting (aligned to +/- strand): ref base T (2/10); new base G (13/0); total (16/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.02e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAG > NZ_CP009273/1153982‑1154147 | aTGCAGATTGGTGATTACAGCAACGCAGGAGGGCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa < 1:992764/90‑1 (MQ=255) aTGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCCACTGCGGGTgaa > 2:827274/1‑90 (MQ=38) gCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACt < 1:629778/90‑1 (MQ=255) gCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGaaaat > 2:92868/1‑88 (MQ=255) aTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGGTGGCtt > 1:1459036/1‑90 (MQ=255) aTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCCACTGCGGGAAAACTGTCTGGTGGCtt > 2:1379366/1‑90 (MQ=255) aTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAACTGGCTGGGGGCtt > 1:776143/1‑90 (MQ=255) tACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGGGGGCTTc > 1:152948/1‑90 (MQ=255) tACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCCACCGCGGGGAAACTGTCTGGGGGCTTc > 2:591436/1‑90 (MQ=255) acacCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTcc < 2:1751589/90‑1 (MQ=255) cacCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCt < 1:827274/90‑1 (MQ=255) ccAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCCACTGCGGGGAAACTATCTGGTGGCTTccctt > 2:588519/1‑87 (MQ=255) aCGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAAATGTCTGGTGGCTTCCTTTTCa > 2:343322/1‑90 (MQ=255) cGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaa < 1:1707927/90‑1 (MQ=255) gcagGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCCACTGCGGGGAAAATGTCTGGGGGGTTCCTGTTTAAAATg > 2:1707927/1‑90 (MQ=255) aggtcaggtTTTCCACGGCGACATTACGCGTTATATGGCGGGGGACCCCACTGCGGGGAAACTGTCTTGGGGCTTCCCGTTCAAAATGTa > 1:1011870/1‑90 (MQ=255) ggtcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAACTGTCTGGGGGGTTCCTGTTTAAAATGTaa > 2:1509343/1‑89 (MQ=255) tcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGGGGGCTTCCTGTTCAAAATGTACgg > 1:1311828/1‑90 (MQ=255) gCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCg < 1:567158/90‑1 (MQ=255) cGGGTGACCCGACTGCGGGGAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTg > 2:463203/1‑90 (MQ=255) tGACCCGACTGCGGGTAAACTGTCTGTTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa < 1:656806/90‑1 (MQ=255) tGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa < 2:67809/90‑1 (MQ=255) tGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa < 1:21834/90‑1 (MQ=255) tGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa < 1:1639385/90‑1 (MQ=255) gACCCCACTGCGGGGAAAATGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaaa > 2:1045931/1‑90 (MQ=255) cGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAg > 2:144156/1‑90 (MQ=255) | ATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAG > NZ_CP009273/1153982‑1154147 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |