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.8% 170.3 / ‑2.5 45coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (17/27);  total (17/28)
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.52e-01

AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948752‑3948841
                                                |                                          
aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                             <  1:1904836/48‑1 (MQ=255)
 aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                           <  1:581062/48‑1 (MQ=255)
   ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggg                                        >  1:1590516/1‑49 (MQ=255)
             ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                               <  1:2010947/49‑1 (MQ=255)
             ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                               <  1:1425079/49‑1 (MQ=255)
             ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                               <  1:1163484/49‑1 (MQ=255)
              tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                              >  1:1486110/1‑49 (MQ=255)
              tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                              >  1:963622/1‑49 (MQ=255)
                tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                             <  1:672726/48‑1 (MQ=255)
                 ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                            <  1:969194/48‑1 (MQ=255)
                    ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                        <  1:594416/49‑1 (MQ=255)
                    ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                         <  1:486959/48‑1 (MQ=255)
                    ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                          <  1:1920003/47‑1 (MQ=255)
                    ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                          <  1:2320709/47‑1 (MQ=255)
                     cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       >  1:477604/1‑49 (MQ=255)
                     cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       <  1:848/49‑1 (MQ=255)
                      ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                       >  1:2260958/1‑48 (MQ=255)
                         aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   <  1:1857499/49‑1 (MQ=255)
                         aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   <  1:67645/49‑1 (MQ=255)
                         aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:198982/1‑49 (MQ=255)
                           tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                   >  1:368094/1‑47 (MQ=255)
                             ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga               <  1:1069782/49‑1 (MQ=255)
                             ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga               <  1:200983/49‑1 (MQ=255)
                             ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 <  1:664939/47‑1 (MQ=255)
                             ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:1600770/1‑47 (MQ=255)
                               gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:626411/1‑45 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:1191918/1‑44 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                 >  1:752362/1‑44 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt            <  1:2021396/49‑1 (MQ=255)
                                ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt            <  1:1581428/49‑1 (MQ=255)
                                  ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           >  1:1843333/1‑48 (MQ=255)
                                  ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc           <  1:1736259/48‑1 (MQ=255)
                                   tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:2000968/1‑48 (MQ=255)
                                    aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg          >  1:832651/1‑47 (MQ=255)
                                     cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa        <  1:389564/48‑1 (MQ=255)
                                     cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa        <  1:2288019/48‑1 (MQ=255)
                                     cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       <  1:1158594/49‑1 (MQ=255)
                                     cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg       <  1:328867/49‑1 (MQ=255)
                                      tCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      >  1:1721063/1‑49 (MQ=255)
                                        aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      <  1:47809/47‑1 (MQ=255)
                                        aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc      <  1:88014/47‑1 (MQ=255)
                                        aaaaTGCCGGGGGGGTGGAAATGGGATCTCCAGTGAGAGATCGAAGCgg    >  1:1381612/1‑49 (MQ=255)
                                         aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:1269592/49‑1 (MQ=255)
                                         aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:2142473/49‑1 (MQ=255)
                                         aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   <  1:1707984/49‑1 (MQ=255)
                                         aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt   >  1:101261/1‑49 (MQ=255)
                                         aaaTGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  <  1:750315/49‑1 (MQ=255)
                                           aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc  >  1:2283511/1‑48 (MQ=255)
                                                |                                          
AACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTC  >  NC_000913/3948752‑3948841

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