Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,077,282 G→T 12.7% intergenic (+29/+79) PP_4470 → / ← mgtE Arc domain‑containing transcriptional regulator/magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,077,2820GT12.7% 106.1 / 4.7 39intergenic (+29/+79)PP_4470/mgtEArc domain‑containing transcriptional regulator/magnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (16/18);  new base T (2/3);  total (18/21)
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.64e-01

CTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTG  >  NC_002947/5077143‑5077427
                                                                                                                                           |                                                                                                                                                 
cTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCGTAAGCTGGc                                                                                                                                          <  2:351046/149‑1 (MQ=255)
cTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCGTAAGCTGGc                                                                                                                                          >  1:351046/1‑149 (MQ=255)
   cATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCtttt                                                                                                                                      <  2:669126/150‑1 (MQ=255)
   cATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCtttt                                                                                                                                      >  1:454668/1‑150 (MQ=255)
    aTGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCTTACGCTGGCttttt                                                                                                                                     <  2:59584/150‑1 (MQ=255)
     tGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATTCGGAAATGGCCGCAAACAACTGGCGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCtttttt                                                                                                                                    <  2:120621/150‑1 (MQ=255)
       aGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTtgt                                                                                                                                  <  1:511242/150‑1 (MQ=255)
                              cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCTGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                           >  2:417820/1‑150 (MQ=255)
                              cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGACCACGAAAAAGTCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                           <  2:760151/150‑1 (MQ=255)
                              cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGTGTCTGACTGGCAGCAAGGc                                                                                                           >  2:561931/1‑150 (MQ=255)
                              cGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGc                                                                                                           >  2:505078/1‑150 (MQ=255)
                               gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttcgt                                                                                                                                  >  2:346311/1‑123 (MQ=39)
                               gCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCTTTTTcgt                                                                                                                                  <  1:346311/126‑1 (MQ=39)
                                       tCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCTGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCt                                                                                                  <  1:417820/150‑1 (MQ=255)
                                       tCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCCGCAAGGCCATTTCCCt                                                                                                  >  1:288708/1‑150 (MQ=255)
                                         ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGc                                                                                                >  1:807573/1‑150 (MQ=255)
                                         ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttc                                                                                                                                    <  2:812084/114‑2 (MQ=39)
                                         ccACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACAAAAAAGCCAGCTTACGCTGGCtttttc                                                                                                                                    >  1:812084/1‑113 (MQ=39)
                                           aCCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCgg                                                                                              >  1:212401/1‑150 (MQ=255)
                                                  gcagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGtt                                                                                       <  1:561931/150‑1 (MQ=255)
                                                   cagAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTc                                                                                      <  2:807573/150‑1 (MQ=255)
                                                         gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          >  2:331885/1‑140 (MQ=255)
                                                         gCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGAtt                                                                                          <  1:331885/140‑1 (MQ=255)
                                                           cTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAggg                                                                              >  1:490919/1‑150 (MQ=255)
                                                               gTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTg                                                                                                                     <  2:347687/107‑1 (MQ=255)
                                                               gTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTg                                                                                                                     >  1:347687/1‑107 (MQ=255)
                                                                 cGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGAtt                                                                        <  1:791445/150‑1 (MQ=255)
                                                                           gCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTacaac                                                              <  1:505078/150‑1 (MQ=255)
                                                                             tCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGt                                                                                                                            >  2:353079/1‑86 (MQ=255)
                                                                             tCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGt                                                                                                                            <  1:353079/86‑1 (MQ=255)
                                                                                     tgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCgg                                                                                              <  1:510640/108‑1 (MQ=255)
                                                                                     tgGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCgg                                                                                              >  2:510640/1‑108 (MQ=255)
                                                                                      ggAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGt                                                   <  2:454668/150‑1 (MQ=255)
                                                                                        aaaTGGCCTCAGACGCCTGCTGAGGGCTGCAGCGAACGCGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGTCATTTCCCTGCGGCATTGTTCCCACAGGGAGGATTGGCTTACAACAGGTACAGCGTgg                                                 <  2:212401/150‑1 (MQ=255)
                                                                                                 cAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCaa                                        <  2:490919/150‑1 (MQ=255)
                                                                                                       cctcctGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGa                                  >  1:13473/1‑150 (MQ=255)
                                                                                                                                    aGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCAcagca     >  1:477103/1‑150 (MQ=255)
                                                                                                                                      ccAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGt   <  2:13473/150‑1 (MQ=255)
                                                                                                                                       cAGCGTAAGCTGGCTTGTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTg  >  2:53905/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CTGCATGAGCGCGAGTTGCTGCAGCGTTTCCGCCAACTGTCCCACCGCCAGCAGAACGCCTTGGTCGCCCTCATCGCTCACGATGTGGAAATGGCCGCAGACGCCTCCTGAGGGCTGCAGCGAACACGAAAAAGCCAGCGTAAGCTGGCTTTTTTGTGGGCGTCTGACTGGCAGCAAGGCCATTTCCCTGCGGGATTGTTCCCACAGGGAGGATTGGCTTACAACAGGAACAGCGTGGCCAGCCCAAGGAAGATGAAGAAACCACCGCTGTCGGTCACAGCAGTG  >  NC_002947/5077143‑5077427

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