Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,000,456 A→C 11.1% K118N (AAA→AAC PP_1786 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,000,4560AC11.1% 131.3 / 10.7 54K118N (AAA→AACPP_1786glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base C (2/4);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

GGCCGTTATGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGC  >  NC_002947/2000307‑2000602
                                                                                                                                                     |                                                                                                                                                  
ggCCGTTATGACTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCaaa                                                                                                                                                    <  2:1007125/150‑1 (MQ=255)
           cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGt                                                                                                                                         <  1:1055996/150‑1 (MQ=255)
              ggCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGa                                                                                                                                      <  1:802737/150‑1 (MQ=255)
               gCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGaa                                                                                                                                     <  1:1317674/150‑1 (MQ=255)
                 ttCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGaggt                                                                                                                                   <  1:1160804/150‑4 (MQ=255)
                     cACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCaa                                                                                                                               >  2:238644/1‑150 (MQ=255)
                      aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGtt                                                                                                                                        <  1:294274/140‑1 (MQ=255)
                      aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAg                                                                                                                              >  2:335432/1‑150 (MQ=255)
                      aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAg                                                                                                                              <  1:335432/150‑1 (MQ=255)
                      aCGAAAGTGGTGATGCGGAGATTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGtt                                                                                                                                        >  2:294274/1‑140 (MQ=255)
                        gAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCaaa                                                                                                                                                    <  1:698984/126‑1 (MQ=255)
                        gAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCaaa                                                                                                                                                    >  2:698984/1‑126 (MQ=255)
                                     cGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGcac                                                                                                               <  1:549116/150‑1 (MQ=255)
                                      ggAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGcaca                                                                                                              >  2:1041251/1‑150 (MQ=255)
                                             aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGag                                                                                                                                     <  2:835806/120‑2 (MQ=255)
                                             aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGag                                                                                                                                     >  1:835806/1‑119 (MQ=255)
                                                      gCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                   <  1:244239/129‑1 (MQ=255)
                                                      gCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                   >  2:244239/1‑129 (MQ=255)
                                                      gCAGTTCTCAATGCGGCCGACAGCGCCGTTAAGCTCAAGACGGATGTTGCGAAACAGATCCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGa                                                                                              <  1:945362/150‑1 (MQ=255)
                                                       cAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGag                                                                                                                                     <  2:124024/110‑2 (MQ=255)
                                                       cAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGag                                                                                                                                     >  1:124024/1‑109 (MQ=255)
                                                         gTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGTGGAAGCACATTTCCAACGCTTCGGACgg                                                                                           >  1:1266867/1‑150 (MQ=255)
                                                         gTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACgg                                                                                           >  2:1217492/1‑150 (MQ=255)
                                                              cAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGa                                                                                      <  2:468875/150‑1 (MQ=255)
                                                               aaTGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAg                                                                                     <  1:379519/150‑1 (MQ=255)
                                                                    ggCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCg                                                                                >  2:460830/1‑150 (MQ=255)
                                                                    ggCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGAAGGGCGGAGGGGCg                                                                                >  1:298227/1‑150 (MQ=255)
                                                                      ccGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCa                                                                                                        >  2:1164845/1‑124 (MQ=255)
                                                                      ccGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCa                                                                                                        <  1:1164845/124‑1 (MQ=255)
                                                                      ccGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGAc                                                                              >  2:39483/1‑150 (MQ=255)
                                                                      ccGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGAc                                                                              <  2:430242/150‑1 (MQ=255)
                                                                         aCAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTc                                                                           <  1:940053/150‑1 (MQ=255)
                                                                         aCAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTc                                                                           <  1:39483/150‑1 (MQ=255)
                                                                         aCAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGACTTTTCACTTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCCAGCTGTGGCGGTAGCACATTTCTAACGCTTCGGACGGGCGGAGGGGCGACTTc                                                                           >  1:747538/1‑150 (MQ=255)
                                                                              gcCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTAct                                                                      >  2:481773/1‑150 (MQ=255)
                                                                                 gTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACtctc                                                                   <  1:238644/150‑1 (MQ=255)
                                                                                          aaGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTgg                                                          >  2:973527/1‑150 (MQ=255)
                                                                                             ccGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGCCCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGc                                                                                 <  1:914045/124‑1 (MQ=255)
                                                                                             ccGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGc                                                                                 >  2:914045/1‑124 (MQ=255)
                                                                                              cGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTa                                                      >  1:972830/1‑150 (MQ=255)
                                                                                                               gACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTGCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCa                                     <  1:1041251/150‑1 (MQ=255)
                                                                                                                  cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGa                                                                               >  1:440946/1‑105 (MQ=255)
                                                                                                                  cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGa                                                                               <  2:440946/105‑1 (MQ=255)
                                                                                                                  cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTa                                  >  1:1414101/1‑150 (MQ=255)
                                                                                                                   tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGc                                                                                 <  1:583284/102‑1 (MQ=255)
                                                                                                                   tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGc                                                                                 >  2:583284/1‑102 (MQ=255)
                                                                                                                            aCCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTAc                        <  2:534808/150‑1 (MQ=255)
                                                                                                                               tGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                     >  1:99342/1‑150 (MQ=255)
                                                                                                                               tGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                     >  1:318041/1‑150 (MQ=255)
                                                                                                                                gACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCa                                     >  2:853015/1‑133 (MQ=255)
                                                                                                                                gACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCa                                     <  1:853015/133‑1 (MQ=255)
                                                                                                                                gACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCt                             >  2:1396689/1‑141 (MQ=255)
                                                                                                                                gACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCt                             <  1:1396689/141‑1 (MQ=255)
                                                                                                                                 aCCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAAc                   >  1:813817/1‑150 (MQ=255)
                                                                                                                                  ccTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGTGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCTATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACg                  <  2:1266867/150‑1 (MQ=255)
                                                                                                                                         cccTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAg           >  2:171674/1‑150 (MQ=255)
                                                                                                                                                  cAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGACGCAGCCTTCTCACGATTGGGACGGGCGGAGGGGCGACGTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCAACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagc  >  1:247917/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCCGTTATGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGC  >  NC_002947/2000307‑2000602

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: