Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,918 C→T 5.5% I212I (ATC→ATT 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,9180CT 5.5% 237.7 / 2.8 74I212I (ATC→ATTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base T (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

ACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAAC  >  minE/835848‑835987
                                                                      |                                                                     
aCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                       >  1:3624603/1‑71 (MQ=255)
 cccGGCATGTGTGGCGACCTGCGCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                       >  1:4018316/1‑70 (MQ=255)
 cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                       >  1:3070801/1‑70 (MQ=255)
   cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                                    <  1:4018509/71‑1 (MQ=255)
        tgtgtgGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGa                                                               <  1:1152876/71‑1 (MQ=255)
          tgtgGCGACCTGCCCGAGCGGTTCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                             <  1:3843982/71‑1 (MQ=255)
          tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                              <  1:3039922/70‑1 (MQ=255)
          cgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                                               <  1:201995/68‑1 (MQ=255)
              gCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAAt                                                          <  1:1605920/70‑1 (MQ=255)
                gACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                                                      <  1:2757191/56‑1 (MQ=255)
                  ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                                                     >  1:1070637/1‑55 (MQ=255)
                  ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                                    >  1:2224999/1‑56 (MQ=255)
                   cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                                                     >  1:1070854/1‑54 (MQ=255)
                   cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                                    >  1:2985414/1‑71 (MQ=255)
                       ccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                                                <  1:1017280/71‑1 (MQ=255)
                        cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                                               >  1:853566/1‑71 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                                    <  1:2819306/48‑1 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                                    <  1:1439225/48‑1 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                           >  1:2295813/1‑57 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                           >  1:697619/1‑57 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                                                >  1:727784/1‑68 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                             <  1:3693529/71‑1 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                             >  1:1832239/1‑71 (MQ=255)
                          aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                             >  1:1429452/1‑71 (MQ=255)
                            cGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                           <  1:1031246/71‑1 (MQ=255)
                            cGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt                                           <  1:3105730/71‑2 (MQ=255)
                               tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                             <  1:2986992/50‑1 (MQ=255)
                               tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                             <  1:3485686/66‑1 (MQ=255)
                                gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                       >  1:3744958/1‑71 (MQ=255)
                                 cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                       <  1:3394863/70‑1 (MQ=255)
                                    tttACAAGCGTGAAGAAGATGGCATCGTCCTGATGGACCAGGATAAATGCCGt                                                     <  1:4032851/53‑1 (MQ=255)
                                    tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                                    <  1:2983292/54‑1 (MQ=255)
                                    tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                                   <  1:1977388/71‑1 (MQ=255)
                                      tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                                                     >  1:4123587/1‑35 (MQ=255)
                                      tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                                 >  1:1857657/1‑71 (MQ=255)
                                      tACAAGAGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa                                  >  1:539242/1‑70 (MQ=255)
                                      tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                                                  >  1:2352836/1‑54 (MQ=255)
                                      tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                                                  >  1:1667466/1‑54 (MQ=255)
                                       aCAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt                                >  1:2751332/1‑71 (MQ=255)
                                        cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                               <  1:2122629/71‑1 (MQ=255)
                                        cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                               <  1:1445361/71‑1 (MQ=255)
                                        cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                               <  1:2085022/71‑1 (MQ=255)
                                        cAAGCGCGAGGAGGACGGCATCGTCCTGATCGACCAGGACAAGTGCCGCGGCTGGCGGATGTGCa                                     <  1:1434330/65‑1 (MQ=255)
                                         aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                             <  1:1203018/40‑1 (MQ=255)
                                         aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                              <  1:3258100/71‑1 (MQ=255)
                                         aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                              <  1:770933/71‑1 (MQ=255)
                                           gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                           >  1:3891022/1‑40 (MQ=255)
                                           gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                           >  1:3903465/1‑40 (MQ=255)
                                           gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                             <  1:3126887/54‑1 (MQ=255)
                                            cGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                                                  <  1:1464355/48‑1 (MQ=255)
                                            cGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                             <  1:3506649/53‑1 (MQ=255)
                                              tGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGt                                            <  1:917889/52‑1 (MQ=255)
                                                aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGc                                      <  1:1409670/56‑1 (MQ=255)
                                                aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc                        >  1:836903/1‑70 (MQ=255)
                                                aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                       >  1:1113775/1‑71 (MQ=255)
                                                aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                       >  1:2578843/1‑71 (MQ=255)
                                                  gaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                                               >  1:328541/1‑45 (MQ=255)
                                                     gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                              >  1:166444/1‑59 (MQ=255)
                                                     gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACa                   >  1:3400660/1‑70 (MQ=255)
                                                      aTGGGATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                               <  1:1931453/57‑1 (MQ=255)
                                                        ggCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                               <  1:2478135/55‑1 (MQ=255)
                                                         gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                                 >  1:1487657/1‑52 (MQ=255)
                                                          cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaaa             >  1:3841752/1‑71 (MQ=255)
                                                            tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc           <  1:1790183/71‑1 (MQ=255)
                                                              gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTa         <  1:608344/71‑1 (MQ=255)
                                                              gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTa         <  1:1253443/71‑1 (MQ=255)
                                                               tCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                                   <  1:342487/44‑1 (MQ=255)
                                                               tCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAc        <  1:949667/71‑1 (MQ=255)
                                                                ccTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                           >  1:909908/1‑35 (MQ=255)
                                                                   gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGcc                         >  1:3569854/1‑50 (MQ=255)
                                                                   gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                       <  1:3551819/52‑1 (MQ=255)
                                                                   gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTc     >  1:3344781/1‑70 (MQ=255)
                                                                    aTCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCaa   >  1:1176971/1‑71 (MQ=255)
                                                                     tCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGa                             >  1:1196998/1‑44 (MQ=255)
                                                                     tCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGa                             >  1:1558340/1‑44 (MQ=255)
                                                                     tCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGa                             >  1:1642629/1‑44 (MQ=255)
                                                                     tCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAAc  >  1:3172560/1‑71 (MQ=255)
                                                                      |                                                                     
ACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAAC  >  minE/835848‑835987

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