Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 595,035 | 0 | T | G | 56.0% | 1.1 / 18.5 | 25 | V289G (GTG→GGG) | cusA | Cu(+)/Ag(+) efflux RND transporter permease subunit CusA |
Reads supporting (aligned to +/- strand): ref base T (0/11); new base G (14/0); total (14/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.24e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAA > NZ_CP009273/594953‑595120 | gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTg > 2:717712/1‑90 (MQ=255) gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg < 2:1724225/81‑1 (MQ=255) aTTGGCCTGAAGATGCGCCGGAGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc < 2:592401/90‑1 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGGAAAAACGcc > 2:1508638/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGcc > 1:320842/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAACAcc > 2:262515/1‑90 (MQ=255) ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGATCCCGGAAAAAAGcc > 2:712076/1‑90 (MQ=255) gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGTCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCCGGCAAAAAAGccc > 2:1255406/1‑90 (MQ=255) agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAGGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg < 1:652048/90‑1 (MQ=255) tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGGTCCGGCAAAAACGCCCGAGGAGGGa > 2:846444/1‑90 (MQ=255) gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCCGCCAAAAAGCCCCAGAAGGGAt > 1:1761831/1‑90 (MQ=255) gggggAATTGGCGAGCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGAAAAAACGCCCGAGAAGTGATcgcc < 1:1581222/89‑1 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATCCGGCAAAAAAACCCGGGAAGGGATggccc > 1:285466/1‑89 (MQ=255) gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc < 1:1056163/90‑1 (MQ=255) ttGCCGAACTAACCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa < 1:1032429/90‑1 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCCGCAAAAAAACCCCAGAAGGGGTCGGCGCCGGGa > 1:579170/1‑90 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCCGCAAAAAAACCCGGGAAGGGAGCGCCGCCGTGGAggga > 1:1085474/1‑88 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGAGGAGGGATCGCCGCCGGGAAGGAc > 2:1457388/1‑90 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc < 2:1085475/90‑1 (MQ=255) aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCCAAAAAGCCCGAGAAGGGATCGCCCCCCTGGAGGACAAACt > 2:1097267/1‑90 (MQ=255) ggcgaaggcgaaGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa < 1:808444/90‑1 (MQ=255) cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc < 1:283860/90‑1 (MQ=255) aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc < 1:623192/90‑1 (MQ=255) ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg < 2:1761832/90‑1 (MQ=255) ggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCCAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAAAGCTg > 2:1311044/1‑90 (MQ=255) gcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa < 1:1311043/90‑1 (MQ=255) ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 1:1444628/90‑1 (MQ=255) ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa < 1:1212601/90‑1 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa > 2:367872/1‑90 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa < 2:579170/90‑1 (MQ=255) | GTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAA > NZ_CP009273/594953‑595120 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |