Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,765 A→C 6.9% V125G (GTA→GGA)  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,7650AC 6.9% 179.0 / 4.5 58V125G (GTA→GGA) PP_0162ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base C (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCAGATAGGCCCTGTGCAGTT  >  NC_002947/185619‑185897
                                                                                                                                                  |                                                                                                                                    
agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                                                    >  2:178019/1‑149 (MQ=255)
       gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                                                             >  7:142352/1‑149 (MQ=255)
        tcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGttt                                                                                                                            >  5:61961/1‑149 (MQ=255)
          ttCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGttt                                                                                                                            >  6:344121/1‑147 (MQ=255)
            cATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTg                                                                                                                           <  2:297244/146‑1 (MQ=255)
            cATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                                                         <  6:221297/148‑1 (MQ=255)
            cATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                                                         <  6:61961/148‑1 (MQ=255)
                   cagcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGa                                                                                                                 <  3:214640/149‑1 (MQ=255)
                          tAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCg                                                                                                          >  4:10735/1‑149 (MQ=255)
                            gCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                                               <  4:63807/142‑1 (MQ=255)
                            gCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                                               >  3:63807/1‑142 (MQ=255)
                             cAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCACGCATTGAGCCAGCAGGAATGCGCGCCGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTg                                                                                                        >  6:225345/1‑148 (MQ=255)
                              aGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGc                                                                                                       >  8:100704/1‑148 (MQ=255)
                                tGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGc                                                                                                           >  7:163164/1‑142 (MQ=255)
                                tGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGc                                                                                                           <  8:163164/142‑1 (MQ=255)
                                tGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCaa                                                                                                     >  7:197469/1‑148 (MQ=255)
                                  aTAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATAc                                                                                                  >  8:25847/1‑149 (MQ=255)
                                        cgcATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCaa                                                                                            <  1:178019/149‑1 (MQ=255)
                                               gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTc                                                                                     <  7:242603/149‑1 (MQ=255)
                                                       gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCcagcag                                                                               >  5:31055/1‑147 (MQ=255)
                                                          cTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCt                                                                          >  3:103011/1‑149 (MQ=255)
                                                           tCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                                                   >  8:263972/1‑139 (MQ=255)
                                                           tCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                                                   <  7:263972/139‑1 (MQ=255)
                                                           tCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTc                                                                         >  5:137009/1‑149 (MQ=255)
                                                                         ggAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCaa                                                                                            >  8:205269/1‑116 (MQ=255)
                                                                         ggAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCaa                                                                                            <  7:205269/116‑1 (MQ=255)
                                                                           agaTCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                                         <  6:47066/149‑1 (MQ=255)
                                                                           agaTCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                                         <  1:136806/149‑1 (MQ=255)
                                                                           agaTCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                                                         <  3:10735/149‑1 (MQ=255)
                                                                            gaTCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAgg                                                        <  6:282775/149‑1 (MQ=255)
                                                                                cccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                                                    <  8:197469/149‑1 (MQ=255)
                                                                                 ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGc                                                                                                <  3:93671/104‑1 (MQ=255)
                                                                                 ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGc                                                                                                >  4:93671/1‑104 (MQ=255)
                                                                                 ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATACGTGCAATACGCTCAAGCAATTcc                                                                                    >  6:229605/1‑116 (MQ=255)
                                                                                 ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATACGTGCAATACGCTCAAGCAATTcc                                                                                    <  5:229605/116‑1 (MQ=255)
                                                                                        ttCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGt                                                <  8:324860/145‑1 (MQ=255)
                                                                                                  gCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTCCCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATg                                  >  5:203686/1‑149 (MQ=255)
                                                                                                  gCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTCCCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATg                                  >  5:208699/1‑149 (MQ=255)
                                                                                                  gCCTGATGCTTCAGCCCCTCGCATGGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACAt                                   <  7:25847/148‑1 (MQ=255)
                                                                                                       aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATgg                                 >  8:28611/1‑145 (MQ=255)
                                                                                                       aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATgg                                 <  7:28611/145‑1 (MQ=255)
                                                                                                       aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGcc                              >  6:107834/1‑148 (MQ=255)
                                                                                                         gCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATggg                                >  2:260935/1‑144 (MQ=255)
                                                                                                         gCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATggg                                <  1:260935/144‑1 (MQ=255)
                                                                                                                ccccTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCa                    <  6:31055/149‑1 (MQ=255)
                                                                                                                ccccTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCa                    <  8:142352/149‑1 (MQ=255)
                                                                                                                  ccTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAggg                                                       <  7:257301/112‑1 (MQ=255)
                                                                                                                  ccTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAggg                                                       >  8:257301/1‑112 (MQ=255)
                                                                                                                  ccTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTc                      >  2:109518/1‑145 (MQ=255)
                                                                                                                   cTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTcc                     >  8:302555/1‑145 (MQ=255)
                                                                                                                   cTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTcc                     <  7:302555/145‑1 (MQ=255)
                                                                                                                               cagcagGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCAGATAGGCCCTGTGCa     <  1:109518/149‑1 (MQ=255)
                                                                                                                                  cagGAATGCGCGACGTCCCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCAGATAGGCCCTGTGCAGtt  <  6:203686/149‑1 (MQ=255)
                                                                                                                                  cagGAATGCGCGACGTCCCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCAGATAGGCCCTGTGCAGtt  <  6:208699/149‑1 (MQ=255)
                                                                                                                                       aTGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGt                                                <  6:183450/98‑1 (MQ=255)
                                                                                                                                       aTGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGt                                                >  5:183450/1‑98 (MQ=255)
                                                                                                                                           gcgACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                                                     <  6:252196/89‑1 (MQ=255)
                                                                                                                                           gcgACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                                                     >  5:252196/1‑89 (MQ=255)
                                                                                                                                                  |                                                                                                                                    
AGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGGGCCAAAGCCTCCAGATAGGCCCTGTGCAGTT  >  NC_002947/185619‑185897

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