Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,096 C→T 11.7% intergenic (+64/‑66) tdcG‑I → / → PP_5441 L‑serine dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947358,0960CT11.7% 132.3 / 7.5 51intergenic (+64/‑66)tdcG‑I/PP_5441L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base T (3/3);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-01

CCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACAC  >  NC_002947/357959‑358243
                                                                                                                                         |                                                                                                                                                   
ccGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGc                                                                                                                                          <  3:298201/149‑1 (MQ=255)
ccGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTCAAGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGc                                                                                                                                          <  1:209025/149‑1 (MQ=25)
   gcgcCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCAGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCggg                                                                                                                                       <  5:320592/149‑1 (MQ=25)
      ccGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCtt                                                                                                                                    <  2:89349/149‑1 (MQ=255)
        gATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTCCAATCCCTGTAGGAGCGGGCTTGc                                                                                                                                  <  1:35645/149‑1 (MQ=255)
             gcccGACATGTACAAAGAGACGGCGCGCGGCGGCCTGGCGGTCAGCGCTATGGTGTGCTGATTTTCTGAGCTGACCTTCTGGGCGTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTTGGAGCGGGCTTGCCcgc                                                                                                                              <  4:309203/145‑1 (MQ=38)
              cACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCcgcg                                                                                                                             >  8:91115/1‑148 (MQ=255)
                 gACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                                           >  5:312267/1‑147 (MQ=255)
                    aaGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                  >  1:272322/1‑121 (MQ=255)
                    aaGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                                  <  2:272322/121‑1 (MQ=255)
                      gTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGccc                                                                                                                    >  6:260809/1‑149 (MQ=255)
                       tACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCAACTACAGGAACTGTGCAATCCCGGTAGGAGCGGGCTTGCCCGCGAAGAGGccca                                                                                                                   >  7:19416/1‑148 (MQ=255)
                         aaaagttaCCTCGCGCGGGGGCCTGGCGTGCTGCGCTATCGAGGGCTGTTTTGCTTAGCTGACCTTCTGGCCCTCTTCTCGTGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTATGTGCGGGCTTGCCCGCGAAGAGGCCCTTg                                                                                                                 <  6:241613/142‑1 (MQ=25)
                            cgagcCCTCCCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTCCTGATTTTCTGACCTCCCCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGCATCCGTGCAACCCCTCTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAgg                                                                                                              <  1:76472/148‑1 (MQ=37)
                            agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                                       >  3:157551/1‑124 (MQ=255)
                            agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                                       <  4:157551/124‑1 (MQ=255)
                               gaCCTCGCGCGGCGGCCGGGCGCTCAGCGCTAGCGCGTGCTGATTTTCTGAGCGGACCTGCTGGCCCTCTGCGCGGGCAAGCCCGCTCCTACAGCACCTGTGCAATCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTg                                                                                                           <  5:144783/149‑1 (MQ=37)
                                    cgcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCGGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATg                                                                                                        >  8:63418/1‑147 (MQ=255)
                                     gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGCTTTTCTGAGCTGCCCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAg                                                                                                     <  5:307896/149‑1 (MQ=255)
                                     gcgcgGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGCTTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAg                                                                                                     <  5:312603/149‑1 (MQ=255)
                                          gcggcCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGc                                                                                                 >  6:304916/1‑148 (MQ=255)
                                                 ggCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGAGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAAc                                                                                         >  6:134049/1‑149 (MQ=255)
                                                        aGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTa                                                                                   <  5:304916/148‑1 (MQ=255)
                                                        aGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTAt                                                                                  <  5:260809/149‑1 (MQ=255)
                                                          cgcTATCGTGTGCTTTTTTTGTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTAAAGCAACTGTGCAATCCCTGTTGGAGGGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGc                                                                                <  8:25219/149‑1 (MQ=37)
                                                            cTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCct                                                                              >  2:93549/1‑149 (MQ=255)
                                                               tCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCtt                                                                           >  2:230239/1‑149 (MQ=255)
                                                                 gAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCaa                                                                                               <  4:263206/127‑1 (MQ=255)
                                                                 gAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCaa                                                                                               >  3:263206/1‑127 (MQ=255)
                                                                      cTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTACTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTg                                                                    >  8:93763/1‑149 (MQ=255)
                                                                                        cTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGc                                                  <  7:91115/149‑1 (MQ=255)
                                                                                        atgcTGGCCCCCTTCGCGGCCACCCCCGCCCAAACACCACCTGTCCACCCCCTGTAGGAGCGGGCTTCCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGc                                                  <  5:15488/146‑1 (MQ=255)
                                                                                              ggCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTAc                                              <  7:63418/147‑1 (MQ=255)
                                                                                                  tctTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcg                                        >  5:93235/1‑149 (MQ=255)
                                                                                                     tCGCGGGCAAGCCCGCTCCTACATGAACTGTGCAATCCCTGTTGGAGGGGGCTTGCCCGCGAATAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGAtcgtc                                      <  4:222216/148‑1 (MQ=255)
                                                                                                       gcgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCgg                                    <  6:312267/148‑1 (MQ=255)
                                                                                                        cgGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCCACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGAt                                  >  2:201536/1‑149 (MQ=255)
                                                                                                                cccGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt                          >  5:51891/1‑149 (MQ=37)
                                                                                                                cccGCTCCAACCGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGAt                          <  7:11845/149‑1 (MQ=255)
                                                                                                                  cGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCgg                                                                 <  2:42527/108‑1 (MQ=21)
                                                                                                                  cGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCgg                                                                 >  1:42527/1‑108 (MQ=21)
                                                                                                                  cGCTCCTACAGGAACTGTGCAATCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa                        <  6:93235/149‑1 (MQ=255)
                                                                                                                       cTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg                   >  7:255420/1‑149 (MQ=37)
                                                                                                                       cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg                   >  3:162088/1‑149 (MQ=255)
                                                                                                                       cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg                   >  1:116454/1‑149 (MQ=255)
                                                                                                                          cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg                  <  5:134049/147‑1 (MQ=255)
                                                                                                                          cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGCTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGc                >  6:264340/1‑149 (MQ=255)
                                                                                                                                   tgCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCgg        <  7:93763/148‑1 (MQ=255)
                                                                                                                                     cAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGa     <  1:93549/149‑1 (MQ=255)
                                                                                                                                       aTCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGaca   <  4:162088/149‑1 (MQ=255)
                                                                                                                                        tCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGacac  <  8:291855/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
CCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGCGTTGTCGGTGACAC  >  NC_002947/357959‑358243

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