Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,399 G→A 11.3% A39A (GCG→GCA 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,3990GA11.3% 150.1 / 9.0 53A39A (GCG→GCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (30/17);  new base A (3/3);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACA  >  minE/835330‑835469
                                                                     |                                                                      
cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt                                                                       >  1:947244/1‑71 (MQ=255)
    ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt                                                                   >  1:3936105/1‑71 (MQ=255)
    ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt                                                                   >  1:3887534/1‑71 (MQ=255)
      cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                                 <  1:796317/71‑1 (MQ=255)
         ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                               <  1:3895725/70‑1 (MQ=255)
         ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                               <  1:815090/70‑1 (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)
              ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                         <  1:3461557/71‑1 (MQ=255)
                tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt                                                                    <  1:2760596/58‑1 (MQ=255)
                tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                             >  1:2640198/1‑65 (MQ=255)
                 gTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                             >  1:3949634/1‑64 (MQ=255)
                  ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt                                                                   <  1:3319294/57‑1 (MQ=255)
                  ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                                                  >  1:197828/1‑58 (MQ=255)
                  ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                                                  >  1:2777811/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)
                        tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTgg                                                                     >  1:236573/1‑49 (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:2804611/1‑65 (MQ=255)
                                aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCg                                             >  1:3292367/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:818788/44‑1 (MQ=255)
                                         ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                         <  1:1358102/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)
                                                     aaGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                                 >  1:1061547/1‑40 (MQ=255)
                                                     aaGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaa                  <  1:4137788/71‑1 (MQ=255)
                                                       ggCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaaa                >  1:1071799/1‑71 (MQ=255)
                                                             gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                               >  1:2147852/1‑50 (MQ=255)
                                                                tatGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaa       <  1:2394829/71‑1 (MQ=255)
                                                                tatGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaa        <  1:2669502/70‑1 (MQ=255)
                                                                  tGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAAt     <  1:958741/71‑1 (MQ=255)
                                                                  tGCGTGGTTCAACAACGTGGAAACCAAGCCGGGACAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaaa      >  1:1278114/1‑70 (MQ=255)
                                                                    cGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                               >  1:2774606/1‑43 (MQ=255)
                                                                    cGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATa    <  1:2095320/70‑1 (MQ=255)
                                                                     gTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACa  <  1:3024812/71‑1 (MQ=255)
                                                                     |                                                                      
CATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACA  >  minE/835330‑835469

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: