Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,909 C→T 5.1% 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,9090CT 5.1% 262.3 / 2.7 78I209I (ATC→ATTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base T (2/2);  total (37/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 = 8.19e-01

TGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATC  >  minE/835841‑835978
                                                                    |                                                                     
tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt                                                                     >  1:3699780/1‑71 (MQ=255)
  ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                                   >  1:4021566/1‑71 (MQ=255)
  ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                                                                   >  1:911336/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)
               tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                                                                 <  1:2435676/60‑1 (MQ=255)
               tgtgtgGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGa                                                      <  1:1152876/71‑1 (MQ=255)
                gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                  <  1:35998/58‑1 (MQ=255)
                gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                  <  1:1035948/58‑1 (MQ=255)
                gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                  <  1:3808084/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:1832239/1‑71 (MQ=255)
                                 aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                    <  1:3693529/71‑1 (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)
                                        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:2122629/71‑1 (MQ=255)
                                               cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:2085022/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: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:2578843/1‑71 (MQ=255)
                                                       aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg              >  1:1113775/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)
                                                                    |                                                                     
TGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATC  >  minE/835841‑835978

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: