| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 822,963 | 0 | T | G | 52.8% | 38.7 / 80.9 | 53 | D492A (GAC→GCC) | ybhF | ATP‑binding cassette domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base T (2/23); new base G (28/0); total (30/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-09 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA > NZ_CP009273/822886‑823045 | gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGCCGGAAg > 1:1080170/1‑90 (MQ=255)gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAg > 2:1950631/1‑90 (MQ=255)gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAg > 2:301441/1‑90 (MQ=255)gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAAAGCCGGGAg > 2:1323849/1‑90 (MQ=255)gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAg > 2:582399/1‑90 (MQ=255) ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAgg > 1:238124/1‑89 (MQ=255) ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGTCgg > 2:906337/1‑90 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc < 2:1519878/90‑1 (MQ=255) cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGgtcgtc > 2:1479503/1‑90 (MQ=255) aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtcg < 1:162832/90‑1 (MQ=255) cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaaa < 1:1564284/90‑1 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGAAGGCGGGTCGGCGAGaaaa > 2:2018066/1‑89 (MQ=255) gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGAGAAAAAg > 1:969501/1‑90 (MQ=255) aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGa < 1:2018066/90‑1 (MQ=255) acagtttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa < 2:1253544/86‑1 (MQ=255) tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTc > 2:314172/1‑90 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 2:21407/90‑1 (MQ=255) ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt < 2:1080170/90‑1 (MQ=255) aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:809007/84‑1 (MQ=255) aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc > 2:809007/1‑84 (MQ=255) ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCAt > 2:1386289/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc < 1:1911404/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc < 2:1841428/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGGCCGGGTTATGc > 1:472924/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGGCGGGTCGTCGGGAAACAGAAAGTTCGGGTTATGc > 1:2080423/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGc > 1:726369/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGc > 1:116227/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGc > 2:605477/1‑90 (MQ=255) cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAATCGGGTCGTCGGGGAACAAAAAGGCCGGGTCATGCATc > 2:1080898/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCCCGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:1564725/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGGAGGGTTGTCGGGAAAAAAAAAGTCCGGGTCATGGATCAg > 1:2085710/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAAAATGTCCGGGTCATGCATCAg > 2:700160/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1180639/90‑1 (MQ=255) gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAgg > 2:499726/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGGGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGGGAAAAgg > 1:491408/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAAAGTGCGGGTCATGCATCAGCGCACAgg > 1:1155566/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAAAAAAAAGTTCGGGTCATGCATCAGCGAACaag > 1:251720/1‑88 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTCATGGATCAGCGAACAgg > 1:969669/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGGCGTCGGGAAAAAAAAAGTCCGGGTCATGCATCAGCGAACAgg > 1:929311/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGTCAAAGCCGGAAGTCGGGTCGTCGAGAAACAAAAAATCCGGGTCATTCATCAGCGAACaag > 2:1254491/1‑88 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAgg > 2:200351/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc < 1:1689422/90‑1 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAGCGAACAGgcc > 1:236833/1‑90 (MQ=255) gCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcca < 2:1588221/90‑1 (MQ=255) aTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcca < 1:1254132/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:969501/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:2085710/90‑1 (MQ=255) tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa < 2:929311/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 1:1386289/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:116227/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 1:499726/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:1437286/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:969669/90‑1 (MQ=255) tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGt > 2:1074024/1‑37 (MQ=255) tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGt < 1:1074024/37‑1 (MQ=255) | GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA > NZ_CP009273/822886‑823045 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |