Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,846 G→C 42.1% L188L (CTG→CTC)  narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,8460GC42.1% 56.6 / 120.2 102L188L (CTG→CTC) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base C (29/14);  total (55/47)
Fisher's exact test for biased strand distribution p-value = 2.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 107 total aligned reads displayed.

AGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCT  >  minE/835776‑835914
                                                                      |                                                                     
aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                                                       <  1:1663632/71‑1 (MQ=255)
   aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat                                                                    >  1:3935238/1‑70 (MQ=255)
   aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat                                                                    >  1:3079316/1‑70 (MQ=255)
     gagATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAA‑tcc                                                                 <  1:5372620/69‑1 (MQ=255)
     gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑ccc                                                                 <  1:560082/71‑1 (MQ=255)
     gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑cc                                                                  <  1:691822/70‑1 (MQ=255)
       gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                                               <  1:1687524/71‑1 (MQ=255)
       gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                                               <  1:2563743/71‑1 (MQ=255)
       gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                                               <  1:5643594/71‑1 (MQ=255)
       gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                                               <  1:5176742/71‑1 (MQ=255)
              aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtgt                                                        >  1:1436486/1‑71 (MQ=255)
                  aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                                                     >  1:2746859/1‑70 (MQ=255)
                   gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGa                                                   >  1:4997169/1‑71 (MQ=255)
                   gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCa                                                             <  1:370947/61‑1 (MQ=255)
                   gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCa                                                             <  1:2105302/61‑1 (MQ=255)
                   gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg                                                    <  1:4811165/70‑1 (MQ=255)
                   gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg                                                    <  1:1012600/70‑1 (MQ=255)
                   gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg                                                    <  1:4585217/70‑1 (MQ=255)
                    ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGc                                                              <  1:4370896/59‑1 (MQ=255)
                     ttgAAAACAACTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGAcc                                                 >  1:3261578/3‑71 (MQ=255)
                     tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                                                                     <  1:149834/52‑1 (MQ=255)
                     tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                                                                     <  1:3212023/52‑1 (MQ=255)
                     tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                                                                     <  1:5537493/52‑1 (MQ=255)
                     tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                                                         >  1:2505592/1‑63 (MQ=255)
                      tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGCACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                                                >  1:2191295/2‑71 (MQ=255)
                      tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                                                >  1:664389/2‑71 (MQ=255)
                      tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                                                >  1:727956/2‑71 (MQ=255)
                      tgAAAACACCTTCATGATGTACATGCCGCGCCTGTGCGAACACTGCCTCaa                                                                     >  1:1353259/2‑51 (MQ=255)
                      cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                                                         >  1:5626662/1‑62 (MQ=255)
                      cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgt                                                          <  1:2833587/61‑1 (MQ=255)
                      cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                                                >  1:139118/1‑71 (MQ=255)
                      cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                                                >  1:83831/1‑71 (MQ=255)
                      cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                                                >  1:4604916/1‑71 (MQ=255)
                       gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑cc                                                                  >  1:2313407/1‑52 (MQ=255)
                       gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑cc                                                                  >  1:4653343/1‑52 (MQ=255)
                       gAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTg                                               >  1:4134574/1‑71 (MQ=255)
                        aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGc                                              <  1:2847320/71‑1 (MQ=255)
                        aaaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGc                                              >  1:559860/1‑71 (MQ=255)
                         aaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                                                     >  1:1943985/1‑48 (MQ=38)
                         aaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                                                     >  1:5518181/1‑48 (MQ=38)
                          aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGc                                              >  1:3584047/1‑69 (MQ=255)
                           acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCg                                           <  1:3126183/71‑1 (MQ=255)
                           acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCg                                           <  1:1157043/71‑1 (MQ=255)
                           acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGccca                                           <  1:3879190/71‑2 (MQ=255)
                           acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGccca                                           <  1:3838918/71‑2 (MQ=255)
                            cacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGccc                                            >  1:2062836/1‑69 (MQ=255)
                             acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg                                         <  1:2812730/71‑1 (MQ=255)
                             acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg                                         <  1:2239291/71‑1 (MQ=255)
                             acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg                                         <  1:3395560/71‑1 (MQ=255)
                             acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg                                         <  1:1122157/71‑1 (MQ=255)
                              ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                                                       <  1:1034980/41‑1 (MQ=255)
                              ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg                                         <  1:5009832/70‑1 (MQ=255)
                              ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGc                                        >  1:1982968/1‑71 (MQ=255)
                                 tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGt                                     >  1:4422061/1‑71 (MQ=255)
                                 tCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggc                                     >  1:1641382/1‑70 (MQ=255)
                                  cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                                                     >  1:3309811/1‑39 (MQ=25)
                                  cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCgg                                      <  1:860058/69‑1 (MQ=255)
                                   atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                                                     >  1:2073067/1‑38 (MQ=25)
                                   atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCgg                                      <  1:73443/68‑1 (MQ=255)
                                   atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgc                                   >  1:3207135/1‑71 (MQ=255)
                                    tgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg                                  <  1:5125726/71‑1 (MQ=255)
                                    tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgct                                  <  1:5492776/71‑2 (MQ=255)
                                    tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCACAGC‑TGCGTGGCGACCTGCCCAAGCGgcgc                                   >  1:3424952/1‑70 (MQ=255)
                                      atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                <  1:4593942/71‑1 (MQ=255)
                                      atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                <  1:5454100/71‑1 (MQ=255)
                                      atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgcta                                 <  1:1889141/70‑3 (MQ=255)
                                           cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc                            >  1:2959457/3‑70 (MQ=255)
                                             tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAg                         >  1:4295982/1‑71 (MQ=255)
                                              gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac                        >  1:5474769/1‑69 (MQ=255)
                                              gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac                        >  1:1914269/1‑69 (MQ=255)
                                               ccGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                       >  1:4131440/1‑71 (MQ=255)
                                               ccGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGT‑GCGACCTGCCCGAGCGGTGCGATTTACAAGCGt                      >  1:4347390/1‑71 (MQ=255)
                                                 gcgcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                                                     >  1:3401348/1‑39 (MQ=255)
                                                 gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                          >  1:4792923/1‑66 (MQ=255)
                                                 gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                          >  1:2057523/1‑66 (MQ=255)
                                                 gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                          >  1:3909760/1‑66 (MQ=255)
                                                   gcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTga                    <  1:1424948/70‑1 (MQ=255)
                                                   gcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaa                   <  1:4008683/71‑1 (MQ=255)
                                                     cTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgct                                  <  1:1982515/54‑2 (MQ=255)
                                                      tgtgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                                                >  1:863524/1‑39 (MQ=255)
                                                      tgtgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaag                  <  1:5114094/69‑1 (MQ=255)
                                                       gtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                   >  1:3075944/1‑67 (MQ=255)
                                                        tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                            <  1:2743499/57‑1 (MQ=255)
                                                        tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                <  1:5674081/69‑1 (MQ=255)
                                                        tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCg                                       >  1:2189811/1‑46 (MQ=12)
                                                        tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                <  1:5420271/69‑1 (MQ=255)
                                                        tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                <  1:1690358/69‑1 (MQ=255)
                                                         gCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg                                  >  1:2820836/1‑50 (MQ=255)
                                                         gCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCC‑AGCGGTGCg                                  >  1:1934623/1‑49 (MQ=255)
                                                         gCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa                         >  1:4783294/1‑58 (MQ=255)
                                                         gCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa                         >  1:298998/1‑58 (MQ=255)
                                                          cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                     >  1:1109417/1‑62 (MQ=255)
                                                           gAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga              >  1:5406577/1‑68 (MQ=255)
                                                            aaCACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg           >  1:2365420/1‑70 (MQ=255)
                                                            aaCACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAagaagacg            <  1:809713/69‑3 (MQ=255)
                                                                cTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                <  1:4544790/61‑1 (MQ=255)
                                                                cTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCg      <  1:3433784/71‑1 (MQ=255)
                                                                     tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATTGCATCGTCCt  >  1:854048/1‑70 (MQ=255)
                                                                     tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc       >  1:1163641/1‑65 (MQ=255)
                                                                     tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt  >  1:5067144/1‑70 (MQ=255)
                                                                      |                                                                     
AGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCT  >  minE/835776‑835914

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