Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,031,221 T→A 13.1% 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,2210TA13.1% 180.2 / 18.4 61intergenic (‑290/‑146)ygfT/uacTputative oxidoreductase, Fe‑S subunit/nucleotide‑binding subunit/uric acid permease
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base A (4/4);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGA  >  NC_000913/3031091‑3031357
                                                                                                                                  |                                                                                                                                        
tCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGAt                                                                                                                         <  2:1072749/148‑1 (MQ=255)
     agagCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCt                                                                                                                    <  3:1033751/148‑1 (MQ=255)
                   aTAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtc                                                                                                      >  1:451899/1‑148 (MQ=255)
                    tAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                     <  2:451899/148‑1 (MQ=255)
                     aaaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTg                                                                                                        >  4:385559/1‑144 (MQ=255)
                     aaaaaCGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTg                                                                                                        <  3:385559/144‑1 (MQ=255)
                                              aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                           <  4:517275/148‑1 (MQ=255)
                                              aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaa                                                                                                                                 <  4:559093/94‑2 (MQ=255)
                                              aCATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTAAACCACGaa                                                                                                                                 >  3:559093/1‑93 (MQ=255)
                                                        aTGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACg                                                                 >  3:963920/1‑148 (MQ=255)
                                                                tGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                    <  4:837325/105‑1 (MQ=255)
                                                                tGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtctc                                                                                                    >  3:837325/1‑105 (MQ=255)
                                                                   aaTATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                      <  4:963920/148‑1 (MQ=255)
                                                                            atGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGAc                                                  <  3:9815/143‑1 (MQ=255)
                                                                            atGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGAc                                                  >  4:9815/1‑143 (MQ=255)
                                                                                         aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                     >  3:765079/1‑79 (MQ=255)
                                                                                         aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                     <  4:765079/79‑1 (MQ=255)
                                                                                         aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                    >  1:213542/1‑128 (MQ=255)
                                                                                         aGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                    <  2:213542/128‑1 (MQ=255)
                                                                                          gATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACg                               >  3:398934/1‑148 (MQ=255)
                                                                                            tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAt                                                                                                                                 <  1:790400/48‑1 (MQ=255)
                                                                                            tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGAt                                                                                                                                 >  2:790400/1‑48 (MQ=255)
                                                                                            tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                      <  4:839783/123‑1 (MQ=255)
                                                                                            tAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGttt                                                      >  3:839783/1‑123 (MQ=255)
                                                                                               ttCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGc                          <  4:398934/148‑1 (MQ=255)
                                                                                                tCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACg                               >  3:256571/1‑142 (MQ=255)
                                                                                                tCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACg                               <  4:256571/142‑1 (MQ=255)
                                                                                                   aTTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtttt                      >  2:825223/1‑148 (MQ=255)
                                                                                                    ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                     <  1:766468/68‑1 (MQ=255)
                                                                                                    ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGtct                                                                                                     >  2:766468/1‑68 (MQ=255)
                                                                                                    ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                           <  1:788475/94‑1 (MQ=255)
                                                                                                    ttGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                           >  2:788475/1‑94 (MQ=255)
                                                                                                        cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCa                                                                        <  2:987309/93‑1 (MQ=255)
                                                                                                        cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCa                                                                        >  1:987309/1‑93 (MQ=255)
                                                                                                        cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACg                                                                 >  4:1034654/1‑100 (MQ=255)
                                                                                                        cAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACg                                                                 <  3:1034654/100‑1 (MQ=255)
                                                                                                           gTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                           <  3:498495/87‑1 (MQ=255)
                                                                                                           gTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTg                                                                           >  4:498495/1‑87 (MQ=255)
                                                                                                            ttctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGtt                                              <  3:1003483/115‑1 (MQ=255)
                                                                                                            ttctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGtt                                              >  4:1003483/1‑115 (MQ=255)
                                                                                                             tctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTgt              >  4:211447/1‑145 (MQ=255)
                                                                                                             gctcttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTgg              <  3:211447/145‑1 (MQ=255)
                                                                                                                cttctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTggcgggc         >  2:112317/1‑148 (MQ=255)
                                                                                                                   ctcATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGaa      <  1:825223/148‑1 (MQ=255)
                                                                                                                     cATTTTCGTTTTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTccc                                                                                            <  4:3733/60‑1 (MQ=255)
                                                                                                                     cATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAg    >  1:870886/1‑148 (MQ=255)
                                                                                                                     cATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTATAACTccc                                                                                            >  3:3733/1‑60 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCGCGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGt                             <  2:840783/121‑1 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCGCGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGt                             >  1:840783/1‑121 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTc                                            <  1:835556/106‑1 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTc                                            >  2:835556/1‑106 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                        <  1:489326/126‑1 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                        >  2:489326/1‑126 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGga  <  1:112317/148‑1 (MQ=255)
                                                                                                                       ttttCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGga  >  3:1016555/1‑148 (MQ=255)
                                                                                                                        aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                        >  1:733581/3‑125 (MQ=255)
                                                                                                                        aatCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCtt                        <  2:733581/123‑1 (MQ=255)
                                                                                                                         atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                         >  1:694589/2‑91 (MQ=255)
                                                                                                                         atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                         <  2:694589/90‑1 (MQ=255)
                                                                                                                         atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                         >  1:453284/2‑91 (MQ=255)
                                                                                                                         atCGTGGTTAAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCAtgtg                                                         <  2:453284/90‑1 (MQ=255)
                                                                                                                                 ttAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                    >  4:1020508/1‑88 (MQ=255)
                                                                                                                                 ttAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAg                                                    <  3:1020508/88‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
TCAAAAGAGCGCTATCCAGATAAAAACGATTAAATGCGAGAGTGCGACATGCCAGAATGATTAATGAAATATCAATATGAGAAAATTACAGATAATTCCATTGTCAGGTTCTCTTCTCATTTTCGTGGTTTAACCACGATTATGTGATTGCCTCCGCAATTCCTGTCTCTAACTCCCCTTCCTCGCAAAAACTGGCACTCCACGAGCATGTGTTTAGACAGTTTCATTAACGTAAACGGTTGCTTTTTACTCTGGCGGGCGAAAGGA  >  NC_000913/3031091‑3031357

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