Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,020,755 Δ1 bp 7.0% coding (394/468 nt) PP_4424 → HTH‑type transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,020,7550A. 7.0% 166.9 / 7.1 57coding (394/468 nt)PP_4424HTH‑type transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base . (2/2);  total (29/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 = 9.55e-01

GTTCAAGCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTAATTTGCCTCTCTAAC  >  NC_002947/5020610‑5020897
                                                                                                                                                 |                                                                                                                                              
gtgcAAGCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCaa                                                                                                                                            >  1:1282845/4‑150 (MQ=255)
    aaGCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGt                                                                                                                                        >  1:89398/1‑150 (MQ=255)
      gCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACAta                                                                                                                                              <  2:18582/142‑2 (MQ=255)
      gCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATc                                                                                                                                              >  1:18582/1‑142 (MQ=255)
      gCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGttt                                                                                                                                      >  2:846902/1‑150 (MQ=255)
      gCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGttt                                                                                                                                      <  1:846902/150‑1 (MQ=255)
         ggCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCaa                                                                                                                                   >  1:456550/1‑150 (MQ=255)
         ggCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCaa                                                                                                                                   <  1:455423/150‑1 (MQ=255)
          gCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGAGGCTGACTTTGTACTCCTCGTAACATAAATGGATATGGAGGACTACGAGAACTTCACCCGCGGGGTCTTCTACGAAAACAATGACATCAAGAATTTCaga                                                                                                                                  >  1:605775/1‑148 (MQ=255)
                   aaTACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTc                                                                                                                         <  2:456550/150‑1 (MQ=255)
                       cTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGa                                                                                                                     <  2:89398/150‑1 (MQ=255)
                        tCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGAt                                                                                                                    >  2:1163055/1‑150 (MQ=255)
                          cTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATgg                                                                                                                  >  2:603651/1‑150 (MQ=255)
                             aGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATc                                                                                                               <  2:1282845/150‑1 (MQ=255)
                               gTCATGAGCGGATCATATGTCACGGGTGATGCTGACTTTATACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGCGTTTCAAAACAATGGTCGTGATGGATCgg                                                                                                             <  2:605775/150‑1 (MQ=255)
                                    gAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTaaa                                                                                                        <  1:794231/150‑1 (MQ=255)
                                        ggATCCCATGTCACAGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTc                                                                                                                                     >  1:808837/1‑117 (MQ=255)
                                        ggATCCCATGTCACAGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTc                                                                                                                                     <  2:808837/117‑1 (MQ=255)
                                            ccTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGAt                                                                                                >  2:166389/1‑150 (MQ=255)
                                                    aCGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACtt                                                                                        <  1:166389/150‑1 (MQ=255)
                                                       ggTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCa                                                                                     <  1:603651/150‑1 (MQ=255)
                                                         tGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGCATCTGTTAAAGGCCGGATTCACAATTCCAAt                                                                                   >  2:802458/1‑150 (MQ=255)
                                                               tGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATtgatga                                                                             <  1:752555/150‑1 (MQ=255)
                                                                   tttGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCt                                                                         >  2:31567/1‑150 (MQ=255)
                                                                   tttGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCt                                                                         >  1:328094/1‑150 (MQ=255)
                                                                   tttGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCt                                                                         >  1:1272088/1‑150 (MQ=255)
                                                                          tcctcGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCAcc                                                                  >  1:1160356/1‑150 (MQ=255)
                                                                          tcctcGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAAT‑‑‑‑‑‑‑‑‑AGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAAt                                                                                   <  2:767982/124‑1 (MQ=255)
                                                                          tcctcGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAAT‑‑‑‑‑‑‑‑‑AGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAAt                                                                                   >  1:767982/1‑124 (MQ=255)
                                                                            ctcGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACGTCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACctc                                                                >  1:1247373/1‑150 (MQ=255)
                                                                             tcGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATtgat                                                                               <  2:59382/134‑1 (MQ=255)
                                                                             tcGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATtgat                                                                               >  1:59382/1‑134 (MQ=255)
                                                                                     tCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCtt                                                       >  1:185043/1‑150 (MQ=255)
                                                                                           gATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCaggag                                                 <  2:1272088/150‑1 (MQ=255)
                                                                                             tatGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAggg                                               <  1:1163055/150‑1 (MQ=255)
                                                                                                 gaggaCTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTAc                                           >  2:1117781/1‑150 (MQ=255)
                                                                                                  aggaCTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGAC‑‑CAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACt                                                                                         <  1:666059/101‑1 (MQ=255)
                                                                                                  aggaCTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGAC‑‑CAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACt                                                                                         >  2:666059/1‑101 (MQ=255)
                                                                                                   ggaCTAAGAGAACTTCACCCGAAGGTTCTTCTCCGAGAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGAGCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTc                                         <  1:648899/150‑1 (MQ=255)
                                                                                                     aCTACGAGAACTTCACACGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCt                                       <  2:979092/150‑1 (MQ=255)
                                                                                                          gagaACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAaga                                  >  1:170014/1‑150 (MQ=255)
                                                                                                          gagaACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAaga                                  <  1:1117781/150‑1 (MQ=255)
                                                                                                                ttCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCGCCTCTTGTCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGc                            >  2:1041983/1‑150 (MQ=255)
                                                                                                                    cccGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGaa                        <  1:802458/150‑1 (MQ=255)
                                                                                                                     ccGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAAtt                                                                                  >  1:195315/1‑91 (MQ=255)
                                                                                                                     ccGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAAtt                                                                                  <  2:195315/91‑1 (MQ=255)
                                                                                                                        aaGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATcc                    <  2:328094/150‑1 (MQ=255)
                                                                                                                          ggTTCTTCTACGAAAACAATGACGTCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCtt                  <  2:1247373/150‑1 (MQ=255)
                                                                                                                           gTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTa                 >  2:434327/1‑150 (MQ=255)
                                                                                                                               ttctACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTc                                                                                       <  2:770840/76‑1 (MQ=255)
                                                                                                                               ttctACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTc                                                                                       >  1:770840/1‑76 (MQ=255)
                                                                                                                               ttctACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCagga                                                  >  2:342390/1‑113 (MQ=255)
                                                                                                                               ttctACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCagga                                                  <  1:342390/113‑1 (MQ=255)
                                                                                                                                   aCGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTAATTTGCct         >  2:109598/1‑150 (MQ=255)
                                                                                                                                          cAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGTTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTAATTTGCCTCTCTAAc  <  1:31567/150‑1 (MQ=255)
                                                                                                                                          cAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTAAttt             <  2:82217/139‑1 (MQ=255)
                                                                                                                                          cAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTAAttt             >  1:82217/1‑139 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTTCAAGCAGGCCATTAAAAATACTCCTGAGGTCATGAGCGGATCCTATGTCACGGGTGATGCTGACTTTGTACTCCTCGTAACCTCAAAGGATATGGAGGACTACGAGAACTTCACCCGAAGGTTCTTCTACGAAAACAATGACATCAAGAGTTTCAAAACAATGGTCGTGATGGATCGGGTAAAGGCCGGATTCACACTTCCAATTGATGATGCTTAGCACCTCTTGGCTCTTCAGGAGGGCTACTCCTGAAGAGCTAGCAGAAATCCTTAATTTGCCTCTCTAAC  >  NC_002947/5020610‑5020897

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: ‑