| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 172,755 | 0 | A | G | 56.4% | 17.7 / 39.2 | 39 | T388A (ACG→GCG) | clcA | H(+)/Cl(‑) exchange transporter ClcA |
| Reads supporting (aligned to +/- strand): ref base A (4/13); new base G (22/0); total (26/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.93e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACC > NZ_CP009273/172676‑172841 | ggCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TAt < 1:20153/90‑1 (MQ=255) cTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGcc < 2:460070/90‑1 (MQ=255) aCTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATg < 1:429164/90‑1 (MQ=255) tgtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATggg < 2:206044/90‑1 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTGC‑TTTTGCCGGGATgggg > 1:819152/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTC‑TTTTGCCGGGATgggg > 2:894145/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTC‑TATTTCCGGAATgggg > 2:112686/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTGCCGGGATtggg > 2:923484/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGC‑TATTTCCGGGATgggg > 1:285832/1‑90 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTC‑TATTTCCGGGATgggggga > 2:264782/1‑88 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGTTTGC‑TATTTCCGGGATGGGGGGAtt > 2:886607/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTC‑TTTTTCCCGGATgggggggtt > 2:1036392/1‑86 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTC‑TATTGCCGGGATGGGGGGAtt > 1:283682/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTTCCGGAATGGGGGCAtt > 2:971112/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTGCCGGGATGGGGGCAtt > 1:913628/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTTTTTGCCGG‑ATGGGGGGAtt > 2:494842/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGC‑TATTTCCGGGATGGGGGGAtt > 1:470198/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGC‑TATTTCCGGGATGGGGGCAtt > 1:451119/1‑90 (MQ=255) aCCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGGATTACTg > 2:544148/1‑90 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TTTTGCCGGGAAGGGGGCATTACTgg > 2:742777/1‑90 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTTCCGGGATGGGGGGATTACTgg > 2:1102240/1‑90 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTGCCGGGATGGGGGGATTACTgg > 1:126912/1‑90 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGC‑TATTGCCGGGATGGGGGCATTACttg > 2:19877/1‑88 (MQ=255) gCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTGCCGGGATGGGGGCATTACTggcg > 2:36016/1‑90 (MQ=255) cGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGC‑TATTTCCGGAATGGGGGCATTACTGGCGGCATc > 1:92681/1‑90 (MQ=255) ttGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGGATGGGGGGATTACTGGCGGCATCTATTcgcg > 1:471840/1‑90 (MQ=255) tGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGT‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTcgcgc < 1:804413/90‑1 (MQ=255) cGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGcc < 1:682605/90‑1 (MQ=255) ttGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGt < 1:1036392/90‑1 (MQ=255) cTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACg < 1:313498/90‑1 (MQ=255) gTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACggg > 2:286642/1‑90 (MQ=255) tCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGAt < 2:819152/90‑1 (MQ=255) atatCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct < 2:66468/90‑1 (MQ=255) tatCACCTTGAGGCGGGGGCGTTTGC‑TATTGCCGGGATGGGGGCATTACTTGCGGCATCTATTTGCGCGCCGTTAACGGGGGGCAttctg > 1:899229/1‑90 (MQ=255) gAGGCGGGGGCGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGgag > 2:41010/1‑90 (MQ=255) ggCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGAt < 2:126912/90‑1 (MQ=255) ggggACGTTTTC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTTGCGCGCCGTTAACGGGGGTCATTCTGGTTTTGGGGGTGAc > 1:284387/1‑90 (MQ=255) gggACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAcc < 1:203011/90‑1 (MQ=255) gggACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAcc < 1:923484/90‑1 (MQ=255) | GGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGC‑TATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACC > NZ_CP009273/172676‑172841 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |