Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,558,329 | 0 | C | A | 55.6% | 9.2 / 39.3 | 27 | R3L (CGT→CTT) | mdtM | multidrug efflux MFS transporter MdtM |
Reads supporting (aligned to +/- strand): ref base C (2/10); new base A (15/0); total (17/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.82e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAA > NZ_CP009273/4558241‑4558418 | 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 < 1:508845/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGa < 1:400979/90‑1 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGa < 2:269118/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) cGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGaa < 1:322841/90‑1 (MQ=255) cGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGaa < 2:79397/90‑1 (MQ=255) | GATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGA‑GAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAA > NZ_CP009273/4558241‑4558418 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |