Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,951,444 T→C 24.7% Q170Q (CAA→CAG fhuF ← ferric iron reductase involved in ferric hydroximate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,951,4440TC24.7% 138.5 / 41.0 113Q170Q (CAA→CAGfhuFferric iron reductase involved in ferric hydroximate transport
Reads supporting (aligned to +/- strand):  ref base T (51/34);  new base C (25/3);  total (76/37)
Fisher's exact test for biased strand distribution p-value = 4.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.24e-01

Only 100 of 113 total aligned reads displayed.

GATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTG  >  minE/2951381‑2951505
                                                               |                                                             
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTTTAGTGCTTGcac                                                          >  1:672389/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:2133545/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:1747178/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:194906/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:690389/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:3155695/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:479191/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:854750/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:909085/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:1199276/1‑69 (MQ=255)
gATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcac                                                          >  1:230608/1‑69 (MQ=255)
   cAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAAc                                                       >  1:1538575/1‑69 (MQ=255)
   cAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAAc                                                       <  1:2481635/69‑1 (MQ=255)
   cAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAAc                                                       <  1:3124557/69‑1 (MQ=255)
               aTTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGa                                                   >  1:667189/1‑61 (MQ=255)
               aTTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAgcgc                                           >  1:1195730/1‑69 (MQ=255)
               aTTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAgcgc                                           >  1:2044942/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCGCAACCGGAACCAGCGCCTg                                        >  1:813632/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCgg                                                    >  1:528748/1‑57 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGa                                                   >  1:1577717/1‑58 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                <  1:1753213/61‑1 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                <  1:2882900/61‑1 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGcc                                          >  1:1351390/1‑67 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCt                                         >  1:1568557/1‑68 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTg                                        >  1:585901/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTg                                        >  1:2553467/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTg                                        >  1:2686237/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTg                                        >  1:1928466/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTg                                        >  1:796267/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTg                                        >  1:918892/1‑69 (MQ=255)
                  aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCGTCTAGTGCTTGCACAACCGGAACCa                                               >  1:633862/1‑62 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGcaca                                                         >  1:1572349/1‑51 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCg                                                     >  1:3043877/1‑55 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                >  1:366181/1‑60 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                >  1:2133789/1‑60 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                >  1:3004114/1‑60 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                >  1:1380353/1‑60 (MQ=255)
                   cTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAAcc                                                >  1:1740141/1‑60 (MQ=255)
                          tAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTa                                >  1:449411/1‑69 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               <  1:2184935/69‑1 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               >  1:1872049/1‑69 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               >  1:2462026/1‑69 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               >  1:441505/1‑69 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               <  1:598413/69‑1 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               >  1:1517488/1‑69 (MQ=255)
                           aaGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTa                                >  1:258505/1‑68 (MQ=255)
                            aGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTaa                               <  1:2579235/68‑1 (MQ=255)
                            aGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAAc                              <  1:3850/69‑1 (MQ=255)
                            aGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAAc                              <  1:687053/69‑1 (MQ=255)
                             gTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAAc                                                       <  1:3286193/43‑1 (MQ=255)
                              tttGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGt                            <  1:2229139/69‑1 (MQ=255)
                              tttGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGt                            <  1:485128/69‑1 (MQ=255)
                                       gagctctcCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGtt                           >  1:1464785/4‑61 (MQ=255)
                                       gATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCg                   >  1:944731/1‑69 (MQ=255)
                                       gATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCAtt                     >  1:118019/1‑67 (MQ=255)
                                       gATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCg                   >  1:1619914/1‑69 (MQ=255)
                                       gATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCg                   >  1:2687350/1‑69 (MQ=255)
                                       gATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCg                   >  1:3366761/1‑69 (MQ=255)
                                       gATCTCTCCGGTCGCTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCg                   >  1:3003378/1‑69 (MQ=255)
                                                   cGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGa       >  1:1158427/1‑69 (MQ=255)
                                                   cGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGa       >  1:1572354/1‑69 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACGAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:670115/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:1138418/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:2033586/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:503903/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:351684/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:601609/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:1666838/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:3144271/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:1341708/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:801556/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:3075890/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:2689273/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:131704/69‑1 (MQ=255)
                                                    gCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGaa      <  1:1279786/69‑1 (MQ=255)
                                                     cTTCTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAAc                              <  1:1396382/44‑1 (MQ=255)
                                                      ttCTAGTGCTTCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgt   <  1:1715789/69‑1 (MQ=255)
                                                       tCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgt   <  1:1743350/69‑1 (MQ=255)
                                                       tCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgt   <  1:2838831/69‑1 (MQ=255)
                                                       tCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgt   <  1:2109732/69‑1 (MQ=255)
                                                       tCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgt   <  1:2724361/69‑1 (MQ=255)
                                                        cTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgtg  >  1:1190807/1‑69 (MQ=255)
                                                        cTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgtg  >  1:2735503/1‑69 (MQ=255)
                                                        cTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgtg  >  1:3100235/1‑69 (MQ=255)
                                                        cTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAAtgtg  >  1:2041667/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACTTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:1362539/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGa       >  1:3407375/1‑64 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2400276/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2851267/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:128895/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2502376/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2451538/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:1487136/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:629302/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:3429277/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:1686338/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2148879/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2117357/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtgtg  >  1:2000681/1‑69 (MQ=255)
                                                        cTAGTGCCTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGCTGCTGCGGCGAAtg    >  1:2239587/1‑67 (MQ=255)
                                                               |                                                             
GATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTG  >  minE/2951381‑2951505

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