Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,396 T→C 12.5% Y38Y (TAT→TAC 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,3960TC12.5% 125.1 / 10.1 48Y38Y (TAT→TACnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (28/14);  new base C (3/3);  total (31/17)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAAT  >  minE/835327‑835466
                                                                     |                                                                      
gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                                                                       <  1:372211/71‑1 (MQ=255)
   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: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)
                           tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTgg                                                                  >  1:236573/1‑49 (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: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: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)
                                                        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)
                                                                     |                                                                      
GTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAAT  >  minE/835327‑835466

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: