| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 595,035 | 0 | T | G | 56.1% | 15.8 / 49.1 | 41 | V289G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
| Reads supporting (aligned to +/- strand): ref base T (2/16); new base G (23/0); total (25/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.48e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAA > NZ_CP009273/594946‑595120 | gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 2:1261116/90‑1 (MQ=255)gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 1:690095/90‑1 (MQ=255) gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa < 1:617942/90‑1 (MQ=255) aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc > 2:1274760/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 2:808472/75‑1 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCgg > 1:174606/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCCGGGGGGGGACCCTGCGGTcccg > 1:1328218/1‑88 (MQ=255) aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaa < 2:772727/90‑1 (MQ=255) gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaac < 2:735627/90‑2 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCCGCaaaaaaacc > 1:858616/1‑86 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGATCCGGCAAAAAAGcc > 2:1188099/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGAGCCGGCcaaaaaggc > 2:626967/1‑86 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGAACCGGCAAAAAAGcc > 2:824601/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAACGcc > 1:1128221/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAACGcc > 2:996935/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGccc > 1:417275/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGccc > 1:330859/1‑90 (MQ=255) cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg < 1:1334150/90‑1 (MQ=255) ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgaga < 2:863875/90‑1 (MQ=255) tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGAGAAGGGa > 1:1085859/1‑90 (MQ=255) tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGGGa > 1:1122996/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCCGGAAGGGATcgccg > 2:290536/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATcgccc > 1:995952/1‑89 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCGCGGAGGGGTCGccgccgggg > 1:771851/1‑87 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa < 2:1085859/90‑1 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAACCCGAGAAGGGATCGCCGCCGGGa > 1:1265964/1‑90 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCCAAAAAGCCCGAGAAGTGGTTGGCGCCGTGGAggga > 2:627225/1‑88 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc < 1:577479/90‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc < 2:1174221/90‑1 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATCCCGCAAAAAAGCCCCAGAAGGGGTCGCCCCCGGGGAGGGCaaaaa > 2:323680/1‑88 (MQ=255) gcgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa < 2:32894/90‑1 (MQ=255) gcgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCCAGAAGTGATCGCCGCCGTGAAGGGCAAAATGGaaa > 1:1120834/1‑90 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc < 2:1328218/90‑1 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc < 1:1295320/90‑1 (MQ=255) gaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGGTCGCCGCCGGGGAGGACAAAATGGgaaag > 1:786269/1‑88 (MQ=255) aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc < 2:192946/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc < 2:881220/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGc > 1:32894/1‑90 (MQ=255) gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa < 2:79309/90‑1 (MQ=255) ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag < 1:503161/90‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCTGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 2:1170533/90‑1 (MQ=255) cgggcgggGGGGGGGTCCTTCGCTCCGGGAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGAAGGGCAAACGGGGAAAGCTGTAAAGGAGTc > 2:327161/1‑90 (MQ=255) gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGTACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg > 2:1157722/1‑90 (MQ=255) gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg > 1:716754/1‑90 (MQ=255) ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 1:290536/90‑1 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa < 1:287138/90‑1 (MQ=255) | GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAA > NZ_CP009273/594946‑595120 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |