Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,016 G→C 22.5% G1158G (GGG→GGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,0160GC22.5% 238.7 / 64.5 120G1158G (GGG→GGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (46/47);  new base C (16/11);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 125 total aligned reads displayed.

GATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCG  >  minE/834947‑835084
                                                                     |                                                                    
gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGa                                                                     >  1:5364934/1‑71 (MQ=255)
gATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGa                                                                     >  1:3396998/1‑71 (MQ=255)
    gAAGTATTCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGAc                                                                 >  1:2012149/1‑71 (MQ=255)
     aaGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGAc                                                                 <  1:388643/70‑1 (MQ=255)
      aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGa                                                                  >  1:1540180/1‑68 (MQ=255)
      aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                               <  1:114199/71‑1 (MQ=255)
      aGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                               <  1:5112215/71‑1 (MQ=255)
       gtattCAACGCCAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCa                                                               >  1:4559720/4‑70 (MQ=255)
       gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                                                              >  1:1337076/1‑71 (MQ=255)
       gTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCat                                                              >  1:277424/1‑71 (MQ=255)
          ttCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCatgat                                                           >  1:1140239/1‑71 (MQ=255)
            cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGAt                                                                    >  1:1074713/1‑60 (MQ=255)
            cAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGAt                                                                    >  1:783057/1‑60 (MQ=255)
               cAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTAcc                                                      <  1:3832943/71‑1 (MQ=255)
                 ccAACGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAc                                                    >  1:4609626/2‑71 (MQ=255)
                   aaCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgc                                                  >  1:1958381/1‑71 (MQ=255)
                   aaCGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgc                                                  >  1:906405/1‑71 (MQ=255)
                   aaCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgc                                                  >  1:828027/1‑71 (MQ=255)
                    aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTat                                                       >  1:2708417/1‑64 (MQ=255)
                    aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTat                                                       >  1:1890195/1‑64 (MQ=255)
                    aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTat                                                       >  1:2962100/1‑64 (MQ=255)
                    aCGGCGCGCTGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCAcgcc                                                 >  1:3788152/1‑70 (MQ=255)
                      ggTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCAcgcg                                                 >  1:4350576/1‑69 (MQ=255)
                       gTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAg                                              <  1:1291937/71‑1 (MQ=255)
                        tGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAgg                                             <  1:859508/71‑1 (MQ=255)
                          cgctGACTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGaa                                           <  1:412202/69‑1 (MQ=255)
                            ctGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                         >  1:447976/1‑71 (MQ=255)
                             tGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTTCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                         >  1:4915988/1‑70 (MQ=255)
                              gACTGCCCGTGCGGTTGTCAGGCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTa                                       >  1:4388736/1‑71 (MQ=255)
                              gACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTa                                       >  1:2540680/1‑71 (MQ=255)
                               aCTGCCCGCGCGGTGGTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAt                                      >  1:91440/1‑71 (MQ=255)
                                   cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCg                                    <  1:661469/69‑1 (MQ=255)
                                   cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                  <  1:4116101/71‑1 (MQ=255)
                                   cccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                   >  1:2246713/1‑70 (MQ=255)
                                   cccGCGCGGTGGTCAACCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAttat                                   >  1:45183/1‑67 (MQ=255)
                                    ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                        <  1:1621651/64‑1 (MQ=255)
                                    ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                        <  1:3707454/64‑1 (MQ=255)
                                    ccGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGt                                        <  1:3891855/64‑1 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                               <  1:2566865/38‑1 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCa                                                               <  1:2461041/38‑1 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTa                                                        >  1:1371684/1‑45 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTa                                                        >  1:293482/1‑45 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTa                                                        >  1:641543/1‑45 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                  <  1:2860641/67‑1 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                  <  1:3040612/67‑1 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                  <  1:679706/67‑1 (MQ=255)
                                       tGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcc                              >  1:2168996/1‑71 (MQ=255)
                                        gCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcct                             >  1:1216563/1‑71 (MQ=255)
                                         cGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctg                            >  1:3045363/1‑71 (MQ=255)
                                            ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                         <  1:1044911/71‑1 (MQ=255)
                                            ttGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAAcctgcct                         <  1:2602551/71‑1 (MQ=255)
                                              gTCAGCCAGCGTGTTCCGGCAGGGATGACCatgatg                                                          >  1:1405320/1‑36 (MQ=255)
                                              gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACg                                         >  1:1033878/1‑53 (MQ=255)
                                              gTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTgg                       >  1:4585703/1‑71 (MQ=255)
                                              gTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCa                                               <  1:3277177/47‑1 (MQ=16)
                                              gTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCa                                               <  1:2518405/47‑1 (MQ=16)
                                              gTCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCa                                               <  1:4638057/47‑1 (MQ=16)
                                               tCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGtt                                  <  1:745313/59‑1 (MQ=255)
                                               tCAGCCAACGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCAt                                      <  1:3224699/55‑1 (MQ=255)
                                                cagccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGtt                     >  1:4710332/1‑71 (MQ=255)
                                                cagccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAgg                                             >  1:2456129/1‑47 (MQ=16)
                                                 agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                    >  1:2082186/1‑71 (MQ=255)
                                                 agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                    >  1:168575/1‑71 (MQ=255)
                                                 agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                    >  1:839341/1‑71 (MQ=255)
                                                 agccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTc                    >  1:832669/1‑71 (MQ=255)
                                                  gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                 <  1:4065311/57‑1 (MQ=255)
                                                  gccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                 <  1:2997136/57‑1 (MQ=255)
                                                   ccagcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                  <  1:4492449/71‑1 (MQ=255)
                                                   ccaacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCgg                  <  1:3500068/71‑1 (MQ=255)
                                                     agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                 >  1:1111952/1‑70 (MQ=255)
                                                     agcGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                 <  1:5320507/70‑1 (MQ=255)
                                                     aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAgg                                             <  1:442422/40‑1 (MQ=12)
                                                     aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCg                   <  1:1771948/66‑1 (MQ=255)
                                                     aacGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaa                <  1:4544682/69‑1 (MQ=255)
                                                      acGTGTACCGCCGGGCATGACCATGATGTATCACGCCCAGGaa                                           >  1:4853511/2‑43 (MQ=12)
                                                       cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                  <  1:2805172/67‑1 (MQ=255)
                                                       cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                  <  1:3043868/67‑1 (MQ=255)
                                                       cGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCgg                  <  1:5510869/67‑1 (MQ=255)
                                                        gtgtACCGCCGGGCATGACCATGATGTATCACGCCCAGGAACGCATTATGAATATTCCTGGTTCGGaagt              >  1:540723/1‑68 (MQ=255)
                                                         tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                   <  1:4901642/48‑1 (MQ=255)
                                                         tgtTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGt                                   <  1:5269924/48‑1 (MQ=255)
                                                           ttCCGGCAGGGATGACCATGATGTACCACGCGCAg                                              >  1:5210523/1‑35 (MQ=255)
                                                           ttCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGa                 <  1:3109862/64‑1 (MQ=255)
                                                             ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa               >  1:5110952/1‑64 (MQ=255)
                                                             ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGaaa               >  1:630580/1‑64 (MQ=255)
                                                             ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc         >  1:5149989/1‑70 (MQ=255)
                                                             ccGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTAccc         >  1:1956105/1‑70 (MQ=255)
                                                              cGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa       <  1:100192/71‑1 (MQ=255)
                                                              cGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCa        <  1:2256105/70‑1 (MQ=255)
                                                                 cAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAg    <  1:5439251/71‑1 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCaa       >  1:3655850/1‑67 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc   >  1:5220430/1‑71 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc   >  1:2105595/1‑71 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc   <  1:711814/71‑1 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc   <  1:2919289/71‑1 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc   <  1:4290125/71‑1 (MQ=255)
                                                                  aGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGc   <  1:3324322/71‑1 (MQ=255)
                                                                   gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTa                                 <  1:4034070/40‑1 (MQ=255)
                                                                   gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg  <  1:1873344/71‑1 (MQ=255)
                                                                   gggATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCg  >  1:32498/1‑71 (MQ=255)
                                                                     |                                                                    
GATTGAAGTCTTCAACAGCAACGGTGCTCTGACTGCCCGTGCGGTTGTCAGCCAGCGTGTTCCGGCAGGGATGACCATGATGTACCACGCGCAGGAACGTATCGTTAACCTGCCTGGTTCGGAAATTACCCAACAGCG  >  minE/834947‑835084

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