Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,187,215 C→G D622E (GAC→GAG rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,187,2150CG98.2% 172.1 / NA 56D622E (GAC→GAGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base G (22/32);  minor base T (0/1);  total (23/33)
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.97e-01

CTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGT  >  NC_000913_3_bme_pgi/4187067‑4187357
                                                                                                                                                    |                                                                                                                                              
cTGTGGGCGGTGCAATAGTGTGGATGATTGGGCCGAAAGGTCGGCCTTGCTCCATCGTCAGCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCAGTCTCGGTCGGAAACCGACCGTTATTTTTGCGGAc                                                                                                                                                <  2:622957/149‑1 (MQ=255)
    tGGCCGGGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGa                                                                                                                                            <  2:310237/149‑1 (MQ=255)
      gCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTCCTCCGTCGTTGGTCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATc                                                                                                                                          <  1:23231/149‑1 (MQ=255)
                  tgtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCAAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGc                                                                                                                              <  2:443300/149‑1 (MQ=255)
                     ggATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTc                                                                                                                           >  2:50486/1‑149 (MQ=255)
                      gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGTCTTCg                                                                                                                          >  2:50404/1‑149 (MQ=255)
                      gatgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCCCATTCTCGATCTGAAACCGACCGTTATTTTTGCGGCGCAGATCATGTACACCGGCTTcc                                                                                                                          >  2:87339/1‑148 (MQ=255)
                         gatTGTAGCGAAAGGTCCGCATTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACGCCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCt                                                                                                                       <  2:749865/149‑1 (MQ=255)
                              tACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCa                                                                                                                  <  1:477376/149‑1 (MQ=255)
                                ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGGCAGATCATGTACACCGGCTTCGCCTATGCAgc                                                                                                                >  1:53498/1‑149 (MQ=255)
                                ccGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCAGGGTAAAAAAGCAATCTCCAAAATGCTGAAAACCTGCTACCGCATTCCCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCCTCGCCTATGccgg                                                                                                                >  2:33667/1‑146 (MQ=255)
                                  gAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGACTTCGCCTATGCAgcgc                                                                                                              >  2:667884/1‑149 (MQ=255)
                                  gAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTAAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAgcgc                                                                                                              >  2:44908/1‑149 (MQ=255)
                                  gAAAGGGGTGCCTTACTCCAACGCTAAACAAGCGCGGGGTAAAAAAGCGATTTCCAAAAGGCTGAACACCTGCTACCGCATGCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAgcgc                                                                                                              <  2:427148/149‑1 (MQ=255)
                                            ccTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGc                                                                                                    >  1:717294/1‑149 (MQ=255)
                                                   ccacacacAACCAGGCGCTGGGTAAAAAAGCAATCCCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                             <  2:280716/142‑1 (MQ=255)
                                                    cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTAGACCGGCTTCGCCTATGCGGGCCGTTCTGGTTCATCTGTTg                                                                                            >  1:172971/1‑149 (MQ=255)
                                                    cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCGGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTg                                                                                            <  2:509605/149‑1 (MQ=255)
                                                       cGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTa                                                                                         <  1:681373/149‑1 (MQ=255)
                                                       cGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTa                                                                                         <  1:778605/149‑1 (MQ=255)
                                                            acgacgcgcTGGGTAAAAAAGCAACCTCCAAACTGCTGAACAGCTGCTACCGCATTCTCTGTCTGAAACCGACCGTTATTTTAGCGAAGCTGGTGAGTTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGag                                                                                    <  2:490722/144‑2 (MQ=255)
                                                              cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGGGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGa                                                                                  <  2:547165/149‑1 (MQ=255)
                                                              cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGt                                                                                                      <  1:313706/129‑1 (MQ=255)
                                                              cAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGt                                                                                                      >  2:313706/1‑129 (MQ=255)
                                                               aGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGAc                                                                                 >  1:775865/1‑149 (MQ=255)
                                                                   gcgtggTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATgg                                                                             <  1:87583/145‑1 (MQ=255)
                                                                                gCAATCTCCAAAATGCTGCACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaa                                                                <  2:306315/149‑1 (MQ=255)
                                                                                 cAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGagaag                                                               <  1:696265/149‑1 (MQ=255)
                                                                                   aTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAgaagaa                                                             <  1:44908/149‑1 (MQ=255)
                                                                                      tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCAGGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAaca                                                          >  2:647526/1‑149 (MQ=255)
                                                                                      tcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTAAAACCGACCGTTATACTTGCGGAGCAGAGCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGTTATCGAGGACATGGTCATCCCGGAGAAGAAaca                                                          >  2:291210/1‑149 (MQ=255)
                                                                                       ccAAAATGCTGAACACCTGCTACCGCCTTCTCGGTCTGAAACCGACCGTTATTTTTGCGGATCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                         <  1:107611/149‑1 (MQ=255)
                                                                                       ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATTGTCGTCCCGGAGAAGAAacac                                                         <  2:452362/149‑1 (MQ=255)
                                                                                       ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                         <  2:796230/149‑1 (MQ=255)
                                                                                       ccAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCGATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                         <  1:352510/149‑1 (MQ=255)
                                                                                               cTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAatcat                                                 >  1:542011/1‑149 (MQ=255)
                                                                                                      ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAg                                          <  1:729706/149‑1 (MQ=255)
                                                                                                       cttcTACCGCACTCTCGGTCTGACACCGATCGTTATTTTTGCGGCGCAGGCCACGTACACAGGCTTGGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                         <  1:506861/146‑1 (MQ=255)
                                                                                                       cTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACAGGGTCATCCCGGAAAAGAAAACCCAAATCACCTcccggg                                         >  2:375832/1‑145 (MQ=255)
                                                                                                          cTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcag                                      >  2:152399/1‑149 (MQ=255)
                                                                                                             ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaag                                   <  1:283278/149‑1 (MQ=255)
                                                                                                              cgaaTTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagc                                  <  2:660783/146‑1 (MQ=255)
                                                                                                              cGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagc                                  >  1:213606/1‑149 (MQ=255)
                                                                                                                  ttCTCGGTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTAATCCCGGAGAAGTAACACGAATTCATCTCCGAGgcagaagcagaa                              <  1:635727/149‑1 (MQ=255)
                                                                                                                       ggTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGc                         <  2:572601/149‑1 (MQ=255)
                                                                                                                        gTCTGAAACCGAGCGGTATGTATGCGGAGCAGATCATGTAGCACGGCTTCGCCTATGCATCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAAGCATCTCCGAGGCAGAAGCAGAAGTTGCt                        <  2:453980/149‑1 (MQ=255)
                                                                                                                        gTCTGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCTAGGCAGAAGCAGAAGTTGCt                        <  2:53498/149‑1 (MQ=255)
                                                                                                                        gTCTGAAACCGACCGGTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCt                        <  1:132220/149‑1 (MQ=255)
                                                                                                                           tGAAACCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaa                     >  1:358451/1‑149 (MQ=255)
                                                                                                                             aaaCCGACCGTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                   >  1:321359/1‑149 (MQ=255)
                                                                                                                                 cggccgGTATTGTTGGGGAGCAGGTCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAg               <  2:290435/146‑1 (MQ=255)
                                                                                                                                   accgTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCAGCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCGTATCCGAGGCTGATGCGGAGGTTGCTGAAGttgcgga             >  1:123523/1‑144 (MQ=255)
                                                                                                                                      gTTATTTTTGCGGAGCAGATCATGTACACTGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCTCGGAGAAGAAACACGAAATCATCTCTTATGCGGAAGCAGAAGTTGCTGAAATTCAGGAGCa          <  1:390873/149‑1 (MQ=255)
                                                                                                                                      gTTATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGACCa          >  2:374063/1‑149 (MQ=255)
                                                                                                                                       ttATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTAAGGAGCAg         >  1:35305/1‑149 (MQ=255)
                                                                                                                                       ttATTTTTGCGGAGCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCGCCGAGGCAGAAGCAGCAGTTGCTGAAATTCAGGAGcgg         >  1:215876/1‑147 (MQ=255)
                                                                                                                                              tGCGGAGCAGATCATGTTCTCCGCCTTCGCCTACGCAGCGCGTTCTGGTGCATCTGTTGGTATCGAGGACAAGGCCATCCCGGAGAAGAAACACGAAATGATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgt  <  2:401853/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGT  >  NC_000913_3_bme_pgi/4187067‑4187357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: