Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,999 T→C 22.8% G239G (GGT→GGC 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,9990TC22.8% 129.3 / 37.1 70G239G (GGT→GGCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (15/39);  new base C (5/11);  total (20/50)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

TGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGA  >  minE/835931‑836064
                                                                    |                                                                 
tGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCaa                                                                 <  1:838168/71‑1 (MQ=255)
  ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt                                                               <  1:604585/71‑1 (MQ=255)
  ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt                                                               <  1:2965818/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:3246154/69‑1 (MQ=255)
           gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa                                                        <  1:1843920/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)
                    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: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: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:39260/1‑71 (MQ=255)
                                                     cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg            >  1:2590286/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:364190/45‑1 (MQ=255)
                                                              gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa                             <  1:359474/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)
                                                                    |                                                                 
TGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGA  >  minE/835931‑836064

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: