Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 13.0% intergenic (‑150/+17) ygfI ← / ← yggE putative LysR‑type DNA‑binding transcriptional regulator YgfI/uncharacterized protein YggE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3230AT13.0% 140.3 / 12.3 54intergenic (‑150/+17)ygfI/yggEputative LysR‑type DNA‑binding transcriptional regulator YgfI/uncharacterized protein YggE
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base T (4/3);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

CCCGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGC  >  NC_000913/3067182‑3067471
                                                                                                                                             |                                                                                                                                                    
cccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                                                                              <  2:396807/150‑1 (MQ=255)
 ccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                                                                             <  2:1000220/150‑1 (MQ=255)
 ccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCt                                                                                                                                             >  1:1396044/1‑150 (MQ=255)
 ccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCAGTAGCCCGACGTATTCGTCGGGCt                                                                                                                                             <  1:1907302/150‑1 (MQ=255)
   gCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCTaa                                                                                                                                           <  1:955428/150‑2 (MQ=255)
               gCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAATATAGGTGATTACGTCATTGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgt                                                                                                                               <  2:1045130/150‑1 (MQ=255)
                 ttAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                                                                                                             <  1:145541/150‑1 (MQ=255)
                      cGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAg                                                                                                                        >  2:182517/1‑150 (MQ=255)
                       gTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAgg                                                                                                                       >  1:511189/1‑150 (MQ=255)
                            agtcagCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGttt                                                                                                                  >  1:1831169/1‑150 (MQ=255)
                                 gCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                                                                             >  1:1023583/1‑150 (MQ=255)
                                 gCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                                                                             >  1:1016737/1‑150 (MQ=255)
                                   cTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                                                                           >  2:564925/1‑150 (MQ=255)
                                         ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                                                                     <  1:1526859/150‑1 (MQ=255)
                                         ttGATGAAACCGCGAGCAAAGATAGGTGATAACGTCATGGTTTTACAGAAAATTACAGAAAAAGGCCGCAATATCGGGTACCGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                                                                     <  2:588065/150‑1 (MQ=255)
                                           gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                                                                   >  2:434636/1‑150 (MQ=255)
                                            attaaaCCGCCAGCAAACTAAGGTCATTACGTCATGGTTTAACATAAAATAACCCAAAACGGAGCAAATATCGGGTAACGGCTTACGCCCGACTATTACGTCGGGCTACAAATATATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                                                                                 <  2:222811/147‑1 (MQ=255)
                                                      gAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCt                                                                                        <  2:1885072/150‑1 (MQ=255)
                                                       aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAAGATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGGTCTa                                                                                       >  2:1728593/1‑150 (MQ=255)
                                                          aaaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACt                                                                                    >  2:671480/1‑150 (MQ=255)
                                                           aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                                                                   <  1:1810413/150‑1 (MQ=255)
                                                           aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                                                                   <  1:564925/150‑1 (MQ=255)
                                                               tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                                                                               <  1:671480/150‑1 (MQ=255)
                                                               tAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGAAGAAAACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGTAAg                                                                               >  1:1902128/1‑150 (MQ=255)
                                                                aGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                                                                              <  1:182517/150‑1 (MQ=255)
                                                                 ggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                                                                             >  2:351226/1‑150 (MQ=255)
                                                                 ggTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGAc                                                                             >  1:1299321/1‑150 (MQ=255)
                                                                                      cataaaaatacataaaaATGAGCCAATATCAGGTAAAGGCATTAAACCGAAGAATAAGTCGGTCTACAAATATTATTGAACTAAAAGTGTTTTAGAGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa                                                        >  2:1650974/4‑150 (MQ=255)
                                                                                        gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg                                                                 >  1:1388138/1‑139 (MQ=255)
                                                                                        gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg                                                                 <  2:1388138/139‑1 (MQ=255)
                                                                                        gaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACt                                                      >  1:936991/1‑150 (MQ=255)
                                                                                              tacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAg                                                >  2:970723/1‑150 (MQ=255)
                                                                                                cagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCg                                              >  2:1064592/1‑150 (MQ=255)
                                                                                                 agaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCgg                                             <  2:511189/150‑1 (MQ=255)
                                                                                                  gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGc                                            >  2:1472932/1‑150 (MQ=255)
                                                                                                       aggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCt                                       >  1:1185868/1‑150 (MQ=255)
                                                                                                         gaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCg                                     >  1:1762867/1‑150 (MQ=255)
                                                                                                                    gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTg                          >  1:1459760/1‑150 (MQ=255)
                                                                                                                        aaaGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCg                      >  1:1399942/1‑150 (MQ=255)
                                                                                                                          aGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGa                    >  1:267510/1‑150 (MQ=255)
                                                                                                                           ggCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGac                   >  2:626107/1‑150 (MQ=255)
                                                                                                                           ggCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGac                   <  2:1396044/150‑1 (MQ=255)
                                                                                                                            gCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGaca                  <  2:1831169/150‑1 (MQ=255)
                                                                                                                             cATTAGCCCGCCGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                 >  1:1882507/1‑150 (MQ=255)
                                                                                                                             cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                 <  2:1299321/150‑1 (MQ=255)
                                                                                                                             cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                 >  2:1184661/1‑150 (MQ=255)
                                                                                                                              aTTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACAcc                <  2:1762867/150‑1 (MQ=255)
                                                                                                                                tAGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCgg              >  1:1337342/1‑150 (MQ=255)
                                                                                                                                 aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc             >  1:191840/1‑150 (MQ=255)
                                                                                                                                 aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc             >  1:198775/1‑150 (MQ=255)
                                                                                                                                 aGCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc             >  1:579217/1‑150 (MQ=255)
                                                                                                                                 aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc             >  1:964924/1‑150 (MQ=255)
                                                                                                                                     cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGATGTTGAGATCGACCTGATCATCAAACTGTATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACTCCGGCGCgg         <  2:1144811/150‑1 (MQ=255)
                                                                                                                                            tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGc  <  2:267510/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCCGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGGC  >  NC_000913/3067182‑3067471

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