Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 595,032 | 0 | T | G | 63.0% | ‑4.1 / 25.6 | 27 | V288G (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 = 1.37e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGA > NZ_CP009273/594943‑595119 | tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt < 1:446735/90‑1 (MQ=255) tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggg > 2:675222/1‑90 (MQ=255) gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggc > 2:1220580/1‑87 (MQ=255) 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) | TCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGA > NZ_CP009273/594943‑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: ‑ |