Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 884,868 | T→A | 7.0% | L399* (TTA→TAA) | mdfA → | multidrug efflux system protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 884,868 | 0 | T | A | 7.0% | 200.0 / 5.3 | 57 | L399* (TTA→TAA) | mdfA | multidrug efflux system protein |
Reads supporting (aligned to +/- strand): ref base T (26/27); new base A (2/2); total (28/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01 |
TGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAAT > NC_000913/884736‑884966 | tgctgcAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGAt < 3:796837/148‑1 (MQ=255) gcAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGa > 2:203283/1‑148 (MQ=255) gcAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTAAAAAGATaa > 2:652362/1‑138 (MQ=255) gcAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTAAAAAGATaa < 1:652362/138‑1 (MQ=255) aTCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACg > 3:868536/1‑148 (MQ=255) tttACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAg < 4:868536/148‑1 (MQ=255) tttACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAg < 1:203283/148‑1 (MQ=255) aCCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGaa < 3:127156/144‑1 (MQ=255) aCCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGaa > 4:127156/1‑144 (MQ=255) ccGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTGATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATaa > 2:916793/1‑121 (MQ=255) ccGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCGTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATaa < 1:916793/121‑1 (MQ=255) tGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAg < 1:67796/140‑1 (MQ=255) tGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAg > 2:67796/1‑140 (MQ=255) tGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTaaaaa > 3:52448/1‑148 (MQ=255) aaTCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTaaa > 2:160180/1‑102 (MQ=255) aaTCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTaaa < 1:160180/102‑1 (MQ=255) tCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTAAAAAGATa > 4:298111/1‑105 (MQ=255) tCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTAAAAAGATa < 3:298111/105‑1 (MQ=255) cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGAt < 2:1075634/108‑1 (MQ=255) cAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGAt > 1:1075634/1‑108 (MQ=255) aaCATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGAt > 1:968200/1‑98 (MQ=255) aaCATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGAt < 2:968200/98‑1 (MQ=255) aCATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTtgt < 4:52448/148‑1 (MQ=255) aCATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTtgt > 4:607124/1‑148 (MQ=255) cATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTaaaaaaa > 3:880833/1‑131 (MQ=255) cATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTaaaaaaa < 4:880833/131‑1 (MQ=255) cATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTtgtg > 3:469674/1‑148 (MQ=255) ggCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGAt > 1:122886/1‑97 (MQ=255) ggCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGAt < 2:122886/97‑1 (MQ=255) gAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATtctc < 4:537164/104‑1 (MQ=255) gAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATtctc > 3:537164/1‑104 (MQ=255) aaCGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTa < 3:607124/148‑1 (MQ=255) ttAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGa < 1:184111/72‑1 (MQ=255) ttAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGa > 2:184111/1‑72 (MQ=255) tAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTATAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgt < 4:469674/148‑1 (MQ=255) tctcTTCAACCTTGTCAAGGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGc < 3:705625/107‑1 (MQ=255) tctcTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGc > 4:705625/1‑107 (MQ=255) ctctTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATtct < 4:339510/126‑1 (MQ=255) ctctTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATtct > 3:339510/1‑126 (MQ=255) tCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATaaa > 2:135766/1‑59 (MQ=255) tCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATaaa < 1:135766/59‑1 (MQ=255) tCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAAt > 1:216627/1‑147 (MQ=255) tCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAAt < 2:216627/147‑1 (MQ=255) aCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATg > 3:560239/1‑81 (MQ=255) aCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATg < 4:560239/81‑1 (MQ=255) cGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTaa > 4:840847/1‑117 (MQ=255) cGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTaa < 3:840847/117‑1 (MQ=255) gAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATa > 1:193139/1‑42 (MQ=255) gAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATa < 2:193139/42‑1 (MQ=255) ggCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTtg < 4:699485/82‑1 (MQ=255) ggCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTtg > 3:699485/1‑82 (MQ=255) tGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgt < 1:274303/107‑1 (MQ=255) tGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGgt > 2:274303/1‑107 (MQ=255) gATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTa < 1:761573/95‑1 (MQ=255) gATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTa > 2:761573/1‑95 (MQ=255) tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAg < 1:149223/39‑1 (MQ=255) tCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAg > 2:149223/1‑39 (MQ=255) | TGCTGCAAATGCTGATCTTTACCGTTGGTATTGAAATCAGCAAACATGCCTGGCTGAACGGGGGCAACGGACTGTTTAATCTCTTCAACCTTGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAAT > NC_000913/884736‑884966 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |