Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,909 C→T 14.3% I209I (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,9090CT14.3% 223.4 / 23.8 91I209I (ATC→ATTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (30/48);  new base T (6/7);  total (36/55)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATC  >  minE/835840‑835978
                                                                     |                                                                     
cTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCg                                                                      <  1:3433784/71‑1 (MQ=255)
     tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATTGCATCGTCCt                                                                  >  1:854048/1‑70 (MQ=255)
     tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc                                                                       >  1:1163641/1‑65 (MQ=255)
     tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                  >  1:5067144/1‑70 (MQ=255)
       aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGGGAAGAAGATGGCATCGTCCTGAt                                                               <  1:4859008/71‑1 (MQ=255)
         cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                                             <  1:4446315/71‑1 (MQ=255)
           cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                           <  1:2858466/71‑1 (MQ=255)
            ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                                          <  1:4683953/71‑1 (MQ=255)
                tgcgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                                      >  1:2592635/4‑71 (MQ=255)
                 gcgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                                      >  1:4385442/3‑70 (MQ=255)
                  tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                     >  1:534635/1‑70 (MQ=255)
                  tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                     >  1:4396349/1‑70 (MQ=255)
                  tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                     >  1:1775170/1‑70 (MQ=255)
                  tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                     >  1:456002/1‑70 (MQ=255)
                  tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                     >  1:2016248/1‑70 (MQ=255)
                      gCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                  >  1:5052181/1‑69 (MQ=255)
                           cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                                   >  1:2232410/1‑47 (MQ=255)
                           cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                           <  1:4317549/71‑1 (MQ=255)
                             gCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                                                 >  1:3817871/1‑47 (MQ=255)
                              cccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                                         <  1:512664/70‑1 (MQ=255)
                              cccAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                                         <  1:2476477/70‑1 (MQ=255)
                                cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                                      >  1:5396429/1‑71 (MQ=255)
                                  aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                    <  1:4305586/71‑1 (MQ=255)
                                    cGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATg                                                <  1:4212846/57‑1 (MQ=255)
                                    cGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt                                  <  1:3102966/71‑2 (MQ=255)
                                     ggTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAt                                 <  1:2869767/71‑1 (MQ=255)
                                      gTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                                          <  1:3894254/45‑1 (MQ=255)
                                      gTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                                          <  1:3645624/45‑1 (MQ=255)
                                      gTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc                                                          <  1:4924439/45‑1 (MQ=255)
                                       tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                    <  1:4264622/50‑1 (MQ=255)
                                       tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                                         <  1:2644161/61‑1 (MQ=255)
                                       tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                  <  1:1157081/68‑1 (MQ=255)
                                       tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgt                               <  1:3130205/71‑1 (MQ=255)
                                       tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATAGACCAGGATa                                                    <  1:4733414/50‑1 (MQ=255)
                                        gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                              >  1:3141955/1‑71 (MQ=255)
                                        gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                              >  1:3749414/1‑71 (MQ=255)
                                         cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                              >  1:5115664/1‑70 (MQ=255)
                                           aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt                           <  1:4883487/71‑1 (MQ=255)
                                           aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt                           <  1:2886675/71‑1 (MQ=255)
                                           aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt                           <  1:3329174/71‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                              <  1:936548/35‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                              <  1:821844/35‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                              <  1:2013520/35‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                              <  1:1611923/35‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                    <  1:262897/45‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                    <  1:3165576/45‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:371163/71‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:2778101/71‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:564971/71‑1 (MQ=255)
                                            tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:870498/71‑1 (MQ=255)
                                            tctACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt                                  <  1:4137480/61‑2 (MQ=255)
                                             ctACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCAt                           >  1:4756963/2‑69 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:2745206/69‑1 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:4344729/69‑1 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:2259651/69‑1 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                          <  1:4979942/69‑1 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa                         <  1:1514189/70‑1 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa                         >  1:4708042/1‑70 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                        >  1:5483122/1‑71 (MQ=255)
                                              tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                        >  1:5361167/1‑71 (MQ=255)
                                              tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCAtaag                        >  1:3331763/1‑68 (MQ=255)
                                                cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:5134126/71‑1 (MQ=255)
                                                cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:1800628/71‑1 (MQ=255)
                                                cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:1970781/71‑1 (MQ=255)
                                                cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:2055133/71‑1 (MQ=255)
                                                cAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCAtaag                        <  1:1404041/69‑4 (MQ=255)
                                                 aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                     <  1:1942584/71‑1 (MQ=255)
                                                  aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                           >  1:2214513/1‑48 (MQ=255)
                                                  aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                           >  1:1531544/1‑48 (MQ=255)
                                                   gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaa                                                   >  1:755488/1‑39 (MQ=255)
                                                   gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaa                                                   >  1:903774/1‑39 (MQ=255)
                                                   gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg                                >  1:128815/1‑58 (MQ=255)
                                                        aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg              <  1:1540661/71‑1 (MQ=255)
                                                         agaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                                         >  1:64715/1‑43 (MQ=255)
                                                         agaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                                         >  1:3814929/1‑43 (MQ=255)
                                                          gaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                           >  1:1698569/1‑40 (MQ=255)
                                                          gaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa                         >  1:1324929/1‑58 (MQ=255)
                                                          gaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                        >  1:704110/1‑59 (MQ=255)
                                                          gaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTa            >  1:3800940/1‑71 (MQ=255)
                                                          gaagaCGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt                                  >  1:1228482/1‑48 (MQ=255)
                                                          gaagaCGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt                                  >  1:5656719/1‑48 (MQ=255)
                                                              acggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTAaaaa        <  1:5016134/69‑1 (MQ=255)
                                                              aTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgt                               <  1:3741857/48‑1 (MQ=255)
                                                               cggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACaaaa       <  1:309817/70‑1 (MQ=255)
                                                               cggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACaaaa       <  1:1296422/70‑1 (MQ=255)
                                                                 gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa     <  1:386177/71‑1 (MQ=255)
                                                                 gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa     <  1:4788190/71‑1 (MQ=255)
                                                                  cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt                           <  1:3364544/48‑1 (MQ=255)
                                                                   aTCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaaa    <  1:2707284/70‑1 (MQ=255)
                                                                    tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt                       >  1:527778/1‑50 (MQ=255)
                                                                    tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt                       >  1:527203/1‑50 (MQ=255)
                                                                    tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc  <  1:1990543/71‑1 (MQ=255)
                                                                     |                                                                     
CTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATC  >  minE/835840‑835978

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: