Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,358 T→C 6.7% R272R (CGT→CGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3580TC 6.7% 206.1 / 4.5 74R272R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base C (3/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGC  >  minE/832289‑832428
                                                                     |                                                                      
catcatCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt                                                                        <  1:866089/70‑1 (MQ=255)
 atcatcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTa                                                                      <  1:1878824/71‑1 (MQ=255)
 atcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTa                                                                      <  1:3473666/71‑1 (MQ=255)
 atcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTa                                                                      <  1:2476045/71‑1 (MQ=255)
    atcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                                                    <  1:3042733/70‑1 (MQ=255)
      cGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAg                                                                 >  1:1611567/1‑71 (MQ=255)
       gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                                                <  1:3844702/71‑1 (MQ=255)
       gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                                                >  1:2569466/1‑71 (MQ=255)
       gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                                                <  1:531809/71‑1 (MQ=255)
       gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                                                <  1:84732/71‑1 (MQ=255)
        ccTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAggg                                                               <  1:2463365/71‑1 (MQ=255)
          tGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCAc                                                             >  1:3366837/1‑71 (MQ=255)
          tGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCAc                                                             >  1:3325524/1‑71 (MQ=255)
          tGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGAc                                                             >  1:467969/1‑71 (MQ=255)
           gggggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTAAAAGGGAcc                                                            >  1:3904245/1‑70 (MQ=255)
            ggggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCa                                                           <  1:134540/71‑1 (MQ=255)
             gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                                          <  1:353513/71‑1 (MQ=255)
             gggTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCa                                                           >  1:317598/1‑70 (MQ=255)
               gTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAgg                                                                <  1:3732973/63‑1 (MQ=255)
               gTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaaaa                                                        <  1:759853/71‑1 (MQ=255)
                  aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                     >  1:3450809/1‑71 (MQ=255)
                  aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                                                     >  1:838937/1‑71 (MQ=255)
                       tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGCACCAAAACTGTTGcc                                                <  1:1592419/71‑1 (MQ=255)
                         ccGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCg                                               >  1:2477889/1‑70 (MQ=255)
                            cAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                                            >  1:1664343/1‑70 (MQ=255)
                                cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                                       <  1:108177/71‑1 (MQ=255)
                                cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                                       <  1:3465792/71‑1 (MQ=255)
                                 gcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                         >  1:3088918/1‑68 (MQ=255)
                                  cgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGa                                     <  1:479214/71‑1 (MQ=255)
                                   gTACGCCGGACGCCCACTTCTTTACCGAACTACGCTACAAAGGCACAAAAAcc                                                      >  1:1529866/1‑52 (MQ=255)
                                   gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAggg                                                               <  1:3381511/44‑1 (MQ=255)
                                   gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGa                                                              >  1:591485/1‑45 (MQ=255)
                                   gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGa                                                              >  1:1628907/1‑45 (MQ=255)
                                   gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGAcc                                                            <  1:2224803/47‑1 (MQ=255)
                                    tACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaa                                                                   <  1:420838/39‑1 (MQ=255)
                                    tACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACt                                   <  1:217214/71‑1 (MQ=255)
                                     aCCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaa                                                                   >  1:2124047/1‑38 (MQ=255)
                                     aCCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGa                                     >  1:894911/1‑68 (MQ=255)
                                      ccccGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                                                    <  1:3627714/36‑1 (MQ=255)
                                       cccGGATGCTCACTTCTTTACTGAAGTGCGTTACaa                                                                   <  1:3251747/36‑1 (MQ=255)
                                         cGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGt                                                    >  1:1109946/1‑49 (MQ=255)
                                         cGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCt                              <  1:223840/71‑1 (MQ=255)
                                          ggATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGAcc                                                            <  1:3524582/40‑1 (MQ=255)
                                            aTGCTCACTTCTTTACTTAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                           <  1:2710426/71‑1 (MQ=255)
                                            aTGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGa                            <  1:1685507/70‑1 (MQ=255)
                                               cTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt                         <  1:2569078/70‑1 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCaa                                                          <  1:1687722/36‑1 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGc                               >  1:2013815/1‑63 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGc                               >  1:1877977/1‑63 (MQ=255)
                                                tCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATc                        >  1:2709380/1‑70 (MQ=255)
                                                 cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                           >  1:2469899/1‑66 (MQ=255)
                                                 cACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGaa                           >  1:1332984/1‑66 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAAt                         <  1:1206059/67‑1 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                      >  1:945253/1‑70 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                      >  1:975520/1‑70 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                      >  1:1840944/1‑70 (MQ=255)
                                                  aCTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGc                      >  1:280352/1‑70 (MQ=255)
                                                   cttcttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCa                    <  1:3150950/71‑1 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaaa                  >  1:3131431/1‑70 (MQ=255)
                                                      cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                 >  1:493160/1‑71 (MQ=255)
                                                       tttACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                         <  1:2742556/46‑1 (MQ=255)
                                                       tttACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                         <  1:2584679/46‑1 (MQ=255)
                                                       tttACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCaa                   <  1:2380354/68‑1 (MQ=255)
                                                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCacac                                         >  1:1450416/1‑43 (MQ=255)
                                                          aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt              <  1:357423/70‑1 (MQ=255)
                                                           cTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGc            <  1:2708205/71‑1 (MQ=255)
                                                            tGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTcac                                           >  1:851099/1‑39 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCg                       >  1:3872703/1‑57 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCg                       >  1:576091/1‑57 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCg                       >  1:1129232/1‑57 (MQ=255)
                                                              aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg           <  1:2276500/69‑1 (MQ=255)
                                                                 tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCt       <  1:219463/70‑1 (MQ=255)
                                                                   cGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtgt     <  1:1111121/70‑1 (MQ=255)
                                                                    gTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtgg   <  1:2952137/71‑1 (MQ=255)
                                                                    gTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtgg   <  1:3665311/71‑1 (MQ=255)
                                                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGGCTACGCTGAAATCGCCAAACTGTGCg           <  1:1124910/62‑1 (MQ=255)
                                                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc  <  1:2074091/71‑1 (MQ=255)
                                                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc  <  1:1655177/71‑1 (MQ=255)
                                                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc  <  1:3627765/71‑1 (MQ=255)
                                                                     ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc  <  1:884587/71‑1 (MQ=255)
                                                                     |                                                                      
CATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGC  >  minE/832289‑832428

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: