Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,004 A→G 41.9% K821R (AAA→AGA) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,0040AG41.9% 17.1 / 45.9 43K821R (AAA→AGA) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (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

GAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTA  >  minE/833936‑834072
                                                                    |                                                                    
gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa                                                                      >  1:4454508/1‑69 (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)
               gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa                                                                >  1:2439283/1‑60 (MQ=255)
                cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt                                                     <  1:3401193/70‑1 (MQ=255)
                cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                                                    <  1:771430/71‑1 (MQ=255)
                cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa                                                    <  1:4377260/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)
                                                                    |                                                                    
GAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTA  >  minE/833936‑834072

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: