Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,985 A→T 9.0% N62I (AAC→ATC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9850AT 9.0% 145.6 / 5.8 55N62I (AAC→ATC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base T (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

CTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAT  >  NC_002947/3762837‑3763130
                                                                                                                                                    |                                                                                                                                                 
cTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAAc                                                                                                                                                  >  1:622580/1‑150 (MQ=255)
    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                                                                                                                                       <  2:912799/98‑1 (MQ=255)
                                                               gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGATCATGTTCCCGGt                                                                                                                                       >  1:912799/1‑98 (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                                                             <  1:337130/150‑1 (MQ=255)
                                                                                     ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTgg                                                             <  2:649988/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)
                                                                                                                                                    |                                                                                                                                                 
CTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGACGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAT  >  NC_002947/3762837‑3763130

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: ‑