Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,334 T→C 37.4% D264D (GAT→GAC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3340TC37.4% 76.4 / 96.7 99D264D (GAT→GACnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base C (19/18);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 105 total aligned reads displayed.

CTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAA  >  minE/832264‑832404
                                                                      |                                                                      
ctgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGAt                                                                        >  1:2694472/1‑71 (MQ=255)
 tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg                                                                       >  1:4384016/1‑71 (MQ=255)
  gctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGAt                                                                        >  1:1039279/1‑69 (MQ=255)
  gctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCACGGAt                                                                        >  1:2394206/1‑69 (MQ=255)
        tGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCAct                                                                 >  1:3605696/1‑70 (MQ=255)
         ggTACAACTCCAGCTACATCATCGCCTGGGGGGCAAACGTGCCGCAGACGCGTACCCCGGATGCTCAct                                                                 <  1:350138/69‑1 (MQ=255)
          gTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCActtct                                                              >  1:3538926/1‑71 (MQ=255)
          gTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCActtc                                                               <  1:3426758/70‑1 (MQ=255)
          gTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                                                                <  1:1521/69‑1 (MQ=255)
           tACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                             <  1:2389662/71‑1 (MQ=255)
            aCAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                             >  1:5534822/1‑70 (MQ=255)
             taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCttt                                                            <  1:4595201/69‑1 (MQ=255)
             taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCttt                                                            <  1:1688819/69‑1 (MQ=255)
             taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTa                                                           >  1:5046981/2‑71 (MQ=255)
              aaCTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCt                                                                     <  1:2662608/60‑1 (MQ=255)
              aaCTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCt                                                                     <  1:4730690/60‑1 (MQ=255)
               aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCTCAGACACGTACGCCGGACGCCCACTTCTTTAcc                                                         >  1:3817375/1‑70 (MQ=255)
                 tCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCActt                                                                <  1:1938902/62‑1 (MQ=255)
                 tCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTAccg                                                        <  1:2584830/70‑3 (MQ=255)
                  ccAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGAt                                                                        >  1:1726886/1‑53 (MQ=255)
                   cAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGa                                                       <  1:650272/69‑1 (MQ=255)
                   cAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGa                                                       <  1:4799217/69‑1 (MQ=255)
                         cctcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt                                                >  1:3551809/3‑70 (MQ=255)
                         catcatCGCCTGGGGTTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtc                                                               >  1:4994262/1‑55 (MQ=255)
                         catcatCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACgc                                                <  1:4611402/70‑2 (MQ=255)
                         catcatCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt                                                >  1:2484509/1‑70 (MQ=255)
                          atcatcGCCTGGGGGTCTATCGTACCGCAGACACGTACGCCGGACGCCCActt                                                                <  1:134433/53‑1 (MQ=255)
                          atcatcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTa                                              >  1:656725/1‑71 (MQ=255)
                           tcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTa                                                           >  1:2728192/1‑57 (MQ=255)
                            catcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTAc                                             >  1:5535609/1‑70 (MQ=255)
                            catcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTAc                                             >  1:4222716/1‑70 (MQ=255)
                             atcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACa                                            >  1:4172383/1‑70 (MQ=255)
                             atcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaa                                           <  1:2475442/71‑1 (MQ=255)
                             atcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                            >  1:1620070/1‑70 (MQ=255)
                             atcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                            >  1:3839864/1‑70 (MQ=255)
                             atcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                            >  1:2794444/1‑70 (MQ=255)
                              tcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTAccg                                                        <  1:812926/57‑3 (MQ=255)
                              tcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACAGAAGTACGCTACaaa                                          >  1:2735956/1‑71 (MQ=255)
                              tcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtc                                                               <  1:2392652/50‑1 (MQ=255)
                              tcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaaa                                          <  1:3325938/71‑1 (MQ=255)
                               cGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAg                                         >  1:1353362/1‑71 (MQ=255)
                                gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAg                                         <  1:5040542/70‑1 (MQ=255)
                                     ggggTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCActa                                   <  1:2941623/71‑3 (MQ=255)
                                     ggggTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCActa                                   <  1:3623406/71‑3 (MQ=255)
                                      gggTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTaa                                  <  1:5529506/71‑1 (MQ=255)
                                      gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCAc                                                                  >  1:4414876/1‑39 (MQ=255)
                                      gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                  <  1:243501/71‑1 (MQ=255)
                                      gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                  <  1:3847510/71‑1 (MQ=255)
                                      gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                  <  1:1378252/71‑1 (MQ=255)
                                       ggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaaa                                 <  1:2029313/71‑1 (MQ=255)
                                         tctaaCGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAAc                               <  1:5613126/68‑1 (MQ=255)
                                          cAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                        >  1:113967/1‑61 (MQ=255)
                                           aaaCGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                             <  1:5332486/39‑1 (MQ=255)
                                           aaaCGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                             <  1:3573448/39‑1 (MQ=255)
                                           aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGaa                                                      >  1:1881890/1‑46 (MQ=255)
                                           aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                        <  1:3776140/60‑1 (MQ=255)
                                           aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                             >  1:3052503/1‑71 (MQ=255)
                                            aaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaaa                                 <  1:1482629/66‑1 (MQ=255)
                                             aCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGtc                           <  1:3143133/71‑2 (MQ=255)
                                             aCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGt                            <  1:164303/70‑1 (MQ=255)
                                             aCGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAg                                                     >  1:464371/1‑45 (MQ=255)
                                             aCGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAg                                                     >  1:2466735/1‑45 (MQ=255)
                                             aCGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAAc                               >  1:1477693/1‑67 (MQ=255)
                                              cGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                             <  1:3762885/36‑1 (MQ=255)
                                              cGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                               >  1:270907/1‑66 (MQ=255)
                                              cGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                          >  1:824326/1‑71 (MQ=255)
                                               gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                          >  1:1725713/1‑70 (MQ=255)
                                               gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTg                          >  1:2942999/1‑70 (MQ=255)
                                               gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGc                         >  1:654801/1‑71 (MQ=255)
                                                tGCCGCAGACGCGTACCCCGGATGCTCACTTCtttt                                                           >  1:3518208/1‑35 (MQ=255)
                                                tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTa                                                           >  1:4410367/1‑36 (MQ=255)
                                                tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGcc                        >  1:350784/1‑71 (MQ=255)
                                                 gCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCg                       >  1:2668483/1‑71 (MQ=255)
                                                  ccGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                      <  1:5044771/71‑1 (MQ=255)
                                                  ccGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                      <  1:2411531/71‑1 (MQ=255)
                                                      aGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaaa                                          >  1:5423738/1‑47 (MQ=255)
                                                         cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa               <  1:568557/71‑1 (MQ=255)
                                                         cacgTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTaa                                  >  1:2967376/3‑52 (MQ=255)
                                                          gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                 >  1:4665990/1‑68 (MQ=255)
                                                            gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                             >  1:482381/1‑54 (MQ=255)
                                                            gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                      <  1:4666541/61‑1 (MQ=255)
                                                             tACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                            <  1:2496207/38‑1 (MQ=255)
                                                              aCCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                            <  1:4415845/37‑1 (MQ=255)
                                                               cgccGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTa          >  1:1199696/3‑70 (MQ=255)
                                                               ccccGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                            >  1:2666763/1‑36 (MQ=255)
                                                               ccccGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                               <  1:1543223/49‑1 (MQ=255)
                                                               ccccGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                               <  1:3603495/49‑1 (MQ=255)
                                                                 ccGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAAcca                             >  1:5219047/1‑47 (MQ=255)
                                                                  cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct      >  1:2756726/1‑71 (MQ=255)
                                                                  cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct      >  1:3479817/1‑71 (MQ=255)
                                                                  cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActc       >  1:259288/1‑68 (MQ=255)
                                                                   ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                      <  1:4998232/54‑1 (MQ=255)
                                                                   ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGt                      <  1:284573/54‑1 (MQ=255)
                                                                   ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGc       >  1:1864301/1‑69 (MQ=255)
                                                                   ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg     <  1:2691524/71‑1 (MQ=255)
                                                                   ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg     <  1:5512562/71‑1 (MQ=255)
                                                                    gATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTg     >  1:689125/1‑70 (MQ=255)
                                                                     aTGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa   <  1:4310725/71‑1 (MQ=255)
                                                                      tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa  <  1:2307398/71‑1 (MQ=255)
                                                                      tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa   >  1:3174580/1‑70 (MQ=255)
                                                                      |                                                                      
CTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAA  >  minE/832264‑832404

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: