Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,398 G→T 6.0% A286S (GCT→TCT)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3980GT 6.0% 378.7 / 5.6 116A286S (GCT→TCT) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (55/54);  new base T (4/3);  total (59/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 119 total aligned reads displayed.

CGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGC  >  minE/832330‑832468
                                                                    |                                                                      
cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                                                                      >  1:2756726/1‑71 (MQ=255)
cGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCACTAAAACCATCGCCATTACCCCTGACTActct                                                                      >  1:3479817/1‑71 (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)
         aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                                                              >  1:2624974/1‑70 (MQ=255)
             cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                                          <  1:2724908/70‑1 (MQ=255)
             cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                         >  1:3395995/1‑71 (MQ=255)
             cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                         >  1:270543/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)
                      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)
                            ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGtaa                                                                  <  1:2958244/47‑1 (MQ=255)
                            ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                                  <  1:2893078/47‑1 (MQ=255)
                            ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaaa                                                                  <  1:95249/47‑1 (MQ=255)
                            ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                                      <  1:617241/59‑1 (MQ=255)
                            ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                                      <  1:1713914/59‑1 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                                             >  1:4828489/1‑50 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACTCTGAAATCGcc                                                             >  1:4834499/1‑50 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                                                                >  1:1235976/1‑47 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                                          >  1:1711066/1‑53 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                            <  1:3783935/51‑1 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                            <  1:3993389/51‑1 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                            >  1:5662090/1‑51 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                                                            <  1:638277/51‑1 (MQ=255)
                              aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggct                                         <  1:2914552/70‑1 (MQ=255)
                                aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                       >  1:2249693/1‑70 (MQ=255)
                                aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                       >  1:4894099/1‑70 (MQ=255)
                                 aaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                       <  1:4086787/69‑1 (MQ=255)
                                    ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                   >  1:1196310/1‑70 (MQ=255)
                                    ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                   >  1:4125620/1‑70 (MQ=255)
                                     gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                 <  1:2448142/71‑1 (MQ=255)
                                      aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa                                <  1:4646452/71‑1 (MQ=255)
                                       ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                                                          <  1:63704/44‑1 (MQ=255)
                                       ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                               <  1:4162029/71‑1 (MQ=255)
                                         aaaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                               >  1:1259197/1‑69 (MQ=255)
                                          aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctggc                                      <  1:2482576/61‑1 (MQ=255)
                                          aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACa                             >  1:5118343/1‑70 (MQ=255)
                                           aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                         >  1:2172988/1‑41 (MQ=255)
                                               gTTGCCGTCACACCAGACTACGCTGAAAACGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                       >  1:4448240/1‑71 (MQ=255)
                                                ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                           <  1:3690683/66‑1 (MQ=255)
                                                ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                      <  1:1910169/71‑1 (MQ=255)
                                                   ccGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt                   <  1:3196267/71‑1 (MQ=255)
                                                      tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                             <  1:894360/58‑1 (MQ=255)
                                                      tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGa                    <  1:3866823/67‑1 (MQ=255)
                                                       cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                              >  1:2958692/1‑40 (MQ=255)
                                                       cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                              >  1:2355399/1‑40 (MQ=255)
                                                       cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc                         >  1:2585835/1‑61 (MQ=255)
                                                       cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg               >  1:2149322/1‑71 (MQ=255)
                                                        acacCAGACTACGCTGAAATCGCCAAACTGTGCGAt                                                 >  1:3032130/1‑36 (MQ=255)
                                                        acacCAGACTACGCTGAAATCGCCAAACTGTGCGAt                                                 >  1:4251374/1‑36 (MQ=255)
                                                        acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                               <  1:2837562/54‑1 (MQ=255)
                                                        acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                               <  1:2830678/54‑1 (MQ=255)
                                                        acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                           >  1:592428/1‑58 (MQ=255)
                                                        acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGc              <  1:5383525/71‑1 (MQ=255)
                                                          acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCaa            <  1:3964059/71‑1 (MQ=255)
                                                           ccAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                       <  1:4590889/43‑1 (MQ=255)
                                                           ccAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAAc                              <  1:5030950/52‑1 (MQ=255)
                                                            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                          >  1:169297/1‑39 (MQ=255)
                                                            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGTACCGAAAc                              >  1:532683/1‑51 (MQ=255)
                                                            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAAc                              >  1:3821236/1‑51 (MQ=255)
                                                            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAAc                              >  1:4357905/1‑51 (MQ=255)
                                                            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg          >  1:4531944/1‑71 (MQ=255)
                                                            cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg          >  1:5291713/1‑71 (MQ=255)
                                                             aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg         >  1:1044184/1‑71 (MQ=255)
                                                             aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg         <  1:3169226/71‑1 (MQ=255)
                                                             aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg         <  1:4082993/71‑1 (MQ=255)
                                                              gACTACTCTGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGaaa                               >  1:5658483/1‑48 (MQ=19)
                                                              gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc        <  1:2195511/71‑1 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                      <  1:2641041/56‑1 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                      <  1:2056672/56‑1 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                      <  1:2506091/56‑1 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                      <  1:4579405/56‑1 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg       >  1:3504627/1‑71 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg       >  1:4369683/1‑71 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg       >  1:3295111/1‑71 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg       >  1:2937726/1‑71 (MQ=255)
                                                               aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc        >  1:5640309/1‑70 (MQ=255)
                                                                 tACTCTGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCAct                      <  1:2780895/54‑2 (MQ=19)
                                                                 tACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                       <  1:2699977/37‑1 (MQ=255)
                                                                 tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                 <  1:667677/43‑1 (MQ=255)
                                                                 tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                 <  1:712533/43‑1 (MQ=255)
                                                                 tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt     >  1:1171134/1‑71 (MQ=255)
                                                                   cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg   >  1:4875849/1‑71 (MQ=255)
                                                                   cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg   <  1:2331287/71‑1 (MQ=255)
                                                                    |                                                                      
CGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGC  >  minE/832330‑832468

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: