Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,282,323 Δ1 bp 8.6% coding (86/300 nt) PP_2883 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,282,3230G. 8.6% 179.4 / 11.2 58coding (86/300 nt)PP_2883hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base . (2/3);  total (26/32)
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.71e-01

ATAGCGCTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTC  >  NC_002947/3282178‑3282468
                                                                                                                                                 |                                                                                                                                                  
aTAGCGCTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCAc                                                                                                                                                >  1:72138/1‑150 (MQ=255)
aTAGCGCTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCAGCCATCGCGGCGTGCAGCCAc                                                                                                                                                >  1:73273/1‑150 (MQ=255)
         ccAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                       <  1:1069690/150‑1 (MQ=255)
           aGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtggt                                                                                                                                    >  1:868102/1‑150 (MQ=255)
            gtgtTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCACCACCTTGAggtggctg                                                                                                                                   <  1:829935/150‑4 (MQ=255)
               ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCAc                                                                                                                                >  2:682189/1‑150 (MQ=255)
               ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCAc                                                                                                                                <  2:396483/150‑1 (MQ=255)
                   aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCACCACCTTGAGGTGGCTGCACGc                                                                                                                              <  1:1058945/148‑1 (MQ=255)
                   aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCATTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCACCACCTTGAGGTGGCTGCACGc                                                                                                                              >  2:1058945/1‑148 (MQ=255)
                    ggccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCaa                                                                                                                           >  2:1418529/1‑150 (MQ=255)
                     gccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAg                                                                                                                          >  1:956092/1‑150 (MQ=255)
                      ccggccTGTTTGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGCGGAGCGTCAGTTGCCCGCCCTGGATGCCAACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAgg                                                                                                                         <  1:461596/150‑1 (MQ=255)
                      ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAgg                                                                                                                         <  2:1158894/150‑1 (MQ=255)
                      ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAgg                                                                                                                         <  1:682189/150‑1 (MQ=255)
                           cTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGTACGCCAAGGGCcac                                                                                                                    <  2:1028951/150‑1 (MQ=255)
                            tGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCcaca                                                                                                                   <  2:72138/150‑1 (MQ=255)
                            tGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCcaca                                                                                                                   <  2:73273/150‑1 (MQ=255)
                                 gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCAc                                                                                                                                <  2:1252080/132‑1 (MQ=255)
                                 gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCAc                                                                                                                                >  1:1252080/1‑132 (MQ=255)
                                     tCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAgg                                                                                                          <  1:560008/150‑1 (MQ=255)
                                     tCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAgg                                                                                                          <  1:518113/150‑1 (MQ=255)
                                     tCGGGGGCTTGGCAACGAAAGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAgg                                                                                                          >  1:877298/1‑150 (MQ=255)
                                             ttGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                          >  1:314352/1‑111 (MQ=255)
                                             ttGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                          <  2:314352/111‑1 (MQ=255)
                                                  aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                             >  1:1123214/1‑150 (MQ=255)
                                                  aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                             >  2:1010511/1‑150 (MQ=255)
                                                         gCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGTGc                                                                                      <  2:1123214/150‑1 (MQ=255)
                                                              gcgcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCg                                                                                 <  1:1418529/150‑1 (MQ=255)
                                                               cgcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCgg                                                                                <  1:703809/150‑1 (MQ=255)
                                                                 ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTg                                                                              >  1:1274981/1‑150 (MQ=255)
                                                                    ggCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAAc                                                                           <  2:956092/150‑1 (MQ=255)
                                                                     gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGc                                                                                                                       >  2:426701/1‑105 (MQ=255)
                                                                     gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGc                                                                                                                       <  1:426701/105‑1 (MQ=255)
                                                                           ggTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtgtg                                                                    <  1:823337/150‑1 (MQ=255)
                                                                                ccAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCa                                                               >  1:1183690/1‑150 (MQ=255)
                                                                                  aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                      >  1:243688/1‑78 (MQ=255)
                                                                                  aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                      <  2:243688/78‑1 (MQ=255)
                                                                                  aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAg                                                             <  2:868102/150‑1 (MQ=255)
                                                                                  aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAg                                                             >  1:632767/1‑150 (MQ=255)
                                                                                          cGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGaaa                                                     <  1:436621/150‑1 (MQ=255)
                                                                                                ttGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGcc                                               <  1:1100245/150‑1 (MQ=255)
                                                                                                     aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                         <  2:105833/135‑1 (MQ=255)
                                                                                                     aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                         >  1:105833/1‑135 (MQ=255)
                                                                                                     aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgt                                          >  2:323755/1‑150 (MQ=255)
                                                                                                       cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTg                                                                                                     <  1:594334/89‑1 (MQ=255)
                                                                                                       cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTg                                                                                                     >  2:594334/1‑89 (MQ=255)
                                                                                                         cTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTg                                      >  2:1169345/1‑150 (MQ=255)
                                                                                                          tGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTgg                                     >  2:1391598/1‑150 (MQ=255)
                                                                                                                 cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                              >  1:348187/1‑150 (MQ=255)
                                                                                                                          tCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGc                     <  2:1183690/150‑1 (MQ=255)
                                                                                                                            aGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCt                   <  2:632767/150‑1 (MQ=255)
                                                                                                                              gccgccATCGCGGCGTGCAGCCACCTCAAGGAGGTGCACGCCAAGCGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCa                 >  2:630915/1‑150 (MQ=255)
                                                                                                                                     tCGCGGCTTGCAGCCGCCTCAAGGTGGTGCACGCCAGGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTGGGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGCTCTCGGCCTCATCCGGGc          >  2:1122434/1‑150 (MQ=255)
                                                                                                                                     ccttggcgTGCACCACCTTGAGGTGGCTGCACGCCGCGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGa                                                            >  2:377190/5‑100 (MQ=21)
                                                                                                                                     ccttggcgTGCACCACCTTGAGGTGGCTGCACGCCGCGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGa                                                            <  1:377190/96‑1 (MQ=21)
                                                                                                                                         ggcgTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt      <  1:64991/150‑1 (MQ=255)
                                                                                                                                         ggcgTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt      <  2:877298/150‑1 (MQ=255)
                                                                                                                                             tGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTc  <  2:348187/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATAGCGCTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTC  >  NC_002947/3282178‑3282468

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: