Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 595,032 | 0 | T | G | 54.2% | ‑4.5 / 16.9 | 24 | V288G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
Reads supporting (aligned to +/- strand): ref base T (1/10); new base G (13/0); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.61e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGA > NZ_CP009273/594946‑595119 | gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggg > 2:1989512/1‑89 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 2:1846567/90‑1 (MQ=255) cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGGCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg < 1:549741/90‑1 (MQ=255) cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtg > 2:523483/1‑90 (MQ=255) aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc > 1:805248/1‑90 (MQ=255) agggTCCAGATTGGCCCGGAGATTCGCCGGGGCATTGCCTAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATcc < 1:1372383/88‑1 (MQ=255) tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGTCGGGCGGGGGGGTGATCCTGCGATCCGGCa > 2:1852015/1‑90 (MQ=255) aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc < 1:1283004/90‑1 (MQ=255) aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGAAAAAAc < 1:523483/90‑1 (MQ=255) ggCCCGGAGATTCGCCGGGGCATTGCCTAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGGCCGGCCAAAAAGcc > 2:1372383/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCCAAAAAGcc > 1:1596426/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCGGCGCGCCGGCaaaaaaacc > 1:1195574/1‑86 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGAt > 2:818855/1‑90 (MQ=255) cgcgGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:1195493/90‑1 (MQ=255) cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 2:1542802/90‑1 (MQ=255) cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc < 1:1603763/90‑1 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCCGCGATCCGGCAAAAAAGCCCGAGAAGGGAGtgccg > 1:1986265/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCCGCGGACCCGCCAAAAAACCCCAGAAGGGGTcgccg > 1:1180633/1‑90 (MQ=255) gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGGGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc < 1:1742810/90‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc < 1:647528/90‑1 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACaaaat > 2:1122177/1‑88 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc < 1:219664/69‑1 (MQ=255) cgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATcgccgc > 2:219664/1‑69 (MQ=255) gaaggcgaagTGGCGGGGGGGGGGGGGGTCCTTCGATCCGGCAAAAAAGGCCCCGAAGGGGACGCCCCCGTGGAGGAAAAAATTGaaaaa > 1:51479/1‑88 (MQ=255) tGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCCAGAAGGGATCGCCGCCCTGAAGGACAAACTGGAAAAGCTGAAAagta > 2:2002316/1‑90 (MQ=255) gcgggcggGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGGTCGCCGCCGGGAAGGAAAAAATGGAAAAGCTGAAAagtagt > 1:1382495/1‑90 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 1:1538485/90‑1 (MQ=255) gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGAAGTCTGCCgg > 2:1496366/1‑90 (MQ=255) ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 1:818855/90‑1 (MQ=255) ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 2:500020/90‑1 (MQ=255) | GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGA > NZ_CP009273/594946‑595119 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |