Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,990 A→T 8.8% I64F (ATC→TTC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9900AT 8.8% 153.1 / 7.0 57I64F (ATC→TTC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base T (3/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 = 4.22e-01

TGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGAT  >  NC_002947/3762841‑3763139
                                                                                                                                                     |                                                                                                                                                     
tGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGa                                                                                                                                                       <  1:293715/150‑1 (MQ=255)
  gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGAAGTGGCACTAACGAGGCCGAGTACGACCACGGCGTCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTTGACGTGACCAAGCACATGAta                                                                                                                                                     >  1:326453/1‑149 (MQ=255)
        gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGATCTACGTGACCGGGAACATGATCCCGGTc                                                                                                                                               >  1:1081133/1‑150 (MQ=255)
          cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                             <  2:628699/150‑1 (MQ=255)
              gCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTAACCGGGAACATGATCCCGGTCACGGTc                                                                                                                                         <  1:1046318/150‑1 (MQ=255)
                aGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                                                                                                       <  1:678671/150‑1 (MQ=255)
                aGGCGCCCACGGCCTGTGGAACGAGGAACGACGTGACACACTCGAGGTCGTGTACGGCCACGGGGCAGAGAGTAACGCGTACTATGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                                                                                                       >  2:971860/1‑150 (MQ=255)
                  gcgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCACGGTCACGGTCGAgc                                                                                                                                     >  2:998768/1‑150 (MQ=255)
                     cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGcc                                                                                                                                  <  1:498151/150‑1 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGCGGATAACGCCTACAAGGCGGACAAGATCAAGCGCGCCTAGGCGTTCGACGTGACCGGGAACATGATCCCGGTCCCCGTCGAGCGccc                                                                                                                                 >  2:148281/1‑149 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                                >  1:362362/1‑135 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGt                                                                                                                                                <  2:362362/135‑1 (MQ=255)
                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCt                                                                                                                                 >  2:1108639/1‑150 (MQ=255)
                        aCGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACGGGGAACATGATCCCGGTCGCGGTCGATCTCCttg                                                                                                                               >  2:797256/1‑148 (MQ=255)
                             cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGAt                                                                                                                          >  2:163568/1‑150 (MQ=255)
                                  gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGc                                                                                                                     >  2:599351/1‑150 (MQ=255)
                                  gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGc                                                                                                                     >  1:767281/1‑150 (MQ=255)
                                      gaggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAg                                                                                                                                      >  1:119828/1‑129 (MQ=255)
                                      gaggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAg                                                                                                                                      <  2:119828/129‑1 (MQ=255)
                                          ggCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgctgc                                                                                                             <  2:326453/149‑1 (MQ=255)
                                                   gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                                                    >  2:66873/1‑150 (MQ=255)
                                                           gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGt                                                                                                                                                >  1:912799/1‑98 (MQ=255)
                                                           gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGt                                                                                                                                                <  2:912799/98‑1 (MQ=255)
                                                            aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                                           <  1:222026/150‑1 (MQ=255)
                                                              gCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAACATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGCCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCACCTGCg                                                                                         <  1:577616/150‑1 (MQ=255)
                                                                    tACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACtt                                                                                   <  2:748601/150‑1 (MQ=255)
                                                                                 ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                             >  1:1054938/1‑79 (MQ=255)
                                                                                 ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                             <  2:1054938/79‑1 (MQ=255)
                                                                                 ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTgg                                                                      <  2:649988/150‑1 (MQ=255)
                                                                                 ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTgg                                                                      <  1:337130/150‑1 (MQ=255)
                                                                                  cGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGa                                                                     >  1:604839/1‑150 (MQ=255)
                                                                                        tAACGCCTACATGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCACGCCGTGCTGCAACCACTTACCCCGCCTACGGCACTTGCCGTGCTGCTGGACAACgg                                                               <  2:584622/150‑1 (MQ=255)
                                                                                             ccTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCtt                                                          <  2:1062754/150‑1 (MQ=255)
                                                                                              cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                         <  1:66873/150‑1 (MQ=255)
                                                                                                aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcccc                                                                                                >  1:942289/1‑109 (MQ=255)
                                                                                                aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcccc                                                                                                <  2:942289/109‑1 (MQ=255)
                                                                                                   aGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcac                                                    >  1:430127/1‑150 (MQ=255)
                                                                                                    ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCTGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                   <  2:604839/150‑1 (MQ=255)
                                                                                                    ggCGAAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                                   <  1:797256/150‑1 (MQ=255)
                                                                                                           aaGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                                            >  2:979399/1‑150 (MQ=255)
                                                                                                                  acctcgcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTGGGTCACAGCGCCCGGACAAGc                                     <  1:496662/146‑1 (MQ=255)
                                                                                                                       gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                         <  2:185677/93‑1 (MQ=255)
                                                                                                                       gcCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                         >  1:185677/1‑93 (MQ=255)
                                                                                                                         cTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                              <  2:450068/150‑1 (MQ=255)
                                                                                                                             gCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAAt                          <  1:599351/150‑1 (MQ=255)
                                                                                                                             gCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGGGCTGCTGGACAACGGGCCTTCGTCACAACCGCCGGCCAAGCCATGGGGCaag                          >  2:63493/1‑149 (MQ=255)
                                                                                                                                     cGTGACCGGGATCATGTTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCg                  >  1:1085824/1‑150 (MQ=255)
                                                                                                                                       tGACCGGGATCATGTTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                                 >  1:306331/1‑85 (MQ=255)
                                                                                                                                       tGACCGGGATCATGTTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                                 <  2:306331/85‑1 (MQ=255)
                                                                                                                                        gACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc               <  1:971860/150‑1 (MQ=255)
                                                                                                                                        gACCGGGAACATGATCAAAGTGACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc               >  1:933512/1‑150 (MQ=255)
                                                                                                                                            gggAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAt           >  2:1102142/1‑150 (MQ=255)
                                                                                                                                             ggAACATGATCCCGGTCACGGTCGAGCGTCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                                <  2:583674/128‑1 (MQ=255)
                                                                                                                                             ggAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                                >  1:583674/1‑128 (MQ=255)
                                                                                                                                                  aTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTTCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg     >  1:529837/1‑150 (MQ=255)
                                                                                                                                                  aTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg     >  1:529787/1‑150 (MQ=255)
                                                                                                                                                    gATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGa   >  2:124960/1‑150 (MQ=255)
                                                                                                                                                    gATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGa   >  2:1045613/1‑150 (MQ=255)
                                                                                                                                                     aTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAANGAt  <  2:699680/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGAT  >  NC_002947/3762841‑3763139

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