Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,306 G→C 8.8% I13M (ATC→ATG)  yhfZ ← putative DNA‑binding transcriptional regulator YhfZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,3060GC 8.8% 155.8 / 5.3 54I13M (ATC→ATG) yhfZputative DNA‑binding transcriptional regulator YhfZ
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base C (3/2);  total (29/25)
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.93e-01

CGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGG  >  NC_000913/3512166‑3512444
                                                                                                                                            |                                                                                                                                          
cGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtgat                                                                                                                                          <  1:749287/143‑1 (MQ=255)
     cGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgacga                                                                                                                                     <  1:680807/143‑1 (MQ=255)
       cTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGacgacg                                                                                                                                   >  2:441055/1‑143 (MQ=255)
       cTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGacgacg                                                                                                                                   >  1:528423/1‑143 (MQ=255)
          cgcgTGATTCCAGCGTTTGCAGCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctg                                                                                                                                <  1:784044/143‑3 (MQ=255)
              tGATTCCAGCGTTTTCAACGCGGCCTGCGTCGGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttctt                                                                                                                            <  1:545948/143‑1 (MQ=255)
                 ttCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                         <  2:29955/143‑1 (MQ=255)
                 ttCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACAGCCCTTTTCACCCAACAAATAACGGGCCAGCGTCCTGATGACGACGCCTGCttttt                                                                                                                         >  2:108312/1‑143 (MQ=255)
                  tCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                        <  1:644747/143‑1 (MQ=255)
                     aGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGAt                                                                                                                     >  2:290667/1‑143 (MQ=255)
                     aGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACCCATTCTTTTTTCAt                                                                                                                     >  2:298731/1‑143 (MQ=255)
                        gTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                                                                                  >  1:264443/1‑143 (MQ=255)
                               tcgcgGCCTGCGTCATGCCATCGGCTGAACGGCATTCATTTGCCAGCGCATCTATGGTTTTCAATCGATTGCCGCACCTTTCTCCCTACTTATAACGGTCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGa                                                                                                           <  2:12751/142‑1 (MQ=255)
                               aCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATCTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGa                                                                                                           <  1:256876/143‑1 (MQ=255)
                               aCGCGGCCCGAGTCTGGCCAACGGATGAACGGCATTCAGTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACGAACGACGGACCAGCGTCGTGGTGACGACGCCTTCTTTTTTGATAAACGTTCGa                                                                                                           <  1:225589/143‑1 (MQ=255)
                                    gCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                      <  1:238741/143‑1 (MQ=255)
                                              ggCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAtttt                                                                                            <  1:91723/143‑1 (MQ=255)
                                               tccAACGGATGTACGGCATGCATTTGCCAGCGCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTc                                                                                           <  2:725736/142‑1 (MQ=255)
                                                  aaCGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGCTAAACGTTCGACGAATAATACAtttcacgt                                                                                        >  2:252084/1‑138 (MQ=255)
                                                       aTGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                   <  2:240140/143‑1 (MQ=255)
                                                       aTGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                   >  1:659991/1‑143 (MQ=255)
                                                         gAACGTCATTCATTTGCCAGCACACCTATGGTTTTCGATCGATTACCGCACTTTTCACCCAACAGATAACGGGCCAGCGTCGGCATGCCGATACCTTCTTTTTGGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGa                                                                                 <  2:360199/143‑1 (MQ=255)
                                                             ggCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGCTAACCGTTCGACGCATAATCAATTTTCAGTAAATTCAatat                                                                             >  2:363195/1‑143 (MQ=255)
                                                              gCATTTATTTGCGTGCTCATCTATGGTTTTCAATCGATTACCGCACGTTTCCCCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAtt                                                                            <  1:286408/143‑1 (MQ=255)
                                                                   cattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTatat                                                                       <  1:73296/143‑1 (MQ=255)
                                                                       tGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTc                                                                   <  1:78653/143‑1 (MQ=255)
                                                                            gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTAAGTAAATTGAATATTTATATCTTCAGGaa                                                              >  2:25403/1‑143 (MQ=255)
                                                                                aTCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTTTTTATTGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATGTTCAGGAATTTg                                                          >  2:200563/1‑143 (MQ=255)
                                                                                  cTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                        <  1:609461/143‑1 (MQ=255)
                                                                                  cTATGCTTTTCAACCGATTACCGCTCTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                        <  1:67139/143‑1 (MQ=255)
                                                                                     tGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTGAGGATTGTTATTATTTATATCTTCAGGAATTTGAAGAt                                                     >  1:545874/1‑143 (MQ=255)
                                                                                      ggTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATa                                                    >  1:11991/1‑143 (MQ=255)
                                                                                        tttaCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGACGATAAg                                                  >  1:246074/1‑143 (MQ=255)
                                                                                               tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATCATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAAc                                           >  1:105815/1‑143 (MQ=255)
                                                                                               tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAAc                                           >  1:328260/1‑143 (MQ=255)
                                                                                                cGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACa                                          >  1:466440/1‑143 (MQ=255)
                                                                                                        gCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                  >  2:561400/1‑143 (MQ=255)
                                                                                                          aCTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGGTATACGGGCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgtg                                <  1:587801/143‑1 (MQ=255)
                                                                                                           cTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTgg                               >  1:643189/1‑143 (MQ=255)
                                                                                                            ttttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCGTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGa                              >  1:58890/1‑143 (MQ=255)
                                                                                                                cACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACg                          >  2:580012/1‑143 (MQ=255)
                                                                                                                 aCCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACgg                         <  1:739897/143‑1 (MQ=255)
                                                                                                                    caacatATAACGGGCCCGCGTCGTGATGACGGCGCCTTCTTTTTTGATAAACGTTGTACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGa                      <  2:92428/143‑1 (MQ=255)
                                                                                                                     aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATCTCTTCAGGACTTTGAAGCTACGGCGGAACAGCGAGACGTGAATACGACGAg                     >  2:760844/1‑143 (MQ=255)
                                                                                                                         aaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAGACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCAc                 >  2:512916/1‑143 (MQ=255)
                                                                                                                          atgaCGGGCCAGCGTCGTGATGACGGCGCCTTCTTTTTTGATAAACGTTCGACGCATAAGAAATGTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACt                <  2:264717/140‑1 (MQ=255)
                                                                                                                            aaCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTc              <  1:434132/143‑1 (MQ=255)
                                                                                                                               gtgCCGGCGTCGTGATGACGAGGCCGGCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTcaca           <  2:72144/141‑1 (MQ=255)
                                                                                                                                ggCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGAAAAGTCTCAACAGCGAGATGTGGAAACGGCGAGGAACttaccat          >  2:318324/1‑138 (MQ=255)
                                                                                                                                    ggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAttttt      >  2:767613/2‑143 (MQ=255)
                                                                                                                                    ggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtcttc      >  2:208762/2‑142 (MQ=255)
                                                                                                                                    ggCGTCGTCATCACGACGCCTTCTTTATTGATAAGCGTTCGGCGCAGAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACGTCACAttttt      >  2:612709/2‑143 (MQ=255)
                                                                                                                                    aGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtgttt      >  1:533699/1‑143 (MQ=255)
                                                                                                                                      cgtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTc    <  1:55017/143‑1 (MQ=255)
                                                                                                                                        tcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTCTTCCgg  >  2:146242/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGG  >  NC_000913/3512166‑3512444

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