Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,346,840 (T)8→7 coding (45/105 nt) ymiB → putative protein YmiB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,346,8330T.93.3% 49.6 / ‑1.7 15coding (38/105 nt)ymiBputative protein YmiB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (10/4);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATC  >  NC_000913/1346753‑1346906
                                                                                |                                                                         
ctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGt                                                                           <  2:234999/81‑1 (MQ=255)
    ggTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttt                                                                       <  1:235000/81‑1 (MQ=255)
         aCAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATTACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttata                                                                  <  1:235001/81‑4 (MQ=255)
                       ggTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAAt                                                   >  1:235015/1‑81 (MQ=255)
                       ggTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAAt                                                   <  1:235002/81‑1 (MQ=255)
                            aataaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGt                                              >  1:235021/1‑81 (MQ=255)
                                   ggAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTCTTATATATTTATTTGTAATCCAGTTTTGGa                                        >  2:235027/1‑80 (MQ=255)
                                    gAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGaa                                       >  2:235018/1‑80 (MQ=255)
                                         acaTGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGc                                 <  1:235003/81‑1 (MQ=255)
                                           aTGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATACTTATTTGTAATCCAGTTTTGGAAAAACGc                                 >  2:235014/1‑79 (MQ=255)
                                                 tACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGtttt                          >  1:235022/1‑80 (MQ=255)
                                                       aaaTTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGa                   <  2:235004/81‑1 (MQ=255)
                                                         aTTAACGGCGGCTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGaaa                 <  1:235005/81‑1 (MQ=255)
                                                                 cggcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGtt          >  1:235020/1‑80 (MQ=255)
                                                                  ggcTAAACGCATTGTTTTTTTTTATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAGGTCAGTTa         >  2:235026/1‑80 (MQ=255)
                                                                  ggcTAAACGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTa         >  2:235019/1‑80 (MQ=255)
                                                                      aaaCGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTaaaaaa    >  2:235023/1‑81 (MQ=255)
                                                                        aCGCATTGTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATc  >  2:235017/1‑81 (MQ=255)
                                                                                |                                                                         
CTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATC  >  NC_000913/1346753‑1346906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: