Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,974,658 G→A 8.2% R183H (CGC→CAC)  PP_2601 → IclR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,974,6580GA 8.2% 131.6 / 3.2 49R183H (CGC→CAC) PP_2601IclR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/16);  new base A (2/2);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CCTCCGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCC  >  NC_002947/2974520‑2974807
                                                                                                                                          |                                                                                                                                                     
ccTCCGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGc                                                                                                                                            >  2:250113/1‑150 (MQ=255)
 cTCCGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCg                                                                                                                                           >  2:490860/1‑150 (MQ=255)
    tcagCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTgg                                                                                                                                        >  2:129349/4‑150 (MQ=255)
    cGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTCCACCGCGCCCGCGAGCGTgg                                                                                                                                        >  1:837394/1‑150 (MQ=255)
     gggCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGt                                                                                                                                       <  2:193822/150‑1 (MQ=255)
     gggCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCACGCTCGCGGGCGCGGt                                                                                                                                       >  2:921336/1‑150 (MQ=255)
         aCGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCACGCTCGCGGGCGCGGtgcag                                                                                                                                   <  1:921336/150‑5 (MQ=255)
           gCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGcc                                                                                                                                 >  1:249645/1‑150 (MQ=255)
            ccTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGCTGTTGCGCACAGACAATTCTTACTCCGACCAGTTCGGCCCGAACGCGCCACGGTCCACCAACGAGCTGCTCGCCTACCTGCACTGCGCCCGCGAACGTGGCTATCCCt                                                                                                                                >  2:539236/1‑150 (MQ=255)
                  tggcCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCg                                                                                                                          >  2:620318/1‑150 (MQ=255)
                   ggcCAGCCTGGATGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGa                                                                                                                         >  2:19636/1‑150 (MQ=255)
                    gcCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAg                                                                                                                        <  1:626204/150‑1 (MQ=255)
                    gcCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAg                                                                                                                        <  1:539236/150‑1 (MQ=255)
                           tGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCt                                                                                                                 >  1:141823/1‑150 (MQ=255)
                                        ggCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCAc                                                                                                    <  1:620318/150‑1 (MQ=255)
                                            cTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGAGTGGTTGTGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCa                                                                                                >  1:622923/1‑150 (MQ=255)
                                               cAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGcc                                                                                             >  1:160341/1‑150 (MQ=255)
                                                 gATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCACGCTCGCGGGCGCGGtgc                                                                                                                                     >  1:369387/1‑106 (MQ=255)
                                                 gATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCACGCTCGCGGGCGCGGtgc                                                                                                                                     <  2:369387/108‑3 (MQ=255)
                                                  aTGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTg                                                                                          >  2:95459/1‑150 (MQ=255)
                                                                 ggTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTgcgc                                                                           >  2:821315/1‑150 (MQ=255)
                                                                  gTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGcgcc                                                                          <  2:160341/150‑1 (MQ=255)
                                                                  gTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACATGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGcgcc                                                                          >  2:263818/1‑150 (MQ=255)
                                                                     ttGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCgccgcc                                                                       >  2:265538/1‑150 (MQ=255)
                                                                          gaccccgaccAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCaa                                                                  >  1:347378/1‑150 (MQ=255)
                                                                                gaccAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGcgccgc                                                                        <  1:505958/138‑1 (MQ=255)
                                                                                gaccAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGcgccgc                                                                        >  2:505958/1‑138 (MQ=255)
                                                                                gaccAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGAc                                                            >  1:419120/1‑150 (MQ=255)
                                                                                    aGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCCCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAgg                                                        <  1:19636/150‑1 (MQ=255)
                                                                                          gCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTg                                                  <  1:250113/150‑1 (MQ=255)
                                                                                            ccGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAAACCCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGACCGCAACGAGACGCGGTGATTGGc                                                >  1:347054/1‑150 (MQ=255)
                                                                                                   cgcCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGa                                         <  1:129349/150‑1 (MQ=255)
                                                                                                   cgcCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGa                                         <  2:790974/150‑1 (MQ=255)
                                                                                                   cgcCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGAGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGa                                         <  1:490860/150‑1 (MQ=255)
                                                                                                    gcCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAg                                        >  2:589310/1‑150 (MQ=255)
                                                                                                             cACCTACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCgg                               <  1:122475/150‑1 (MQ=255)
                                                                                                              aCCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGAtt                                                   <  2:440060/129‑1 (MQ=255)
                                                                                                              aCCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGAtt                                                   >  1:440060/1‑129 (MQ=255)
                                                                                                                  acgaGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCa                          <  2:347378/150‑1 (MQ=255)
                                                                                                                    gaGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAgc                        >  2:454560/1‑150 (MQ=255)
                                                                                                                           tCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGc                 >  1:424776/1‑150 (MQ=255)
                                                                                                                             gCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTg               >  2:1036461/1‑150 (MQ=255)
                                                                                                                               cTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGc                                                >  1:837947/1‑115 (MQ=255)
                                                                                                                               cTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGc                                                <  2:837947/115‑1 (MQ=255)
                                                                                                                               cTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGc             >  2:209487/1‑150 (MQ=255)
                                                                                                                               cTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGCTGCGCCGCCCGCAAAGTGACGAGGCGATTGGCGTGCTGAGCATCGCCGCGCCAAGCGCGCGGCTGGc             >  2:346042/1‑150 (MQ=255)
                                                                                                                                   cTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGACACCTCCTCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGa                                         <  1:981074/118‑1 (MQ=255)
                                                                                                                                   cTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGa                                         >  2:981074/1‑118 (MQ=255)
                                                                                                                                        ccGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAgccg    >  1:224563/1‑150 (MQ=255)
                                                                                                                                          gcgcCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAgccgcc  >  1:539159/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
CCTCCGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCC  >  NC_002947/2974520‑2974807

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