Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,321 | 0 | T | G | 55.2% | 25.2 / 31.2 | 29 | T6P (ACC→CCC) | mdtM | multidrug efflux MFS transporter MdtM |
Reads supporting (aligned to +/- strand): ref base T (4/9); new base G (16/0); total (20/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCA > NZ_CP009273/4558236‑4558410 | ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 2:370416/1‑90 (MQ=255) ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa > 1:641277/1‑90 (MQ=255) ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa > 1:85494/1‑89 (MQ=255) gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACg < 1:397017/90‑1 (MQ=255) cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATaaa > 2:642639/1‑89 (MQ=255) aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca < 2:331281/90‑1 (MQ=255) aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca < 1:654773/90‑1 (MQ=255) cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAacaaacaa < 1:164441/90‑1 (MQ=255) gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAAGGGCATAAAAAAAAACCCCCTTTcc > 2:578592/1‑90 (MQ=255) gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACAGGGGATAACAAAAAACTCCTTGTcc > 2:49661/1‑90 (MQ=255) gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAAAAACTCCTTTTcc > 2:164441/1‑90 (MQ=255) aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGAAAAAAAAAAAAACCTTGTCCgg > 1:79397/1‑90 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACCCCTTTTCCGGG‑GAAcc > 2:340243/1‑90 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCAAAAAAAAAAAAACCCTTTCCGGA‑Gaaac > 1:127559/1‑88 (MQ=255) aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAAAACCCTTTCCGGAGGAACCCCaa > 1:112485/1‑90 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAACAAAAAAATCCTTGTCCGGA‑GAAACCCAACgg > 2:16155/1‑90 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGAAAAAAAAAAAAACCTTTTCCGGG‑GAAACCCCACgg > 1:157375/1‑90 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAAAACCCTGTCCGGG‑GAACCCCAACgg > 1:269118/1‑90 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGa < 2:269118/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGa < 1:508845/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGa < 1:400979/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAAACCCTTTCCGGG‑GAAACCCCACGGGGGAAAAGa > 1:209454/1‑90 (MQ=255) aaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATa < 2:403568/90‑1 (MQ=255) gCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACCCCCTGTCCGGA‑GGACCCCAACGGGGAAAAAGAAAATCAATAACGt > 2:341663/1‑90 (MQ=255) cGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTc > 2:60612/1‑90 (MQ=255) cGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCGGA‑GAAACCCAACGGGGGAAAAGAAAATCAATAcccgc > 2:421145/1‑87 (MQ=255) cGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCTTTTCCGGG‑GGACCCCAACGGGGGAAAAGAAATTCAATACCCTc > 1:332288/1‑90 (MQ=255) gTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGGAAAAGAAAATCAATACCGTCg > 2:144960/1‑90 (MQ=255) gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc < 2:239032/90‑1 (MQ=255) tAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCa < 2:85494/90‑1 (MQ=255) | GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCA > NZ_CP009273/4558236‑4558410 |
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: ‑ |