Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 595,035 | 0 | T | G | 71.9% | ‑0.1 / 23.2 | 33 | V289G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
Reads supporting (aligned to +/- strand): ref base T (0/9); new base G (23/0); total (24/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.57e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGC > NZ_CP009273/594948‑595123 | gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa < 2:922030/90‑1 (MQ=255) gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCgg > 1:977175/1‑89 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCCGGGATCCgg > 1:615053/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTcccg > 2:1122375/1‑88 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCgg > 2:52624/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCgg > 2:466803/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGGCCgg > 1:314581/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTCCgg > 2:1018105/1‑90 (MQ=255) ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATcccg > 2:1065208/1‑88 (MQ=255) tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCa > 2:13983/1‑90 (MQ=255) ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCaa > 2:305488/1‑90 (MQ=255) aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaa < 2:166867/90‑1 (MQ=255) aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc < 2:237569/90‑1 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGcc > 2:1011688/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTCCGGTCCGGCAAAAACAcc > 2:803162/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAACGcc > 2:230880/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGccc > 1:402653/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGccc > 2:409820/1‑90 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGAt > 2:54475/1‑90 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGAt > 2:45268/1‑90 (MQ=255) gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc < 2:1166220/90‑1 (MQ=255) cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa < 2:977175/90‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCAGGGAAAAAGCCCCAGCAGTGGGCGCCGCCGTGGAggtc > 1:809118/1‑88 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGGAAAAAAGCCCGAGAAGGGATCGCCGCCGCGAAGGAc > 2:987210/1‑90 (MQ=255) gAACTAAACAGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCCCCGGGGAgggc > 2:862589/1‑88 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGAGATCCCGCGATCCGGCCAAAAACCACGAGAAGTGATCGCCGCCGCGAAGGACAAACt > 2:134833/1‑90 (MQ=255) ggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCCAGAAGGGATCGCCGCCGGGAAGGGAAAACTGGAAAAGCTg > 1:445691/1‑90 (MQ=255) aagTGGCGGGCGGGGGGGGGATCCTTGGATCCCGCGAAAACGCCCTTGGATTGGGCGCCGGCGGGGAGGGCGAAATGGGAAAGCTGaaaa > 1:431934/1‑90 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 2:402653/90‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 2:255095/90‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 2:724076/90‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc < 1:999249/90‑1 (MQ=255) gggcgggGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCCAGAAGGGATCGCCGCCGGGAAGGGCAAAATGGAAAAGCTGAAAAGTAGTCt > 2:1202431/1‑90 (MQ=255) ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 2:809118/90‑1 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa < 2:384188/90‑1 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa < 2:144069/90‑1 (MQ=255) gtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc < 1:276395/90‑1 (MQ=255) | GCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGC > NZ_CP009273/594948‑595123 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |