Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,031,102 A→T 8.3% M137L (ATG→TTG)  mgtA → ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,031,1020AT 8.3% 137.0 / 3.4 48M137L (ATG→TTG) mgtAATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base T (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAGTGCCC  >  NC_002947/3030957‑3031248
                                                                                                                                                 |                                                                                                                                                  
gCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATgg                                                                                                                                                 >  7:297630/1‑149 (MQ=255)
        tCAAGGCCCTGAAAAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGAAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCt                                                                                                                                          >  5:312319/1‑148 (MQ=255)
          aaGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCATGGTCATGATGATGATGACCTTGGTCAggtcg                                                                                                                                         <  5:138691/147‑5 (MQ=38)
          aaGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCa                                                                                                                                                    <  4:123099/136‑1 (MQ=255)
              cccTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCGGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCagcagc                                                                                                                                   >  5:313210/1‑149 (MQ=255)
                 tGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCctg                                                                                                                                >  3:270540/1‑149 (MQ=255)
                       acCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGACTCAGCAGCCTGCTGCGa                                                                                                                          >  5:307623/1‑149 (MQ=255)
                        cccGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGAt                                                                                                                         >  4:310496/1‑149 (MQ=255)
                         ccGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGAtt                                                                                                                        <  4:133320/149‑1 (MQ=255)
                           gTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCt                                                                                                                      >  6:20352/1‑149 (MQ=255)
                                  tATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                                 >  6:256263/1‑147 (MQ=255)
                                   aTGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTcagca                                                                                                                                     >  8:207234/1‑126 (MQ=255)
                                    tGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAg                                                                                                                 <  3:310496/145‑1 (MQ=255)
                                           ctgACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCGGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTc                                                                                                      <  6:313210/149‑1 (MQ=255)
                                                       gggATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGc                                                                                          <  3:209669/149‑1 (MQ=255)
                                                              gCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAg                                                                                               <  6:177384/137‑1 (MQ=255)
                                                              gCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAg                                                                                               >  5:177384/1‑137 (MQ=255)
                                                              gCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGAcgcg                                                                                   <  6:21686/149‑1 (MQ=255)
                                                               cTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACgcgc                                                                                  >  4:21648/1‑149 (MQ=255)
                                                               cTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACgcgc                                                                                  <  3:21648/149‑1 (MQ=255)
                                                                    tCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACGAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCaaa                                                                             <  4:57612/149‑1 (MQ=255)
                                                                    tCACAGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACAGCTCGAACAAGGCTGCCGACGCGCTCaaa                                                                             <  4:53314/149‑1 (MQ=255)
                                                                      aCCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAg                                                                            >  8:64855/1‑148 (MQ=255)
                                                                       ccGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGc                                                                           >  1:286551/1‑148 (MQ=255)
                                                                         gACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAAt                                                                        <  5:20352/149‑1 (MQ=255)
                                                                            tactGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGt                                                                     >  2:82732/1‑149 (MQ=255)
                                                                              ctGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTgc                                                                   <  2:286551/149‑1 (MQ=255)
                                                                                 gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAgg                                                                                              >  5:159415/1‑119 (MQ=255)
                                                                                 gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAgg                                                                                              <  6:159415/119‑1 (MQ=255)
                                                                                 gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCa                                                                <  5:104844/149‑1 (MQ=255)
                                                                                 gCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCa                                                                <  4:270540/149‑1 (MQ=255)
                                                                                     ccGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAAt                                                                        >  5:174230/1‑137 (MQ=255)
                                                                                     ccGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAAt                                                                        <  6:174230/137‑1 (MQ=255)
                                                                                     ccGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTa                                                             <  8:297630/148‑1 (MQ=255)
                                                                                      cGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATTGTCAGCCTCAGCATCCTGCTGCGATTCTGGCAGGAGTACCGCTCGTACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCTCACTACg                                                           <  6:312319/149‑1 (MQ=255)
                                                                                              gcAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGgt                                                   <  3:100110/149‑1 (MQ=255)
                                                                                                  gCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTgc                                                                   >  1:297847/1‑129 (MQ=255)
                                                                                                  gCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTgc                                                                   <  2:297847/129‑1 (MQ=255)
                                                                                                    gAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTgcg                                             >  8:162851/1‑149 (MQ=255)
                                                                                                        cAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGcgccgc                                         <  7:64855/149‑1 (MQ=255)
                                                                                                         aGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCcgcg                                        >  7:89045/1‑149 (MQ=255)
                                                                                                           acgacgCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAg                                      >  5:42455/1‑149 (MQ=255)
                                                                                                                gCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAg                                      >  7:23811/1‑144 (MQ=255)
                                                                                                                       tGACCATGGTCATGATGATGATGACCTTGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATCGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCaa                          >  4:221809/1‑149 (MQ=255)
                                                                                                                        gACCATGGTCATGATGATGATGACCTTGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAg                                   <  6:203732/139‑1 (MQ=39)
                                                                                                                        gACCATGGTCATGATGATGATGACCTTGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAg                                   >  5:203732/1‑139 (MQ=39)
                                                                                                                                        atcatGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGcgcg          >  5:68011/1‑148 (MQ=255)
                                                                                                                                              aCCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAGTGc    <  7:162851/148‑1 (MQ=255)
                                                                                                                                               ccATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAGTGccc  <  6:68011/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAACTGCTCAAGGCCCTGAACAACCCGTTCATCTATGTACTGCTGACCCTGGCCGGGATCAGCTTCGTCACCGACTACTGGCTGCCGGTCAGCGCAGGCGAAGCAGACGACGCCGACCTGACCAAGGTCATCATCATCATGACCATGGTCAGCCTCAGCAGCCTGCTGCGATTCTGGCAGGAGTACCGCTCGAACAAGGCTGCCGACGCGCTCAAAGCAATGGTGCGCACTACGGCCACGGTGTTGCGCCGCGAGCAGATCGGGCAAGCCCCGCGCCTGCGCGAAGTGCCC  >  NC_002947/3030957‑3031248

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