| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,081,789 | 0 | T | G | 57.9% | 10.4 / 19.6 | 19 | H476P (CAC→CCC) | frvR | putative frv operon regulatory protein |
| Reads supporting (aligned to +/- strand): ref base T (2/6); new base G (11/0); total (13/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCTTTCACCGCGTACATTCGAATAATGAAAACTTCCGGGGGC > NZ_CP009273/4081725‑4081877 | aaTCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAAAATGGGGTGGTGCGGCTTATTGTGCA‑CAAAt > 1:343442/1‑90 (MQ=255) aGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGc < 1:432026/90‑1 (MQ=255) tttCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATTGGGTTGTGCGGCTAATTTCGCA‑CAAATTTTCCggg > 1:73728/1‑90 (MQ=255) tCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTa < 2:224411/90‑1 (MQ=255) ttcgccttcgcGGGCGATGATGCGTTGTGCTTCATCGGCGGGAATATTGGGTTGTGGGACTAATTTCGCC‑CAAATTTTCCGGGGGATATg > 2:277626/1‑90 (MQ=255) tcgccttcgcGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATTGGGGTGTGCGGCTAATTTGGCC‑CAAAATTGCCGGGGGatatct > 2:432026/1‑88 (MQ=255) ttcgcGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCa < 1:425492/90‑1 (MQ=255) tcgcGGGCGATGATGCGTTGTGCTTCATCGGCGGGAATATGGGGTTGTGCGGCTAATTGCGCA‑CCAATTTTCCCGGGGATATTTTTCCAg > 2:96196/1‑90 (MQ=255) gggCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCttt < 1:370522/90‑1 (MQ=255) atgCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCTTTCACCGCGt < 1:301533/90‑1 (MQ=255) atgCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCTTTCACCGCGt < 2:378521/90‑1 (MQ=255) ttGTGCTTCATCGGCGGTAATATGGGGTTGTGGGGCTAATTGCGCA‑CAAAATTTTCGGGGAATATTTTTCCCGCTTTTACCCCGTACAtt > 2:176557/1‑90 (MQ=255) gtgCTTCATCGGCGGTAATATGGGGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGGAATATGTTGCCAGCTTTCACCGCGTACATTCg > 1:97778/1‑90 (MQ=255) gtgCTTCATCGGCGGTAATATGGGGGTTTGGGGCTAATTGCGCC‑CAAAATTTCCGGGGAAAATTTTTCCAGCTTTTCCCGCGGACAtttg > 2:28998/1‑88 (MQ=255) gtgCTTCATCGGCGGTAATATGGGGGTGTGCGGCTAATTTCGCC‑CAAATTTTCCGGGGAATATGTTTCCCGCTTTTCCCCCGGACATTCg > 2:52168/1‑90 (MQ=255) gtgCTTCATCGGCGGGAATATTGGGGTGTGGGGCTAATTGCGCACCAAATTTCCCGGG‑AATATGTTGCCAGCTTTCACCGCGGACATTCg > 2:247300/1‑90 (MQ=255) ggcggTAATATGGTGTTGTGCGACTAATTTCGCA‑CAAATTTGCCGGGGAAAATGTTGCCCGCTTTCACCCCGTACATTCGAATAATGaaa > 1:371114/1‑90 (MQ=255) ggcggTAATATGGGGTTGTGCGACTAATTGCGCA‑CAAATTTTCCGGGGAATATGTTGCCCGCTTTCACCCCGTACATTCGAATAATGaaa > 1:208444/1‑90 (MQ=255) ggTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCTTTCACCGCGTACATTCGAATAATGAAAACTTCCggggg > 1:32400/1‑90 (MQ=255) gtgtTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCTTTCACCGCGTACATTCGAATAATGAAAACTTCCGGGGGc > 2:402483/1‑90 (MQ=255) | AATCAGGTTTTCACCTTCGCCTTCGCGGGCGATGATGCGTTGTGCTTCATCGGCGGTAATATGGTGTTGTGCGACTAATTGCGCA‑CAAATTTGCCGGGTAATATGTTGCCAGCTTTCACCGCGTACATTCGAATAATGAAAACTTCCGGGGGC > NZ_CP009273/4081725‑4081877 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |