Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,972 G→T 41.3% L810L (CTG→CTT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,9720GT41.3% 45.9 / 40.5 46L810L (CTG→CTTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (14/13);  new base T (6/13);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

ACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAG  >  minE/833907‑834038
                                                                 |                                                                  
aCGGCGACCTGGTACGAAAAAGACGACATGAATACTTAGGATATGCATCCGTTTATTCACCCGCTGTCTgc                                                               >  1:184044/1‑71 (MQ=255)
    cGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTgcggcg                                                           <  1:1316571/71‑1 (MQ=255)
      aCCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGt                                                         <  1:3340782/71‑1 (MQ=255)
       ccTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTc                                                        >  1:1082343/1‑71 (MQ=255)
       ccTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc                                                        >  1:645282/1‑71 (MQ=255)
       ccTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc                                                        >  1:523411/1‑71 (MQ=255)
        cTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCg                                                       <  1:1147388/71‑1 (MQ=255)
          ggTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGAt                                                     >  1:539908/1‑71 (MQ=255)
           gTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGAt                                                     <  1:2894652/70‑1 (MQ=255)
             aCGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCg                                                  >  1:1174350/1‑71 (MQ=255)
                 aaaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGcgt                                              <  1:2898435/71‑3 (MQ=255)
                 aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt                                              >  1:3946166/1‑71 (MQ=255)
                 aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt                                              <  1:2567963/71‑1 (MQ=255)
                  aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg                                             <  1:3714343/71‑1 (MQ=255)
                  aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTc                                                        >  1:1896657/1‑60 (MQ=255)
                  aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg                                             <  1:942542/71‑1 (MQ=255)
                  aaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTg                                             >  1:2116689/1‑71 (MQ=255)
                       cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa                                         <  1:3126815/70‑1 (MQ=255)
                       cgacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg                                        <  1:404217/71‑1 (MQ=255)
                         ataTGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa                                         <  1:2987246/66‑1 (MQ=255)
                               aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTgg                                            <  1:1288601/59‑1 (MQ=255)
                                 aCTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGa                                          <  1:3113469/59‑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:2271006/70‑1 (MQ=255)
                                                 cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa               <  1:429712/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)
                                                     tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGaa                                         >  1:2600892/1‑40 (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:921612/70‑3 (MQ=255)
                                                           tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca    <  1:781429/70‑3 (MQ=255)
                                                           tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca    <  1:292052/70‑3 (MQ=255)
                                                             cGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca    >  1:3470978/1‑67 (MQ=255)
                                                             cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg     >  1:2037592/1‑68 (MQ=255)
                                                             cGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATcgcg     >  1:1838794/1‑68 (MQ=255)
                                                              gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag  >  1:1536389/1‑70 (MQ=255)
                                                                 gTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATc                    <  1:2554677/49‑1 (MQ=255)
                                                                 |                                                                  
ACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAG  >  minE/833907‑834038

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: