Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,974 T→G 29.3% A1144A (GCT→GCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,9740TG29.3% 85.6 / 68.4 82A1144A (GCT→GCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (33/25);  new base G (18/6);  total (51/31)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGG  >  minE/834904‑835041
                                                                      |                                                                   
gTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGct                                                                     <  1:3521631/71‑1 (MQ=255)
 tGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGct                                                                     <  1:1526458/70‑1 (MQ=255)
   aaGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   >  1:4154054/1‑70 (MQ=255)
   aaGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   >  1:4765087/1‑70 (MQ=255)
            ccAAAGATCTGGGTATCGCCGTTAACGACTGGATTGAAGTCTTCAACAGCAACGGTTCTCTGACTGCCCGt                                                         <  1:618572/71‑1 (MQ=255)
            ccAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGt                                                         >  1:2396364/1‑71 (MQ=255)
             cAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGt                                                         >  1:158524/1‑70 (MQ=255)
              aaaGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGc                                                       >  1:1479942/1‑71 (MQ=255)
                aGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCgg                                                     <  1:3748388/71‑1 (MQ=255)
                aGATCTGGGTATCGCCGAAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGt                                                    >  1:133434/1‑71 (MQ=255)
                 gATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGt                                                    <  1:3713258/71‑1 (MQ=255)
                    cTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTg                                                  >  1:1986201/1‑70 (MQ=255)
                    cTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTg                                                  >  1:1525425/1‑70 (MQ=255)
                    cTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTg                                                  >  1:4659357/1‑70 (MQ=255)
                      gggTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGt                                                 >  1:3106948/1‑69 (MQ=255)
                        gTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTcagc                                             <  1:2836867/71‑1 (MQ=255)
                        gTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTcagc                                             <  1:4109995/71‑1 (MQ=255)
                         tATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGtt                                                   <  1:5235198/64‑1 (MQ=255)
                           tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:946040/46‑1 (MQ=255)
                           tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:4329922/46‑1 (MQ=255)
                           tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:2095486/46‑1 (MQ=255)
                           tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:5548519/46‑1 (MQ=255)
                           tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:5433227/46‑1 (MQ=255)
                           tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTcagccag                                          >  1:3045687/1‑71 (MQ=255)
                             gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACt                                                               >  1:3904438/3‑48 (MQ=12)
                             gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACt                                                               >  1:433713/3‑48 (MQ=12)
                             gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTca                                               >  1:2057095/3‑64 (MQ=255)
                             gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACg                                        >  1:3982122/3‑71 (MQ=255)
                             gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACg                                        >  1:5386486/3‑71 (MQ=255)
                             gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACg                                        >  1:1991539/3‑71 (MQ=255)
                             gCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGGCAGCg                                        >  1:1210870/1‑71 (MQ=255)
                             gCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCg                                        >  1:1421242/1‑71 (MQ=255)
                             gCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCg                                        >  1:4128563/1‑71 (MQ=255)
                             gCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCg                                        >  1:5591984/1‑71 (MQ=255)
                                gATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtg                                      <  1:3609679/70‑1 (MQ=255)
                                 aTAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtgta                                    <  1:1640830/71‑2 (MQ=255)
                                 aTAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGtt                                    >  1:1871631/1‑71 (MQ=255)
                                  tAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtgtac                                   >  1:25384/1‑69 (MQ=255)
                                  tAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtgtac                                   >  1:1578511/1‑69 (MQ=255)
                                  tAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTc                                   <  1:4348712/71‑1 (MQ=255)
                                   aaCGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTAcc                                  >  1:5064802/1‑71 (MQ=255)
                                   aaCGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTcc                                  >  1:4437493/1‑71 (MQ=255)
                                   aaCGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTcc                                  >  1:4115595/1‑71 (MQ=255)
                                      gACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtg                                      <  1:2911574/64‑1 (MQ=255)
                                      gACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtg                                      <  1:4456176/64‑1 (MQ=255)
                                      gACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACgtg                                      <  1:4463087/64‑1 (MQ=255)
                                      gACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGc                               <  1:3507533/71‑1 (MQ=255)
                                       aCTGGATTGAAGTCTTCAACAGAAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCa                              <  1:4126060/71‑1 (MQ=255)
                                        cTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAg                             >  1:1206532/1‑71 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGt                                                 >  1:5522397/1‑48 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:4852113/69‑1 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:4414455/69‑1 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:3770936/69‑1 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:234117/69‑1 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:3777645/69‑1 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:1166388/69‑1 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGa                          >  1:5364934/1‑71 (MQ=255)
                                           gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGa                          >  1:3396998/1‑71 (MQ=255)
                                            aTTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGt                                                    >  1:2849532/1‑44 (MQ=255)
                                               gAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGAc                      >  1:2012149/1‑71 (MQ=255)
                                                aaGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGAc                      <  1:388643/70‑1 (MQ=255)
                                                 aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGa                       >  1:1540180/1‑68 (MQ=255)
                                                 aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                    <  1:114199/71‑1 (MQ=255)
                                                 aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                    <  1:5112215/71‑1 (MQ=255)
                                                  gtattCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCa                    >  1:4559720/4‑70 (MQ=255)
                                                  gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                   >  1:1337076/1‑71 (MQ=255)
                                                  gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                   >  1:277424/1‑71 (MQ=255)
                                                  gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                   >  1:4921961/1‑71 (MQ=255)
                                                  gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                   >  1:4671006/1‑71 (MQ=255)
                                                     ttCAACGCCAACGGCGCGCTGACTGCCCGCGCGGGGGTCAGCCAACgtgt                                     <  1:4139325/50‑1 (MQ=12)
                                                     ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:3374290/59‑1 (MQ=255)
                                                     ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:828517/59‑1 (MQ=255)
                                                     ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgat                >  1:1140239/1‑71 (MQ=255)
                                                       cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGAt                         >  1:783057/1‑60 (MQ=255)
                                                       cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGAt                         >  1:1074713/1‑60 (MQ=255)
                                                        aaCAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTa             >  1:930638/1‑71 (MQ=255)
                                                          cAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTAcc           <  1:3832943/71‑1 (MQ=255)
                                                           gccAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCa          >  1:1251324/3‑71 (MQ=255)
                                                            ccAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAc         >  1:4609626/2‑71 (MQ=255)
                                                              aaCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgc       >  1:1958381/1‑71 (MQ=255)
                                                              aaCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgc       >  1:906405/1‑71 (MQ=255)
                                                              aaCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgc       >  1:828027/1‑71 (MQ=255)
                                                               aCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAg                             <  1:3711934/48‑1 (MQ=255)
                                                               aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTat            >  1:2962100/1‑64 (MQ=255)
                                                               aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTat            >  1:2708417/1‑64 (MQ=255)
                                                               aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTat            >  1:1890195/1‑64 (MQ=255)
                                                               aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgcc      >  1:3788152/1‑70 (MQ=255)
                                                                cGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATg                        <  1:2474331/52‑1 (MQ=255)
                                                                 ggTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg      >  1:4350576/1‑69 (MQ=255)
                                                                  gTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAg   <  1:1291937/71‑1 (MQ=255)
                                                                   tGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAgg  <  1:859508/71‑1 (MQ=255)
                                                                      tctGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGc                               <  1:5557475/39‑1 (MQ=255)
                                                                      |                                                                   
GTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGG  >  minE/834904‑835041

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