Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,221,722 T→C 6.9% V267A (GTC→GCC)  PP_4601 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,221,7220TC 6.9% 156.6 / 2.0 58V267A (GTC→GCC) PP_4601LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base C (2/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GCGACTCGAGCTTCCAGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGA  >  NC_002947/5221584‑5221862
                                                                                                                                           |                                                                                                                                            
gCGACTCGAGCTTCCAGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGgctgctgc                                                                                                                                   >  1:925721/1‑150 (MQ=255)
   aCTCGAGCTTCCAGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCtgctgctgc                                                                                                                                >  2:589551/1‑150 (MQ=255)
     tCGAGCTTCCAGCTTATCCAGGCCATCGACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCAGCGAACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTgcgc                                                                                                                              >  1:617319/1‑150 (MQ=255)
     tCGAGCTTCCAGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTgcgc                                                                                                                              <  1:871965/150‑1 (MQ=255)
     tCGAGCTTACAGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTgcgc                                                                                                                              >  1:997181/1‑150 (MQ=255)
               aGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGc                                                                                                                    >  1:495733/1‑150 (MQ=255)
                ggttAACCAGCCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGcc                                                                                                                   <  1:1036677/148‑1 (MQ=255)
                gCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGcc                                                                                                                   >  1:77013/1‑150 (MQ=255)
                     tCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCt                                                                                                              >  2:26812/1‑150 (MQ=255)
                     tCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGACCGGgctgt                                                                                                                                      >  2:715194/1‑125 (MQ=255)
                     tCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGACCGGgctgt                                                                                                                                      <  1:715194/126‑2 (MQ=255)
                        aGGNCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTAACGGGGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCgg                                                                                                           >  2:700047/1‑150 (MQ=255)
                        aGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCtgctgctgc                                                                                                                                <  1:334414/129‑1 (MQ=255)
                        aGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCtgctgctgc                                                                                                                                >  2:334414/1‑129 (MQ=255)
                             aTCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGAtt                                                                                                      >  2:673629/1‑150 (MQ=255)
                                  cgccGGCATATGCTGCGCGATCATGTCACTGGACTGCGGCCTGGAAGAACTCAGCGAGCACCTGCGCATGGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGa                                                                                                 <  2:418771/148‑1 (MQ=255)
                                    ccGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGc                                                                                               >  2:33889/1‑150 (MQ=255)
                                     cGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGc                                                                                                                        >  2:817476/1‑124 (MQ=255)
                                     cGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGc                                                                                                                        <  1:817476/124‑1 (MQ=255)
                                      ggCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCaa                                                                                             <  1:26812/150‑1 (MQ=255)
                                      ggCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCaa                                                                                             >  2:48069/1‑150 (MQ=255)
                                      ggCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACATGCGCATCGTAACGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAACaa                                                                                             >  2:48103/1‑150 (MQ=255)
                                       gCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAAt                                                                                            <  2:495733/150‑1 (MQ=255)
                                                         tGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCggg                                                                          >  1:757743/1‑150 (MQ=255)
                                                               tGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCGCGGCGATTGCCTAGCAATGCTTTGCCGAGCCGCGGGGGCTGt                                                                    >  2:68143/1‑150 (MQ=255)
                                                                    tGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAg                                                               <  1:700047/150‑1 (MQ=255)
                                                                       ggCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGcc                                                            >  2:598540/1‑150 (MQ=255)
                                                                           tGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCt                                                        <  2:418839/150‑1 (MQ=255)
                                                                             gAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGa                                                      >  2:981463/1‑150 (MQ=255)
                                                                                    tCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATcggcgg                                               >  1:538509/1‑150 (MQ=255)
                                                                                    tCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATcggcgg                                               >  1:211807/1‑150 (MQ=255)
                                                                                     cAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCAACAGTACTTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGt                                              <  1:68143/150‑1 (MQ=255)
                                                                                            cacgcgcgcATCGTTCGGTTGGCAGACGCGGGACTCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGACGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTgg                                       <  2:286646/145‑1 (MQ=255)
                                                                                               cTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTg                                                                                                     <  1:395863/85‑1 (MQ=255)
                                                                                               cTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTg                                                                                                     >  2:395863/1‑85 (MQ=255)
                                                                                               cTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCggg                                                                          >  1:1105472/1‑112 (MQ=255)
                                                                                               cTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCggg                                                                          <  2:1105472/112‑1 (MQ=255)
                                                                                                     aTCGTACCGGTGGCAGACGCCGCCATCCACAGCCCAGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTgg                              >  1:262415/1‑150 (MQ=255)
                                                                                                           ccGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCgcggc                        <  1:33889/150‑1 (MQ=255)
                                                                                                           ccGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCgcggc                        <  1:673629/150‑1 (MQ=255)
                                                                                                                ggCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTaaa                   <  1:598540/150‑1 (MQ=255)
                                                                                                                  cAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGc                 >  1:160249/1‑150 (MQ=255)
                                                                                                                   aGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGcc                >  1:942869/1‑150 (MQ=255)
                                                                                                                       gccgccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTc            <  2:160249/150‑1 (MQ=255)
                                                                                                                       gccgccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTc            >  1:40660/1‑150 (MQ=255)
                                                                                                                         cgccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAg                                                                                 <  2:662610/79‑1 (MQ=255)
                                                                                                                         cgccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAg                                                                                 >  1:662610/1‑79 (MQ=255)
                                                                                                                          gccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAgc         <  2:757743/150‑1 (MQ=255)
                                                                                                                           ccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGc                                      <  2:435530/120‑1 (MQ=255)
                                                                                                                           ccATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGc                                      >  1:435530/1‑120 (MQ=255)
                                                                                                                            cATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAgctg       >  1:212296/1‑150 (MQ=255)
                                                                                                                              tCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAgctgct     <  2:77013/150‑1 (MQ=255)
                                                                                                                               ccACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTGTCTATGATCGGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCtt    <  2:997181/150‑1 (MQ=255)
                                                                                                                                  gcaGCCCGACCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGcc                <  1:935764/135‑1 (MQ=255)
                                                                                                                                  gcaGCCCGACCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGcc                >  2:935764/2‑136 (MQ=255)
                                                                                                                                  acaGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGa  >  1:793642/1‑150 (MQ=255)
                                                                                                                                       ccGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCaa                                                                                             <  2:526964/54‑1 (MQ=255)
                                                                                                                                       ccGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCaa                                                                                             >  1:526964/1‑54 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GCGACTCGAGCTTCCAGCTTATCCAGGCCATCAACACCGGCATATGCTGCGCGATCATGCCACTGGACTGCGGCCTGGAAGACCTCAGCGAGCACCTGCGCATCGTACCGGTGGCAGACGCCGCCATCCACAGCCCGGTCGGGCTGCTGCTGCGCCGCAGCGAGCCGCGCTCGGCGATTGCCGAGCAATGCTTTGCCGAGGCGCGGGGGCTGTTTCAGGCCAGCTGATCGGCGGTAGCCTGGCGGTACTGGCGCGGCGTAAAGCCGGTCAGCTGCTTGA  >  NC_002947/5221584‑5221862

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: