Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,005 A→C 41.9% K821N (AAA→AAC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0050AC41.9% 3.2 / 53.0 43K821N (AAA→AAC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base C (8/10);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAG  >  minE/833938‑834075
                                                                   |                                                                      
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)
             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)
                    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:1893072/1‑71 (MQ=255)
                                                 tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg                    >  1:1222129/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)
                                                               cGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTggg       >  1:4726129/1‑70 (MQ=255)
                                                                gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGa                            <  1:5235413/48‑1 (MQ=255)
                                                                gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGt      >  1:2200177/1‑70 (MQ=255)
                                                                  aaaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATtctc                               >  1:3297926/1‑43 (MQ=255)
                                                                  aaaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa     >  1:1998257/1‑69 (MQ=255)
                                                                   aaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa     >  1:3336904/1‑68 (MQ=255)
                                                                   aaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAg  <  1:3145946/71‑1 (MQ=255)
                                                                   aaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAg  <  1:2508654/71‑1 (MQ=255)
                                                                   |                                                                      
ATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAAG  >  minE/833938‑834075

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: