Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,426,109 A→T 21.0% D54E (GAT→GAA napB ← nitrate reductase, small, cytochrome C550 subunit, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,426,1090AT21.0% 137.3 / 28.7 76D54E (GAT→GAAnapBnitrate reductase, small, cytochrome C550 subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (4/56);  new base T (8/8);  total (12/64)
Fisher's exact test for biased strand distribution p-value = 2.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGATACTTCCGGC  >  minE/1426047‑1426163
                                                              |                                                      
tGATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGc                                                  <  1:3139154/69‑1 (MQ=255)
tGATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGc                                                  <  1:1094534/69‑1 (MQ=255)
tGATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGc                                                  <  1:559757/69‑1 (MQ=255)
tGATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGc                                                  <  1:7716/69‑1 (MQ=255)
 gATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGc                                                  <  1:2060730/68‑1 (MQ=255)
    aaCCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCtt                                              >  1:1930579/1‑69 (MQ=255)
    aaCCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCtt                                              >  1:2550075/1‑69 (MQ=255)
    aaCCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCtt                                              >  1:1565985/1‑69 (MQ=255)
        ttCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTgg                                          <  1:2679513/69‑1 (MQ=255)
        ttCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTgg                                          <  1:686482/69‑1 (MQ=255)
        ttCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTgg                                          <  1:359329/69‑1 (MQ=255)
         tCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCTGCATCCGATCCTGCTCTTTTGGc                                         <  1:2609431/69‑1 (MQ=255)
          cAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCa                                        >  1:282029/1‑69 (MQ=255)
          cAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCa                                        >  1:1679293/1‑69 (MQ=255)
          cAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCa                                        >  1:1877244/1‑69 (MQ=255)
          cAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCa                                        >  1:1046817/1‑69 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGCTTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:3308014/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:496626/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:3246939/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:2886524/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:1004985/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:376764/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:2492840/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:2901939/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:2367886/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:2168954/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:3185773/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:1043251/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:527909/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:1912939/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:1371907/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:1413034/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:643493/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:1524752/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:163121/69‑1 (MQ=255)
              gCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCAGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:728794/69‑1 (MQ=255)
               cTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:2890437/68‑1 (MQ=255)
                 aTGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAAt                                 <  1:2136994/69‑1 (MQ=255)
                 aTGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAAt                                 <  1:1484911/69‑1 (MQ=255)
                 aTGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAAt                                 <  1:2680607/69‑1 (MQ=255)
                    cGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGc                              <  1:1617740/69‑1 (MQ=255)
                    cGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGc                              <  1:642696/69‑1 (MQ=255)
                     gggATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCg                                    <  1:428330/62‑1 (MQ=255)
                         tCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCtt                         <  1:1648214/69‑1 (MQ=255)
                         tCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCtt                         <  1:385350/69‑1 (MQ=255)
                         tCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCtt                         <  1:2894429/69‑1 (MQ=255)
                         tCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCtt                         <  1:1168810/69‑1 (MQ=255)
                         tCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCta                         <  1:1584927/69‑2 (MQ=255)
                           aTTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCtt                         <  1:2012251/67‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:2541644/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:545303/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:394981/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:320239/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:811213/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:2821284/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:1436990/43‑1 (MQ=255)
                             tGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:2335458/43‑1 (MQ=255)
                              ggcggcTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCtgc                    >  1:873543/1‑68 (MQ=255)
                              ggcggcTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCtgc                    >  1:657285/1‑68 (MQ=255)
                              ggcggcTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCtgc                    >  1:1861260/1‑68 (MQ=255)
                              ggcggcTGATTCACATAGTTCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCtgc                    >  1:24996/1‑68 (MQ=255)
                                 ggcTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTc                 >  1:2976752/1‑69 (MQ=255)
                                     gATTCACATAGTTCAGCGGCATCCGATCCTGctct                                               <  1:1902733/35‑1 (MQ=255)
                                        tCACATAGTTCAGCGGCATCCGATCCTGCTTTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGa            <  1:2681268/67‑1 (MQ=255)
                                        tCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGa            <  1:25232/67‑1 (MQ=255)
                                        tCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGa            <  1:54516/67‑1 (MQ=255)
                                        tCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGa            <  1:1726362/67‑1 (MQ=255)
                                        tCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGa            <  1:2727984/67‑1 (MQ=255)
                                                 tCATCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:1959971/68‑1 (MQ=255)
                                                 tCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:2087407/68‑1 (MQ=255)
                                                 tCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:1904078/68‑1 (MQ=255)
                                                 tCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:3426727/68‑1 (MQ=255)
                                                 tCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:1349728/68‑1 (MQ=255)
                                                 tCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:108656/68‑1 (MQ=255)
                                                 tCAGCGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCACCTTCCTGCGTCCCGGATACTTCCGGc  <  1:1616757/68‑1 (MQ=255)
                                                     cGGCATCCGTTCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGCGTCCCGGATACTTCCGGc  <  1:1958147/64‑1 (MQ=255)
                                                              |                                                      
TGATAACCTTCAACGCTATGCGGGATCATTGGCGGCTGATTCACATAGTTCAGCGGCATCCGATCCTGCTCTTTTGGCATCCGAATGGCCCCTTCCTGTGTCCCGGATACTTCCGGC  >  minE/1426047‑1426163

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: