Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,908,070 G→C 8.5% intergenic (‑1/‑281) PP_4317 ← / → PP_4318 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,908,0700GC 8.5% 129.6 / 5.5 47intergenic (‑1/‑281)PP_4317/PP_4318hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base C (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGA  >  NC_002947/4907921‑4908217
                                                                                                                                                     |                                                                                                                                                   
cTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATg                                                                                                                                                     >  1:336627/1‑150 (MQ=255)
cTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATg                                                                                                                                                     <  2:1157659/150‑1 (MQ=255)
 tgttgtTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATgg                                                                                                                                                    >  2:706167/1‑150 (MQ=255)
     gttgctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATATTGATCTTCGCGCCTTCAGGAGGCCGAGCGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGc                                                                                                                                                >  2:512274/1‑150 (MQ=255)
         ctgttgctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCGTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAAGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGt                                                                                                                                            <  2:1010105/150‑1 (MQ=14)
               ctgttgcCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                                                      >  2:597793/1‑150 (MQ=39)
                     ccgttgccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAgg                                                                                                                                >  2:772232/1‑150 (MQ=39)
                           ccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                  >  1:224099/1‑142 (MQ=39)
                           ccgttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                  <  2:224099/142‑1 (MQ=38)
                             gttgccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCa                                                                                                                           <  1:805952/147‑1 (MQ=39)
                             gttgccggtgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCa                                                                                                                           >  2:805952/1‑147 (MQ=38)
                                 ccgttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATggg                                                                                                                    <  1:628954/150‑1 (MQ=39)
                                   gttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                  <  1:597793/150‑1 (MQ=37)
                                    ttgccgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGcc                                                                                                                 >  1:738544/1‑150 (MQ=39)
                                        cgttgccgttgATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                             <  2:833296/150‑1 (MQ=35)
                                            gCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgcgc                                                                                                         >  1:512274‑M2/1‑146 (MQ=255)
                                                                      tcttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCACAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGt                                                                               >  2:4635/1‑150 (MQ=17)
                                                                        ttgatcttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGa                                                                             <  1:730043/150‑1 (MQ=32)
                                                                              ttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCCGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGAc                                                                       >  1:147174/1‑150 (MQ=17)
                                                                                    ttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGa                                                                 <  2:1122160/150‑1 (MQ=32)
                                                                                    ttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGACGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGa                                                                 >  1:519523/1‑150 (MQ=17)
                                                                                       atcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGAcccctgcgct                                                              >  1:359465/1‑144 (MQ=17)
                                                                                          ttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGTGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgcag                                                           >  2:832574/1‑150 (MQ=17)
                                                                                          ttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgcag                                                           >  1:792924/1‑150 (MQ=32)
                                                                                            gatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAggg                                                                          >  2:307109/1‑133 (MQ=33)
                                                                                            gatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAggg                                                                          <  1:307109/133‑1 (MQ=33)
                                                                                                ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgc                                                                                                           >  1:140495‑M2/1‑94 (MQ=255)
                                                                                                ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgc                                                                                                           <  2:140495‑M2/96‑3 (MQ=255)
                                                                                                       cTTCAGGAGGCCGAACGCAGGCCTTGGGGCGGCGGGTGAGGGGGATGGATGCCCGTCAAGCGTTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGga                                              <  2:864921/150‑1 (MQ=11)
                                                                                                       cTTCAGGAGGCCGAACGCAGGCCTTGCTTAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGga                                              <  1:907183/150‑1 (MQ=17)
                                                                                                            ggaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGatgatg                                         <  1:706167/150‑1 (MQ=32)
                                                                                                             gaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATgg                                        <  1:966539/150‑1 (MQ=32)
                                                                                                             gaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCACAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATgg                                        <  1:4635/150‑1 (MQ=17)
                                                                                                                gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                     >  1:340507/1‑150 (MQ=32)
                                                                                                                      cGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGc                               <  2:307801/150‑1 (MQ=32)
                                                                                                                             cTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAGGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                        <  2:267877/150‑1 (MQ=17)
                                                                                                                              ttGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                       >  2:169403/1‑150 (MQ=32)
                                                                                                                              ttGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                       <  1:307765/150‑1 (MQ=32)
                                                                                                                               tGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTgg                      <  2:738544/150‑1 (MQ=32)
                                                                                                                                  ggAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTa                   >  1:837499/1‑150 (MQ=32)
                                                                                                                                        aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTTGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt             <  1:865974/150‑1 (MQ=17)
                                                                                                                                        aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt             >  2:392631/1‑150 (MQ=32)
                                                                                                                                        aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt             <  2:837499/150‑1 (MQ=32)
                                                                                                                                        aGGTCACCGCCATGCACGCCCGTCAAGCGCTGAGTCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt             <  2:671720/150‑1 (MQ=11)
                                                                                                                                          gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGt           >  1:948873/1‑150 (MQ=32)
                                                                                                                                            gACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt             <  2:534921/146‑1 (MQ=12)
                                                                                                                                            gACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACttttt             >  1:534921/1‑146 (MQ=12)
                                                                                                                                            gACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg         >  2:76443/1‑150 (MQ=12)
                                                                                                                                             aCGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGCTGGGGCCTTATGCGCGGTCCTGCCGGTAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc        >  1:329100/1‑150 (MQ=255)
                                                                                                                                                   aTGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa  >  2:745807/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                   
CTGTTGTTGCTGTTGCTGTTGCCGTTGCCGTTGCCGTTGCCGTTGCCGTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGA  >  NC_002947/4907921‑4908217

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