Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,000,466 G→T 11.1% V122F (GTT→TTT)  PP_1786 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,000,4660GT11.1% 150.2 / 8.2 54V122F (GTT→TTT) PP_1786glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base T (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 = 7.97e-01

CTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCAGCCATTGA  >  NC_002947/2000318‑2000611
                                                                                                                                                    |                                                                                                                                                 
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)
                          cGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGcac                                                                                                                        <  1:549116/150‑1 (MQ=255)
                           ggAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGcaca                                                                                                                       >  2:1041251/1‑150 (MQ=255)
                                  aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGag                                                                                                                                              >  1:835806/1‑119 (MQ=255)
                                  aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACGCCGCCGGCTTTGag                                                                                                                                              <  2:835806/120‑2 (MQ=255)
                                           gCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                            >  2:244239/1‑129 (MQ=255)
                                           gCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                            <  1:244239/129‑1 (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                                                                                        <  2:440946/105‑1 (MQ=255)
                                                                                                       cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGa                                                                                        >  1:440946/1‑105 (MQ=255)
                                                                                                       cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTa                                           >  1:1414101/1‑150 (MQ=255)
                                                                                                        tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGc                                                                                          >  2:583284/1‑102 (MQ=255)
                                                                                                        tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGc                                                                                          <  1:583284/102‑1 (MQ=255)
                                                                                                                 aCCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTAc                                 <  2:534808/150‑1 (MQ=255)
                                                                                                                    tGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                              >  1:318041/1‑150 (MQ=255)
                                                                                                                    tGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                              >  1:99342/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)
                                                                                                                                                cggcgTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCAGCCATTGa  <  1:973527/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAGCCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGCACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCAGCCATTGA  >  NC_002947/2000318‑2000611

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