Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,819,795:1 +A 8.9% intergenic (+32/‑443) rpoS → / → PP_1624 RNA polymerase sigma S factor/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,819,7951.A 8.9% 175.2 / 9.3 45intergenic (+32/‑443)rpoS/PP_1624RNA polymerase sigma S factor/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base . (19/22);  new base A (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

GGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGG  >  NC_002947/1819647‑1819943
                                                                                                                                                     |                                                                                                                                                    
gggCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGAtc                                                                                                                                                       >  6:315672/1‑149 (MQ=255)
            ggCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGtttt                                                                                                                                          >  6:169760/1‑149 (MQ=255)
             gCCTTACCCGTGAACGGGTGCGGCAAATCCATGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGttttt                                                                                                                                         <  7:87173/149‑1 (MQ=255)
               cTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTtgt                                                                                                                                       <  6:119047/149‑1 (MQ=255)
               cTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTtgt                                                                                                                                       <  6:274230/149‑1 (MQ=255)
                ttACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTtgtg                                                                                                                                      >  7:250295/1‑149 (MQ=255)
                    ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGAGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTgcg                                                                                                                                  >  4:312810/1‑149 (MQ=255)
                            gggTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCAt                                                                                                                          <  2:16813/149‑1 (MQ=255)
                                  ggCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTATCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCggg                                                                                                                                                 <  5:65264/120‑1 (MQ=255)
                                  ggCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAACTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCggg                                                                                                                                                 >  6:65264/1‑120 (MQ=255)
                                   gCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGt                                                                                                                   >  6:260794/1‑149 (MQ=255)
                                        tCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTgcgc                                                                                                              <  1:76978/149‑1 (MQ=255)
                                                gAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGTGATCACCGGGGCGtttttt                                                                                                                                        >  1:264643/1‑114 (MQ=39)
                                                gAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGAGATCACCGGGGCGtttttt                                                                                                                                        <  2:264643/115‑2 (MQ=38)
                                                 aaGGGCTCAAGCGCTTGCGCGAGATCCTTGAAATGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGAGCTGTAACAAAAAACGCCCAGGTGATCTCGGGTCGTTTTTTTGTGCGCCAGGCATgg                                                                                                                        >  7:232172/1‑130 (MQ=255)
                                                       tCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGc                                                                                               >  7:90720/1‑149 (MQ=255)
                                                         aaGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                             <  5:260794/149‑1 (MQ=255)
                                                         aaGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                             >  4:149526/1‑149 (MQ=255)
                                                           gcgcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                           <  3:149526/149‑1 (MQ=255)
                                                             gcTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGtt                                                                                         <  5:315672/149‑1 (MQ=255)
                                                               ttGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                                       >  3:109532/1‑149 (MQ=255)
                                                                 gcgcGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                      <  5:169760/148‑1 (MQ=255)
                                                                           cTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                            <  7:30231/148‑1 (MQ=255)
                                                                               aaaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                        <  6:34472/148‑1 (MQ=255)
                                                                                aaaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCACGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTTCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa                                                                      >  4:221706/1‑149 (MQ=255)
                                                                                 aaGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                                     >  3:67851/1‑149 (MQ=255)
                                                                                      tGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                                <  4:67851/149‑1 (MQ=255)
                                                                                        gTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                               <  8:90720/148‑1 (MQ=255)
                                                                                                   gAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                   <  4:109532/149‑1 (MQ=255)
                                                                                                        aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                               >  4:35811/1‑148 (MQ=255)
                                                                                                        aTTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                               >  2:68210/1‑148 (MQ=255)
                                                                                                          tGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGAGCAACCCGCTTGGTTGTGGAGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacaca                                             >  3:248843/1‑148 (MQ=255)
                                                                                                               cAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                    <  4:107813/104‑1 (MQ=255)
                                                                                                               cAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                    >  3:107813/1‑104 (MQ=255)
                                                                                                                 gTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCaa                                     <  3:312810/149‑1 (MQ=255)
                                                                                                                        gcgcTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                               <  1:68210/148‑1 (MQ=255)
                                                                                                                        gcgcTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg                              >  7:290423/1‑149 (MQ=255)
                                                                                                                         cgcTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                             <  3:35811/149‑1 (MQ=255)
                                                                                                                              tAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa                         <  8:250295/148‑1 (MQ=255)
                                                                                                                                   aaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                          <  1:286747/94‑1 (MQ=11)
                                                                                                                                   aaaaaCGCCCCGAGATCACCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                          >  2:286747/1‑94 (MQ=11)
                                                                                                                                           cccGGTGATCTCGTGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                             <  3:115216/131‑1 (MQ=37)
                                                                                                                                           cccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                             >  4:115216/1‑131 (MQ=255)
                                                                                                                                           cccGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAggg           >  1:321870/1‑149 (MQ=255)
                                                                                                                                                  aTCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                            <  1:61589/109‑1 (MQ=38)
                                                                                                                                                  aTCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                            >  2:61589/1‑109 (MQ=38)
                                                                                                                                                    ctcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCggg  >  2:101884/1‑149 (MQ=34)
                                                                                                                                                    ctcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg   <  2:321870/148‑1 (MQ=34)
                                                                                                                                                    ctcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg   <  4:248843/148‑1 (MQ=34)
                                                                                                                                                    ctcGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTAGCggg  >  7:284361/1‑149 (MQ=21)
                                                                                                                                                     |                                                                                                                                                    
GGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATCTCGGGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGG  >  NC_002947/1819647‑1819943

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: