Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,156,907 C→A 6.8% intergenic (‑73/+128) PP_3659 ← / ← PP_3660 transcriptional regulator, AraC family/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,156,9070CA 6.8% 183.5 / 2.2 59intergenic (‑73/+128)PP_3659/PP_3660transcriptional regulator, AraC family/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base A (2/2);  total (28/31)
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.19e-01

ACGTAGTTCTGCTCCTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACAGTCAACCCACCTTTTGC  >  NC_002947/4156760‑4157050
                                                                                                                                                   |                                                                                                                                               
aCGTAGTTCTGCTCCTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATAccc                                                                                                                                               >  2:1017079/1‑150 (MQ=255)
  gTAGTTCTGCTCCTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTa                                                                                                                                             >  1:346464/1‑150 (MQ=255)
           cTCCTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGaa                                                                                                                                    <  2:1568535/150‑1 (MQ=255)
             ccTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGc                                                                                                                                  >  1:751814/1‑150 (MQ=255)
              cTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGcc                                                                                                                                 >  2:337795/1‑150 (MQ=255)
              cTTGAAGTGAGTCCCGCACGAGACTGATCGCCCTGCTCAATCTGTCCGCACAGAACTCCATCTGGGGCGCGGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGACGATACATATAGCCAATACCCTATGGTCAGAAgag                                                                                                                                 >  1:823437/1‑148 (MQ=255)
                tGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGccac                                                                                                                               <  1:419723/150‑1 (MQ=255)
                     tGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAgg                                                                                                                          <  2:346464/150‑1 (MQ=255)
                     tGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAgg                                                                                                                          <  2:297366/150‑1 (MQ=255)
                            cGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGAt                                                                                                                   >  1:191060/1‑150 (MQ=255)
                             gCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGAtt                                                                                                                  >  1:1167866/1‑150 (MQ=255)
                              cACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCTGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTg                                                                                                                 <  1:484209/150‑1 (MQ=255)
                              cACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTg                                                                                                                 <  1:1017079/150‑1 (MQ=255)
                              cACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTg                                                                                                                 <  2:633542/150‑1 (MQ=255)
                                ggAGGCTGATGGTCGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGt                                                                                                               <  2:138463/149‑1 (MQ=255)
                                      tGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCaa                                                                                                         <  1:1839556/150‑1 (MQ=255)
                                       gATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAAc                                                                                                        <  2:1167866/150‑1 (MQ=255)
                                                 gtgCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTggggg                                                                                              >  2:578178/1‑150 (MQ=255)
                                                           gggCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGa                                                                                     <  2:234886/149‑1 (MQ=255)
                                                           gggCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAAGACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGa                                                                                     >  1:234886/1‑149 (MQ=255)
                                                             gCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATAc                                                                                  >  2:882479/1‑150 (MQ=255)
                                                                 aaaaCTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCa                                                                              >  1:1122238/1‑150 (MQ=255)
                                                                  aaaCTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCAc                                                                             <  1:882479/150‑1 (MQ=255)
                                                                      tCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACaaacagc                                                                                                                        >  1:1088154/1‑99 (MQ=255)
                                                                      tCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGc                                                                                                                        <  2:1088154/103‑1 (MQ=255)
                                                                            tGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGTAGTCCAGATACAGCACCGTCTTCaaa                                                                   <  1:278441/150‑1 (MQ=255)
                                                                                        gCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGt                                                       <  1:649382/150‑1 (MQ=255)
                                                                                          tAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATAACCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGAt                                                                                    >  1:1553955/1‑119 (MQ=255)
                                                                                          tAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATAACCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGAt                                                                                    <  2:1553955/119‑1 (MQ=255)
                                                                                              tCGATACCGAAGGCTGCCTTGCAAATCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCtgt                                                 <  2:1899860/150‑1 (MQ=255)
                                                                                              tCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCtgt                                                 >  1:188384/1‑150 (MQ=255)
                                                                                              tCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATAACCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAatat                                                                <  1:1033613/135‑1 (MQ=255)
                                                                                              tCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATAACCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAatat                                                                >  2:1033613/1‑135 (MQ=255)
                                                                                               cGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCtgtg                                                <  2:191060/150‑1 (MQ=255)
                                                                                                 aTACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGc                                              >  2:1396110/1‑150 (MQ=255)
                                                                                                     cGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGa                                                                                     >  2:1134224/1‑107 (MQ=255)
                                                                                                     cGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGa                                                                                     <  1:1134224/107‑1 (MQ=255)
                                                                                                     cGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTACAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGc                                          >  1:1582623/1‑150 (MQ=255)
                                                                                                      gAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCa                                         >  1:249002/1‑150 (MQ=255)
                                                                                                         ggCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGt                                                                          >  2:1353237/1‑114 (MQ=255)
                                                                                                         ggCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGt                                                                          <  1:1353237/114‑1 (MQ=255)
                                                                                                               cTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGa                                >  2:1222688/1‑150 (MQ=255)
                                                                                                                  gCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAaga                             >  1:849394/1‑150 (MQ=255)
                                                                                                                  gCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGCCAGTTCAGATTGAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAaga                             <  2:751814/150‑1 (MQ=255)
                                                                                                                    aaaGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTACCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCCGGCACGCCCGCGAagaaa                           >  2:156680/1‑148 (MQ=255)
                                                                                                                         aTACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCt                                                                                                                       <  1:1037921/53‑1 (MQ=255)
                                                                                                                         aTACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCt                                                                                                                       >  2:1037921/1‑53 (MQ=255)
                                                                                                                            cacGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACa                   <  2:249002/150‑1 (MQ=255)
                                                                                                                              cGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACAGt                 <  1:1222688/150‑1 (MQ=255)
                                                                                                                                 tgaaccATTTAGCCGATCCCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACAGTCaa              <  2:527688/146‑1 (MQ=255)
                                                                                                                                 cGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACAGTCaa              <  1:1275495/150‑1 (MQ=255)
                                                                                                                                       atatAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCtgt                                                            <  1:27244/98‑1 (MQ=255)
                                                                                                                                       atatAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCtgt                                                            >  2:27244/1‑98 (MQ=255)
                                                                                                                                       atatAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGaa                               >  1:346080/1‑127 (MQ=255)
                                                                                                                                       atatAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGaa                               <  2:346080/127‑1 (MQ=255)
                                                                                                                                           aGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTAc                                                     <  1:33928/101‑1 (MQ=255)
                                                                                                                                           aGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTAc                                                     >  2:33928/1‑101 (MQ=255)
                                                                                                                                             ccAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACAGTCAACCCACCTTTTGc  >  1:888127/1‑150 (MQ=255)
                                                                                                                                               aaTACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACa                   >  2:409284/1‑131 (MQ=255)
                                                                                                                                               aaTACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACa                   <  1:409284/131‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACGTAGTTCTGCTCCTTGAAGTGCGTCCCGCACGAGGCTGATCGACGTGGTGCCTTTGTGGGCAGAAAACTCCATCTGGGGCTCCGTTGCTAACTCGATACCGAAGGCTGCCTTGCAAAGCATACACGCCGATCCATATAGCCAATACCCTATGTTCAGAAGCCACCAAGGGCTGGATTGGTCTCCAACTTCCTGGGGGCAGTCCAGATACAGCACCGTCTTCAAATATCTGTGCTGTACCTGTGGCGGGCACGCCCGCGAAGAGACCAGTACAGTCAACCCACCTTTTGC  >  NC_002947/4156760‑4157050

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: