Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,903 T→C 7.5% D207D (GAT→GAC 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,9030TC 7.5% 228.1 / 8.7 80D207D (GAT→GACnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base C (3/3);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACA  >  minE/835834‑835970
                                                                     |                                                                   
cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                                                    >  1:1010523/1‑71 (MQ=255)
cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                                                    >  1:3143601/1‑71 (MQ=255)
  aaCACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                                                                   >  1:1651483/1‑70 (MQ=255)
  aaCACTGCCTCAATCCCAGATGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACgg                                                                   >  1:2350301/1‑70 (MQ=255)
   acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCa                                                                 <  1:71421/71‑1 (MQ=255)
       tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt                                                             >  1:3699780/1‑71 (MQ=255)
         ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                           >  1:911336/1‑71 (MQ=255)
         ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                           >  1:4021566/1‑71 (MQ=255)
          cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                          <  1:2695292/71‑1 (MQ=255)
          cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                          <  1:1308139/71‑1 (MQ=255)
             aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                          >  1:757825/1‑68 (MQ=255)
              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)
                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGc                                                                  <  1:1662162/51‑1 (MQ=255)
                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                                         <  1:2435676/60‑1 (MQ=255)
                      tgtgtgGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGa                                              <  1:1152876/71‑1 (MQ=255)
                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                          <  1:3808084/58‑1 (MQ=255)
                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                          <  1:35998/58‑1 (MQ=255)
                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                          <  1:1035948/58‑1 (MQ=255)
                        tgtgGCGACCTGCCCGAGCGGTTCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                            <  1:3843982/71‑1 (MQ=255)
                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                             <  1:3039922/70‑1 (MQ=255)
                        cgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                                              <  1:201995/68‑1 (MQ=255)
                           ggCGCCCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcgt                                                          >  1:2612787/1‑52 (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)
                                       gAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                                         <  1:325348/43‑1 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                   <  1:1439225/48‑1 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                   <  1:2819306/48‑1 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                          >  1:697619/1‑57 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                          >  1:2295813/1‑57 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                               >  1:727784/1‑68 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                            >  1:1429452/1‑71 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                            <  1:3693529/71‑1 (MQ=255)
                                        aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                            >  1:1832239/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)
                                               cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                                       >  1:578630/1‑37 (MQ=255)
                                               cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                                       >  1:3754100/1‑37 (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:1667466/1‑54 (MQ=255)
                                                    tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                                 >  1:2352836/1‑54 (MQ=255)
                                                     aCAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt               >  1:2751332/1‑71 (MQ=255)
                                                      cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg              <  1:2085022/71‑1 (MQ=255)
                                                      cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg              <  1:2122629/71‑1 (MQ=255)
                                                      cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg              <  1:1445361/71‑1 (MQ=255)
                                                      cAAGCGCGAGGAGGACGGCATCGTCCTGATCGACCAGGACAAGTGCCGCGGCTGGCGGATGTGCa                    <  1:1434330/65‑1 (MQ=255)
                                                       aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                            <  1:1203018/40‑1 (MQ=255)
                                                       aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg             <  1:770933/71‑1 (MQ=255)
                                                       aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg             <  1:3258100/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)
                                                                     |                                                                   
CGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACA  >  minE/835834‑835970

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: