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.G77.7% 306.1 / 94.2 103coding (1171/1521 nt)yifBputative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base . (11/12);  new base G (35/45);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 123 total aligned reads displayed.

TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGA  >  NC_000913/3948750‑3948845
                                                  |                                               
ttAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGcc                                                  >  2:1769645/1‑50 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                               >  1:417148/1‑50 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                               >  1:4098456/1‑50 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgg                                               >  2:3727751/1‑50 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                                <  2:917064/49‑1 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                                <  2:2949050/49‑1 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                                <  1:2564947/49‑1 (MQ=255)
  aaCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCg                                                <  2:711089/49‑1 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggg                                              >  2:1571561/1‑50 (MQ=255)
   aCGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggg                                              <  1:2309638/50‑1 (MQ=255)
     ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggg                                            >  1:565295/1‑50 (MQ=255)
     ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggg                                            <  2:2668610/50‑1 (MQ=255)
     ggTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggg                                            >  1:1383940/1‑50 (MQ=255)
      gTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggg                                           <  2:3802956/50‑1 (MQ=255)
       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCggggggg                                          >  1:3389414/1‑49 (MQ=255)
       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGggt                                         >  2:2771731/1‑50 (MQ=255)
       tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGggt                                         >  1:2870104/1‑50 (MQ=255)
         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCgggggggt                                         <  2:4149929/49‑3 (MQ=255)
         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGggtg                                        <  2:5153927/49‑1 (MQ=255)
         gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGggtg                                        <  1:4393914/49‑1 (MQ=255)
         gcgcTGCTTTCGCCCGGCGCTACCGTTTTACTCAAAATGCCgggggggt                                         <  1:4336489/49‑3 (MQ=38)
           gctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTa                                     <  1:5206296/50‑1 (MQ=255)
            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                     <  1:4226434/50‑1 (MQ=255)
            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                     <  1:5024242/50‑1 (MQ=255)
            ctgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTa                                     <  1:2217769/50‑1 (MQ=255)
             tgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                    >  1:1108662/1‑50 (MQ=255)
             cgctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTaa                                    >  2:1383507/2‑50 (MQ=255)
              gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                   >  1:3441474/1‑50 (MQ=255)
              gctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATg                                  >  2:5386025/1‑50 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAAt                                   <  1:1702871/49‑1 (MQ=255)
               ctTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATgg                                 >  1:3351553/1‑50 (MQ=255)
                tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                 >  1:241279/1‑50 (MQ=255)
                tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                 >  1:3936186/1‑50 (MQ=255)
                tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                 >  2:3476269/1‑50 (MQ=255)
                tttCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATgg                                 >  2:205561/1‑50 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATggg                                <  1:4567673/49‑1 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                               >  2:2092583/1‑50 (MQ=255)
                  tctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTTGGTAATggg                                <  2:4469892/49‑1 (MQ=255)
                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGa                               <  1:2918146/49‑1 (MQ=255)
                   ctcCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGAt                              <  1:1207155/49‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                           <  2:3336378/50‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtctc                           <  1:156728/50‑1 (MQ=255)
                      ccGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGAtc                             <  1:2323700/48‑1 (MQ=255)
                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTcc                          <  2:1681675/50‑1 (MQ=255)
                       cGGCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTcc                          <  2:4470498/50‑1 (MQ=255)
                        ggCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCa                         <  2:2535263/50‑1 (MQ=255)
                        ggCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAg                        <  2:3647970/50‑1 (MQ=255)
                        ggCACTACCGTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                        >  1:1246702/1‑50 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                        <  2:2698719/50‑1 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                        <  2:4464605/50‑1 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAg                        >  2:718748/1‑50 (MQ=255)
                         gCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGt                       <  2:4468009/50‑1 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                       <  2:3523927/50‑1 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                       <  1:4811584/50‑1 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                       >  2:431683/1‑50 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGt                       >  2:1100497/1‑49 (MQ=255)
                          cACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGt                       >  1:3854965/1‑49 (MQ=255)
                           aCTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTg                      <  1:1293085/50‑1 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGt                       >  2:5011525/1‑48 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                     <  1:1868484/50‑1 (MQ=255)
                            cTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGt                       >  1:2561183/1‑47 (MQ=255)
                             tACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTga                     >  1:3146251/1‑48 (MQ=255)
                              aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                     >  1:3575/1‑48 (MQ=255)
                              aCCGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                     >  1:3850097/1‑48 (MQ=255)
                              aCCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTga                     >  2:1403136/1‑47 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgaga                   <  1:1391470/49‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTgag                    <  2:3903613/48‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTgagaga                 <  1:5207535/50‑1 (MQ=255)
                               ccGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTgaga                   <  1:578281/48‑1 (MQ=255)
                                cGTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                     >  1:290266/1‑46 (MQ=255)
                                 gTTTTACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                     >  1:577946/1‑45 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTga                     >  2:3464275/1‑44 (MQ=255)
                                  ttttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGGga                     >  1:213216/1‑44 (MQ=255)
                                   tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc               >  2:4729087/1‑49 (MQ=255)
                                   tttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATc               >  2:3577088/1‑49 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg              <  2:5019196/49‑1 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg              <  1:1664222/49‑1 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg              <  2:3725546/49‑1 (MQ=255)
                                    ttACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg              >  1:118315/1‑48 (MQ=255)
                                     tACTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCg              >  1:2837238/1‑48 (MQ=255)
                                      aCTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGaa            <  2:5260043/49‑1 (MQ=255)
                                      aCTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg           <  1:4636934/49‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAg           <  2:4908134/49‑1 (MQ=255)
                                       cTCAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGc          >  1:626033/1‑50 (MQ=255)
                                        tCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCgg        >  1:4037430/1‑50 (MQ=255)
                                         cAAAATGCCGGGGGGGTGGTAATGGGATCTCCAGCGAGAGATCGAAGCgg        >  1:2577883/1‑50 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCg         <  2:4968417/48‑1 (MQ=255)
                                          aaaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGt       <  1:5372551/50‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      >  1:1796328/1‑50 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      <  2:519565/50‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      <  2:91245/50‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      <  1:1137525/50‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      <  2:4363473/50‑1 (MQ=255)
                                           aaaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      <  2:3029682/50‑1 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      >  1:2724787/1‑49 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      >  2:316203/1‑49 (MQ=255)
                                            aaTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCg     <  1:4948317/50‑1 (MQ=255)
                                             aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTc      >  1:3916907/1‑48 (MQ=255)
                                             aTGCCGGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCga    <  2:2605053/50‑1 (MQ=255)
                                              tGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCgaga  <  2:2489290/50‑1 (MQ=255)
                                                  |                                               
TTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGCCGGGGGGTGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGA  >  NC_000913/3948750‑3948845

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