Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,911,777 G→A 14.0% intergenic (+72/+134) mgtA → / ← yjgF magnesium transporter/ketoacid‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,911,7770GA14.0% 156.9 / 13.9 64intergenic (+72/+134)mgtA/yjgFmagnesium transporter/ketoacid‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (4/5);  total (28/36)
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.61e-01

ATAAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCA  >  minE/2911709‑2911846
                                                                    |                                                                     
aTAAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:3290706/71‑1 (MQ=255)
  aaaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                   >  1:100364/1‑71 (MQ=255)
  aaaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                   >  1:4000063/1‑71 (MQ=255)
   aaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaa                                                                    >  1:1414919/1‑69 (MQ=255)
   aaCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAAt                                                                  >  1:2340116/1‑71 (MQ=255)
    aCTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATg                                                                 <  1:4279937/71‑1 (MQ=255)
       gggCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:1312315/64‑1 (MQ=255)
        ggCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCtt                                                             <  1:2150770/71‑1 (MQ=255)
         gCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGt                                                                      <  1:550927/61‑1 (MQ=255)
           aCGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaa                                                                    >  1:532283/1‑61 (MQ=255)
           aCGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATc                                                          <  1:1393271/71‑1 (MQ=255)
            cGATAGCCCCGATTCAGATTGTAGACAAAGTGAGATTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:4294125/59‑1 (MQ=255)
             gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGt                                                                      <  1:4837651/57‑1 (MQ=255)
             gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                         <  1:5329835/70‑1 (MQ=255)
             gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                         >  1:5592202/1‑70 (MQ=255)
             gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATc                                                          <  1:1027106/69‑1 (MQ=255)
             gATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCg                                                        >  1:2046015/1‑71 (MQ=255)
              aTAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCgg                                                       >  1:2608133/1‑71 (MQ=255)
               tAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGc                                                      >  1:5080350/1‑71 (MQ=255)
                aGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                   >  1:1647753/1‑57 (MQ=255)
                 cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                     >  1:5580942‑M2/49‑68 (MQ=255)
                 cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                     >  1:1035075‑M2/49‑68 (MQ=255)
                 cttatcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATccggat                                                     >  1:2375641‑M2/49‑68 (MQ=255)
                    ccGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCcaaca                                                 <  1:3010212/71‑1 (MQ=255)
                    ccGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATACGGCGTAAATGCCTTATCCGGCcaac                                                  >  1:664712/1‑70 (MQ=255)
                     tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                <  1:1471091‑M2/27‑1 (MQ=255)
                     tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                <  1:548195‑M2/27‑1 (MQ=255)
                     tcgtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGGAtaacaa                                                <  1:591027‑M2/27‑1 (MQ=255)
                      gATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACaaa                                               <  1:3773129/71‑1 (MQ=255)
                       gtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCgg                                                       <  1:2777510‑M2/20‑1 (MQ=255)
                       gtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGgataacaata                                              >  1:179359‑M2/43‑69 (MQ=255)
                       gtgccggaatacatcccatatttcggcacgctcaccggatgcGGCATCAGTGCATTATCCGgataacaata                                              <  1:1359697‑M2/29‑3 (MQ=255)
                       aTTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:2748961/48‑1 (MQ=255)
                       aTTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:3686778/48‑1 (MQ=255)
                        ttCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:2204925/47‑1 (MQ=255)
                         tCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:396003/46‑1 (MQ=255)
                         tCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGc                                                                <  1:406676/51‑1 (MQ=255)
                             aTTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAAt                                         >  1:1179855/1‑70 (MQ=255)
                                 tAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:816051/38‑1 (MQ=255)
                                 tAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:2207174/38‑1 (MQ=255)
                                   gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                   >  1:5045404/1‑38 (MQ=255)
                                   gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTaaa                                                                   >  1:4444367/1‑38 (MQ=255)
                                   gACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTg                                   <  1:3545159/70‑1 (MQ=255)
                                    aCAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTa                                                                     <  1:4234065/35‑1 (MQ=255)
                                     cAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAg                                >  1:1643894/1‑71 (MQ=255)
                                         gTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCgg                                                       >  1:1000726/1‑44 (MQ=255)
                                         gTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAAt                            <  1:4649130/71‑1 (MQ=255)
                                           gAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATcc                                                         >  1:1803372/1‑40 (MQ=255)
                                            aGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCCGCCAACAAATTCAATAAATTGCAGGAATCAc                         >  1:2271938/1‑71 (MQ=255)
                                             gTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGcc                                                     >  1:3416890/1‑42 (MQ=255)
                                             gTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCa                          >  1:2281783/1‑69 (MQ=255)
                                                ttGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAg                     >  1:812365/1‑71 (MQ=255)
                                                    tCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCt                 >  1:3438460/1‑71 (MQ=255)
                                                     cATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCt                 >  1:3526411/1‑70 (MQ=255)
                                                     cATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCt                 <  1:4446270/70‑1 (MQ=255)
                                                      aTGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGa               <  1:2773408/71‑1 (MQ=255)
                                                      aTGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGa               <  1:5099439/71‑1 (MQ=255)
                                                           ggatgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGc                                  <  1:2635218‑M2/41‑1 (MQ=255)
                                                           ggatgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGc                                  <  1:3166951‑M2/41‑1 (MQ=255)
                                                              tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCg         >  1:5076195‑M2/4‑69 (MQ=255)
                                                              tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa       <  1:4959696‑M2/68‑1 (MQ=255)
                                                              tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa       <  1:1182871‑M2/68‑1 (MQ=255)
                                                              tgcGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTa       <  1:1901803‑M2/68‑1 (MQ=255)
                                                                cGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGc     <  1:5184338‑M2/70‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAAt                            <  1:4642394/47‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg    <  1:4635045/71‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg    <  1:3069909/71‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg    <  1:4971433/71‑1 (MQ=255)
                                                                 ggcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcg    <  1:4256496/71‑1 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                      >  1:5058795/1‑52 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTa                      >  1:1219835/1‑52 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcgc   <  1:2051417/71‑1 (MQ=255)
                                                                  gcgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAgcgc   <  1:3327339/71‑1 (MQ=255)
                                                                   cgTAAATGCCTTATCCTGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa  <  1:1377405/71‑1 (MQ=255)
                                                                   cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa  >  1:3306600/1‑71 (MQ=255)
                                                                   cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa  <  1:3934798/71‑1 (MQ=255)
                                                                   cgTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCa  <  1:4497950/71‑1 (MQ=255)
                                                                    |                                                                     
ATAAACTGGGCACGATAGCCCCGATTCAGATTGTAGACAAAGTGAGTTTTGTTCATGCCGGATGCGGCGTAAATGCCTTATCCGGCCAACAAATTCAATAAATTGCAGGAATCACGTAGGCCTGATAAGCGTAGCGCA  >  minE/2911709‑2911846

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