Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,337 T→C 20.3% A265A (GCT→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3370TC20.3% 138.6 / 30.4 69A265A (GCT→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (21/34);  new base C (7/7);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATC  >  minE/832271‑832406
                                                                  |                                                                     
cTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCAc                                                                    <  1:971018/70‑1 (MQ=255)
 tGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                                                                  >  1:2942217/1‑71 (MQ=255)
   gTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGCCACGTACGCCGGACGCCCActtct                                                                >  1:2614210/1‑71 (MQ=255)
   gTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtct                                                                >  1:2291785/1‑71 (MQ=255)
   gTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACACGTACCCCGGATGCTCActtct                                                                <  1:1890624/71‑1 (MQ=255)
   gTACAACTCCAGCTAAATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtc                                                                 <  1:3567884/70‑1 (MQ=255)
    tACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                               <  1:1598700/71‑1 (MQ=255)
     aCAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCAGGATGATCActtctt                                                               >  1:3898303/1‑70 (MQ=255)
      taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCttt                                                              <  1:2280594/69‑1 (MQ=255)
      cAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCttt                                                              <  1:599472/70‑1 (MQ=255)
        aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTAcc                                                           >  1:2859635/1‑70 (MQ=255)
        aCTCCAGCTACATAATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTAc                                                            >  1:1291364/1‑70 (MQ=255)
         gtCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCActtc                                                                 <  1:1262830/63‑1 (MQ=255)
          tCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTg                                                          >  1:900394/1‑70 (MQ=255)
             aGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGt                                                      >  1:459677/1‑71 (MQ=255)
               cTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCAct                                                                   <  1:463447/56‑1 (MQ=255)
                 aCATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                               <  1:2290592/58‑1 (MQ=255)
                  catcatCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGt                                                  <  1:866089/70‑1 (MQ=255)
                   atcatcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCAc                                                                    >  1:954430/1‑51 (MQ=255)
                   atcatcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTa                                                <  1:1878824/71‑1 (MQ=255)
                   atcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt                                                               <  1:714420/56‑1 (MQ=255)
                   atcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTa                                                <  1:3473666/71‑1 (MQ=255)
                   atcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTa                                                <  1:2476045/71‑1 (MQ=255)
                     catcGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCActt                                                                  <  1:277374/51‑1 (MQ=255)
                      atcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACa                                              <  1:3042733/70‑1 (MQ=255)
                        cGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAg                                           >  1:1611567/1‑71 (MQ=255)
                         gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                          <  1:531809/71‑1 (MQ=255)
                         gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                          <  1:3844702/71‑1 (MQ=255)
                         gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                          <  1:84732/71‑1 (MQ=255)
                         gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                                          >  1:2569466/1‑71 (MQ=255)
                          ccTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAggg                                         <  1:2463365/71‑1 (MQ=255)
                            tGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCAc                                       >  1:3325524/1‑71 (MQ=255)
                            tGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCTTTACCGAAGTACGCTACAAAGGCAc                                       >  1:3366837/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)
                                    aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGa                                                         <  1:1595823/45‑1 (MQ=255)
                                    aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                               >  1:838937/1‑71 (MQ=255)
                                    aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                               >  1:3450809/1‑71 (MQ=255)
                                        gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTa                                                             >  1:2319136/1‑37 (MQ=255)
                                         tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGCACCAAAACTGTTGcc                          <  1:1592419/71‑1 (MQ=255)
                                           ccGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCg                         >  1:2477889/1‑70 (MQ=255)
                                           ccGCAGACACGTACGCCGGACGCCCACTTCTTTAc                                                            <  1:1460777/35‑1 (MQ=19)
                                              cAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTca                      >  1:1664343/1‑70 (MQ=255)
                                                  cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                 <  1:3465792/71‑1 (MQ=255)
                                                  cgcgTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCa                 <  1:108177/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:1628907/1‑45 (MQ=255)
                                                     gTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGa                                        >  1:591485/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)
                                                                  |                                                                     
CTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATC  >  minE/832271‑832406

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