Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,881,128 C→A E228* (GAA→TAA)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,1280CA93.8% 252.4 / 9.4 81E228* (GAA→TAA) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (39/37);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTT  >  NC_000913/1881031‑1881223
                                                                                                 |                                                                                               
aCTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTtatc                                                                                               <  2:1213920/100‑1 (MQ=255)
     cACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGa                                                                                          >  1:1544170/1‑100 (MQ=255)
      aCTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGAt                                                                                         >  2:1594907/1‑100 (MQ=255)
      aCTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGAt                                                                                         >  2:1638387/1‑100 (MQ=255)
       cTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGAtt                                                                                        >  1:295035/1‑100 (MQ=255)
         aaaTGACCACTGGCAATGTTTTCGCTAACGTCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAg                                                                                      >  2:1683373/1‑100 (MQ=255)
         aaaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAg                                                                                      >  2:970053/1‑100 (MQ=255)
          aaTGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCTATGATTAg                                                                                      >  2:1330306/1‑99 (MQ=255)
           aTGACCACTGGCAATGTTTTCGCTAACATCGCACCAGGAACTCAGAGCAATCCCCTACCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgac                                                                                    <  1:750203/100‑1 (MQ=255)
              aCCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAgacgac                                                                                 <  2:642664/100‑1 (MQ=255)
               ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACa                                                                                >  2:127194/1‑100 (MQ=255)
               ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACa                                                                                >  1:871184/1‑100 (MQ=255)
               ccACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACa                                                                                >  1:1574660/1‑100 (MQ=255)
                  cTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTg                                                                             >  2:1342302/1‑100 (MQ=255)
                      naaTGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                         >  2:1577261/2‑100 (MQ=255)
                      cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                         >  1:51561/1‑100 (MQ=255)
                      cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                         <  1:1476673/100‑1 (MQ=255)
                      cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                         <  2:364436/100‑1 (MQ=255)
                      cAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGcc                                                                         <  2:1253457/100‑1 (MQ=255)
                        aTGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCgg                                                                       <  1:1297772/100‑1 (MQ=255)
                        aTGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCgg                                                                       <  2:1013898/100‑1 (MQ=255)
                         tGTTTTCGCTAACATCCCACCAGGAACGCGGAGCAATCCCCTGCCCATCGAGCCACCACTGATGGACAATTTATCCATCATTAGACGACaaggg                                                                            >  1:1158383/1‑90 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                     >  2:1729710/1‑100 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                     >  1:636088/1‑100 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTg                                                                     <  2:1414852/100‑1 (MQ=255)
                          gTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCCAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGACCg                                                                        >  1:337261/1‑97 (MQ=255)
                            tttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGAc                                                                   <  1:35432/100‑1 (MQ=255)
                             ttCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACt                                                                  <  2:871185/100‑1 (MQ=255)
                              tCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACtt                                                                 >  2:793037/1‑100 (MQ=255)
                               cGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACttt                                                                <  1:187223/100‑1 (MQ=255)
                                   aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                            <  2:1250053/100‑1 (MQ=255)
                                   aaCATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATc                                                            <  1:1649062/100‑1 (MQ=255)
                                      aTCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCa                                                         <  1:1289656/100‑1 (MQ=255)
                                      aTCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCa                                                         <  2:382984/100‑1 (MQ=255)
                                       tCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGAGGACAGTGGGCAGGTGACTTTAATCg                                                           >  1:1203/1‑97 (MQ=255)
                                       tCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCAt                                                        <  1:47682/100‑1 (MQ=255)
                                          caccaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATggg                                                     >  2:510718/1‑100 (MQ=255)
                                           accaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACATTGGGCCGGTGACTTTAATCGCATgggg                                                    >  1:1022382/1‑100 (MQ=255)
                                           accaGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATgggg                                                    >  2:642173/1‑100 (MQ=255)
                                              aGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGcc                                                 <  1:1469878/100‑1 (MQ=255)
                                              aGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGcc                                                 <  1:420826/100‑1 (MQ=255)
                                              aGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGcc                                                 <  1:1661570/100‑1 (MQ=255)
                                               ggAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCct                                                <  2:51561/100‑1 (MQ=255)
                                                   cGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCcttctt                                            >  2:1302727/1‑100 (MQ=255)
                                                     cAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTa                                          <  1:644072/100‑1 (MQ=255)
                                                     cAGAGCAATCCCCTGCCCATCGAGACACCACTGAGGGGCAATTTATCCATGATTAGGCGCCAGTGGGCCGGTTACTTTACTCGCATGGGAGCct                                                >  1:510450/1‑94 (MQ=255)
                                                      agagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGTCGGGTGACTTTGATGGCATGGGGGCCTTCTTTAt                                         >  1:1608566/1‑100 (MQ=255)
                                                      agagCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAt                                         <  2:530100/100‑1 (MQ=255)
                                                       nagCAATCTCCAGCCCATCGAGACACCACTGAGGGACAATTTATCCATGATTAGGCGACAGTGGGACGGAGACAATAATCACATGGCGGActtct                                             >  1:163568/2‑95 (MQ=255)
                                                        agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                       <  1:970052/100‑1 (MQ=255)
                                                        agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                       <  2:114379/100‑1 (MQ=255)
                                                        agCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTg                                       <  1:1479648/100‑1 (MQ=255)
                                                         gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTAtt                                        >  1:280365/1‑98 (MQ=255)
                                                         gCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGCGGCCCGGTGACTTTAATCGCATAGGGGCCTTCGTTATTGc                                      >  1:112779/1‑100 (MQ=255)
                                                            aTCCCCTGCCCATCCAAACACCACTGATTGACAATTTATCCATGATTAGACTACAGTGGGCCGGTTACTTTAATCGCATGGGCGCCTTCTTTATTGCGTa                                   <  1:1531849/100‑1 (MQ=255)
                                                               cccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAAtt                                >  1:1461590/1‑100 (MQ=255)
                                                                ccTGCCCATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTg                               <  1:297906/100‑1 (MQ=255)
                                                                    cccATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAg                           <  2:295037/100‑1 (MQ=255)
                                                                    cccATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAg                           <  1:436352/100‑1 (MQ=255)
                                                                     ccATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGa                          <  1:1191530/100‑1 (MQ=255)
                                                                      cATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGac                         <  1:112358/100‑1 (MQ=255)
                                                                      cATCGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGac                         <  1:29711/100‑1 (MQ=255)
                                                                         cGAGACACCACTGATGGACAATTTATCCATGATTAGACGACAGGGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGaca                        >  2:301998/1‑98 (MQ=255)
                                                                          gagaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCa                     >  2:844324/1‑100 (MQ=255)
                                                                            gaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaaa                   >  1:1692262/1‑100 (MQ=255)
                                                                            gaCACCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaaa                   >  1:774109/1‑100 (MQ=255)
                                                                             acacCACTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCaaaa                  <  2:1303855/100‑1 (MQ=255)
                                                                              caccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAg                 <  1:1209589/100‑1 (MQ=255)
                                                                               accacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAgg                <  2:1749091/100‑1 (MQ=255)
                                                                                ccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGga               >  1:1383799/1‑100 (MQ=255)
                                                                                ccacTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTACCAGACACCaaaa                  >  1:137663/1‑97 (MQ=255)
                                                                                 cacTGATGGACAATTTCTCCATGATTAGACGCCAGTGGGCCGGTGACTTTAAGCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgat              <  1:1594896/100‑1 (MQ=255)
                                                                                  acTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgatg             <  2:292870/100‑1 (MQ=255)
                                                                                   cTGATGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGgatg             >  1:666845/1‑99 (MQ=255)
                                                                                      aTGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcg         <  1:844323/100‑1 (MQ=255)
                                                                                       tGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgc        >  1:1912765/1‑100 (MQ=255)
                                                                                       tGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgc        >  2:66267/1‑100 (MQ=255)
                                                                                       tGGACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATcgc        <  2:1574671/100‑1 (MQ=255)
                                                                                         gACAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCgcgc      >  2:708695/1‑100 (MQ=255)
                                                                                           cAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc    >  2:182845/1‑100 (MQ=255)
                                                                                           cAATTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGctc    >  1:1868108/1‑100 (MQ=255)
                                                                                             aTTTATCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGAGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCtt  >  1:662652/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ACTTGCACTAAATGACCACTGGCAATGTTTTCGCTAACATCCCACCAGGAACGCAGAGCAATCCCCTGCCCATCGAGACACCACTGATGGACAATTTCTCCATGATTAGACGACAGTGGGCCGGTGACTTTAATCGCATGGGGGCCTTCTTTATTGCGTAATTGCCAGACACCAAAAGGATGATCGCGCTCTT  >  NC_000913/1881031‑1881223

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