Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,031,221 T→A 9.5% intergenic (‑290/‑146) ygfT ← / → uacT putative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit/uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,031,2210TA 9.5% 172.2 / 6.2 52intergenic (‑290/‑146)ygfT/uacTputative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit/uric acid permease
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base A (3/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGA  >  NC_000913/3031073‑3031367
                                                                                                                                                    |                                                                                                                                                  
ttattTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGttt                                                                                                                                                    >  2:265656/1‑149 (MQ=255)
  attTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTaa                                                                                                                                                  >  1:103730/1‑149 (MQ=255)
          aTTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAtt                                                                                                                                          >  4:540777/1‑149 (MQ=255)
             gCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTAtg                                                                                                                                       <  3:68203/149‑1 (MQ=255)
              ccGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTAtgt                                                                                                                                      <  4:227900/149‑1 (MQ=255)
                gATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATCCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGa                                                                                                                                    <  1:456423/149‑1 (MQ=255)
                 aTCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGAt                                                                                                                                   >  4:90347/1‑149 (MQ=255)
                  tCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGAtt                                                                                                                                  <  1:265656/149‑1 (MQ=255)
                       agagCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTc                                                                                                                             >  3:125362/1‑149 (MQ=255)
                       agagCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTc                                                                                                                             >  2:481120/1‑149 (MQ=255)
                            gcTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCaa                                                                                                                        >  4:240896/1‑149 (MQ=255)
                            gcTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCaa                                                                                                                        <  1:481120/149‑1 (MQ=255)
                                   aGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGt                                                                                                                 <  3:33756/149‑1 (MQ=255)
                                   aGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGt                                                                                                                 <  3:43351/149‑1 (MQ=255)
                                    gATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                                <  1:584908/149‑1 (MQ=255)
                                    gATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                                >  4:408409/1‑149 (MQ=255)
                                         aaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTaa                                                                                                           <  1:25868/149‑1 (MQ=255)
                                         aaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTaa                                                                                                           <  1:162740/149‑1 (MQ=255)
                                            cGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTc                                                                                                        <  4:130888/149‑1 (MQ=255)
                                             gATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGATATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTcc                                                                                                       <  1:458312/149‑1 (MQ=255)
                                             gATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTcc                                                                                                       >  2:283098/1‑149 (MQ=255)
                                                  aaTGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTc                                                                                                  >  3:581753/1‑149 (MQ=255)
                                                     gCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTc                                                                                               <  2:213551/149‑1 (MQ=255)
                                                       gagaGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGc                                                                                             >  3:544635/1‑149 (MQ=255)
                                                         gagTGTGACATCCCAGAATGATTAATTAAATATCAATATGAGAGAATTACAGATAATTCCATTGTCAGGTTCTCTACTGATTTTCGTGGTTTAACTACGGTTATGTGATTGCCTCCGCACTTCCTGTCTCTAACTCCCCTTCCTCCCaa                                                                                           >  4:342360/1‑149 (MQ=255)
                                                               gACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                     <  3:27062/149‑1 (MQ=255)
                                                               gACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                                     <  3:408409/149‑1 (MQ=255)
                                                                              ttAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAt                                                                      <  3:430171/149‑1 (MQ=255)
                                                                                    aaaTATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                                >  2:118487/1‑149 (MQ=255)
                                                                                      atatCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCGTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                              >  4:48025/1‑149 (MQ=255)
                                                                                       tatCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGa                                                             <  1:569708/149‑1 (MQ=255)
                                                                                          cAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGGGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAg                                                          >  4:508173/1‑149 (MQ=255)
                                                                                                    aaaTTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAAc                                                <  3:511714/149‑1 (MQ=255)
                                                                                                       ttACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTa                                             >  3:271165/1‑149 (MQ=255)
                                                                                                              tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGtt                                      <  2:104943/149‑1 (MQ=255)
                                                                                                                  tCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                                  >  3:342270/1‑149 (MQ=255)
                                                                                                                     aTTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCttttt                               >  4:365723/1‑149 (MQ=255)
                                                                                                                         tCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTActc                           >  4:515193/1‑149 (MQ=255)
                                                                                                                              ttctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTggcg                      >  4:128845/1‑149 (MQ=255)
                                                                                                                               tctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTggcgg                     <  3:508173/149‑1 (MQ=255)
                                                                                                                                      tcATTTTCGTGGTTTAACCACGATTATGTGATTGCTACCGCAATTCCTGCCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAg              <  3:342360/149‑1 (MQ=255)
                                                                                                                                        aTTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGga            <  4:125362/149‑1 (MQ=255)
                                                                                                                                         ttttCGTGGTCTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgag           >  4:402801/1‑149 (MQ=255)
                                                                                                                                          aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgt                                                                    <  2:250957/89‑1 (MQ=255)
                                                                                                                                          aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgt                                                                    >  1:250957/3‑91 (MQ=255)
                                                                                                                                          aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgaga          >  3:224533/3‑149 (MQ=255)
                                                                                                                                            tCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGaaa        >  4:232023/1‑149 (MQ=255)
                                                                                                                                            tCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGaaa        <  3:540777/149‑1 (MQ=255)
                                                                                                                                            tCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGaaa        <  3:240896/149‑1 (MQ=255)
                                                                                                                                            tCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGaaa        >  3:312845/1‑149 (MQ=255)
                                                                                                                                               tGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgaga          >  1:414186/1‑144 (MQ=255)
                                                                                                                                               tGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGgaga          <  2:414186/144‑1 (MQ=255)
                                                                                                                                                 gTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACtg   >  1:151883/1‑149 (MQ=255)
                                                                                                                                                  tttAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACtga  >  2:579868/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTATTTCAACATTGCCGATCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGAGAAACACTGA  >  NC_000913/3031073‑3031367

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