Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,414 G→C 19.1% V44V (GTG→GTC 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,4140GC19.1% 97.6 / 17.4 47V44V (GTG→GTCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base C (6/3);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGAT  >  minE/835348‑835482
                                                                  |                                                                    
ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                                  >  1:1289810/1‑71 (MQ=255)
ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                                  >  1:2403902/1‑71 (MQ=255)
ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                    <  1:1053184/69‑1 (MQ=255)
    gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaa                                                              <  1:3125305/68‑1 (MQ=255)
    gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaa                                                              <  1:2572197/68‑1 (MQ=255)
     ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                                             >  1:3842573/1‑71 (MQ=255)
     ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                                             >  1:5082360/1‑71 (MQ=255)
     ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                                             >  1:1586002/1‑71 (MQ=255)
      gaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCaaac                                                            >  1:3144337/2‑69 (MQ=255)
        ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                                                               <  1:314561/66‑1 (MQ=255)
          tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                                        >  1:5021954/1‑71 (MQ=255)
           gTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCggg                                                        >  1:3579405/1‑70 (MQ=255)
           gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                    >  1:2549367/1‑58 (MQ=255)
           gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                    >  1:1070581/1‑58 (MQ=255)
           gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                                        >  1:65983/1‑70 (MQ=255)
            tAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCggg                                                        <  1:5453935/69‑1 (MQ=255)
            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                    <  1:2172006/57‑1 (MQ=255)
            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                    <  1:2614785/57‑1 (MQ=255)
            tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                                                       <  1:5391119/70‑1 (MQ=255)
             aaaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                                                       >  1:763680/1‑69 (MQ=255)
              aaaaCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCGGGcatt                                                    >  1:3831433/1‑68 (MQ=255)
              aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                                              >  1:4597067/1‑61 (MQ=255)
              aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                                              >  1:2311803/1‑61 (MQ=255)
              aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAg                                                    <  1:5291826/71‑1 (MQ=255)
               aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                                      <  1:2925466/52‑1 (MQ=255)
               aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                                      <  1:3860438/52‑1 (MQ=255)
               aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                                      <  1:3225637/52‑1 (MQ=255)
               aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                                      <  1:3669423/52‑1 (MQ=255)
               aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCGGGGCCa                                                     <  1:7934/69‑1 (MQ=255)
               aaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCa                                                     <  1:1014431/69‑1 (MQ=255)
                 aCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                                  >  1:1029114/1‑54 (MQ=255)
                   gTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                                               >  1:1468629/1‑71 (MQ=255)
                      tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                                            >  1:2553376/1‑71 (MQ=255)
                        gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                                            <  1:116334/69‑1 (MQ=255)
                        gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                                            <  1:5587888/69‑1 (MQ=255)
                             gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                                     <  1:4536268/71‑1 (MQ=255)
                             gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                                     <  1:5170359/71‑1 (MQ=255)
                             gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                                     <  1:2235549/71‑1 (MQ=255)
                              ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact                                    <  1:2013259/71‑1 (MQ=255)
                                gTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTgg                                  <  1:1987430/71‑1 (MQ=255)
                                   aaGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGa                                <  1:591039/70‑1 (MQ=255)
                                    aGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaa                              <  1:4420474/71‑1 (MQ=255)
                                       cGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                                               <  1:4216180/51‑1 (MQ=255)
                                         tGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc                                                            <  1:3736977/36‑1 (MQ=255)
                                         tGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGcc                                                      <  1:2714899/42‑1 (MQ=255)
                                          ggAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAg                         >  1:719327/1‑70 (MQ=255)
                                           gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                                            >  1:2942501/1‑50 (MQ=255)
                                           gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTcc                                             >  1:3956612/1‑49 (MQ=255)
                                                tGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaaa                   <  1:165773/70‑1 (MQ=255)
                                                 gCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATa                 >  1:2757157/1‑71 (MQ=255)
                                                                gTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGAt  >  1:810920/1‑71 (MQ=255)
                                                                  |                                                                    
TTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTGGAT  >  minE/835348‑835482

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: