Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,307 A→T 38.8% S255S (TCA→TCT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3070AT38.8% 67.3 / 98.5 85S255S (TCA→TCTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (19/33);  new base T (21/12);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTG  >  minE/832238‑832377
                                                                     |                                                                      
gggCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                                                        <  1:3053453/70‑1 (MQ=255)
gggCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                                                        <  1:3753141/70‑1 (MQ=255)
   cGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACg                                                                    <  1:259794/71‑1 (MQ=255)
    gAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACg                                                                    >  1:98572/1‑70 (MQ=255)
      gctaaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                                                 <  1:1355083/68‑1 (MQ=255)
      gCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                                                 <  1:2662689/71‑1 (MQ=255)
       cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                                                <  1:1191578/71‑1 (MQ=255)
       cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                                                <  1:618696/71‑1 (MQ=255)
       cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                                                >  1:1545580/1‑71 (MQ=255)
       cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                                                 >  1:1217330/1‑70 (MQ=255)
        aaaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCaa                                                                       <  1:2915844/63‑1 (MQ=255)
         aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGc                                                              <  1:3461661/71‑1 (MQ=255)
         aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCACATCATCGCCTGGGGGTCAAAc                                                                     >  1:2568647/1‑63 (MQ=255)
         aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCACATCATCGCCTGGGGGTCAAAc                                                                     >  1:2409408/1‑63 (MQ=255)
          accgACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCg                                                               >  1:2536065/4‑69 (MQ=255)
          aCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCg                                                               >  1:2502871/1‑69 (MQ=255)
              aCGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa                                                           <  1:3652509/69‑1 (MQ=255)
              aCGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAc                                                          <  1:3415441/70‑1 (MQ=255)
                gTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacg                                                       <  1:2666827/71‑1 (MQ=255)
                 tACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGt                                                      <  1:225416/71‑1 (MQ=255)
                 tACCGGAATCTGCCGACTGGTATAACTCCAGCTACTTCATCGCCTGGGGGTCTAACGTACCGCAGAcac                                                        <  1:228291/69‑3 (MQ=255)
                 tACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacg                                                       <  1:1049604/70‑1 (MQ=255)
                  aCCGGAATCTGCCGACTGGTATAACTCCAGCTAGATCATCGCCTGGGGGTCTAACGTACCGCAGACACGt                                                      <  1:1323966/70‑1 (MQ=255)
                  aCCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTa                                                     >  1:2193418/1‑71 (MQ=255)
                  aCCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTa                                                     >  1:2342308/1‑71 (MQ=255)
                   ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa                                                           <  1:3919279/64‑1 (MQ=255)
                   ccGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacgaa                                                     >  1:1522485/1‑68 (MQ=255)
                    cGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAcc                                                   <  1:2679312/71‑1 (MQ=255)
                    cGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAca                                                         >  1:1701875/1‑64 (MQ=255)
                     ggAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAccc                                                  >  1:3743391/1‑71 (MQ=255)
                      gAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTAcgcc                                                 <  1:390411/71‑1 (MQ=255)
                       aaTCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGTTCAAACGTGCCGCAGACGCGTACCCCg                                                <  1:485721/71‑1 (MQ=255)
                       aaTCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCg                                                <  1:3692813/71‑1 (MQ=255)
                       aaTCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGa                                                           >  1:2412859/1‑60 (MQ=255)
                        aactgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCa                                                             <  1:211184/55‑1 (MQ=255)
                           tgccgACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGacg                                            >  1:170482/1‑69 (MQ=255)
                            gccgACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGc                                           >  1:1510936/4‑71 (MQ=255)
                              tgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg                                            <  1:887414/68‑1 (MQ=255)
                                 cTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCAc                                       <  1:971018/70‑1 (MQ=255)
                                  tGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                                                        <  1:2318042/36‑1 (MQ=255)
                                  tGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                                     >  1:2942217/1‑71 (MQ=255)
                                   ggTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                                                        <  1:2889815/35‑1 (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)
                                       taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAg                                                            <  1:1068020/42‑1 (MQ=255)
                                       taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCttt                                 <  1:2280594/69‑1 (MQ=255)
                                       cAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCttt                                 <  1:599472/70‑1 (MQ=255)
                                        aaCTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAc                                                    >  1:1188893/1‑50 (MQ=255)
                                         aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAca                                                         >  1:1200998/1‑43 (MQ=255)
                                         aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAca                                                         >  1:2186408/1‑43 (MQ=255)
                                         aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAc                                                          >  1:637283/1‑43 (MQ=255)
                                         aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAc                                                          >  1:641960/1‑43 (MQ=255)
                                         aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAc                                                          >  1:3615971/1‑43 (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)
                                            ccAGCTACATCATCGCCTGGGGGTCTAACGTACCg                                                               >  1:1608210/1‑35 (MQ=38)
                                            ccAGCTACATCATCGCCTGGGGGTCAAACGTGCCGc                                                              >  1:1116364/1‑36 (MQ=255)
                                              aGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGt                         >  1:459677/1‑71 (MQ=255)
                                                cTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCAct                                      <  1:463447/56‑1 (MQ=255)
                                                cTACATCATCGCCTGGGGGTAAAACGTGCCGCAGACGCGTACCCCGGa                                              >  1:3713743/1‑48 (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:2569466/1‑71 (MQ=255)
                                                          gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg             <  1:84732/71‑1 (MQ=255)
                                                          gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg             <  1:531809/71‑1 (MQ=255)
                                                          gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg             <  1:3844702/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)
                                                                     aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGa                            <  1:1595823/45‑1 (MQ=255)
                                                                     aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg  >  1:838937/1‑71 (MQ=255)
                                                                     aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg  >  1:3450809/1‑71 (MQ=255)
                                                                     |                                                                      
GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTG  >  minE/832238‑832377

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: