Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,322 G→A 32.2% T260T (ACG→ACA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3220GA32.2% 76.6 / 49.0 62T260T (ACG→ACAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (18/24);  new base A (10/10);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCA  >  minE/832254‑832391
                                                                    |                                                                     
gTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacg                                                                     <  1:2666827/71‑1 (MQ=255)
 tACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGt                                                                    <  1:225416/71‑1 (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)
   ccGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacgaa                                                                   >  1:1522485/1‑68 (MQ=255)
    cGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAcc                                                                 <  1:2679312/71‑1 (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)
           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)
                       taaCTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGCCCACTTCttt                                               <  1:2280594/69‑1 (MQ=255)
                       cAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCttt                                               <  1:599472/70‑1 (MQ=255)
                        aaCTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAc                                                                  >  1:1188893/1‑50 (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)
                                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:3844702/71‑1 (MQ=255)
                                          gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                           <  1:531809/71‑1 (MQ=255)
                                          gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                           >  1:2569466/1‑71 (MQ=255)
                                          gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                           <  1:84732/71‑1 (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:3450809/1‑71 (MQ=255)
                                                     aaaCGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTg                >  1:838937/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)
                                                                    |                                                                     
GTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCA  >  minE/832254‑832391

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: