| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 595,035 | 0 | T | G | 63.0% | 3.8 / 23.1 | 27 | V289G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
| Reads supporting (aligned to +/- strand): ref base T (0/10); new base G (17/0); total (17/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGC > NZ_CP009273/594952‑595123 | tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt < 1:106925/90‑1 (MQ=255) gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGC‑gg > 2:774358/1‑89 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGC‑GGTTCgg > 1:170357/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGC‑GGTCCgg > 2:415252/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCgg > 1:220757/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTC‑GAAcccg > 1:674123/1‑88 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGGTGCGGATCCGGGaaaaagca > 2:1010317/1‑87 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGG‑GGTCCGGCCAAAAAGcc > 2:469214/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGcc < 2:122596/90‑1 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GGTCCGGGCAAAAAGcc > 2:1263266/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGcc > 1:118518/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGC‑GATCCGGCAAAAAAGccc > 1:627000/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGccc < 2:305780/90‑1 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGC‑GATCCCGCAAAAAAGccc > 1:691349/1‑90 (MQ=255) cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCgag < 2:989503/90‑1 (MQ=255) agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAg < 1:648194/90‑1 (MQ=255) gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATcgccgc < 1:1263266/90‑1 (MQ=255) cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGGAggg > 1:271769/1‑89 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGC‑GATCCGGCAAAAAAGCCCGAGAAGAGATCGCCGCCGGGAAGGaa > 2:269368/1‑89 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGTCGTGAAGGAc > 1:85067/1‑90 (MQ=255) aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GGTCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGAGAAAATgg > 1:1261853/1‑90 (MQ=255) aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGGCAAACTgg > 1:158874/1‑90 (MQ=255) gcgaaggcgaagTGGCGGGCGGGGGGGGGGTCCTGG‑GCTCCGGCAAAAAAACCCCAAAAGGGATCCCCGCCGGGGAGGGCAAAATGGgaa > 2:478288/1‑90 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc < 1:1189992/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc < 2:224562/90‑1 (MQ=255) ggcgggcgGGGTGGTGATCCTGG‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag < 2:1215968/90‑1 (MQ=255) ggcgggcgGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag < 2:519192/90‑1 (MQ=255) ggtggtGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 1:788856/90‑1 (MQ=255) gtgATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc > 1:1086023/1‑90 (MQ=255) | TGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGC > NZ_CP009273/594952‑595123 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |