Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,554 Δ1 bp 6.9% coding (1087/1116 nt) PP_1104 ← AotO‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,5540A. 6.9% 173.0 / 7.1 58coding (1087/1116 nt)PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base . (2/2);  total (30/28)
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.09e-01

CAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTTGGCC  >  NC_002947/1261407‑1261701
                                                                                                                                                    |                                                                                                                                                   
cAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATa                                                                                                                                                   <  1:993745/150‑1 (MQ=255)
       gcTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTATGCCGCAGCGGCTACCTGttg                                                                                                                                            <  1:482669/150‑2 (MQ=255)
         ttGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGCTTCCGATGTTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGc                                                                                                                                          <  2:904371/150‑1 (MQ=255)
          tGGCAAACCAGTCATTTGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCAGTTGACGGCTGCAATGGTTTTACCGGATCAATTTGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGcc                                                                                                                                         <  1:421974/150‑1 (MQ=255)
           ggCAATCCAGGGATTGGCGAGGTCGGTAGCAGCAGGCATGCGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACACAACACTTGCCGGCTGCAAGGGGTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCg                                                                                                                                        <  1:320176/150‑1 (MQ=255)
           ggCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCg                                                                                                                                        <  1:1266765/150‑1 (MQ=255)
             cAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCACGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGc                                                                                                                                      <  1:300560/150‑1 (MQ=255)
              aaaCCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGcc                                                                                                                                     >  2:1197362/1‑150 (MQ=255)
               aaCCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCca                                                                                                                                    >  2:7011/1‑150 (MQ=255)
                   aGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATg                                                                                                                                >  2:958922/1‑150 (MQ=255)
                           gCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGa                                                                                                                        <  2:807181/150‑1 (MQ=255)
                                gTCGGTAGCAGCAGGCATGAGGGCAGGAGGTCCTTCATTGCGTGTGCTGCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATg                                                                                                                   <  2:727053/150‑1 (MQ=255)
                                             ggCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCa                                                                                                      <  1:1188115/150‑1 (MQ=255)
                                                 tGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTATGCCGCAGCGGCTACCt                                                                                                                                                <  1:379047/104‑1 (MQ=39)
                                                 tGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTATGCCGCAGCGGCTACCt                                                                                                                                                >  2:379047/1‑104 (MQ=255)
                                                     ggCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCgg                                                                                              >  2:319685/1‑150 (MQ=255)
                                                              tCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                                     <  1:7011/150‑1 (MQ=255)
                                                                    gAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAc                                                                               <  2:782602/150‑1 (MQ=255)
                                                                     aaGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAca                                                                              >  2:271390/1‑150 (MQ=255)
                                                                      aGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcag                                                                             <  1:671636/150‑1 (MQ=255)
                                                                             gCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGACCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATCCAGCAGcccc                                                                      >  1:61998/1‑150 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                                     >  2:688540/1‑133 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                                     <  1:688540/133‑1 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCgg                                                                    >  1:898415/1‑150 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCgg                                                                    >  1:490611/1‑150 (MQ=255)
                                                                                   aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                                  <  1:1074361/100‑1 (MQ=255)
                                                                                   aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                                  >  2:1074361/1‑100 (MQ=255)
                                                                                   aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                                <  1:474455/150‑1 (MQ=255)
                                                                                   aCAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                                >  2:474455/1‑150 (MQ=255)
                                                                                     aGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctggct                                                              >  1:955952/1‑150 (MQ=255)
                                                                                     aGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctggct                                                              <  2:864454/150‑1 (MQ=255)
                                                                                        ccACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                           <  1:257996/150‑1 (MQ=255)
                                                                                        ccACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                           <  1:966452/150‑1 (MQ=255)
                                                                                         cACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTg                                                          >  2:99373/1‑150 (MQ=255)
                                                                                          aCTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagcag                                                                          <  2:898590/133‑1 (MQ=255)
                                                                                          aCTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagcag                                                                          >  1:898590/1‑133 (MQ=255)
                                                                                          aCTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTgc                                                         <  1:99373/150‑1 (MQ=255)
                                                                                            ttGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGAcc                                                                                                                      >  1:393114/1‑87 (MQ=255)
                                                                                            ttGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGAcc                                                                                                                      <  2:393114/87‑1 (MQ=255)
                                                                                            ttGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTgcgc                                                       >  1:444838/1‑150 (MQ=255)
                                                                                                    tGCAATGGTTTTACCGGATCAATTGGCCACCAAGCCACTGGCGTGGCCCCGACGTGCACCACCCACAGGGTGGTGGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGt                                               >  1:27410/1‑150 (MQ=255)
                                                                                                    tGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGt                                               <  1:439592/150‑1 (MQ=255)
                                                                                                      cAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGAGATGCTCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                             >  1:1085560/1‑150 (MQ=255)
                                                                                                       aaTGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGAc                                            >  2:317754/1‑150 (MQ=255)
                                                                                                                  cGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGa                                 >  2:145105/1‑150 (MQ=255)
                                                                                                                   ggATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCgg                                                                    >  2:1304273/1‑114 (MQ=255)
                                                                                                                   ggATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCgg                                                                    <  1:1304273/114‑1 (MQ=255)
                                                                                                                    gATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                               <  1:1226491/150‑1 (MQ=255)
                                                                                                                     aTCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTg                              >  1:272489/1‑150 (MQ=255)
                                                                                                                           tGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCGGCCCGGGCTAGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGa                        <  1:958922/150‑1 (MQ=255)
                                                                                                                             gcgaaCAGGTATGCCGCAGCGGCTACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATg                      >  1:838125/4‑150 (MQ=255)
                                                                                                                                aaCAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACt                   >  1:213508/1‑150 (MQ=255)
                                                                                                                                aaCAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACt                   >  1:1346984/1‑150 (MQ=255)
                                                                                                                                       aGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGAATTCGGCGa            >  2:278858/1‑150 (MQ=255)
                                                                                                                                          ccGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATc          >  2:773741/1‑150 (MQ=255)
                                                                                                                                            gCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAg        >  1:577622/1‑150 (MQ=255)
                                                                                                                                            gCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAg        >  2:350432/1‑150 (MQ=255)
                                                                                                                                                  gCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTTggcc  <  2:955952/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTAGCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTTGGCC  >  NC_002947/1261407‑1261701

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