Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,004 A→G 28.6% 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,0040AG28.6% 66.5 / 33.9 49K821R (AAA→AGA) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base G (4/10);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAA  >  minE/833943‑834074
                                                             |                                                                      
tCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTg                                                               >  1:2156777/1‑71 (MQ=255)
 cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                              >  1:681132/1‑71 (MQ=255)
 cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg                                                              >  1:3257707/1‑71 (MQ=255)
  ggATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTgg                                                              <  1:1921612/70‑1 (MQ=255)
   gATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa                                                            <  1:1009121/71‑1 (MQ=255)
    atatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaa                                                           <  1:110971/71‑1 (MQ=255)
     tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa                                                          <  1:2340387/71‑1 (MQ=255)
     tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaaa                                                          <  1:2129344/71‑1 (MQ=255)
             cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa                                                   <  1:429712/70‑1 (MQ=255)
             cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa                                                   <  1:2271006/70‑1 (MQ=255)
                ttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc                                                >  1:583890/1‑70 (MQ=255)
                 tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta                                              >  1:1205943/1‑68 (MQ=255)
                 tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggt                                               <  1:1102742/70‑3 (MQ=255)
                   ttCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATc                                            >  1:1281837/1‑71 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggt                                               <  1:1713221/63‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                        <  1:921612/70‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                        <  1:781429/70‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                        <  1:2048209/70‑3 (MQ=255)
                       tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                        <  1:292052/70‑3 (MQ=255)
                         cGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                        >  1:3470978/1‑67 (MQ=255)
                         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                         >  1:1838794/1‑68 (MQ=255)
                         cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                         >  1:2037592/1‑68 (MQ=255)
                          gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag                                      >  1:1536389/1‑70 (MQ=255)
                             gTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATc                                                        <  1:2554677/49‑1 (MQ=255)
                              tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                         <  1:2575807/63‑1 (MQ=255)
                              tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                         <  1:2063502/63‑1 (MQ=255)
                                     cggTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa                          <  1:1211621/71‑1 (MQ=255)
                                       gTCGATACGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAg                         >  1:210192/1‑70 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAACGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAAt                                  <  1:1941501/60‑1 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGaaa                                                          <  1:3122900/36‑1 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCg                            <  1:1616871/66‑1 (MQ=255)
                                        tCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtg                       >  1:2782389/1‑71 (MQ=255)
                                         cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                         <  1:2342961/52‑1 (MQ=255)
                                           aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaaa                                      <  1:3020136/53‑2 (MQ=255)
                                            tCCGGCCTTGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg                   <  1:945227/71‑1 (MQ=255)
                                             ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                  >  1:3054163/1‑71 (MQ=255)
                                             ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt                  >  1:1433009/1‑71 (MQ=255)
                                                  cTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATtct                               >  1:3820388/1‑53 (MQ=255)
                                                   ctggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc                    >  1:2525992/3‑63 (MQ=255)
                                                     ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc          >  1:1586483/1‑71 (MQ=255)
                                                     ggAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATc          >  1:597351/1‑71 (MQ=255)
                                                      gAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCt         <  1:2591893/71‑1 (MQ=255)
                                                       aaGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTg        <  1:1974148/71‑1 (MQ=255)
                                                        aGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGaa                          >  1:1634938/1‑52 (MQ=255)
                                                        aGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg       >  1:2398472/1‑71 (MQ=255)
                                                         gCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTggg      <  1:1754002/71‑1 (MQ=255)
                                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTgg       <  1:2087617/68‑1 (MQ=255)
                                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa    <  1:3754192/71‑1 (MQ=255)
                                                           gAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTa    <  1:1517705/71‑1 (MQ=255)
                                                             aaaGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTaaa  <  1:3922874/71‑1 (MQ=255)
                                                             |                                                                      
TCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGGTAAA  >  minE/833943‑834074

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: