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.G91.6% 478.5 / 36.9 131coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (4/7);  new base G (42/78);  total (46/85)
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.99e-01

Only 100 of 152 total aligned reads displayed.

TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAG  >  NC_000913/3948750‑3948848
                                                   |                                                 
 tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                                    <  1:2034900/49‑1 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggg                                                >  1:1202510/1‑50 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggg                                                <  2:887503/50‑1 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                                 <  2:328323/49‑1 (MQ=255)
     ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggg                                              <  2:2916858/50‑1 (MQ=255)
      gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggg                                             <  2:1382656/50‑1 (MQ=255)
      gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggg                                             <  1:4987313/50‑1 (MQ=255)
         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggggt                                           <  2:1963824/49‑3 (MQ=255)
            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                        <  1:4111689/50‑1 (MQ=255)
             tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                       >  2:2373627/1‑50 (MQ=255)
             tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAAt                                      >  1:2991994/1‑50 (MQ=255)
              gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                      >  1:5314331/1‑50 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                       <  2:4112627/48‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                       <  1:4272547/48‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                      <  1:1879964/49‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                     >  2:3167632/1‑50 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATgg                                    >  1:4788142/1‑50 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATg                                     <  1:5829637/49‑1 (MQ=255)
                tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                    >  1:3909202/1‑50 (MQ=255)
                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                                   >  1:4833716/1‑50 (MQ=255)
                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                    >  1:5111363/1‑49 (MQ=255)
                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                    >  1:1084003/1‑49 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                    <  2:4661952/48‑1 (MQ=255)
                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                  <  2:2891380/49‑1 (MQ=255)
                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                                  <  1:1526095/49‑1 (MQ=255)
                    tcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                                >  1:5504343/1‑50 (MQ=255)
                     cccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                               >  2:4044937/1‑50 (MQ=255)
                     cccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                               >  2:4984061/1‑50 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                              <  1:3493536/50‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                              >  1:1590278/1‑50 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                              >  2:2416101/1‑50 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                              <  2:423618/50‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                               <  2:3491455/49‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtct                               <  2:4019688/49‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                                 <  1:2079391/47‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGAtctc                              <  1:5666182/49‑1 (MQ=255)
                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                              >  1:5039951/1‑49 (MQ=255)
                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTcc                             <  2:268170/50‑1 (MQ=255)
                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGAtctc                              >  1:1420900/1‑48 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                           <  1:5371777/50‑1 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                           >  1:4309059/1‑50 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  1:2582083/1‑49 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  1:1282111/1‑49 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                         <  2:1462891/50‑1 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                         <  2:618065/50‑1 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  2:2488227/1‑48 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  2:5327389/1‑48 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  1:5588941/1‑48 (MQ=255)
                             tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  1:5677298/1‑47 (MQ=255)
                             tACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                          >  1:2301982/1‑47 (MQ=255)
                              aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                        >  2:2128454/1‑48 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                      <  1:2980275/49‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                       <  1:853796/48‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                       <  1:445142/48‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGACCTCCAGTgaga                      <  1:4155586/49‑1 (MQ=255)
                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga                    <  2:4853860/50‑1 (MQ=255)
                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                        >  1:2311792/1‑46 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                        >  1:934457/1‑45 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt                   <  2:1737025/50‑1 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGAt                   <  2:1823747/50‑1 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                  <  2:4833312/50‑1 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                  <  1:5411702/50‑1 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                  <  2:1565040/50‑1 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                  <  2:3185580/50‑1 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                  <  2:2915563/50‑1 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc                  <  2:2743143/50‑1 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGGGGTAATGGGATCTCCAGTga                        >  2:2359685/1‑44 (MQ=37)
                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg                 <  1:1776766/49‑1 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa                <  1:1402052/49‑1 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa                <  1:2683031/49‑1 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa                <  1:5201580/49‑1 (MQ=255)
                                     tACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa               <  2:3917599/49‑1 (MQ=255)
                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa               <  2:152411/49‑1 (MQ=255)
                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa               <  1:4103801/49‑1 (MQ=255)
                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg              >  2:5280643/1‑50 (MQ=255)
                                       cTCAAAATGCCCAGGGGGTGGTAATAGGATATCCAGTGAGAAATCGAAg              <  2:5207736/49‑1 (MQ=19)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg              <  2:1018941/49‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc             <  2:2205340/50‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc             <  2:54329/50‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc             <  1:3323/50‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc             >  2:3957857/1‑50 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc             <  1:4903743/50‑1 (MQ=255)
                                         cAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt          >  1:1267678/1‑50 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg           <  2:4123719/49‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg           <  1:3718033/49‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg           <  2:2796824/49‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg            <  1:4540034/48‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg            <  2:3662800/48‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg            <  2:5260244/48‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt          <  2:425614/50‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt          <  2:351709/50‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt          <  1:1814533/50‑1 (MQ=255)
                                          aaaaTGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt          <  1:548547/49‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc         <  2:354571/50‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc         >  1:78296/1‑50 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc         <  2:198857/50‑1 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc         >  1:1748068/1‑49 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCg        <  2:4644228/50‑1 (MQ=255)
                                                 ccGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa   <  2:5148015/50‑1 (MQ=255)
                                                  cGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAg  <  2:1148016/50‑1 (MQ=255)
                                                   |                                                 
TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAG  >  NC_000913/3948750‑3948848

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