Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,427,002 T→C 40.8% P43P (CCA→CCG napH ← ferredoxin‑type protein essential for electron transfer from ubiquinol to periplasmic nitrate reductase (NapAB)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,427,0020TC40.8% 47.0 / 88.5 120P43P (CCA→CCGnapHferredoxin‑type protein essential for electron transfer from ubiquinol to periplasmic nitrate reductase (NapAB)
Reads supporting (aligned to +/- strand):  ref base T (41/30);  new base C (36/13);  total (77/43)
Fisher's exact test for biased strand distribution p-value = 8.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 120 total aligned reads displayed.

CAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTG  >  minE/1426936‑1427067
                                                                  |                                                                 
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATgg                                                                 <  1:3414711/69‑1 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATgg                                                                 <  1:12005/69‑1 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:396357/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:1618660/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:1858616/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:1571212/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:2058534/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:2195516/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:2447567/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:3296939/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:855048/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:1206254/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACgg                                                                 >  1:2988252/1‑69 (MQ=255)
cAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGAT‑‑ACACACCAAACCACGGAc                                                               >  1:499244/1‑69 (MQ=255)
    ggAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCg                                                             >  1:2116412/1‑69 (MQ=255)
    ggAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCg                                                             >  1:1802801/1‑69 (MQ=255)
    ggAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCg                                                             >  1:3131140/1‑69 (MQ=255)
    ggAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCg                                                             >  1:841984/1‑69 (MQ=255)
    ggAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCg                                                             >  1:2925460/1‑69 (MQ=255)
    ggAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCg                                                             >  1:2751928/1‑69 (MQ=255)
          gtgtCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaa                                                       <  1:603468/69‑1 (MQ=255)
          gtgtCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaa                                                       <  1:2839182/69‑1 (MQ=255)
                 aTAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGGTCaaaaa                                                     >  1:3372778/1‑64 (MQ=255)
                 aTAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaaaa                                                     >  1:2102315/1‑64 (MQ=255)
                 aTAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaaaa                                                     >  1:2163519/1‑64 (MQ=255)
                 aTAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaaaa                                                     >  1:574883/1‑64 (MQ=255)
                 aTAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaaaa                                                     >  1:2858439/1‑64 (MQ=255)
                 aTAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCaaaaa                                                     >  1:1016566/1‑64 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCTAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:2262109/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:2902253/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:100765/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:1367264/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:2423308/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:1360793/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:2643150/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:1405542/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:1359936/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:1296887/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:2975016/59‑1 (MQ=255)
                        gCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:1172381/59‑1 (MQ=255)
                        gCTACTCCTGTAGTTGCCGTGCAAGATCCACACACCAAACCACGGACCGCTCAAAAACa                                                   <  1:2187376/59‑1 (MQ=255)
                          tACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCa                                          >  1:3097697/1‑66 (MQ=255)
                              gctgTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaagaa                                   >  1:1032624/1‑69 (MQ=255)
                              gctgTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaagaa                                   >  1:3177635/1‑69 (MQ=255)
                              gctgTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:3135057/67‑1 (MQ=255)
                              gctgTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:883134/67‑1 (MQ=255)
                                    gTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     >  1:1220856/1‑61 (MQ=255)
                                    gTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     >  1:1581693/1‑61 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:1663183/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:2544453/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:812097/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:234239/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:2227735/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:623348/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:1957289/1‑69 (MQ=255)
                                       gCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAg                          >  1:2096536/1‑69 (MQ=255)
                                             cAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:2250371/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGCCCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:2991918/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:847207/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:770200/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:1559395/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:3371997/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:3376058/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:1124429/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:1551537/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:2352664/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:596736/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:2821139/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:1412585/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:43583/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:2569452/1‑68 (MQ=255)
                                             cAAGATCCACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGAc                     >  1:2508537/1‑68 (MQ=255)
                                               aGATCCACACACCAAACGATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:3149329/68‑1 (MQ=255)
                                               aGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:1264179/68‑1 (MQ=255)
                                               aGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:1289544/68‑1 (MQ=255)
                                               aGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:2992638/68‑1 (MQ=255)
                                               aGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:1607539/68‑1 (MQ=255)
                                               aGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:2395253/68‑1 (MQ=255)
                                                gATCCACACACCAAACGATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   <  1:3145478/67‑1 (MQ=255)
                                                  tCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAac                >  1:139928/1‑68 (MQ=255)
                                                   ccACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACa                             >  1:2380955/1‑54 (MQ=255)
                                                   ccACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACa                             >  1:2112996/1‑54 (MQ=255)
                                                   ccACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACa                             >  1:1013008/1‑54 (MQ=255)
                                                   ccACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACa                             >  1:1940734/1‑54 (MQ=255)
                                                   ccACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACa                             >  1:3297865/1‑54 (MQ=255)
                                                   ccACACACCAAACCACGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACa                             >  1:2594787/1‑54 (MQ=255)
                                                    cacacaCCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:3077129/45‑1 (MQ=255)
                                                    cacacaCCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:3082228/45‑1 (MQ=255)
                                                    cacacaCCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:1820844/45‑1 (MQ=255)
                                                    cacacaCCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:161900/45‑1 (MQ=255)
                                                    cacacaCCAAACCATGGACCGCTCAAAAACATCCCCAGCACgaag                                     <  1:1535434/45‑1 (MQ=255)
                                                    cacacaCCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGt                   >  1:1432207/1‑63 (MQ=255)
                                                              aCCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGt   <  1:1671268/69‑1 (MQ=255)
                                                              aCCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGt   <  1:3092577/69‑1 (MQ=255)
                                                               ccATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTTACAAAGGCGACGTAACACCAGCCAACGGTg  <  1:1156867/69‑1 (MQ=255)
                                                               ccATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTg  <  1:2677643/69‑1 (MQ=255)
                                                               ccATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTg  <  1:3054565/69‑1 (MQ=255)
                                                               ccATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTg  <  1:2343287/69‑1 (MQ=255)
                                                               ccATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTg  <  1:2131328/69‑1 (MQ=255)
                                                               ccATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTg  <  1:630292/69‑1 (MQ=255)
                                                                  |                                                                 
CAGCGGAACGGTGTCGAATAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTG  >  minE/1426936‑1427067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: