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.3% 205.6 / 1.1 54coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (25/27);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCG  >  NC_000913/3948751‑3948842
                                                  |                                            
tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                               <  1:5694742/49‑1 (MQ=255)
tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                               <  1:4090091/49‑1 (MQ=255)
tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                               <  1:235139/49‑1 (MQ=255)
tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                             >  1:2907/1‑50 (MQ=255)
 aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                             <  1:3022770/49‑1 (MQ=255)
  aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                            <  1:668832/49‑1 (MQ=255)
   cGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggg                                          >  1:667596/1‑50 (MQ=255)
   cGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggg                                          >  1:2840080/1‑50 (MQ=255)
    ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggg                                         >  1:2534825/1‑50 (MQ=255)
    ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggg                                         <  1:5129604/50‑1 (MQ=255)
     gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggg                                       <  1:4598384/50‑1 (MQ=255)
      tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggggg                                      >  1:5347653/1‑49 (MQ=255)
        gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggggt                                     <  1:4195483/49‑3 (MQ=255)
        gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggggt                                     <  1:4947818/49‑3 (MQ=255)
           ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                  <  1:450033/50‑1 (MQ=255)
            tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                 >  1:2126198/1‑50 (MQ=255)
             gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                >  1:5296142/1‑50 (MQ=255)
             gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                >  1:5027379/1‑50 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                 <  1:5808015/48‑1 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                 <  1:2889518/48‑1 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                 <  1:2786045/48‑1 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                 <  1:5948775/48‑1 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                <  1:4423391/49‑1 (MQ=255)
              ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                               <  1:897549/50‑1 (MQ=255)
                ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGAGGTAATgg                              >  1:3524092/1‑49 (MQ=255)
                 tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                              <  1:1241299/48‑1 (MQ=255)
                 tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                            >  1:2364943/1‑50 (MQ=255)
                  ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                            <  1:3300967/49‑1 (MQ=255)
                    cccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                         >  1:3038187/1‑50 (MQ=255)
                     ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                        >  1:4700671/1‑50 (MQ=255)
                     ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                           <  1:2202448/47‑1 (MQ=255)
                      cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                        >  1:3389287/1‑49 (MQ=255)
                      cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                        >  1:3604375/1‑49 (MQ=255)
                      cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTcc                       <  1:1276160/50‑1 (MQ=255)
                       ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                        >  1:4441500/1‑48 (MQ=255)
                        gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                     <  1:1034597/50‑1 (MQ=255)
                         cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                    >  1:3038528/1‑50 (MQ=255)
                          aCTACCGTTTTACTCAAAATGTCGGGGGGGTGGTAATGGGATCTCCAGTg                   <  1:4859833/50‑1 (MQ=255)
                          aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                    >  1:6029405/1‑49 (MQ=255)
                          aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                    >  1:4951904/1‑49 (MQ=255)
                              ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                <  1:4768883/49‑1 (MQ=255)
                              ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                <  1:5866843/49‑1 (MQ=255)
                               cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga              <  1:2773553/50‑1 (MQ=255)
                               cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                  >  1:527004/1‑46 (MQ=255)
                                gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                  >  1:4809880/1‑45 (MQ=255)
                                gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt             <  1:1691097/50‑1 (MQ=255)
                                 ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc            <  1:5161572/50‑1 (MQ=255)
                                 ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc            <  1:1144747/50‑1 (MQ=255)
                                 ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc            <  1:1332170/50‑1 (MQ=255)
                                 ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc            <  1:2613065/50‑1 (MQ=255)
                                  tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc            >  1:3235385/1‑49 (MQ=255)
                                   ttACTCAAAATGCCGGGGGTGTGGTAATGGGATCTCCAGTGAGAGATCg           <  1:3856094/49‑1 (MQ=255)
                                    tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg           >  1:3489088/1‑48 (MQ=255)
                                     aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa         <  1:976375/49‑1 (MQ=255)
                                       tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg      >  1:5685364/1‑50 (MQ=255)
                                       tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg      >  1:4780861/1‑50 (MQ=255)
                                         aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg     <  1:74195/49‑1 (MQ=255)
                                         aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg      <  1:3502361/48‑1 (MQ=255)
                                         aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg      <  1:4830068/48‑1 (MQ=255)
                                         aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt    <  1:2307312/50‑1 (MQ=255)
                                          aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   <  1:5593337/50‑1 (MQ=255)
                                          aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   >  1:3164803/1‑50 (MQ=255)
                                          aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   >  1:2675096/1‑50 (MQ=255)
                                           aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   >  1:5505234/1‑49 (MQ=255)
                                           aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   >  1:1389626/1‑49 (MQ=255)
                                           aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   >  1:512847/1‑49 (MQ=255)
                                           aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCg  <  1:2828596/50‑1 (MQ=255)
                                            aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc   >  1:555368/1‑48 (MQ=255)
                                                  |                                            
TAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCG  >  NC_000913/3948751‑3948842

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