Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,575 T→C 7.4% V135A (GTC→GCC)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,5750TC 7.4% 157.6 / 2.5 54V135A (GTC→GCC) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAACCAT  >  NC_002947/604432‑604717
                                                                                                                                               |                                                                                                                                              
gcTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCgg                                                                                                                                          <  1:942974/150‑1 (MQ=255)
 cTGGTGCCTTGCATGACCATCGCACTGGAAGAGCCGCACTCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgc                                                                                                                                         <  2:303699/150‑1 (MQ=255)
     tGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                                                                                                                                     >  1:80865/1‑150 (MQ=255)
        cTTGCATGACCATCGCACTGGAAGAGCTCGACCTAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTg                                                                                                                                  >  2:1198549/1‑150 (MQ=255)
              tGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCTTCGACTTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGcc                                                                                                                            <  2:262082/150‑1 (MQ=255)
                        aCTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGt                                                                                                                  >  1:628560/1‑150 (MQ=255)
                           ggAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa                                                                                                               >  2:108265/1‑150 (MQ=255)
                              ccagCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACgg                                                                                                            >  2:922567/3‑150 (MQ=255)
                                       cccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGa                                                                                                   >  2:192106/1‑150 (MQ=255)
                                                 ctggctgTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGct                                                                                         <  2:80865/150‑2 (MQ=255)
                                                  tggctgTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCATCACCGGCGTACGGGACAAGCTTGCAc                                                                                        >  2:66836/1‑150 (MQ=255)
                                                    gcggtgCTGCGCCTGTCCACCTGGGTGCTGATCACTCTTGGCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACct                                                                                      <  1:66836/147‑1 (MQ=255)
                                                      tgtgCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACg                                                                                                                               >  1:1076239/1‑107 (MQ=255)
                                                      tgtgCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACg                                                                                                                               <  2:1076239/107‑1 (MQ=255)
                                                      tgtgCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACctgc                                                                                    <  1:127581/150‑1 (MQ=255)
                                                              gcCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTg                                                                            >  2:451225/1‑150 (MQ=255)
                                                                cTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGaa                                                                          <  2:628560/150‑1 (MQ=255)
                                                                  gTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTAc                                                                                                       >  1:668678/1‑119 (MQ=255)
                                                                  gTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTAc                                                                                                       <  2:668678/119‑1 (MQ=255)
                                                                      aCCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAg                                                                         >  1:552723/1‑145 (MQ=255)
                                                                      aCCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAg                                                                         <  2:552723/145‑1 (MQ=255)
                                                                      aCCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGaa                                                                    <  1:960369/150‑1 (MQ=255)
                                                                        cTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGt                                                                  >  2:338961/1‑150 (MQ=255)
                                                                        cTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGt                                                                  >  1:494260/1‑150 (MQ=255)
                                                                         tGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTc                                                                 <  1:294891/150‑1 (MQ=255)
                                                                            gAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAg                                                              <  1:338961/150‑1 (MQ=255)
                                                                                  aTCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGc                                                                                                           <  1:57991/99‑1 (MQ=255)
                                                                                  aTCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGc                                                                                                           >  2:57991/1‑99 (MQ=255)
                                                                                            tCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAg                                                     <  1:580218/143‑1 (MQ=255)
                                                                                            tCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAg                                                     >  2:580218/1‑143 (MQ=255)
                                                                                                cGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTg                                          <  1:192106/150‑1 (MQ=255)
                                                                                                 gcgCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTTCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTgg                                         >  1:1076225/3‑150 (MQ=255)
                                                                                                  tGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGgc                                        <  2:494260/150‑1 (MQ=255)
                                                                                                          cgcgTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTg                                <  1:108265/150‑1 (MQ=255)
                                                                                                             gtcgtGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGcc                             <  2:513894/150‑1 (MQ=255)
                                                                                                                  gATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATggg                        >  2:449957/1‑150 (MQ=255)
                                                                                                                   aTCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGt                       >  1:1153783/1‑150 (MQ=255)
                                                                                                                    tCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACaa                                                      <  2:1085090/118‑1 (MQ=255)
                                                                                                                    tCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACaa                                                      >  1:1085090/1‑118 (MQ=255)
                                                                                                                     cAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTTCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGa                     <  2:1076225/150‑1 (MQ=255)
                                                                                                                         gAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTc                 <  1:1198549/150‑1 (MQ=255)
                                                                                                                                tgGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGa          >  1:434490/1‑150 (MQ=255)
                                                                                                                                tgGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGa          >  1:1090459/1‑150 (MQ=255)
                                                                                                                                   aaCTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGt                       <  1:1091142/134‑1 (MQ=255)
                                                                                                                                   aaCTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGt                       >  2:1091142/1‑134 (MQ=255)
                                                                                                                                     cTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCg                >  1:546796/1‑139 (MQ=255)
                                                                                                                                     cTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCg                <  2:546796/139‑1 (MQ=255)
                                                                                                                                        gCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAACCAt  >  2:151515/1‑150 (MQ=255)
                                                                                                                                        gCCAAGGCCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAAc     <  1:221582/147‑1 (MQ=255)
                                                                                                                                        gCCAAGGCCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAAc     >  2:221582/1‑147 (MQ=255)
                                                                                                                                        gCCAAGGCCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAAc     >  2:220823/1‑147 (MQ=255)
                                                                                                                                        gCCAAGGCCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAAc     <  1:220823/147‑1 (MQ=255)
                                                                                                                                          cAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATgg                         >  2:1049813/1‑125 (MQ=255)
                                                                                                                                          cAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATgg                         <  1:1049813/125‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAAGCTTGCACCTGCGCTGCGTGAAGCAGAAGTCAAGGCGAACAAGCGCTGATCCTGGCGCTGCCTGGCCATGGGTGAGTTCGAGCTGATCAACCAT  >  NC_002947/604432‑604717

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