Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,304,764 C→G 7.8% intergenic (+17/+15) cysS → / ← PP_2906 cysteine‑‑tRNA ligase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,304,7640CG 7.8% 152.7 / 3.7 51intergenic (+17/+15)cysS/PP_2906cysteine‑‑tRNA ligase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base G (2/2);  total (23/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 = 7.63e-01

ATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCG  >  NC_002947/3304616‑3304911
                                                                                                                                                    |                                                                                                                                                    
aTCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCAc                                                                                                                                                      <  6:122404/149‑1 (MQ=255)
 tCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCAcc                                                                                                                                                     >  3:5458/1‑149 (MQ=255)
      gcgcgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGc                                                                                                                                                <  1:307316/149‑1 (MQ=255)
         cgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGcccc                                                                                                                                             <  4:2253/149‑1 (MQ=255)
             tGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCt                                                                                                                                         >  2:152451/1‑149 (MQ=255)
               aaGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAgt                                                                                                                                          <  1:132998/146‑2 (MQ=255)
               aaGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAgt                                                                                                                                          >  2:132998/1‑145 (MQ=255)
                         cGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCtt                                                                                                                              <  6:2360/148‑1 (MQ=255)
                            aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAc                                                                                                                          <  8:92729/149‑1 (MQ=255)
                              aaGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAccc                                                                                                                        >  2:199432/1‑149 (MQ=255)
                                gAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGc                                                                                                                      <  1:152451/149‑1 (MQ=255)
                                      ggCTGAGTCTGACCGCATTCGCGATCAGCTGAGCGGCATGGGTGTAGTGCTGGAGGACAGCAAGGGGACGACTACCTGGGGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑gg                                                                                                               <  2:107307/149‑1 (MQ=255)
                                      ggCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGc                                                                                                                  <  8:319140/147‑1 (MQ=255)
                                           aGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGt                                                                                                          <  5:177683/149‑1 (MQ=255)
                                               tGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGCCTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCaa                                                                                                      >  5:15726/1‑149 (MQ=255)
                                                gACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGcccc                                                                                                                                             >  3:236187/1‑110 (MQ=255)
                                                gACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGcccc                                                                                                                                             <  4:236187/110‑1 (MQ=255)
                                                gACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCaa                                                                                                      >  7:125380/1‑148 (MQ=255)
                                                 aCCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGATGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCa                                                                                                            <  7:172370/141‑1 (MQ=255)
                                                 aCCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGATGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCa                                                                                                            >  8:172370/1‑141 (MQ=255)
                                                  ccGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATAc                                                                                                   >  7:180420/1‑149 (MQ=255)
                                                   cGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATAcc                                                                                                  <  6:15726/149‑1 (MQ=255)
                                                      aTTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGc                                                                                               >  3:44706/1‑149 (MQ=255)
                                                        tCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGc                                                                                               <  7:49655/147‑1 (MQ=255)
                                                        tCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGc                                                                                               >  8:49655/1‑147 (MQ=255)
                                                         cgcgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                >  3:140982/1‑114 (MQ=255)
                                                         cgcgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                <  4:140982/114‑1 (MQ=255)
                                                            gATCAGCTGACCGCCATGGATGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGcc                                                                                         >  3:329489/1‑149 (MQ=255)
                                                                aGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCccca                                                                                      <  4:44706/148‑1 (MQ=255)
                                                                 gCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCcccacc                                                                                    <  8:125380/149‑1 (MQ=255)
                                                                     aCCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCa                                                                                                            >  6:73048/1‑121 (MQ=255)
                                                                     aCCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCa                                                                                                            <  5:73048/121‑1 (MQ=255)
                                                                         ccATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGc                                                                            <  3:128611/149‑1 (MQ=255)
                                                                             gggTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTcc                                                                        >  6:137142/1‑149 (MQ=255)
                                                                                          gaggaCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAg                                                            <  5:137142/148‑1 (MQ=255)
                                                                                               cAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGaa                                                      <  1:199432/149‑1 (MQ=255)
                                                                                               cAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGaa                                                      <  4:329489/149‑1 (MQ=255)
                                                                                                 gCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAc                                                    >  3:90383/1‑149 (MQ=255)
                                                                                                  cAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAc                                                    >  4:159553/1‑148 (MQ=255)
                                                                                                      ggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGc                                               <  8:180420/149‑1 (MQ=255)
                                                                                                       gCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCACCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACAAGCCTGGACCCCCCCAACTCGCTTCCCCCAATATCCACCCCGAAACAATg                                                >  1:242126/1‑147 (MQ=255)
                                                                                                         atgtcTACCTGGCGGCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCGCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAACCCATGCAgg                                            <  6:292212/145‑1 (MQ=255)
                                                                                                         aCGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAgg                                            >  8:157661/1‑149 (MQ=255)
                                                                                                                ccTGGCTTCTGGCTGACTGAAGTAACTGGGGCTTCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcgacga                                     <  4:5458/149‑1 (MQ=255)
                                                                                                                ccTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCAACCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcgac                                       <  4:90383/147‑1 (MQ=255)
                                                                                                                   ggCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATcg                                  >  3:308596/1‑149 (MQ=255)
                                                                                                                             tgaatgaaGTATCTTGGGCTGCACCGGAGCCCCATCTTTTATTTTCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACTAt                        <  2:242126/149‑1 (MQ=255)
                                                                                                                              gactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                       >  7:126600/1‑149 (MQ=255)
                                                                                                                                 tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGa                     >  1:180096/1‑148 (MQ=255)
                                                                                                                                   aaGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                  <  3:159553/149‑1 (MQ=255)
                                                                                                                                        aCTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGCTCGACGATCGCCCCCACGATCGACAACCCCa             >  4:107378/1‑149 (MQ=255)
                                                                                                                                          tGGGGCTGCGGTGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAg            >  1:89074/1‑149 (MQ=255)
                                                                                                                                                   aCCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCg  <  2:180096/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGC‑GGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCG  >  NC_002947/3304616‑3304911

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