| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,371,397 | 0 | A | C | 56.2% | 33.4 / 43.6 | 32 | K103N (AAA→AAC) | ycjS | Gfo/Idh/MocA family oxidoreductase |
| Reads supporting (aligned to +/- strand): ref base A (12/2); new base C (0/18); total (12/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.03e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTC > NZ_CP009273/1371313‑1371479 | ggTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAAccgcc > 2:503623/1‑90 (MQ=255) tAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAAccgcccgc > 1:241425/1‑90 (MQ=255) cGTCTGCTCACCTAACCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCAt < 2:119304/90‑1 (MQ=255) gTCTCCTCCCCTACCCGTTTTCATTACGAACATCCCCTGAGGGCACTGGAAGGGGGCTGCCTTGTGTTGTGCAAAAACCCCCCCGCCATg < 1:791429/90‑1 (MQ=255) agcTGCTCCCCTACCCGTTTTCATTACGAACATCCCCTGATGGCACTGGAAGCGGGCTCCCATGTGATGTGCGAAAACCCGCCCGCCATg < 2:1047838/88‑1 (MQ=255) gCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGc > 1:889014/1‑90 (MQ=255) aaCCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAAcct > 1:1047838/1‑88 (MQ=255) aaCCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAg > 2:876830/1‑90 (MQ=255) aaCCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAg > 2:1242064/1‑90 (MQ=255) aaCCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAg > 1:352708/1‑90 (MQ=255) tttttCATTACAAACATCCCCTGATGGCACGGGAGGGGGGCTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAGgcgc < 1:3384/89‑1 (MQ=255) aTTACGAACATCCCCTGATGGCACTGGAAGCGGGCTGCCATGGGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGaaa < 2:725908/90‑1 (MQ=255) tACGAACATCCCCTGTTGGCACGGGAAGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAAtg < 2:643558/90‑1 (MQ=255) ggATGGCACGGGAGGGGGGCTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACcgcgc < 1:136882/89‑1 (MQ=255) gCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCaaa > 2:363778/1‑90 (MQ=255) cACTGGAAGGGGGCTGCCATGTGATGTGCAAAAACCCCCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAAc < 1:876830/90‑1 (MQ=255) cTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATg < 2:1130084/45‑1 (MQ=255) cTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATg > 1:1130084/1‑45 (MQ=255) ggggggCGGCCGTGTGATGGGCAAAAACCCCCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCaa < 1:1355477/87‑1 (MQ=255) agggggCTGCCTTGTGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCaa < 2:765736/87‑1 (MQ=255) accGGGCTGCCGTGTGATGTGCAAAACCCCCCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCaa < 1:1103093/88‑1 (MQ=255) aGCGGGCTGCCATGTGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCaa < 2:363608/90‑1 (MQ=255) aGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCaa < 2:352708/90‑1 (MQ=255) gCTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGc < 2:1206417/90‑1 (MQ=255) gCTGCCATGTGATGTGCAAAAACCCCCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGc < 1:346941/90‑1 (MQ=255) cTGCCATGTGATGTGCAAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCt < 2:744332/90‑1 (MQ=255) cTGCCATGGGATGTGCAAAAACCCCCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCt < 2:1211113/90‑1 (MQ=255) cTGCAAGGTGATGTGCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCt < 2:828045/90‑1 (MQ=255) gtgatgtgCGAAAACCCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTAc < 2:241425/90‑1 (MQ=255) gtgCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGGCtt > 1:465945/1‑90 (MQ=255) gtgCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACtt > 1:451971/1‑90 (MQ=255) gCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTc > 1:1315905/1‑90 (MQ=255) | GGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTACGACTTTC > NZ_CP009273/1371313‑1371479 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |