| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 520,722 | 0 | T | G | 30.3% | 53.5 / 10.0 | 33 | V669G (GTA→GGA) | rhsD | rhs element protein RhsD |
| Reads supporting (aligned to +/- strand): ref base T (8/15); new base G (10/0); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.05e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTC > NZ_CP009273/520637‑520810 | tACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGAcag > 1:486491/1‑90 (MQ=255) cGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGgcagca > 2:74173/1‑90 (MQ=255) cGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGAcagca > 1:576093/1‑90 (MQ=255) aTGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑gg > 1:265421/1‑90 (MQ=255) gctgcCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑Gggc < 2:171543/90‑1 (MQ=255) ctgcCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGGCAGCA‑ggggg > 2:182580/1‑88 (MQ=255) ctgcCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGACAGCA‑Gggcg > 1:143306/1‑90 (MQ=255) ctgcCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGACACCG‑Gggcg > 2:492365/1‑90 (MQ=255) tgtCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGGCAGCAGGGGGGGCTTTgcgg > 2:354598/1‑87 (MQ=255) tcatcaGCTTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGc < 1:492365/90‑1 (MQ=255) tcatcaGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGc < 1:182580/90‑1 (MQ=255) cTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCaa < 1:211792/90‑1 (MQ=255) cTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCaa < 2:576093/90‑1 (MQ=255) aaCGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTAt < 1:138869/90‑1 (MQ=255) aaCGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGGACAGCA‑GGGGGGCTTTTTCGGGGGGGCGCAACGGTTTCATTAt > 2:488672/1‑90 (MQ=255) gCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGaa > 1:49597/1‑90 (MQ=255) cGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGaaa < 2:265421/90‑1 (MQ=255) ttCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTc < 2:155946/90‑1 (MQ=255) gTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac < 1:600415/90‑1 (MQ=255) gTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac < 1:588660/90‑1 (MQ=255) gTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac < 2:331125/90‑1 (MQ=255) aCGATGCGCTGGACCGGCTGGGACAGCA‑GGGGGGCTTTGACGGGGGGGCGCAACGTTATCATTATGACCTGGCCGGAAAAATCACACaga > 2:600415/1‑90 (MQ=255) gATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGt > 1:252721/1‑90 (MQ=255) gcgcTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac > 1:349621/1‑82 (MQ=255) gcgcTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac < 2:349621/82‑1 (MQ=255) gACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac < 2:343391/76‑1 (MQ=255) gACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcacac > 1:343391/1‑76 (MQ=255) gACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAg < 2:329048/90‑1 (MQ=255) ggCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGAc > 2:233998/1‑90 (MQ=255) ggCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGAc > 1:153896/1‑90 (MQ=255) ggCTGGGACAGCA‑GGGCGGCTTTTACGGGCGGGCGCAACGTTATCATTATGACCTGACCGGGAAACTCACACAGAGGGAGGATGAggggc > 2:227648/1‑88 (MQ=255) tGGTACAGCA‑GGGCGGCTTTGACGGGCGGGCGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTg > 2:316912/1‑90 (MQ=255) gTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTc < 2:143306/90‑1 (MQ=255) | TACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCA‑GGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTC > NZ_CP009273/520637‑520810 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |