Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,000 G→T 31.0% 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,0000GT31.0% 79.7 / 24.9 45A820S (GCG→TCG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base T (4/10);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGG  >  minE/833930‑834070
                                                                      |                                                                      
cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                                                        <  1:404217/71‑1 (MQ=255)
             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:2048209/70‑3 (MQ=255)
                                    tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                    <  1:292052/70‑3 (MQ=255)
                                    tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                    <  1:921612/70‑3 (MQ=255)
                                    tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca                                    <  1:781429/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:2063502/63‑1 (MQ=255)
                                           tCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg                                     <  1:2575807/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)
                                                                      |                                                                      
CGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGTCGGCCATCTGGG  >  minE/833930‑834070

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: