Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,971 T→C 25.8% G1143G (GGT→GGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,9710TC25.8% 148.6 / 56.7 93G1143G (GGT→GGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (18/6);  total (56/37)
Fisher's exact test for biased strand distribution p-value = 9.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGG  >  minE/834902‑835041
                                                                     |                                                                      
gAGTGAAGCCGATGCCAAAGATCTTGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTg                                                                       >  1:4234490/1‑71 (MQ=255)
gAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTg                                                                       >  1:4294503/1‑71 (MQ=255)
gAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTg                                                                       >  1:2770190/1‑71 (MQ=255)
  gTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGct                                                                     <  1:3521631/71‑1 (MQ=255)
   tGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGct                                                                     <  1:1526458/70‑1 (MQ=255)
     aaGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   >  1:4765087/1‑70 (MQ=255)
     aaGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   >  1:4154054/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:4659357/1‑70 (MQ=255)
                      cTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTg                                                  >  1:1525425/1‑70 (MQ=255)
                      cTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTg                                                  >  1:1986201/1‑70 (MQ=255)
                        gggTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGt                                                 >  1:3106948/1‑69 (MQ=255)
                          gTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTcagc                                             <  1:4109995/71‑1 (MQ=255)
                          gTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTcagc                                             <  1:2836867/71‑1 (MQ=255)
                           tATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGtt                                                   <  1:5235198/64‑1 (MQ=255)
                             tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:4329922/46‑1 (MQ=255)
                             tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:946040/46‑1 (MQ=255)
                             tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:2095486/46‑1 (MQ=255)
                             tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:5433227/46‑1 (MQ=255)
                             tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct                                                                   <  1:5548519/46‑1 (MQ=255)
                             tCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTcagccag                                          >  1:3045687/1‑71 (MQ=255)
                               gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACt                                                               >  1:433713/3‑48 (MQ=12)
                               gtcGATAACGACTGGGTGGAAGTATTCAACGCCAACGGCGCGCTGACt                                                               >  1:3904438/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:4128563/1‑71 (MQ=255)
                               gCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCg                                        >  1:5591984/1‑71 (MQ=255)
                               gCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCg                                        >  1:1421242/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:4115595/1‑71 (MQ=255)
                                     aaCGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTcc                                  >  1:4437493/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:4414455/69‑1 (MQ=255)
                                             gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:3770936/69‑1 (MQ=255)
                                             gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAgg                            <  1:4852113/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)
                                                                     |                                                                      
GAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGG  >  minE/834902‑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: