Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,313 G→A 37.2% V257V (GTG→GTA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3130GA37.2% 76.8 / 84.2 86V257V (GTG→GTAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (20/34);  new base A (20/12);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 2.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCC  >  minE/832244‑832382
                                                                     |                                                                     
gctaaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                                                      <  1:1355083/68‑1 (MQ=255)
gCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                                                      <  1:2662689/71‑1 (MQ=255)
 cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                                                     <  1:618696/71‑1 (MQ=255)
 cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                                                     <  1:1191578/71‑1 (MQ=255)
 cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                                                     >  1:1545580/1‑71 (MQ=255)
 cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                                                      >  1:1217330/1‑70 (MQ=255)
   aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGc                                                                   <  1:3461661/71‑1 (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:2342308/1‑71 (MQ=255)
            aCCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTa                                                          >  1:2193418/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)
                            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)
                                 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:2186408/1‑43 (MQ=255)
                                   aCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAca                                                              >  1:1200998/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:84732/71‑1 (MQ=255)
                                                    gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                  <  1:531809/71‑1 (MQ=255)
                                                    gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                  <  1:3844702/71‑1 (MQ=255)
                                                    gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                  >  1:2569466/1‑71 (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)
                                                                   gTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTa                                     >  1:2319136/1‑37 (MQ=255)
                                                                    tGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGCACCAAAACTGTTGcc  <  1:1592419/71‑1 (MQ=255)
                                                                     |                                                                     
GCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCC  >  minE/832244‑832382

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: