Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,328 C→G 23.9% T262T (ACC→ACG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,3280CG23.9% 118.1 / 35.9 67T262T (ACC→ACGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base G (9/7);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACT  >  minE/832258‑832395
                                                                      |                                                                   
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)
                     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:2476045/71‑1 (MQ=255)
                                atcatcGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTa                                     <  1:3473666/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:3844702/71‑1 (MQ=255)
                                      gCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAgg                               <  1:531809/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: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: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)
                                                                      |                                                                   
CGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACT  >  minE/832258‑832395

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: