Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,031,395 G→C 8.9% S33T (AGC→ACC)  regA → photosynthetic apparatus regulatory protein RegA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,031,3950GC 8.9% 166.1 / 7.3 56S33T (AGC→ACC) regAphotosynthetic apparatus regulatory protein RegA
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base C (2/3);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CACCCTGACCGAACTGCGCCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGAACTGGACC  >  NC_002947/1031250‑1031544
                                                                                                                                                 |                                                                                                                                                     
cACCCTGACCGAACTGCGCCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGt                                                                                                                                                   <  1:1175484/150‑1 (MQ=255)
cACCCTGACCGAACTGCGCCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGCTCATGGCCCGGGCCATGACCCgg                                                                                                                                                   <  2:1005427/150‑2 (MQ=255)
 aCCCTGACCGAACTGCGCCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTc                                                                                                                                                  >  2:269066/1‑150 (MQ=255)
            aCTGCGCCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCgt                                                                                                                                       <  1:182438/150‑1 (MQ=255)
             cctcgcttGCCATAGGGCAGGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCgtg                                                                                                                                      >  2:1339348/8‑150 (MQ=255)
              ctcgcTTGCCATAGGGCAGGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCgtgt                                                                                                                                     >  1:823186/3‑150 (MQ=255)
               tcgcttGCCATAGGGCAGGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTGGCGGTTTCCGtgtg                                                                                                                                    >  2:718507/6‑150 (MQ=255)
                cgcCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGa                                                                                                                                   <  2:212077/150‑1 (MQ=255)
                cgcCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGa                                                                                                                                   <  1:269066/150‑1 (MQ=255)
                       ccTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGCTCATGGCCCGGGCCATGACCCgggt                                                                                                                                                 >  2:973718/1‑126 (MQ=255)
                       ccTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGCTCATGGCCCGGGCCATGACCCgggt                                                                                                                                                 <  1:973718/129‑4 (MQ=255)
                          agggCAGGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAg                                                                                                                         <  1:1339348/148‑1 (MQ=255)
                               aaGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTg                                                                                                                    >  1:1131020/1‑150 (MQ=255)
                                          tgaatgaTGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGc                                                                                                         <  1:1280030/150‑1 (MQ=255)
                                                 tgaGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTc                                                                                                  <  1:917752/150‑1 (MQ=255)
                                                  gaGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCt                                                                                                 >  2:328333/1‑150 (MQ=255)
                                                   aGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTg                                                                                                <  2:1131020/150‑1 (MQ=255)
                                                   aGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGACGAAGAAGGCTTGATTCTg                                                                                                <  1:1404711/150‑1 (MQ=255)
                                                    gCGAAGAAAACCAGGTCGAAAGCGAATAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTgg                                                                                               <  2:293734/150‑1 (MQ=255)
                                                      gaagaaAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTAATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGcc                                                                                             >  1:1153165/1‑150 (MQ=255)
                                                                   tCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGCGACGCCGCGGTTTCCGTTTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGc                                                                                >  1:1032945/1‑150 (MQ=255)
                                                                     gAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTTCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCg                                                                              <  2:823707/150‑1 (MQ=255)
                                                                     gAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCg                                                                              <  2:823186/150‑1 (MQ=255)
                                                                           gAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTAc                                                                        >  1:151770/1‑150 (MQ=255)
                                                                              gagCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGcc                                                                     >  2:689119/1‑150 (MQ=255)
                                                                                    ccGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTg                                                               <  1:328333/150‑1 (MQ=255)
                                                                                              tgctgGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTgg                                                                                               <  2:7648/108‑1 (MQ=255)
                                                                                              tgctgGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTgg                                                                                               >  1:7648/1‑108 (MQ=255)
                                                                                                 tggtggATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTg                                                                                                >  1:40459/1‑104 (MQ=255)
                                                                                                 tggtggATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTg                                                                                                <  2:40459/104‑1 (MQ=255)
                                                                                                       aTGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                <  1:896573/146‑1 (MQ=255)
                                                                                                       aTGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                >  2:896573/1‑146 (MQ=255)
                                                                                                        tGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                >  1:709169/1‑145 (MQ=255)
                                                                                                        tGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                <  2:709169/145‑1 (MQ=255)
                                                                                                            gATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTcc                                       >  2:842515/1‑150 (MQ=255)
                                                                                                                ccACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGcc                                   >  2:1287330/1‑150 (MQ=255)
                                                                                                                  aCCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTg                                 <  2:151770/150‑1 (MQ=255)
                                                                                                                       ggcccgggccatgacccgggccatgaGCCGTCGCGGATTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGtgct                            >  2:534106/3‑150 (MQ=255)
                                                                                                                       cACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGtgct                            <  2:1234188/150‑1 (MQ=255)
                                                                                                                         cccGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCGACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGtgctgc                          <  2:1032945/150‑1 (MQ=255)
                                                                                                                            ggctCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTg                                 >  1:105183/4‑140 (MQ=255)
                                                                                                                            ggctCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTg                                 <  2:105183/137‑1 (MQ=255)
                                                                                                                            ggctCATGGCCCGGGCCATGACCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCca                                                                                           >  2:1167255/4‑82 (MQ=255)
                                                                                                                            ggctCATGGCCCGGGCCATGACCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCca                                                                                           <  1:1167255/79‑1 (MQ=255)
                                                                                                                              gTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGccc                     >  1:640628/1‑150 (MQ=255)
                                                                                                                                 aTGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAg                  >  2:735668/1‑150 (MQ=255)
                                                                                                                                     cccGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAgctgc              >  2:995192/1‑150 (MQ=255)
                                                                                                                                      ccGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAgctgct             >  2:904591/1‑150 (MQ=255)
                                                                                                                                       cGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                >  2:1271761/1‑114 (MQ=255)
                                                                                                                                       cGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                >  2:1271632/1‑114 (MQ=255)
                                                                                                                                       cGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                <  1:1271761/114‑1 (MQ=255)
                                                                                                                                       cGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAgatggat                                                <  1:1271632/114‑1 (MQ=255)
                                                                                                                                          gCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGaa         <  1:1287330/150‑1 (MQ=255)
                                                                                                                                           ccATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGAAc        <  2:1153165/150‑1 (MQ=255)
                                                                                                                                            cATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGAACt       >  1:1251904/1‑150 (MQ=255)
                                                                                                                                             aTGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGAACTg      >  1:1354205/1‑150 (MQ=255)
                                                                                                                                                 gCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGAACTGGAcc  <  1:995192/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CACCCTGACCGAACTGCGCCTGCCCTATGGCAAGCGAGGAGATGAATGATGAGCGAAGAAAACCAGGTCGAAAGCGAAGAGCTGCCGCACCTGCTGCTGGTGGATGACGATGCCACCTTCACCCGGGTCATGGCCCGGGCCATGAGCCGTCGCGGTTTCCGTGTGAGCACTGCCAGTTCTGCCGAAGAAGGCTTGATTCTGGCCCAGCAGGACCTGCCGGACTACGCCACGCTGGACCTGAAGATGGATGGCGACTCCGGCCTGGTGCTGCTGCCCAAGCTGCTGGAACTGGACC  >  NC_002947/1031250‑1031544

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