Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,085 A→G 10.9% intergenic (+53/‑77) 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,0850AG10.9% 140.2 / 8.1 55intergenic (+53/‑77)tdcG‑I/PP_5441L‑serine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/29);  new base G (3/3);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-02

TAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357938‑358229
                                                                                                                                                   |                                                                                                                                                
tAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGaa                                                                                                                                                 >  7:56211/1‑149 (MQ=255)
  aTCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAAc                                                                                                                                                <  3:152151/148‑1 (MQ=255)
   tCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtg                                                                                                                                              <  4:82810/149‑1 (MQ=255)
     cGCCCCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                            <  4:264249/149‑1 (MQ=255)
     cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                            <  6:202723/149‑1 (MQ=255)
     cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACtgtg                                                                                                                                            <  8:56211/149‑1 (MQ=255)
     cGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTATTCGCCGGCAAGCCCGCTCCTACAGGAACtgag                                                                                                                                            >  5:241613/1‑147 (MQ=255)
                     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                                                                                                                                    <  2:272322/121‑1 (MQ=255)
                                         aaGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCt                                                                                                                                    >  1:272322/1‑121 (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                                                                                                                         <  4:157551/124‑1 (MQ=255)
                                                 agagACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCggg                                                                                                                         >  3:157551/1‑124 (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                                                                                 >  3:263206/1‑127 (MQ=255)
                                                                                      gAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCaa                                                                                 <  4:263206/127‑1 (MQ=255)
                                                                                           cTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTACTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTg                                                      >  8:93763/1‑149 (MQ=255)
                                                                                                 ttCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGc                                                                                                                                           <  7:198306/58‑1 (MQ=255)
                                                                                                 ttCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGc                                                                                                                                           >  8:198306/1‑58 (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                                                   >  1:42527/1‑108 (MQ=21)
                                                                                                                                       cGCTCCTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCgg                                                   <  2:42527/108‑1 (MQ=21)
                                                                                                                                       cGCTCCTACAGGAACTGTGCAATCCCTGTCGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCa          <  6:93235/149‑1 (MQ=255)
                                                                                                                                            cTACAGGGATTGCACAGTTCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg     >  7:255420/1‑149 (MQ=37)
                                                                                                                                            cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg     >  1:116454/1‑149 (MQ=255)
                                                                                                                                            cTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGtgtg     >  3:162088/1‑149 (MQ=255)
                                                                                                                                               cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTgg    <  5:134049/147‑1 (MQ=255)
                                                                                                                                               cAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGCTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGc  >  6:264340/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TAATCCGCACCATGCGCGATACCGGCGCCGATATGCACGACAAGTACAAAGAGACCTCGCGCGGCGGCCTGGCGGTCAGCGCTATCGAGTGCTGATTTTCTGAGCTGACCTTCTGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGAACTGTGCAATCCCTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAGGCTGATGCAGGTGCAACACAACCGCCTATGCCTCTTTATGGTGCGGTGATGGGCGCCCCGCCTACCGATCGTCGGATGTCGCGATCAGTGTGGGC  >  NC_002947/357938‑358229

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