Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,849,341 | 0 | T | C | 57.1% | 11.8 / 20.7 | 21 | G109G (GGA→GGG) | yidH | YidH family protein |
Reads supporting (aligned to +/- strand): ref base T (8/1); new base C (0/12); total (8/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.42e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATC > NZ_CP009273/3849261‑3849428 | aCGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg > 2:238972/1‑90 (MQ=255) aCGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg > 1:209338/1‑90 (MQ=255) gtgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc > 2:1797/1‑90 (MQ=255) gcgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGa > 1:279479/1‑90 (MQ=255) cgggCTGCAGCCCGGGGTGGGAAACGCGTCTGGTTTGGCGGCTACCCGGCATACAACCCCAGCCCCAAAACGATCACCGCCACGACCATc < 2:75989/88‑1 (MQ=255) ccGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATc < 1:80764/90‑1 (MQ=255) ccGGCTGCAGTCCGGGGTCGGCAATGCGTTTGGTTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc < 2:127302/90‑1 (MQ=255) ccGGCTGAAGTCGGGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACCCCAGCCCAAAACCGATCACCGCCACGACCATc < 1:238972/90‑1 (MQ=255) ccGGCTGAAGCCCGGGGTAGGCAATGCTTCTGTTTTGGCGGCTACCCGGAAAAAAACACCAGCCCAATAACGACCACCGCCACGACCATc < 2:188838/90‑1 (MQ=255) ggCAGTCGGGGGTGGGAAATGCTTTTGGTTTGGCGGCTACCCGGAAAACAACACCAGCCCAAAAACGATCACCGCCACGACCATCAatat < 2:176838/89‑1 (MQ=255) ggTCGGAAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATg < 1:280490/90‑1 (MQ=255) ccGGCAATGCTTCTGTCTTGGCGGCTACCCGGAAAACAACCCCAGCCCAAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt < 2:159787/89‑1 (MQ=255) ggCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta > 1:159787/1‑90 (MQ=255) ggCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta > 2:197025/1‑90 (MQ=255) ggCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta > 2:112323/1‑90 (MQ=255) tGCGTCTGGCTTTGCGGCTACCCGGCATACAACACCAGCCCCAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAg < 2:87076/90‑1 (MQ=255) cgggCTTTGCGGCTACCCGGCATAAAACACCAGCCCCAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGtt < 1:157930/88‑1 (MQ=255) cTGGTTTTGCGGCTACCCGGCAAAAAACCCCAGCCCAAAAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGtt < 2:279479/90‑1 (MQ=255) cTGGCTTTGCGGCTACCCGGCAAACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGtt < 1:283415/90‑1 (MQ=255) ccGGCAAACACCCCCAGCCCCAAAACGATCCCCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTcc < 1:268128/90‑1 (MQ=255) aGTCCCATAACGATCATCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATc > 1:184941/1‑90 (MQ=255) | ACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTTTCAGACGCATC > NZ_CP009273/3849261‑3849428 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |