Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,357 C→A 15.9% H605Q (CAC→CAA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,3570CA15.9% 153.0 / 10.4 50H605Q (CAC→CAAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base A (6/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAG  >  minE/833290‑833426
                                                                   |                                                                     
aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg                                                                    <  1:2975949/71‑1 (MQ=255)
aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg                                                                    <  1:1671413/71‑1 (MQ=255)
     cGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACa                                                                <  1:2215143/70‑1 (MQ=255)
         cGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCAct                                                           >  1:3167630/1‑71 (MQ=255)
            aaaCCGGCTGTCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                        >  1:3302308/1‑71 (MQ=255)
            aaaCCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCAcctc                                                         >  1:2172934/1‑70 (MQ=255)
            aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                        <  1:1320267/71‑1 (MQ=255)
              aCCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCAcctcgtt                                                      >  1:1776245/1‑68 (MQ=255)
                cGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAtt                                                    <  1:152361/71‑1 (MQ=255)
                  gCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAttta                                                  >  1:2680760/1‑70 (MQ=255)
                      gCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                               <  1:1040045/70‑1 (MQ=255)
                       ttgCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACaa                                             <  1:2131717/69‑1 (MQ=255)
                           cGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATaa                                             >  1:2773946/1‑67 (MQ=255)
                            gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACt                                        >  1:775826/1‑71 (MQ=255)
                              tGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTc                                       <  1:3227087/70‑1 (MQ=255)
                                gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag                                    >  1:229897/4‑71 (MQ=255)
                                gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag                                    >  1:3170020/4‑71 (MQ=255)
                                 cGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagc                                   <  1:254711/71‑1 (MQ=255)
                                   tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                               <  1:613701/57‑1 (MQ=255)
                                   tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                               <  1:2309997/57‑1 (MQ=255)
                                     tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAt                                                     >  1:597055/1‑49 (MQ=255)
                                     tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTcca                                     >  1:3372621/1‑65 (MQ=255)
                                     tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                               >  1:1590333/1‑71 (MQ=255)
                                       cccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCa                                          <  1:184800/58‑1 (MQ=255)
                                        ccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCta                                                >  1:188770/1‑51 (MQ=255)
                                         attGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTa                                                      >  1:999774/2‑44 (MQ=255)
                                          ttGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                               <  1:832556/50‑1 (MQ=255)
                                           tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                        >  1:2874469/1‑40 (MQ=255)
                                           tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct                                                        >  1:2879497/1‑40 (MQ=255)
                                           tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCAc                                                            >  1:2821326/1‑36 (MQ=255)
                                            gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat                                               >  1:22360/1‑48 (MQ=255)
                                            gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGtt                         <  1:1186222/70‑1 (MQ=255)
                                            gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa                        >  1:3239879/1‑71 (MQ=255)
                                                 gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                   >  1:617901/1‑71 (MQ=255)
                                                 gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa                   >  1:2589941/1‑71 (MQ=255)
                                                    gCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                               <  1:1024171/56‑1 (MQ=255)
                                                     cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt               <  1:3222984/71‑1 (MQ=255)
                                                     cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt               <  1:1574199/71‑1 (MQ=255)
                                                      gTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTc              <  1:2144655/71‑1 (MQ=255)
                                                         cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                               >  1:2523068/1‑51 (MQ=255)
                                                         cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt                               >  1:685868/1‑51 (MQ=255)
                                                           gcgcGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggc         <  1:3085662/71‑1 (MQ=255)
                                                            cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc                                  <  1:529854/45‑1 (MQ=255)
                                                            cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag                >  1:2353032/1‑60 (MQ=255)
                                                            cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag                >  1:1769060/1‑60 (MQ=255)
                                                             gcgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg       <  1:1408708/71‑1 (MQ=255)
                                                                tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa                    <  1:2789547/55‑1 (MQ=255)
                                                                tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg       >  1:4022660/1‑68 (MQ=255)
                                                                tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg           <  1:593061/64‑1 (MQ=255)
                                                                tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAag    <  1:3654739/71‑1 (MQ=255)
                                                                 cacaTGAACAGCACTTCTTATTTCTATAACCACTccagcc                                  <  1:1225315/40‑1 (MQ=255)
                                                                  acaTGAACAGCACTTCTTATTTCTATAACCACTcc                                      >  1:4035666/1‑35 (MQ=255)
                                                                  acaTGAACAGCACTTCTTATTTCTATAACCACTcc                                      >  1:1666908/1‑35 (MQ=255)
                                                                  acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag  <  1:3111514/71‑1 (MQ=255)
                                                                   |                                                                     
AACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAG  >  minE/833290‑833426

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