Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,355 G→A 20.0% V271V (GTG→GTA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3550GA20.0% 215.9 / 44.2 105V271V (GTG→GTAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (46/38);  new base A (14/7);  total (60/45)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 105 total aligned reads displayed.

CATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTG  >  minE/832289‑832424
                                                                  |                                                                     
cctcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt                                                                    >  1:3551809/3‑70 (MQ=255)
catcatCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACgc                                                                    <  1:4611402/70‑2 (MQ=255)
catcatCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt                                                                    >  1:2484509/1‑70 (MQ=255)
 atcatcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTa                                                                  >  1:656725/1‑71 (MQ=255)
   catcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTAc                                                                 >  1:4222716/1‑70 (MQ=255)
   catcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTAc                                                                 >  1:5535609/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)
     tcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACAGAAGTACGCTACaaa                                                              >  1:2735956/1‑71 (MQ=255)
     tcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaaa                                                              <  1:3325938/71‑1 (MQ=255)
      cGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAg                                                             >  1:1353362/1‑71 (MQ=255)
       gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAg                                                             <  1:2757397/70‑1 (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)
             gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                                      <  1:3847510/71‑1 (MQ=255)
             gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                                      <  1:1378252/71‑1 (MQ=255)
             gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                                      <  1:243501/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)
                  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)
                    aCGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAAc                                                   >  1:1477693/1‑67 (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:2328699/1‑71 (MQ=255)
                      gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGc                                             >  1:654801/1‑71 (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:3479817/1‑71 (MQ=255)
                                         cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                          >  1:2756726/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:5512562/71‑1 (MQ=255)
                                          ggACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTACTCTg                         <  1:2691524/71‑1 (MQ=255)
                                           gATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTg                         >  1:689125/1‑70 (MQ=255)
                                            aTGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                       <  1:4310725/71‑1 (MQ=255)
                                             tGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                       >  1:3174580/1‑70 (MQ=255)
                                              gCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaaa                                                     <  1:4119077/39‑1 (MQ=255)
                                              gCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAAc                                                   <  1:1698369/41‑1 (MQ=255)
                                               cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                        >  1:5423924/1‑51 (MQ=255)
                                               cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                        >  1:3486655/1‑51 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTAGAAAGGGACCAAAAc                                                   >  1:5014385/1‑39 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                 >  1:2005069/1‑41 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                 >  1:3764434/1‑41 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                 >  1:4794281/1‑41 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCg                   >  1:3303168/1‑71 (MQ=255)
                                                 cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                              <  1:5307474/59‑1 (MQ=255)
                                                 cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTa                              <  1:3966909/59‑1 (MQ=255)
                                                 cACTTCTTTACCGAAGTACGCTACAAAGGCACTaaaa                                                    >  1:5687904/1‑37 (MQ=16)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGa                                 >  1:5292437/1‑55 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                  >  1:2640918/1‑70 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                  >  1:2624974/1‑70 (MQ=255)
                                                   cttcttTACTGAAGTGCGTTACAAAGGGACCAAAAc                                                   >  1:4810934/1‑36 (MQ=255)
                                                    ttcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGcc                                            <  1:2940867/42‑1 (MQ=255)
                                                     tcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa              >  1:2717180/1‑71 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa              <  1:2724908/70‑1 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc             >  1:270543/1‑71 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc             >  1:1454583/1‑71 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc             >  1:3395995/1‑71 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACAAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                    >  1:3524381/1‑64 (MQ=255)
                                                         tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtg           <  1:2003354/70‑1 (MQ=255)
                                                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa               <  1:4275270/65‑1 (MQ=255)
                                                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg         >  1:5422094/1‑71 (MQ=255)
                                                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg         >  1:2901675/1‑71 (MQ=255)
                                                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt          >  1:88103/1‑70 (MQ=255)
                                                           cTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGc        <  1:5286533/71‑1 (MQ=255)
                                                             gAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa      >  1:5670754/1‑71 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc             >  1:5057104/1‑63 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg       <  1:1822553/69‑1 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg       <  1:1851614/69‑1 (MQ=255)
                                                               aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc    <  1:1794494/71‑1 (MQ=255)
                                                               aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc    >  1:1294859/1‑71 (MQ=255)
                                                                 tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa      <  1:517313/67‑1 (MQ=255)
                                                                 tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg  >  1:2701369/1‑71 (MQ=255)
                                                                  |                                                                     
CATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTG  >  minE/832289‑832424

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: