Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,354 T→C 16.7% R604R (CGT→CGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,3540TC16.7% 109.7 / 16.2 48R604R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base C (6/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAG  >  minE/833285‑833424
                                                                     |                                                                      
ggAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGt                                                                        >  1:1719091/1‑70 (MQ=255)
 gAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTca                                                                      >  1:3305700/1‑71 (MQ=255)
 gAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTca                                                                      >  1:1790014/1‑71 (MQ=255)
     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:1574199/71‑1 (MQ=255)
                                                          cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt             <  1:3222984/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)
                                                                     |                                                                      
GGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAG  >  minE/833285‑833424

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: