Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 997,447 | 0 | T | G | 48.9% | 43.3 / 12.2 | 45 | V62G (GTA→GGA) | elfG | fimbrial protein |
Reads supporting (aligned to +/- strand): ref base T (6/17); new base G (22/0); total (28/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.15e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCGAT > NZ_CP009273/997360‑997536 | ggCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAg > 1:1477182/1‑90 (MQ=255) ggCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAg > 1:1549761/1‑90 (MQ=255) tactacTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACaa > 1:1327450/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAttt > 1:1244922/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAttt > 1:199571/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAAttt > 1:1268312/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAAAAAAatt > 2:1602704/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGGTTAAGACAAAttt > 2:1649795/1‑90 (MQ=255) tCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTa > 1:658146/1‑90 (MQ=255) ggggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAAAAAAatttttta > 2:1754163/1‑86 (MQ=255) gggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAg < 2:297656/90‑1 (MQ=255) gAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTg > 1:1601065/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAGAAAAAATTTATTAggggg > 2:1212558/1‑89 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa > 1:1235820/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa > 1:1712771/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAATTTATCAGTGGa > 2:41777/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTTTTAGTGGa > 1:1902758/1‑90 (MQ=255) tataTTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaa < 2:1542612/90‑1 (MQ=255) tataTTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAATTTTTTCAGGGGAAcgaa > 2:128949/1‑90 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:199571/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:61618/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:658146/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:126288/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:1074018/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 1:1212558/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 1:128949/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGTGTTTAAAAAAAAATTTTTTAGTGGAAcgaaacggaa > 2:1286560/1‑87 (MQ=255) gCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatc < 2:380209/85‑1 (MQ=255) gCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatc > 1:380209/1‑85 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGat > 1:1682806/1‑90 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGTGGTAAAAAAAAAATTTTTTAGTGGAAAGAATCGAATGat > 2:1224417/1‑90 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGGGTTAAAGAAAAAATTTTTTAGGGGAAAGAATAGAATGat > 2:1138483/1‑90 (MQ=255) gTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGat > 1:979577/1‑90 (MQ=255) aaTTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGattat < 1:1649795/90‑1 (MQ=255) gTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAAAAAAATTTTTCAGGGGGAAGAATCGAATGATTTTTTTGTTTcc > 2:1229411/1‑90 (MQ=255) tAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTTCtg > 2:621542/1‑90 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTg < 1:9762/90‑1 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTCCCTGTGATTg < 1:914164/90‑1 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa < 1:404783/90‑1 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa < 1:833803/90‑1 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa < 1:325359/90‑1 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa < 2:1549761/90‑1 (MQ=255) gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGt > 2:744046/1‑90 (MQ=255) aGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTTCTGTGATTGCGATAAAGACAATGTCAGAAGTgg > 2:1309592/1‑90 (MQ=255) ggTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCg > 1:780632/1‑90 (MQ=255) tAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCGAt < 2:716828/90‑1 (MQ=255) | GGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGTGGCCGAT > NZ_CP009273/997360‑997536 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |