Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,948,805 (G)6→7 coding (1165/1521 nt) yifB ← putative magnesium chelatase YifB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,948,7991.G96.2% 86.6 / ‑2.1 26coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (14/11);  total (15/11)
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.21e-01

TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948751‑3948841
                                                 |                                           
tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                              <  1:653552/49‑1 (MQ=255)
 aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                           >  1:982265/1‑50 (MQ=255)
 aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                           >  1:2697116/1‑50 (MQ=255)
  aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggg                                          >  1:739359/1‑50 (MQ=255)
  aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTCTTACTCAAAATGCCggg                                          >  1:2373675/1‑50 (MQ=255)
     gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggg                                      <  1:2260979/50‑1 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                               <  1:1207617/49‑1 (MQ=255)
               tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                             >  1:2650606/1‑50 (MQ=255)
                 tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                           >  1:2059574/1‑50 (MQ=255)
                  ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                           <  1:848250/49‑1 (MQ=255)
                  ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                          >  1:1354153/1‑50 (MQ=255)
                       ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                     <  1:1047801/50‑1 (MQ=255)
                       ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                     <  1:2607597/50‑1 (MQ=255)
                       ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                     <  1:1739717/50‑1 (MQ=255)
                        gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                    >  1:1750374/1‑50 (MQ=255)
                         cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:844682/1‑50 (MQ=255)
                          aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                  <  1:2080650/50‑1 (MQ=255)
                          aCTACCGTTTTACTCAAAATGCCGGGGTGGTGGTAATGGGATCTCCAGt                   >  1:1036339/1‑49 (MQ=255)
                           cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 <  1:2217548/50‑1 (MQ=255)
                            tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:1221917/1‑47 (MQ=255)
                                  tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           >  1:1944222/1‑49 (MQ=255)
                                   ttACTCAAAATGCCGGGGGGGTGGTAATGGGCTCTCCAGTGAGAGATCg          >  1:1924143/1‑49 (MQ=255)
                                     aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       >  1:2508176/1‑50 (MQ=255)
                                     aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       >  1:996901/1‑50 (MQ=255)
                                      cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      <  1:2697044/50‑1 (MQ=255)
                                      cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      <  1:308026/50‑1 (MQ=255)
                                      cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      <  1:1390370/50‑1 (MQ=255)
                                       tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg     >  1:1486990/1‑50 (MQ=255)
                                       tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg     >  1:1093292/1‑50 (MQ=255)
                                         aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:37927/50‑1 (MQ=255)
                                           aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  >  1:1979712/1‑49 (MQ=255)
                                            aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  >  1:201192/1‑48 (MQ=255)
                                                 |                                           
TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948751‑3948841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: