Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,919 A→G 6.5% K1126R (AAA→AGA) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,9190AG 6.5% 175.8 / 3.7 61K1126R (AAA→AGA) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (19/38);  new base G (2/2);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCG  >  minE/834849‑834989
                                                                      |                                                                      
aGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCaa                                                                        <  1:3696262/71‑1 (MQ=255)
  cGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAg                                                                      <  1:4030723/71‑1 (MQ=255)
  cGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAg                                                                      <  1:2487281/71‑1 (MQ=255)
   gACAACCTGCTGATGCTGACTTTAGGTCGCGGGGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAg                                                                      <  1:604000/70‑1 (MQ=255)
      aaCCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCt                                                                  >  1:55885/1‑71 (MQ=255)
        ccTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTgg                                                                >  1:939944/1‑71 (MQ=255)
            ctgatgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGt                                                              <  1:2698323/69‑1 (MQ=255)
              gatgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCg                                                          >  1:3245974/1‑71 (MQ=255)
               atgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATc                                                           <  1:2356704/69‑1 (MQ=255)
               atgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATc                                                           <  1:660302/69‑1 (MQ=255)
                tgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGcc                                                        >  1:94788/1‑71 (MQ=255)
                tgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGcc                                                        >  1:1517154/1‑71 (MQ=255)
                tgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGcc                                                        <  1:1935477/71‑1 (MQ=255)
                tgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGcc                                                        >  1:3241079/1‑71 (MQ=255)
                tgctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGc                                                         >  1:3506100/1‑70 (MQ=255)
                 gctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGcc                                                        >  1:3800049/1‑70 (MQ=255)
                  ctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGc                                                         >  1:1744017/1‑68 (MQ=255)
                  ctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGa                                                      >  1:101767/1‑71 (MQ=255)
                  ctgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGa                                                      >  1:3017471/1‑71 (MQ=255)
                   tgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAg                                                                      <  1:332618/54‑1 (MQ=255)
                   tgaCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAg                                                                      <  1:532098/54‑1 (MQ=255)
                     aCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGAt                                                     <  1:2692664/69‑1 (MQ=255)
                     aCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGAt                                                     <  1:2235862/69‑1 (MQ=255)
                     aCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGAt                                                     <  1:1084031/69‑1 (MQ=255)
                     aCTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATaa                                                   <  1:2160717/71‑1 (MQ=255)
                      cTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAAc                                                  <  1:792321/71‑1 (MQ=255)
                             tCGCGGTGGACCGATTGTCTGGATCAGCGAAACAGATGCCCGAGAACTGACCATTGTCGATAACGACTggg                                           >  1:4009463/1‑70 (MQ=255)
                              cgcgGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTg                                                                 >  1:1646543/1‑48 (MQ=255)
                                cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTg                                        <  1:3708212/71‑1 (MQ=255)
                                cggtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTg                                        <  1:1986199/71‑1 (MQ=255)
                                  gtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAg                                                                      <  1:2596766/39‑1 (MQ=255)
                                  gtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGAt                                                     <  1:1035404/56‑1 (MQ=255)
                                  gtggtccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGAt                                                     <  1:719589/56‑1 (MQ=255)
                                       ccggtggtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGaa                                      >  1:652050/1‑66 (MQ=255)
                                            tgtcTGGATCAGGGAAACAGATGCCCGAGAACTGACCATTGTCGATAACGACTGGGTGGAAGTATTCAAcg                            >  1:1531574/2‑70 (MQ=255)
                                             gtcTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAg                           <  1:2990809/71‑1 (MQ=255)
                                              tcTGGATCAGCGAAACAGATGCCCGAGAACTGACCATTGTCGATAACGACTGGGTGGAAGTATTCAAcgcc                          <  1:1700512/71‑4 (MQ=255)
                                              tcTGGATCAGCGAAACAGATGCCCGAGAACTGACCATTGTCGATAACGACTGGGTGGAAGTATTCAAcgcc                          <  1:1218425/71‑4 (MQ=255)
                                                tGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATc                                                           <  1:2242557/36‑1 (MQ=255)
                                                tGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACg                                                 <  1:2513791/46‑1 (MQ=255)
                                                tGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTgg                                            <  1:240110/51‑1 (MQ=255)
                                                tGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTg                                        <  1:3865997/55‑1 (MQ=255)
                                                  gTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGc                                                         <  1:3940459/36‑1 (MQ=255)
                                                  gTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGaa                                      >  1:3610945/1‑55 (MQ=255)
                                                  gTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCATAACGACTGGATTGAAGTCTTCAACAGCAAc                       <  1:1318900/69‑1 (MQ=255)
                                                   ttGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCg                                                       <  1:701745/37‑1 (MQ=255)
                                                    tGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGt                    >  1:3302810/1‑71 (MQ=255)
                                                     gAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGAt                                                     >  1:3626790/1‑37 (MQ=255)
                                                     gAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGt                    <  1:993559/70‑1 (MQ=255)
                                                       gTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGc                  >  1:2477386/1‑70 (MQ=255)
                                                        tGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCaa                              <  1:594398/57‑1 (MQ=255)
                                                          acgCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct               <  1:3687096/68‑1 (MQ=255)
                                                          aaGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGc                          <  1:2599645/59‑1 (MQ=255)
                                                          aaGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct               <  1:2021943/70‑1 (MQ=255)
                                                           aGCCGATGCCAAAGATCTGGGTATCGCCGATAACGa                                                <  1:2385670/36‑1 (MQ=255)
                                                            gCCGATGCCAAAGATCTGGGTATCGCCGATAACGAc                                               >  1:1357909/1‑36 (MQ=255)
                                                            gCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGct                 <  1:1768966/66‑1 (MQ=255)
                                                                aTGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct               <  1:4003014/64‑1 (MQ=255)
                                                                aTGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGctct               <  1:3229268/64‑1 (MQ=255)
                                                                aTGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGcc        <  1:414273/71‑1 (MQ=255)
                                                                  gCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCg      >  1:3341898/1‑71 (MQ=255)
                                                                    cAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTg    <  1:4058024/71‑1 (MQ=255)
                                                                     aaaGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGAc            >  1:1402149/1‑62 (MQ=255)
                                                                      aaGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCg  <  1:1399246/71‑1 (MQ=255)
                                                                      |                                                                      
AGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAGATCTGGGTATCGCCGATAACGACTGGATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCG  >  minE/834849‑834989

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: