Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,008 G→A 17.2% E242E (GAG→GAA narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,0080GA17.2% 149.5 / 22.0 64E242E (GAG→GAAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (19/34);  new base A (2/9);  total (21/43)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGT  >  minE/835938‑836077
                                                                      |                                                                     
gctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag                                                                       <  1:203394/71‑1 (MQ=255)
    gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaaat                                                                   >  1:1719179/1‑69 (MQ=255)
    gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                                     <  1:3246154/69‑1 (MQ=255)
    gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                                     <  1:1843920/69‑1 (MQ=255)
    gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                                     <  1:3241472/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:309516/1‑71 (MQ=255)
        aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt                                                               >  1:1890133/1‑71 (MQ=255)
         tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtc                                                              <  1:1561817/71‑1 (MQ=255)
         tgtgCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCAtc                                                              >  1:778851/1‑71 (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)
                      aTGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa                                                                     <  1:2546401/51‑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:2545390/67‑1 (MQ=255)
                                       tCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                                    <  1:3554602/67‑1 (MQ=255)
                                       tCTACTTCAACTGGAAGAGAGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                        >  1:3360921/1‑63 (MQ=255)
                                          aCTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt                                                            <  1:1433275/40‑1 (MQ=255)
                                           cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                              <  1:3729428/69‑1 (MQ=255)
                                           cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGcc                              <  1:2234600/69‑1 (MQ=255)
                                           cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCg                             <  1:1050029/70‑1 (MQ=255)
                                           cTTCAACTCCAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAACCGGGTCAGcc                              <  1:224040/69‑1 (MQ=255)
                                              cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                         >  1:39260/1‑71 (MQ=255)
                                              cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                         >  1:2590286/1‑71 (MQ=255)
                                              cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                         >  1:1174327/1‑71 (MQ=255)
                                               aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                        <  1:1229407/71‑1 (MQ=255)
                                               aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt                        <  1:514625/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:359474/45‑1 (MQ=255)
                                                       gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                                          <  1:364190/45‑1 (MQ=255)
                                                       gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg                <  1:431094/71‑1 (MQ=255)
                                                        agCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAg                >  1:2151562/1‑70 (MQ=255)
                                                        agCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGa               <  1:78065/71‑1 (MQ=255)
                                                         gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                     <  1:1834760/48‑1 (MQ=255)
                                                         gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGa               <  1:2602887/70‑1 (MQ=255)
                                                          cGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgtgt                      <  1:2625560/62‑1 (MQ=255)
                                                             caaGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTc                  <  1:1988407/62‑1 (MQ=255)
                                                              aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                     >  1:1747970/1‑43 (MQ=255)
                                                              aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                     >  1:2661141/1‑43 (MQ=255)
                                                              aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCt          <  1:730330/70‑1 (MQ=255)
                                                              aaGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtg         <  1:1467286/71‑1 (MQ=255)
                                                               aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgt        >  1:1183589/1‑71 (MQ=255)
                                                               aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc        <  1:1451602/71‑2 (MQ=255)
                                                               aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc        <  1:57487/71‑2 (MQ=255)
                                                               aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc        <  1:1001021/71‑2 (MQ=255)
                                                                gTCTGAGAAGTTCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgtg       <  1:3809722/71‑1 (MQ=255)
                                                                 tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgtg       <  1:3616713/70‑1 (MQ=255)
                                                                  cTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtc     >  1:1855850/1‑71 (MQ=255)
                                                                  cTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtc     >  1:2418464/1‑71 (MQ=255)
                                                                   tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcg    >  1:713415/1‑71 (MQ=255)
                                                                     agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtctgt  <  1:3299085/71‑1 (MQ=255)
                                                                      |                                                                     
GCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGT  >  minE/835938‑836077

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: