Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,381 T→C 11.4% R33R (CGT→CGC 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,3810TC11.4% 185.4 / 14.9 70R33R (CGT→CGCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (37/25);  new base C (3/5);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

CTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGA  >  minE/835313‑835446
                                                                    |                                                                 
cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                                                  <  1:2688763/70‑1 (MQ=255)
cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                                                  >  1:696405/1‑70 (MQ=255)
cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                                                  >  1:701179/1‑70 (MQ=255)
  gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg                                                               >  1:2275248/1‑71 (MQ=255)
    acctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc                                                             <  1:3752733/69‑1 (MQ=255)
      ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                                                             <  1:1376149/69‑1 (MQ=255)
      ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                                                             <  1:776056/69‑1 (MQ=255)
         gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                                                          <  1:1161257/69‑1 (MQ=255)
         gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                                                          <  1:612248/69‑1 (MQ=255)
         gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                                                          <  1:687478/69‑1 (MQ=255)
              gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                   <  1:372211/71‑1 (MQ=255)
               tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                                                          >  1:908627/1‑63 (MQ=255)
                 cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt                                                >  1:947244/1‑71 (MQ=255)
                   tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:2318982/1‑64 (MQ=255)
                     ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt                                            >  1:3936105/1‑71 (MQ=255)
                     ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt                                            >  1:3887534/1‑71 (MQ=255)
                       cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                          <  1:796317/71‑1 (MQ=255)
                         gCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                                                  <  1:771015/45‑1 (MQ=255)
                          ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                        <  1:815090/70‑1 (MQ=255)
                          ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                        <  1:3895725/70‑1 (MQ=255)
                          ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa                                                       <  1:1636282/55‑1 (MQ=255)
                            acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                                                         <  1:3749559/51‑1 (MQ=255)
                            acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:247609/1‑55 (MQ=255)
                            acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:3281131/1‑55 (MQ=255)
                            acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac                                     >  1:2683398/1‑71 (MQ=255)
                             tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg                                    >  1:244932/2‑71 (MQ=255)
                             cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                    >  1:2469016/1‑71 (MQ=255)
                               ccTGTTCGGTGACCTGTAAAAACGTCTGTACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                  >  1:2059072/1‑70 (MQ=255)
                               ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:1167491/1‑52 (MQ=255)
                               ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                  <  1:3461557/71‑1 (MQ=255)
                                cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa                                                       <  1:3328891/49‑1 (MQ=255)
                                 tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt                                             <  1:2760596/58‑1 (MQ=255)
                                 tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                      >  1:2640198/1‑65 (MQ=255)
                                  gTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                      >  1:3949634/1‑64 (MQ=255)
                                   ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATg                                                          <  1:940882/43‑1 (MQ=16)
                                   ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt                                            <  1:3319294/57‑1 (MQ=255)
                                   ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                                                  <  1:512008/35‑1 (MQ=255)
                                   ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:2868468/1‑48 (MQ=255)
                                   ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                           >  1:2777811/1‑58 (MQ=255)
                                   ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                           >  1:197828/1‑58 (MQ=255)
                                   ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                               <  1:1224112/70‑1 (MQ=255)
                                    tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac                                     >  1:2321343/1‑63 (MQ=255)
                                     cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                  >  1:2559532/1‑65 (MQ=255)
                                     cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAc                             >  1:1769383/1‑70 (MQ=255)
                                       gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCa                           >  1:3987689/4‑70 (MQ=255)
                                       gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                                                         <  1:1442120/40‑1 (MQ=255)
                                       gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:2612771/1‑44 (MQ=255)
                                         tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAat                                                      >  1:3478885/1‑41 (MQ=255)
                                         tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTgg                                              >  1:236573/1‑49 (MQ=255)
                                          aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                                                     >  1:53289/1‑41 (MQ=255)
                                             tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAatat                                                    <  1:195268/39‑1 (MQ=255)
                                             tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                        <  1:3215261/67‑1 (MQ=255)
                                              gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                         >  1:2514953/1‑49 (MQ=255)
                                               tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                   <  1:3336516/70‑1 (MQ=255)
                                                aaaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                   >  1:3875650/1‑69 (MQ=255)
                                                 aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCg                      >  1:3292367/1‑65 (MQ=255)
                                                 aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCg                      >  1:2804611/1‑65 (MQ=255)
                                                 aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAg                <  1:2807869/71‑1 (MQ=255)
                                                         tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaaa                              >  1:3995356/1‑49 (MQ=255)
                                                         tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                               >  1:3169442/1‑48 (MQ=255)
                                                          ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                  <  1:1358102/44‑1 (MQ=255)
                                                          ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                  <  1:818788/44‑1 (MQ=255)
                                                           gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                       >  1:2079717/1‑54 (MQ=255)
                                                           gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                        >  1:2084506/1‑53 (MQ=255)
                                                            aCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                >  1:233150/1‑44 (MQ=255)
                                                            aCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCga      >  1:3414097/1‑70 (MQ=255)
                                                              cAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAggg              <  1:3489602/60‑1 (MQ=255)
                                                               aGCCGTAAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactga  <  1:2591958/71‑1 (MQ=255)
                                                                  cGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTc          >  1:1178231/1‑60 (MQ=255)
                                                                   gTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaaa                              >  1:3321462/1‑39 (MQ=255)
                                                                    tGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                   <  1:1629368/49‑1 (MQ=255)
                                                                    |                                                                 
CTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGA  >  minE/835313‑835446

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: