Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,190 C→T 14.1% S216S (AGC→AGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1900CT14.1% 250.9 / 25.6 99S216S (AGC→AGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (45/40);  new base T (5/9);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTA  >  minE/832121‑832256
                                                                     |                                                                  
ggCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCt                                                                   <  1:5232913/71‑1 (MQ=255)
 gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                                                  >  1:1552504/1‑71 (MQ=255)
 gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                                                  >  1:1360984/1‑71 (MQ=255)
 gCAATGTCGATGGTTTCTTACGCATCGGGTGAACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCt                                                                   >  1:3361782/1‑70 (MQ=255)
  cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGTCGGTACTTGCTTAAGCTTc                                                                 <  1:5181192/71‑1 (MQ=255)
  cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                                                 <  1:1379888/71‑1 (MQ=255)
  cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                                                 <  1:2256318/71‑1 (MQ=255)
     tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                                                              <  1:4317805/71‑1 (MQ=255)
     tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                                                              >  1:1451138/1‑71 (MQ=255)
       tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                                                            >  1:4296815/1‑71 (MQ=255)
       tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                                                            >  1:1965010/1‑71 (MQ=255)
         gATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                                                          >  1:2149607/1‑71 (MQ=255)
         gATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACt                                                          >  1:2233995/1‑71 (MQ=255)
         gATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACt                                                          >  1:2666366/1‑71 (MQ=255)
           tGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTgg                                                        <  1:2098656/71‑1 (MQ=255)
            ggTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTgg                                                        >  1:4249428/1‑70 (MQ=255)
             gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGt                                                       <  1:5306882/70‑1 (MQ=255)
             gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                      >  1:1920521/1‑71 (MQ=255)
             gTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                      >  1:1816686/1‑71 (MQ=255)
              tttCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTAc                                                     >  1:586300/1‑71 (MQ=255)
               ttCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                                                          <  1:966277/65‑1 (MQ=255)
               ttCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                                                          <  1:1504544/65‑1 (MQ=255)
                tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACt                                                          >  1:5542322/1‑64 (MQ=255)
                tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                    <  1:428408/70‑1 (MQ=255)
                tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                    <  1:4581840/70‑1 (MQ=255)
                tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                                                   >  1:3482011/1‑71 (MQ=255)
                tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                                                   >  1:3910588/1‑71 (MQ=255)
                  ttACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                                                 >  1:291284/1‑55 (MQ=255)
                    aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc                                                  <  1:4645435/68‑1 (MQ=255)
                     cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                                                              <  1:5182964/55‑1 (MQ=255)
                     cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                              <  1:4281585/71‑1 (MQ=255)
                      gcttCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc                                                  <  1:798671/63‑1 (MQ=255)
                      gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                             <  1:3784066/71‑1 (MQ=255)
                        aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc                                               >  1:2123407/1‑67 (MQ=255)
                        aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc                                               >  1:3373802/1‑67 (MQ=255)
                        aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc                                           >  1:3964838/1‑71 (MQ=255)
                         tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                          <  1:4910086/71‑1 (MQ=255)
                         tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                          <  1:121244/71‑1 (MQ=255)
                         tCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                          <  1:2835482/71‑1 (MQ=255)
                          cGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                    <  1:2174563/60‑1 (MQ=255)
                           gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                        >  1:4599508/1‑71 (MQ=255)
                           gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcct                                         <  1:4498985/70‑1 (MQ=255)
                             gTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                      <  1:747071/71‑1 (MQ=255)
                               gCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcc                                       <  1:2590620/68‑1 (MQ=255)
                                  cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                 <  1:168697/71‑1 (MQ=255)
                                   gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                                                 <  1:2787953/38‑1 (MQ=255)
                                    ttATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                                                          <  1:2325871/43‑1 (MQ=12)
                                    ttATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACt                                                          <  1:3373230/43‑1 (MQ=12)
                                    cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                               <  1:1356289/71‑1 (MQ=255)
                                     tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTa                                                               >  1:5105221/1‑38 (MQ=255)
                                     tATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                             >  1:1171821/1‑56 (MQ=255)
                                      aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGAc                                               >  1:3205980/1‑53 (MQ=19)
                                      aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGAc                                               >  1:3175116/1‑53 (MQ=19)
                                      aTCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                                                 <  1:1875241/35‑1 (MQ=255)
                                       tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCg                             <  1:3323946/70‑1 (MQ=255)
                                       tctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                              >  1:5193665/1‑69 (MQ=255)
                                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                             <  1:3387117/53‑1 (MQ=255)
                                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                             <  1:5300695/53‑1 (MQ=255)
                                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                             <  1:1298291/53‑1 (MQ=255)
                                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                            >  1:5165428/1‑70 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                                                            >  1:1631215/1‑37 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                      >  1:1441997/1‑43 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                      >  1:802918/1‑43 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                      >  1:2222665/1‑43 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                      >  1:2146808/1‑43 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                      >  1:687189/1‑59 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                      >  1:1252817/1‑59 (MQ=255)
                                         tctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                      >  1:725684/1‑59 (MQ=255)
                                          ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                              >  1:1235824/1‑66 (MQ=255)
                                          ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                              >  1:1526610/1‑66 (MQ=255)
                                          ctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                         >  1:5052580/1‑71 (MQ=255)
                                           tcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCgatg                          <  1:761100/69‑4 (MQ=255)
                                           tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                        >  1:3511473/1‑71 (MQ=255)
                                             gctgctTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTg                                            >  1:3583570/1‑49 (MQ=16)
                                             gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                       <  1:569964/70‑1 (MQ=255)
                                             gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                       <  1:2597733/70‑1 (MQ=255)
                                                gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAcc                       <  1:3422109/67‑1 (MQ=255)
                                                  ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                   >  1:5620596/1‑69 (MQ=255)
                                                  ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                 <  1:1200017/71‑1 (MQ=255)
                                                    ggCGGTACATGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                <  1:1009608/70‑1 (MQ=255)
                                                        gTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCaa           <  1:3585810/71‑1 (MQ=255)
                                                          aCTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa            <  1:663507/68‑1 (MQ=255)
                                                           cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAAc         <  1:4644879/70‑1 (MQ=255)
                                                           cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt        <  1:163062/71‑1 (MQ=255)
                                                            ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                      >  1:1414845/1‑40 (MQ=255)
                                                            ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctcct                                      >  1:4451249/1‑40 (MQ=255)
                                                            ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                          <  1:5615593/52‑1 (MQ=255)
                                                             tGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                        >  1:631224/1‑37 (MQ=255)
                                                             tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                              >  1:1707352/1‑47 (MQ=255)
                                                             tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                              >  1:4719233/1‑47 (MQ=255)
                                                             tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                              >  1:1723349/1‑47 (MQ=255)
                                                               cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg    <  1:1722190/71‑1 (MQ=255)
                                                                 tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt   <  1:3375919/70‑1 (MQ=255)
                                                                 tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGt   <  1:3382378/70‑1 (MQ=255)
                                                                 tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa               >  1:2239762/1‑58 (MQ=255)
                                                                 tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGa               >  1:5256285/1‑58 (MQ=255)
                                                                  aaGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa  <  1:1006427/70‑1 (MQ=255)
                                                                  aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                            >  1:3515463/1‑44 (MQ=255)
                                                                  aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGt   <  1:4166571/69‑1 (MQ=255)
                                                                     |                                                                  
GGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTA  >  minE/832121‑832256

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: