Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 595,032 | 0 | T | G | 57.5% | 4.0 / 59.1 | 48 | V288G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
Reads supporting (aligned to +/- strand): ref base T (2/18); new base G (27/0); total (29/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.16e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAG > NZ_CP009273/594943‑595121 | tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt < 2:60471/90‑1 (MQ=255) tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt < 1:1762882/90‑1 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 2:748956/90‑1 (MQ=255) gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 2:116324/90‑1 (MQ=255) cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggg > 1:1066969/1‑88 (MQ=255) gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCgg > 1:510430/1‑89 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGGTCC‑gg > 1:1829096/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑gg > 2:321009/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAGCC‑gg > 2:52362/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGGTcccg > 2:959336/1‑88 (MQ=255) tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg > 2:273317/1‑72 (MQ=255) tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg < 1:273317/72‑1 (MQ=255) aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCCGCGGTCC‑GGCgaaa > 2:1090078/1‑90 (MQ=255) aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCaaaa < 1:693569/90‑1 (MQ=255) aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAAc < 1:1851283/90‑1 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCTCCGGTCC‑CGCAAAAAAGcc > 1:1466637/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGATCC‑CGCaaaaaaacc > 1:1182999/1‑86 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGCaaaaaaacc > 2:1883269/1‑86 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGTCTTCGAGCC‑GGGGAAACCGccc > 1:48373/1‑90 (MQ=255) ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCgaga < 1:1301526/90‑1 (MQ=255) ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCgaga < 2:1152422/90‑1 (MQ=255) aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTg < 2:377097/90‑1 (MQ=255) tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCCGCGGTCC‑CGCCAAAAAGCCCCAGAAGGGa > 1:636140/1‑90 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCCCGAAGGGAt > 1:786126/1‑90 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCC‑CGCAAAAAAGCCCCAGAAGGGAt > 2:1808021/1‑90 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGAt > 2:1523472/1‑90 (MQ=255) ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATcgc < 1:321009/90‑1 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCC‑GGCAAAAAAGCCCCAGAAGGGATcgccg > 1:526788/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCC‑GGCAAAAAAGCCCCAGAAGGGATcgccg > 1:409614/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCC‑GGGCAAAAAACCCCAGAAGTGATtgccg > 2:1400355/1‑90 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAAAGCCCGAGAAGGGATcgccg > 2:158719/1‑90 (MQ=255) gggCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGAGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATcgccgc < 1:1400355/90‑1 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGa > 2:448610/1‑90 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATCC‑CGCAAAAAAGCCCGGGAAGGGGTCGccgccgggg > 1:639232/1‑87 (MQ=255) tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa < 2:1448717/90‑1 (MQ=255) gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg < 2:409614/90‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAggga > 1:82016/1‑88 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgggc > 2:52483/1‑88 (MQ=255) aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCC‑GGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGACa > 2:1527645/1‑90 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTTCGAGCC‑GGCAAAAAAGCCCGAGAAGTGAATGCCGCCGCGAAGGGCaaaat > 2:1634673/1‑88 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTTCGCGCCGGGCAAAAAC‑CCCGAGAAGGGATCGCCGCCGCGGAGGAaaaaac > 2:1661544/1‑88 (MQ=255) aCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg > 2:83913/1‑90 (MQ=255) gcgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa < 2:29206/90‑1 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc < 2:401319/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc < 1:272341/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc < 1:1527645/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGc > 2:620089/1‑90 (MQ=255) ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg < 2:510430/90‑1 (MQ=255) ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg < 1:61158/90‑1 (MQ=255) ggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGGAAAGCTg > 2:1642080/1‑90 (MQ=255) cgggcgggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 1:36562/90‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 1:959336/90‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 1:83913/90‑1 (MQ=255) ggggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg < 2:1304375/90‑1 (MQ=255) gggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg > 2:1760927/1‑90 (MQ=255) ggtggtGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 2:742543/90‑1 (MQ=255) tggtgATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg < 1:132625/90‑1 (MQ=255) | TCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAG > NZ_CP009273/594943‑595121 |
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: ‑ |