| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,571,380 | 0 | A | C | 55.6% | 13.9 / 21.6 | 18 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base A (6/2); new base C (0/10); total (6/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAAT > NZ_CP009273/1571304‑1571464 | gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:268913/1‑90 (MQ=255) tCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAg < 2:389196/90‑1 (MQ=255) gAACGCGAAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGa > 1:451903/1‑90 (MQ=255) aaCGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGAc < 1:154220/90‑1 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGc < 2:414003/82‑1 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGc > 1:414003/1‑82 (MQ=255) cGATACCCCAATTCCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 2:17343/90‑1 (MQ=255) cGAAAACCCAATTTCCCCTTTCTTTCAACAGGGTTGAAGCTTTTTATTTTTGCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 1:163368/90‑1 (MQ=255) atcaGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTcac > 2:58878/1‑90 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:331013/90‑1 (MQ=255) ttGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTcc > 2:444843/1‑90 (MQ=255) ttGATTTTTGAATCAATGCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:19670/90‑1 (MQ=255) ttGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:175336/90‑1 (MQ=255) ttGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:268913/90‑1 (MQ=255) ttGATCTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 1:254898/90‑1 (MQ=255) tttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa > 1:140576/1‑90 (MQ=255) tGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAg < 2:182048/90‑1 (MQ=255) cAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat < 2:157757/90‑1 (MQ=255) | GTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAAT > NZ_CP009273/1571304‑1571464 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |