Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,320,842 T→G 7.0% N148H (AAC→CAC)  PP_2041 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,320,8420TG 7.0% 166.1 / 2.5 57N148H (AAC→CAC) PP_2041LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base G (2/2);  total (30/27)
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.76e-01

AGGCTACGGGGTGGCATTACCAGGCGGCTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCGAGGCAGTGGCG  >  NC_002947/2320694‑2320990
                                                                                                                                                    |                                                                                                                                                    
aGGCTACGGGGTGGCATTACCAGGCGGCTCGCTGCGGCTTGAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTg                                                                                                                                                     <  2:1198908/150‑1 (MQ=255)
aGGCTACGGGGTGGCATTACCAGGCGGCTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTg                                                                                                                                                     >  2:599206/1‑150 (MQ=255)
        gggcggCATTACCAGGCGGCTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGcaccac                                                                                                                                             <  1:908403/146‑1 (MQ=255)
               aTTACCAGGCGGCTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGaaa                                                                                                                                      <  2:1469900/150‑1 (MQ=255)
                 tACCAGGCGGCTCGCTGCGGCTTCAGGCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGAGCCAATATAGGTCTTCGGAAAACACGGGCTCGCCGACCACATGGACAAGTTCTACGGCGTGGTGCCCAACGCCGaagaa                                                                                                                                    <  1:129152/150‑4 (MQ=255)
                  aCCAGGCGGCTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATcc                                                                                                                                   >  1:167195/1‑150 (MQ=255)
                          gcTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGcc                                                                                                                           <  2:1282913/150‑1 (MQ=255)
                                tGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTc                                                                                                                     >  1:828695/1‑150 (MQ=255)
                                tGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTc                                                                                                                     <  2:678588/150‑1 (MQ=255)
                                      ttCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACGAGTTCTTCGGCGTTGGGCACCACGCAGAAATCCACCTTGCCCAGATcacacgc                                                                                                               >  1:1572182/1‑144 (MQ=255)
                                         aGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTg                                                                                                            >  2:1450296/1‑150 (MQ=255)
                                           ccAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGc                                                                                                          >  2:1501808/1‑150 (MQ=255)
                                                     gggCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTcgcg                                                                                                >  1:1039384/1‑150 (MQ=255)
                                                            tCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAgg                                                                                         <  2:934020/150‑1 (MQ=255)
                                                            tCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAgg                                                                                         >  1:934020/1‑150 (MQ=255)
                                                              gTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGcc                                                                                       >  2:962799/1‑150 (MQ=255)
                                                               tCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCa                                                                                      >  2:1390846/1‑150 (MQ=255)
                                                               tCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGACTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCATCTTCGCGGCTCAGGccc                                                                                      >  2:419506/1‑149 (MQ=255)
                                                                      cGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCaaa                                                                               <  1:1390846/150‑1 (MQ=255)
                                                                          cccAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGt                                                                           <  2:756399/150‑1 (MQ=255)
                                                                          cccAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGt                                                                           <  1:1501808/150‑1 (MQ=255)
                                                                           ccAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTg                                                                          >  2:159233/1‑150 (MQ=255)
                                                                            caggcaggCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTTGCCACTTTCGCGGCTCAGGCCATCCCAAACGGTGa                                                                         >  1:806963/1‑150 (MQ=255)
                                                                            caggcaggCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGa                                                                         >  2:26548/1‑150 (MQ=255)
                                                                                   gcCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAACCGGTGAGCTTGAt                                                                  >  2:286937/1‑150 (MQ=255)
                                                                                             tACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAg                                                        <  2:1204116/150‑1 (MQ=255)
                                                                                                     ttCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTTTTCGGCGTGGTGCCCAACGCCGaagaa                                                                                                                                    <  2:998266/66‑4 (MQ=21)
                                                                                                     ttCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTTTTCGGCGTGGTGCCCAACGCCGaagaa                                                                                                                                    >  1:998266/1‑63 (MQ=18)
                                                                                                     ttCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGt                                                     <  1:1021236/145‑1 (MQ=255)
                                                                                                     ttCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGt                                                     >  2:1021236/1‑145 (MQ=255)
                                                                                                      tCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTgcctgg                                               <  2:167195/150‑1 (MQ=255)
                                                                                                       cGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTgcctggg                                              >  1:39796/1‑150 (MQ=255)
                                                                                                       cGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCGGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTgcctggg                                              <  2:1123548/150‑1 (MQ=255)
                                                                                                          aaaaCACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGcc                                                                                                         >  1:1481198/1‑88 (MQ=255)
                                                                                                          aaaaCACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGcc                                                                                                         <  2:1481198/88‑1 (MQ=255)
                                                                                                                 gggCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAat                                    <  1:599206/150‑1 (MQ=255)
                                                                                                                   gCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAatga                                  >  2:1285866/1‑150 (MQ=255)
                                                                                                                    cTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAatgat                                 <  2:39796/150‑1 (MQ=255)
                                                                                                                     tCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAatgatg                                <  1:962799/150‑1 (MQ=255)
                                                                                                                      cGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATgg                               <  1:966196/150‑1 (MQ=255)
                                                                                                                       ggCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACATTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGc                              >  1:369072/1‑150 (MQ=255)
                                                                                                                         cGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCg                            <  2:1039384/150‑1 (MQ=255)
                                                                                                                          gACCAAATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCgg                           >  1:606728/1‑150 (MQ=255)
                                                                                                                                    aCAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAAc                                                                              <  1:376316/89‑1 (MQ=255)
                                                                                                                                    aCAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAAc                                                                              >  2:376316/1‑89 (MQ=255)
                                                                                                                                    aCAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGt                 >  1:37524/1‑150 (MQ=255)
                                                                                                                                     cAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTg                <  1:159233/150‑1 (MQ=255)
                                                                                                                                      aaGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCgg                                                            <  1:588037/105‑1 (MQ=255)
                                                                                                                                      aaGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCgg                                                            >  2:588037/1‑105 (MQ=255)
                                                                                                                                        gTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATTGCCTGGCACCCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATTCCCGGGAAGCGTGCGTGGGCCTGGGTAATGAGGGCCGGCAACAACAGTGCCg             >  2:1580924/1‑150 (MQ=255)
                                                                                                                                        gTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTaa                                     <  2:33654/126‑1 (MQ=255)
                                                                                                                                        gTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTaa                                     >  1:33654/1‑126 (MQ=255)
                                                                                                                                        gTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCg             >  1:325642/1‑150 (MQ=255)
                                                                                                                                         gatttcGGCGTGGTGCCCAACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCGa            >  1:1114531/4‑150 (MQ=255)
                                                                                                                                          tcttcGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCGAg           <  2:812239/150‑1 (MQ=255)
                                                                                                                                           cttcGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCGAgg          <  2:606728/150‑1 (MQ=255)
                                                                                                                                            tacGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCGAGGc         >  1:714386/3‑150 (MQ=255)
                                                                                                                                                   ttGGGCACCACGCCGAAATCCACCTTTCCCAGCTTGATCGCTTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACTTTGCCGAGGCAGTGGCg  <  2:1572182/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGCTACGGGGTGGCATTACCAGGCGGCTCGCTGCGGCTTCAGCCAAGGTGATGGGCGACTCGTCAGGGCCGTGCCCAGGCAGGCCGACCCAATACAGGTCTTCGGAAAACACGGGCTCGGCGACCACATGGACAAGTTCTTCGGCGTTGGGCACCACGCCGAAATCCACCTTGCCCAGCTCGATCGCCTGGCCACCTTCGCGGCTCAGGCCATCCCAAACGGTGAGCTTGATCCCCGGGAAGCGTGCCTGGGCCTGGGTAATGATGGCCGGCAACAACAGTGCCGAGGCAGTGGCG  >  NC_002947/2320694‑2320990

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: