Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,993 G→A 25.0% K237K (AAG→AAA narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,9930GA25.0% 114.5 / 36.0 64K237K (AAG→AAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (14/34);  new base A (5/11);  total (19/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAG  >  minE/835925‑836063
                                                                    |                                                                      
gATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAaaag                                                                      >  1:3892619/1‑71 (MQ=255)
gATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAagag                                                                      >  1:3515872/1‑71 (MQ=255)
  tAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCg                                                                    <  1:271886/71‑1 (MQ=255)
  tAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCg                                                                    <  1:3585669/71‑1 (MQ=255)
      tGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCaa                                                                <  1:838168/71‑1 (MQ=255)
        ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt                                                              <  1:2965818/71‑1 (MQ=255)
        ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt                                                              <  1:604585/71‑1 (MQ=255)
          gTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCt                                                            <  1:239832/71‑1 (MQ=255)
           tggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg                                                           <  1:2128510/71‑1 (MQ=255)
            ggctggcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTcag                                                           >  1:1195984/1‑68 (MQ=255)
            ggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCt                                                            >  1:2222021/1‑69 (MQ=255)
             gctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                         <  1:203394/71‑1 (MQ=255)
               tggcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAgaaa                                                        >  1:1152407/1‑68 (MQ=255)
                 gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaaat                                                     >  1:1719179/1‑69 (MQ=255)
                 gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                       <  1:3241472/69‑1 (MQ=255)
                 gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                       <  1:1843920/69‑1 (MQ=255)
                 gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                       <  1:3246154/69‑1 (MQ=255)
                 gcGTATGTTCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAg                                                      >  1:2915481/1‑70 (MQ=255)
                 gcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGt                                                     >  1:3747510/1‑71 (MQ=255)
                  cGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                         <  1:934189/66‑1 (MQ=255)
                    tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                  <  1:1503231/71‑1 (MQ=255)
                    tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa                                                  <  1:1710470/71‑1 (MQ=255)
                     aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                 >  1:1890133/1‑71 (MQ=255)
                     aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                 >  1:309516/1‑71 (MQ=255)
                      tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtc                                                <  1:1561817/71‑1 (MQ=255)
                      tgtgCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCAtc                                                >  1:778851/1‑71 (MQ=255)
                          cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGt                                                              >  1:207690/1‑53 (MQ=255)
                          cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttc                                             <  1:1723363/70‑1 (MQ=255)
                            tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc                                          >  1:241643/1‑71 (MQ=255)
                              aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                        <  1:501911/71‑1 (MQ=255)
                              aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa                                        <  1:2942732/71‑1 (MQ=255)
                               cTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTAt                                       >  1:3086932/1‑71 (MQ=255)
                                 ggATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATcc                                     >  1:2096601/1‑71 (MQ=255)
                                  gATGCCCGTACAAAAAAATCTACTTCAATTGGAAGAGCGGTAAGTc                                                             <  1:2849429/46‑1 (MQ=255)
                                  gATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTc                                                             <  1:583075/46‑1 (MQ=255)
                                   aTGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                       <  1:2546401/51‑1 (MQ=255)
                                    tGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga                                                          <  1:136854/47‑1 (MQ=255)
                                      tccGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCcacgaa                                <  1:2509527/70‑3 (MQ=255)
                                        cGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTAtt                              <  1:1497057/71‑1 (MQ=255)
                                         gTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg                             <  1:960377/71‑1 (MQ=255)
                                          tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga                                                        <  1:1231986/43‑1 (MQ=255)
                                            cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                 >  1:2328939/1‑48 (MQ=255)
                                            cAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGag                           <  1:1155359/70‑2 (MQ=255)
                                              aaaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttc                                             <  1:3307498/50‑1 (MQ=255)
                                              aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaaat                                                     <  1:2405111/42‑3 (MQ=255)
                                               aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTa                                <  1:3062759/62‑1 (MQ=255)
                                               aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg                             <  1:1252982/65‑1 (MQ=255)
                                               aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCCTCTGCTATCCGCGTAtt                              <  1:299188/64‑1 (MQ=255)
                                                aaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTATGCTATCCGCGTATTGAAGCggg                      <  1:3819574/71‑1 (MQ=255)
                                                   aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCa                   <  1:3498115/71‑1 (MQ=255)
                                                    tCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                      <  1:3554602/67‑1 (MQ=255)
                                                    tCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                      <  1:2545390/67‑1 (MQ=255)
                                                    tCTACTTCAACTGGAAGAGAGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                          >  1:3360921/1‑63 (MQ=255)
                                                       aCTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt                                              <  1:1433275/40‑1 (MQ=255)
                                                        cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                <  1:2234600/69‑1 (MQ=255)
                                                        cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                <  1:3729428/69‑1 (MQ=255)
                                                        cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCg               <  1:1050029/70‑1 (MQ=255)
                                                        cTTCAACTCCAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAACCGGGTCAGcc                <  1:224040/69‑1 (MQ=255)
                                                           cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg           >  1:1174327/1‑71 (MQ=255)
                                                           cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg           >  1:2590286/1‑71 (MQ=255)
                                                           cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg           >  1:39260/1‑71 (MQ=255)
                                                            aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt          <  1:514625/71‑1 (MQ=255)
                                                            aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt          <  1:1229407/71‑1 (MQ=255)
                                                             aCTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTac                                       <  1:611707/41‑2 (MQ=12)
                                                               tGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGtgtg       >  1:2997512/1‑71 (MQ=255)
                                                                 gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt     <  1:534380/71‑1 (MQ=255)
                                                                    gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                            <  1:364190/45‑1 (MQ=255)
                                                                    gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                            <  1:359474/45‑1 (MQ=255)
                                                                    gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg  <  1:431094/71‑1 (MQ=255)
                                                                    |                                                                      
GATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAG  >  minE/835925‑836063

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: