Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,299:1 +G 6.0% coding (875/897 nt) ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,2991.G 6.0% 302.9 / 13.1 67coding (875/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (30/33);  new base G (2/2);  total (32/35)
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.11e-01

CATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGC  >  NC_000913/3066152‑3066432
                                                                                                                                                    |                                                                                                                                      
caTCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAt                                                                                                                                         >  1:732539/1‑148 (MQ=255)
  tCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCt                                                                                                                                      >  2:6167/1‑148 (MQ=255)
               gATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                        <  1:971943/148‑1 (MQ=255)
                 ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCATATTGCTGa                                                                                                                      <  2:732539/148‑1 (MQ=255)
                   aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAg                                                                                                                    <  1:6167/148‑1 (MQ=255)
                     cTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                        <  4:455627/142‑1 (MQ=255)
                     cTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                        >  3:455627/1‑142 (MQ=255)
                       gCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAgatc                                                                                                                 >  4:113502/1‑145 (MQ=255)
                           cGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTga                                                                                                            <  3:113502/148‑2 (MQ=38)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAt                                                                                                                                         <  2:178756/113‑1 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAt                                                                                                                                         >  1:178756/1‑113 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                 <  4:287493/135‑1 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                 >  3:287493/1‑135 (MQ=255)
                                     aaaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGa                                                                                                                      >  1:448205/1‑128 (MQ=255)
                                     aaaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGa                                                                                                                      <  2:448205/128‑1 (MQ=255)
                                      aaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                 <  3:395241/148‑1 (MQ=255)
                                       aaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATg                                                                                                                          <  1:667763/122‑1 (MQ=255)
                                       aaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATg                                                                                                                          >  2:667763/1‑122 (MQ=255)
                                                    gACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                 <  4:887255/118‑1 (MQ=255)
                                                    gACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                 >  3:887255/1‑118 (MQ=255)
                                                           tCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                 >  1:771076/1‑127 (MQ=255)
                                                           tCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                 <  2:771076/127‑1 (MQ=255)
                                                           tCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                          <  2:54874/134‑1 (MQ=255)
                                                           tCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                          >  1:54874/1‑134 (MQ=255)
                                                                acGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                                   <  3:620460/136‑1 (MQ=255)
                                                                acGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                                   >  4:620460/1‑136 (MQ=255)
                                                                     cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGa                                                                                                                      <  2:395660/96‑1 (MQ=255)
                                                                     cTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGa                                                                                                                      >  1:395660/1‑96 (MQ=255)
                                                                       tAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                   <  4:385243/82‑1 (MQ=255)
                                                                       tAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                   >  3:385243/1‑82 (MQ=255)
                                                                           cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                              <  4:176458/130‑1 (MQ=255)
                                                                           cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                              >  3:176458/1‑130 (MQ=255)
                                                                             cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCTTGGCGGtata                                                                                                     <  3:917557/105‑1 (MQ=255)
                                                                             cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtata                                                                                                     >  4:917557/1‑105 (MQ=255)
                                                                                  tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATGCTTCAGCATATTGCTGAAg                                                                                                                    >  1:12284/1‑86 (MQ=37)
                                                                                  tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATGCTTCAGCATATTGCTGAAg                                                                                                                    <  2:12284/86‑1 (MQ=37)
                                                                                              gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatgat                                                                                             >  4:73215/1‑96 (MQ=255)
                                                                                              gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatgat                                                                                             <  3:73215/96‑1 (MQ=255)
                                                                                              gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                                        <  2:169100/101‑1 (MQ=255)
                                                                                              gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                                        >  1:169100/1‑101 (MQ=255)
                                                                                                aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                               >  2:669309/1‑108 (MQ=255)
                                                                                                aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                               <  1:669309/108‑1 (MQ=255)
                                                                                                  tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                                      >  4:588337/1‑115 (MQ=255)
                                                                                                  tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                                      <  3:588337/115‑1 (MQ=255)
                                                                                                    taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                     <  4:752280/66‑1 (MQ=255)
                                                                                                    taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                     >  3:752280/1‑66 (MQ=255)
                                                                                                            aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                                   >  3:1024673/1‑92 (MQ=255)
                                                                                                            aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                                   <  4:1024673/92‑1 (MQ=255)
                                                                                                                            gTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAata                                   >  4:525080/1‑124 (MQ=255)
                                                                                                                            gTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATAGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAata                                   <  3:525080/124‑1 (MQ=255)
                                                                                                                             ttAATGACTCGTCCCGTCACGATCTTCAGCAATATg                                                                                                                          >  3:682989/1‑36 (MQ=255)
                                                                                                                             ttAATGACTCGTCCCGTCACGATCTTCAGCAATATg                                                                                                                          <  4:682989/36‑1 (MQ=255)
                                                                                                                                  gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCa     >  3:105336/1‑148 (MQ=255)
                                                                                                                                    cTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                                >  2:832377/1‑103 (MQ=255)
                                                                                                                                    cTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                                <  1:832377/103‑1 (MQ=255)
                                                                                                                                     tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                  >  4:529855/1‑84 (MQ=255)
                                                                                                                                     tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                  <  3:529855/84‑1 (MQ=255)
                                                                                                                                     tCGACCCGTCTAGATCTACAGAATCATGCTGAAGCATCCTGGCGGTATAACGGTCATTTCTGAAAAAATAAAGATTTTCGTGCTTCTGCAACCGCCAGTCCGGTGCATCAAAATATTTACACGTTGCCGCATCCTCCAGTCGCCAgtc  <  4:105336/148‑3 (MQ=255)
                                                                                                                                      cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                                        <  4:1017789/61‑1 (MQ=255)
                                                                                                                                      cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTg                                                                                        >  3:1017789/1‑61 (MQ=255)
                                                                                                                                      cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                                    >  4:73448/1‑113 (MQ=255)
                                                                                                                                      cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                                    <  3:73448/113‑1 (MQ=255)
                                                                                                                                        tCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                               <  2:1100203/68‑1 (MQ=255)
                                                                                                                                        tCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                               >  1:1100203/1‑68 (MQ=255)
                                                                                                                                        tCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                               >  4:1080497/1‑116 (MQ=255)
                                                                                                                                        tCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                               <  3:1080497/116‑1 (MQ=255)
                                                                                                                                              cACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGc       <  1:1035263/134‑1 (MQ=255)
                                                                                                                                              cACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGc       >  2:1035263/1‑134 (MQ=255)
                                                                                                                                                  aTCTTCTGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACAt                                                           <  4:686539/78‑1 (MQ=255)
                                                                                                                                                  aTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACAt                                                           >  3:686539/1‑78 (MQ=255)
                                                                                                                                                   tCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                                <  4:1070842/88‑1 (MQ=255)
                                                                                                                                                   tCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                                >  3:1070842/1‑88 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
CATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGC  >  NC_000913/3066152‑3066432

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