Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,204,350 A→G Q344R (CAG→CGG)  zraR → fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,204,3500AG85.2% 142.5 / ‑0.6 54Q344R (CAG→CGG) zraRfused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/23);  minor base C (8/0);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATCC  >  NC_000913/4204201‑4204493
                                                                                                                                                          |                                                                                                                                               
tCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCa                                                                                                                                                 <  2:862367/150‑1 (MQ=255)
  gCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGc                                                                                                                                              <  2:831023/151‑1 (MQ=255)
    cAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATTCCTCTGCTGGCTGGCCATTTTTTGCAGCGCTTTGCCGAGCGGAATCGCCAGGCGGGAAAAAGGTTTATGCCCCCGGCg                                                                                                                                             >  1:30519/1‑150 (MQ=255)
      ggATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATg                                                                                                                                          <  1:487561/151‑1 (MQ=255)
          ctctACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGAAATTCCTCTGCTTGCTGGCCATTTTTTGCAGCGCCTTTCCGAGGGGAATCGCAAGGGGGGGAAAAGGTTTTAGCCCCCGGGGGAgggt                                                                                                                                        >  2:552935/1‑149 (MQ=255)
            ctactaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGa                                                                                                                                        <  1:703285/147‑1 (MQ=255)
            ctactaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATc                                                                                                                                      <  2:22862/149‑1 (MQ=255)
            ctactaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGAAATCCCCCTTCTGGGCGGCCACTTTTTTCCGCGCTTTTGCGGGGGGTATCTCAAGGAGGGGAAAAGGTTTTCTCCCCCGGCg                                                                                                                                             >  2:1001042/1‑142 (MQ=255)
              actaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGGTGGCCATTTTTTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTggt                                                                                                                                  >  1:440416/1‑149 (MQ=255)
              actaTCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCt                                                                                                                                     <  1:1393/148‑1 (MQ=255)
                  tCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGGTGGGCCTTTTCTGCAGCGCGTTTCCGCGCGGAAACCCAAAGCGGGAAAAGGGTTTACGCCCCGGGCGATGGGTCTGTTGAtt                                                                                                                              >  2:251653/1‑151 (MQ=255)
                   cGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTc                                                                                                                             <  1:271312/151‑1 (MQ=255)
                    gCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCa                                                                                                                            <  1:40073/151‑1 (MQ=255)
                     ccTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGGTATTCCTCTGCTGGGTGGCCCATTTTTGCCGCGCCTTTCCGGGCGGAATCGCAAAGCGGGCAAAGGGTTTTCGCCCCCGGCGGTGGGTCTGGTGATTCAt                                                                                                                           >  1:763389/1‑151 (MQ=255)
                        gAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTa                                                                                                                         <  1:124266/150‑1 (MQ=255)
                        gAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCCCTGCTGGGTGGCCATTTTCTTCAGCGCTTTTGCGAGCGTAATCGCCAGGCGGTAAAAGGGTTTTCGCCCCCGGCGATGGAGCTGTTGTTTCAtt                                                                                                                          >  2:854049/1‑149 (MQ=255)
                        gAATGTGGTGGCGATTGAAGTACCATCGCTGCGGAAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTAc                                                                                                                        <  1:141584/151‑1 (MQ=255)
                            gtggtgGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGAAATCACCCCGCCGGCTGGCCATTTTTTTCAAGGCTTTTGCCAGCGTAAACGCAAAGGGGGAAAAAGGTTTTAGCCCCCGGGGAGGGGTCTGTTTGTTCAtt                                                                                                                          >  1:377898/1‑145 (MQ=255)
                              ggtggCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAAGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTg                                                                                                                   <  1:985678/150‑1 (MQ=255)
                                  gCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATTCCCCCGCTTGCTGGGCCCTTTTTTCCGCGCTTTGGCCGAGGGAATCCCAAAGCGGGAAAAGGGTTTAAGCCCCGGGCGGTGGGTCTTTTGTTTTAATTCGACTTGCCgg                                                                                                              >  2:552569/1‑151 (MQ=255)
                                    gATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGa                                                                                                            <  1:303584/151‑1 (MQ=255)
                                      ttGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCCTTTTTTTCACCGCTTTGCCGAGCGTAATCGCAAAGCGGGAAAAGGTTTTTCGCCCCCGGCGGTGGGTCTGTTGTTTCATTTAGACTGGCCgggga                                                                                                           >  1:605954/1‑148 (MQ=255)
                                          aGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATAtt                                                                                                      <  2:841144/151‑1 (MQ=255)
                                          aGTACCATCGCTGCGGCAACGGCGGAAAGATATACCTCCTGCGGCTGGCCATTCTCTTCTGCACCTTGGCCACCCCAATCCGGACGCGGGGAAAGGATTTTCGCCCCCGGCGGTGGTTCTTTTTGTTTATTACGAcacgcc                                                                                                                >  2:384046/1‑136 (MQ=255)
                                               cATCGCTGCGGCAACGGCGGGAAGACATCCCTCCGCTTGCTGGCCCTTTTTTTCACCGCTTTTCCGAGCGGAAACCCAAGGCGGGAAAAGGGTTTATGCCCCGGGGGATGGGTCTGTTTATTCATTACGACTGGCCGGGGAATATTCGTg                                                                                                  >  1:430443/1‑150 (MQ=255)
                                                aTCGCTGCGGCAACGGCGGGGAGAGAACCCTCTTCCGGCTGGCCCATTTTTGCAGCGGCTTGCCGAGAGGAAACGCAAAGGGGGAAAAGGGTTTACGCCCCGGGCGGTGGGGCTGTTGATTCATTTAGAACGGCCGGGGAATATTtgtggg                                                                                                >  2:96128/1‑149 (MQ=255)
                                                   gCTGCGGCAACGGCGGGAAGAAATCCCCCTCCAGGCCGGCCCTTTTCTGCCGACCTCTGTCCGGCGCAAAACCAAAGCGGGAAAAAGTTTTTAGCCCCCGGCGAGGGGTCTGTTGTTTCATTACTAACGGCCCGGAAAGATTTGGGAGc                                                                                               >  2:972412/1‑149 (MQ=255)
                                                     tGCGGCAACGGCGGGAAGATATCCCTCTGCTGGGTGGCCAATTTTTGCAGCGCTTTGCCGGGCGGAAACGCAAAGCGGGAAAAGGGTTTACGCCCCGGGCGATGGATCTGTTGGTTTATTTAGACTGGGCGGGGAATAATTGTGGGCTgg                                                                                            >  2:358393/1‑150 (MQ=255)
                                                       cGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGAc                                                                                                                     <  1:414407/123‑1 (MQ=255)
                                                        ggCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGaaaa                                                                                        <  1:131787/151‑1 (MQ=255)
                                                             cggcggGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCgg                                                                                   >  1:94605/1‑151 (MQ=255)
                                                                tgggAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTgg                                                                                >  2:116350/2‑151 (MQ=255)
                                                                       atatCCCTCTCCTGGCTGGCCATTTTCTGCACCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGCAATATTCGTGAGCTGGAAACCGCGGTGGAACggg                                                                          >  2:1049548/1‑150 (MQ=255)
                                                                        tatCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGc                                                                         >  2:968249/1‑150 (MQ=255)
                                                                         atCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAg                                                                       <  1:879052/151‑1 (MQ=255)
                                                                                     ggctggcCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTgctgct                                                             <  1:738652/150‑1 (MQ=255)
                                                                                     ggctggcCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGctgctg                                                            >  1:958418/1‑151 (MQ=255)
                                                                                       ctggcCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGAc                                                          >  1:941393/1‑151 (MQ=255)
                                                                                         agcgatttgatgaaGCGCGTTGCCGAGCGGAATCACAAGGCGGGAAACGGCTTTACGCCCAGGGCGATGGATCTGCTGATTCACTACGACTGGCCGGGAAATATTCGCGAGCTGGAAAACGCGGTGGAATGGACAGTGGGGCTGCTGACCg                                                        >  2:1048366/14‑151 (MQ=255)
                                                                                            cATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCgggg                                                     <  1:94445/151‑1 (MQ=255)
                                                                                               tttCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAat                                                  >  1:233553/1‑151 (MQ=255)
                                                                                                ttCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAat                                                  <  2:263097/150‑1 (MQ=255)
                                                                                                  cTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatat                                                >  2:232788/1‑150 (MQ=255)
                                                                                                   tGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatatat                                              >  1:984900/1‑151 (MQ=255)
                                                                                                               gCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAcgc                                    <  2:877814/150‑1 (MQ=255)
                                                                                                                gcGAGCGTAATCGCAAGGCGGTGAAAGGGTTTACGCCCCGGGCGATGGATCTGTGGTTTCATTACGACTGGCAGGGAAATATTCGTGAGCTGGAAAAACCGGTGGAACGGGCAGTGGGGCTGCTGACCGGGGAATATATTTCCGAAcgc                                    >  2:217603/2‑149 (MQ=255)
                                                                                                                ccGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGa                                  <  1:461906/151‑1 (MQ=255)
                                                                                                                    gCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCt                               <  2:457859/150‑1 (MQ=255)
                                                                                                                    gCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTg                              <  1:700773/151‑1 (MQ=255)
                                                                                                                       tAATCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCg                           >  2:560354/1‑151 (MQ=255)
                                                                                                                         aTCGCAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCt                         <  1:25770/151‑1 (MQ=255)
                                                                                                                             cAAGGCGGTAAAAGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCg                     >  1:738910/1‑151 (MQ=255)
                                                                                                                                         aaaaGGTTTTACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGt            >  1:552816/1‑151 (MQ=255)
                                                                                                                                                  tACGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATc   >  2:459923/1‑151 (MQ=255)
                                                                                                                                                   aCGCCCCGGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATcc  >  1:429233/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                               
TCGCCAGGATCTCTACTATCGCCTGAATGTGGTGGCGATTGAAGTACCATCGCTGCGGCAACGGCGGGAAGATATCCCTCTGCTGGCTGGCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGCTGGCGATTGCCAGTACGCCGATCC  >  NC_000913/4204201‑4204493

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