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.G97.1% 125.4 / ‑2.3 34coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (15/18);  total (15/19)
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.37e-01

AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948752‑3948843
                                                |                                             
aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gg                                            >  1:774569/1‑50 (MQ=255)
aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gg                                            >  1:2376904/1‑50 (MQ=255)
aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑g                                             <  1:1256780/49‑1 (MQ=255)
   ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggg                                         >  1:284146/1‑50 (MQ=255)
    gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggg                                        >  1:608333/1‑50 (MQ=255)
      ggCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                      >  1:970581/1‑48 (MQ=255)
       gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                      <  1:2491345/49‑3 (MQ=255)
           tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  >  1:1459266/1‑50 (MQ=255)
            gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 >  1:994940/1‑50 (MQ=255)
             ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 <  1:1627897/49‑1 (MQ=255)
              tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               >  1:1670556/1‑50 (MQ=255)
              tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGGGGTAATgg                               >  1:1617013/1‑50 (MQ=255)
                 ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                             <  1:3896243/49‑1 (MQ=255)
                 ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                             <  1:3857281/49‑1 (MQ=255)
                    ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            <  1:3080960/47‑1 (MQ=255)
                      ggCACTACCGTTTTACTCAAAATGCCGGGGGGGGTGGTAATGGGAtctc                         >  1:3083534/1‑49 (MQ=255)
                      ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                       <  1:2787387/50‑1 (MQ=255)
                       gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                      <  1:2778319/50‑1 (MQ=255)
                        cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     <  1:2042181/50‑1 (MQ=255)
                         aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:2407821/1‑49 (MQ=255)
                          cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:1815065/50‑1 (MQ=255)
                          cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:3337169/50‑1 (MQ=255)
                          cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:1993758/50‑1 (MQ=255)
                           caCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:873373/2‑47 (MQ=255)
                            aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:1971715/1‑48 (MQ=255)
                               gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt              <  1:2990878/50‑1 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:429375/1‑44 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  1:3503519/50‑1 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  1:3828756/50‑1 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  1:3552507/50‑1 (MQ=255)
                                 tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             >  1:1849763/1‑49 (MQ=255)
                                  ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            <  1:1148398/49‑1 (MQ=255)
                                     cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc        <  1:1835175/50‑1 (MQ=255)
                                      tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg       >  1:782431/1‑50 (MQ=255)
                                       cAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg      >  1:1474121/1‑50 (MQ=255)
                                        aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt     >  1:408131/1‑50 (MQ=255)
                                        aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt     >  1:510957/1‑50 (MQ=255)
                                        aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt     <  1:1362566/50‑1 (MQ=255)
                                         aaaTGCCGGGGGGGTGGCAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  1:872555/1‑50 (MQ=255)
                                         aaaTGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCg   <  1:3122890/50‑1 (MQ=255)
                                           aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCga  <  1:2425210/50‑1 (MQ=255)
                                                |                                             
AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948752‑3948843

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