Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,731 C→G 7.3% G136G (GGG→GGC)  PP_0162 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947185,7310CG 7.3% 166.2 / 4.1 55G136G (GGG→GGC) PP_0162ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base G (2/2);  total (26/29)
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.70e-01

TTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTC  >  NC_002947/185588‑185877
                                                                                                                                               |                                                                                                                                                  
ttGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCt                                                                                                                                                <  7:274354/148‑1 (MQ=255)
       aGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGAGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGCTCCCCATCTGTTCGGCTATGGCCTGATGCTTCAGCCCCTCGCATTGa                                                                                                                                        <  5:78449/149‑1 (MQ=255)
          gCGTGCGACGGTCATGGATTCAGCAGTGTTCTTCATAGCGCAGCACGTAGCAGTTATAGAGCGCATTGGCGATGTATCGCTCAAGGGTGCGCAGGGATATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGc                                                                                                                                      <  7:171213/148‑1 (MQ=255)
           cGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCca                                                                                                                                    <  2:183490/149‑1 (MQ=255)
                        tGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGcgcg                                                                                                                       <  5:269803/149‑1 (MQ=255)
                          gATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGAc                                                                                                                     >  6:78802/1‑149 (MQ=255)
                               agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                                >  6:247877/1‑149 (MQ=255)
                                 cagTCGTCTTCATAGCGCAGCACCTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCtt                                                                                                              >  2:63344/1‑149 (MQ=255)
                                  agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                             <  1:133492/149‑1 (MQ=255)
                                  agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                             >  6:46455/1‑149 (MQ=255)
                                      gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                                         >  1:322943/1‑149 (MQ=255)
                                      gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                                         <  4:60362/149‑1 (MQ=255)
                                      gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                                         >  4:332920/1‑149 (MQ=255)
                                       tcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGttt                                                                                                        <  2:322943/149‑1 (MQ=255)
                                         ttCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGAACCTTGGGttt                                                                                                        >  7:54499/1‑147 (MQ=255)
                                            aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAgg                                                                                                    <  5:46455/148‑1 (MQ=255)
                                            aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAgg                                                                                                    <  3:332920/148‑1 (MQ=255)
                                                 gcagcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCACGAATGCGCGACGTACCTAGGGATTGAGGCCGTc                                                                                               >  3:297832/1‑148 (MQ=255)
                                                    gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                           <  5:78802/149‑1 (MQ=255)
                                                    gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                           <  1:219632/149‑1 (MQ=255)
                                                           gCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCGCATCTGTTCGGCGCTGGCCTTATGCTTCAGCCCCTCGCATTGAGCCAGCATGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGt                                                                                     <  1:158109/148‑1 (MQ=255)
                                                                 aTAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATAc                                                                              >  7:134791/1‑149 (MQ=255)
                                                                  tAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACg                                                                             >  2:257864/1‑149 (MQ=255)
                                                                           ttGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATCCTTCAGCCCCACGCATTTAGCCAGCAGGAATCCCCGACGTACCTTGGGTTTGAGCCCGACGAGCATGCCTCCAATACGCTCAAGCaa                                                                    >  3:277777/1‑149 (MQ=255)
                                                                              gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTgg                                                                                                            <  2:257045/106‑1 (MQ=255)
                                                                              gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTc                                                                 >  3:178895/1‑149 (MQ=255)
                                                                              gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGCCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTgg                                                                                                            >  1:257045/1‑106 (MQ=255)
                                                                                gATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                               <  4:178895/149‑1 (MQ=255)
                                                                                   gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                                     >  8:176606/1‑108 (MQ=255)
                                                                                   gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                                     <  7:176606/108‑1 (MQ=255)
                                                                                      gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                               >  7:23314/1‑143 (MQ=255)
                                                                                      gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                               <  8:23314/143‑1 (MQ=255)
                                                                                                    gcACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCTCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGtctc                                            >  8:139020/1‑148 (MQ=255)
                                                                                                    gcACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGtctc                                            >  7:209207/1‑148 (MQ=255)
                                                                                                     cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                                          <  4:297832/149‑1 (MQ=255)
                                                                                                     cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                                          >  3:152935/1‑149 (MQ=255)
                                                                                                              ccccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCCCGACGAACCAGGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                                 >  8:184928/1‑149 (MQ=255)
                                                                                                               cccATCTGTTCGGCGATGGCGTGATGCTTCTGCCCCTCGCATTGAGCCAGCAGGAATTCGCGACGTACCTTGTGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTGCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                                <  4:277777/149‑1 (MQ=255)
                                                                                                                ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAcc                               <  1:63344/149‑1 (MQ=255)
                                                                                                                ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                                >  1:284068/1‑148 (MQ=255)
                                                                                                                 cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                     >  2:219328/1‑142 (MQ=255)
                                                                                                                 cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                     <  1:219328/142‑1 (MQ=255)
                                                                                                                  aTCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAg                             <  8:209207/149‑1 (MQ=255)
                                                                                                                     ctcaatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTg                           <  7:178506/142‑1 (MQ=255)
                                                                                                                     ctcaatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTg                           >  8:178506/7‑148 (MQ=255)
                                                                                                                       ttCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCATGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGt                            <  1:257864/145‑1 (MQ=255)
                                                                                                                         atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc                      >  6:249693/3‑149 (MQ=255)
                                                                                                                         atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc                      <  5:249693/147‑1 (MQ=255)
                                                                                                                             gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                                 <  6:209370/118‑1 (MQ=255)
                                                                                                                             gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                                 >  5:209370/1‑118 (MQ=255)
                                                                                                                                  ccTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCTAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATgg             <  7:78728/149‑1 (MQ=255)
                                                                                                                                      aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc                      >  3:329329/1‑136 (MQ=255)
                                                                                                                                      aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc                      <  4:329329/136‑1 (MQ=255)
                                                                                                                                            cAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGcc    <  4:152935/148‑1 (MQ=255)
                                                                                                                                              gCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTc  >  3:310959/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTC  >  NC_002947/185588‑185877

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