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.G99.0% 373.1 / ‑2.8 96coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (48/47);  total (48/48)
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.08e-01

Only 100 of 119 total aligned reads displayed.

TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948750‑3948843
                                                  |                                             
ttAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                                >  1:3610291/1‑50 (MQ=255)
 tAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                                <  2:5414085/49‑1 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gg                                            >  2:722578/1‑50 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑g                                             <  2:70410/49‑1 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑g                                             <  1:4542484/49‑1 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑g                                             <  2:4730313/49‑1 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggg                                           >  1:3733469/1‑50 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gg                                            <  1:3952539/49‑1 (MQ=255)
    cGGTGGCGCTGCTTTCTCCCGGCACTACCGGTTTACTCAAAATGCC‑‑gggg                                          >  2:3545061/1‑50 (MQ=255)
      gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggg                                        <  2:1229609/50‑1 (MQ=255)
      gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggg                                        >  2:836507/1‑50 (MQ=255)
       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggggg                                       >  2:3004245/1‑49 (MQ=255)
       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggggg                                       >  2:166275/1‑49 (MQ=255)
       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑ggggggg                                       <  1:4255038/50‑2 (MQ=255)
        ggCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                      >  1:1163785/1‑48 (MQ=255)
         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑gggggggt                                      <  1:4268562/49‑3 (MQ=255)
            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                   <  2:3832547/50‑1 (MQ=255)
            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                   <  1:757397/50‑1 (MQ=255)
             tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  >  2:3228304/1‑50 (MQ=255)
              gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 >  2:958266/1‑50 (MQ=255)
              gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 >  2:4827294/1‑50 (MQ=255)
               ctTTCTCCCGGCACTTCCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 <  2:1316866/49‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  <  1:3618003/48‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  <  2:5173418/48‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                  <  1:5562009/48‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 <  1:2520454/49‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 <  2:4792360/49‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                 <  1:2006385/49‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                <  1:4609130/50‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                >  1:3884681/1‑50 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATg                                <  1:5670001/50‑1 (MQ=255)
                tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               >  1:2932735/1‑50 (MQ=255)
                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               >  2:1730736/1‑49 (MQ=255)
                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               >  2:2261755/1‑49 (MQ=255)
                 ttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               >  2:2388499/1‑49 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                              <  2:5280681/49‑1 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               <  2:2817360/48‑1 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                               <  1:3520515/48‑1 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                             >  2:4349782/1‑50 (MQ=255)
                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            >  2:1550180/1‑50 (MQ=255)
                    tcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                           >  1:474549/1‑50 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         >  1:2053928/1‑50 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                           <  1:1054412/48‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            <  1:2383910/47‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAt                            <  1:996857/47‑1 (MQ=255)
                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                         >  2:2751969/1‑49 (MQ=255)
                        ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                       <  1:2612490/50‑1 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTGATGGGATCTCCAg                      >  2:5116728/1‑50 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                      <  1:1758890/50‑1 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  2:1528742/1‑50 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     <  2:1988942/50‑1 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:4956219/1‑50 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:5582095/1‑50 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:2675436/1‑50 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:1859364/1‑50 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     <  2:3113124/50‑1 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTACAGt                     >  2:3327076/1‑50 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCCGt                     >  2:3315990/1‑49 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  2:4696100/1‑49 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:5422432/1‑49 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  2:3326516/1‑49 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:2774668/1‑48 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:3951379/1‑48 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:2469804/50‑1 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   <  1:5553388/50‑1 (MQ=255)
                            cTACCGATTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                     >  1:4937391/1‑48 (MQ=38)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagag                <  1:1649067/50‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                 <  1:3907232/49‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                  <  2:2247946/48‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                  <  1:3254609/48‑1 (MQ=255)
                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgagaga               <  2:222298/50‑1 (MQ=255)
                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  2:3723959/1‑46 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:1973662/1‑45 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  2:2467364/1‑45 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:4945175/1‑44 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                   >  1:1638138/1‑44 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  1:1378866/50‑1 (MQ=255)
                                   tttACTCAAAATGCCGGGGGGGTGGGAATGGGATCTCCAGTGAGAGATc             >  1:3349125/1‑49 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc             <  2:1308749/49‑1 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            <  2:4048971/49‑1 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            <  2:4070531/49‑1 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            >  1:2755971/1‑48 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg            >  1:4010523/1‑48 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGa           <  1:5522362/49‑1 (MQ=255)
                                     tACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa          <  1:351620/49‑1 (MQ=255)
                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa          <  1:809695/49‑1 (MQ=255)
                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa          <  1:331479/49‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg         <  2:2668353/49‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg         <  1:4440934/49‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc        <  1:772152/50‑1 (MQ=255)
                                         cAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg      >  1:1146076/1‑50 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg      <  2:1949046/49‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg       <  2:1625694/48‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt     >  2:5336729/1‑50 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  2:2897127/1‑50 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    <  1:4508522/50‑1 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  2:3522246/1‑49 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  2:2134404/1‑49 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  2:5565980/1‑49 (MQ=255)
                                             aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc    >  1:4413876/1‑48 (MQ=255)
                                                  |                                             
TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCC‑‑GGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGA  >  NC_000913/3948750‑3948843

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