Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,974,671 T→C 7.9% R187R (CGT→CGC)  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,6710TC 7.9% 129.2 / 2.7 51R187R (CGT→CGC) PP_2601IclR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (29/18);  new base C (2/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCCTCGCAGCGCTG  >  NC_002947/2974524‑2974818
                                                                                                                                                   |                                                                                                                                                   
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                                                    <  2:790974/150‑1 (MQ=255)
                                                                                               cgcCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGa                                                    <  1:129349/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                                                              >  1:440060/1‑129 (MQ=255)
                                                                                                          aCCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGAtt                                                              <  2:440060/129‑1 (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)
                                                                                                                                           cgcgAGCGTGGTTATGCCTGGGTCGAGGAAACCCCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTTACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCCTCGCa        <  1:383789/150‑1 (MQ=255)
                                                                                                                                              gAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCCTCGCAgcg     >  1:135820/1‑150 (MQ=255)
                                                                                                                                                 cGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCCTCGCAGCGCTg  <  2:141823/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGGCACGCCTGGCTGGCCAGCCTGGACGATGAGCAGGCGCTGCAGATGGTGTTGCGCCAAGGTATTGCCGACCCCGACCAGTTCGGCCCGAACGCGCCACGTTCCACCGACGAGCTGCTCGCCTACCTGCACCGCGCCCGCGAGCGTGGTTATGCCTGGGTCGAGGAAACCTCCGCCATCGGCACCTCAGCCCTGGCCGCAGTGGTGCGCCGCCCGCAAAGTGACGAGGTGATTGGCGTGCTGAGCATCGCCGGGCCCAGCGCGCGGCTGGCACAGAGCCGCCTCGCAGCGCTG  >  NC_002947/2974524‑2974818

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