Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,000 G→T 40.8% A820S (GCG→TCG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0000GT40.8% 43.1 / 37.9 44A820S (GCG→TCG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base T (8/10);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGT  >  minE/833931‑834058
                                                                     |                                                          
gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                                            <  1:673802/70‑1 (MQ=255)
gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                                            <  1:2757281/70‑1 (MQ=255)
gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                                            <  1:5605469/70‑1 (MQ=255)
    tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                                                         <  1:2070046/69‑1 (MQ=255)
     gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa                                                        >  1:4454508/1‑69 (MQ=255)
      aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                                            >  1:1726496/1‑64 (MQ=255)
       aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                                    <  1:3035158/71‑1 (MQ=255)
             cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                              >  1:5187583/1‑71 (MQ=255)
              ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTggg                                             <  1:1846389/71‑1 (MQ=255)
              ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                              >  1:3689816/1‑70 (MQ=255)
                 tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGTGaaa                                          <  1:1481207/71‑1 (MQ=255)
                 tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa                                          <  1:2915590/71‑1 (MQ=255)
                 tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa                                            <  1:3322576/69‑1 (MQ=255)
                   tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGc                                                           <  1:2124322/52‑1 (MQ=255)
                    gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa                                                  >  1:2439283/1‑60 (MQ=255)
                     cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt                                       <  1:3401193/70‑1 (MQ=255)
                     cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                                      <  1:4377260/71‑1 (MQ=255)
                     cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                                      <  1:771430/71‑1 (MQ=255)
                       tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                                    <  1:78578/55‑1 (MQ=255)
                         cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaaa                                  <  1:1469837/71‑1 (MQ=255)
                         cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACa                                    >  1:4287315/1‑69 (MQ=255)
                         cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa                                    >  1:3522696/1‑69 (MQ=255)
                          gTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt                                         <  1:4112380/63‑1 (MQ=255)
                           tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa                                                         <  1:2210751/46‑1 (MQ=255)
                           tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc                                                    <  1:5515961/51‑1 (MQ=255)
                            ttATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAgg                                >  1:4407565/1‑69 (MQ=255)
                             tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                              >  1:5369519/1‑55 (MQ=255)
                             tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                              >  1:5665279/1‑55 (MQ=255)
                             tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta                              >  1:4662936/1‑68 (MQ=255)
                             tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa                              >  1:3288003/1‑71 (MQ=255)
                             tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa                              >  1:5064259/1‑71 (MQ=255)
                                   tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                        <  1:5645485/70‑3 (MQ=255)
                                      gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag                      >  1:4232274/1‑70 (MQ=255)
                                       ctttCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaa                       >  1:2376203/4‑68 (MQ=255)
                                                   gTCGATCCGGCGGGGGAATCACGCAGCGACTGGGAAATCTACaa                                   <  1:2095785/44‑1 (MQ=255)
                                                   gTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgt        <  1:593408/71‑1 (MQ=255)
                                                     cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGaaa                   >  1:3149122/1‑58 (MQ=255)
                                                      gATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGtgtg     <  1:4110679/71‑1 (MQ=255)
                                                       aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGc    <  1:1316784/71‑1 (MQ=255)
                                                       aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc                                <  1:3609045/43‑1 (MQ=255)
                                                       aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt      <  1:2696619/69‑1 (MQ=255)
                                                       aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc    >  1:217339/1‑71 (MQ=255)
                                                        tCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa                                   <  1:2132342/39‑1 (MQ=16)
                                                        tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg   >  1:1222129/1‑71 (MQ=255)
                                                        tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg   >  1:1893072/1‑71 (MQ=255)
                                                         ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt  >  1:3690035/1‑71 (MQ=255)
                                                          cGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt  <  1:5064572/70‑1 (MQ=255)
                                                            gCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc                               <  1:349508/39‑1 (MQ=255)
                                                                     |                                                          
GACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGT  >  minE/833931‑834058

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: