Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,039 G→T 19.5% A253S (GCG→TCG) ‡ 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,0390GT19.5% 150.6 / 17.6 61A253S (GCG→TCG) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (20/29);  new base T (5/7);  total (25/36)
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.92e-01

CAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGA  >  minE/835969‑836109
                                                                      |                                                                      
cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                                        >  1:975136/1‑71 (MQ=255)
    aaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg                                                                    >  1:4372422/1‑71 (MQ=255)
     aaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGt                                                                   <  1:1879659/71‑1 (MQ=255)
            cTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAAcc                                                              <  1:3325358/69‑1 (MQ=255)
               cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                                                                >  1:688396/1‑64 (MQ=255)
               cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg                                                                >  1:1748580/1‑64 (MQ=255)
               cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGAGGGTCAg                                                                >  1:5459720/1‑64 (MQ=255)
                aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                                                         <  1:1274678/70‑1 (MQ=255)
                    ggAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGc                                                    >  1:1181434/1‑71 (MQ=255)
                     gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                                        >  1:1831097/1‑50 (MQ=255)
                     gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                                        >  1:3881858/1‑50 (MQ=255)
                     gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg                                                                        >  1:4733433/1‑50 (MQ=255)
                     gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                                    <  1:2373275/70‑1 (MQ=255)
                          gCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGaa                                              <  1:3124529/71‑1 (MQ=255)
                          gCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGaa                                              <  1:2603101/71‑1 (MQ=255)
                             gTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg                                                         <  1:3512480/57‑1 (MQ=255)
                              tAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGTGTCAGCCGACCGTGTGCTCAGAAACCt                                          <  1:3052700/71‑1 (MQ=255)
                                aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCg                                                                      >  1:807190/1‑41 (MQ=255)
                                aGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtgt                                        >  1:1209816/1‑71 (MQ=255)
                                aGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCtgc                                        >  1:1711302/1‑70 (MQ=255)
                                  tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                                    <  1:4415156/57‑1 (MQ=255)
                                  tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                                    <  1:3714398/57‑1 (MQ=255)
                                  tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc                                                    <  1:3871235/57‑1 (MQ=255)
                                    tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                                                     >  1:3008409/1‑38 (MQ=255)
                                    tGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCgg                                                                     >  1:4145982/1‑38 (MQ=255)
                                      agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggt                                  <  1:4389993/71‑1 (MQ=255)
                                      agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcgg                                   <  1:508103/70‑1 (MQ=255)
                                        aaatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtc                                 >  1:777092/4‑70 (MQ=255)
                                         aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGaa                                              <  1:1557606/54‑1 (MQ=255)
                                         aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGaa                                              <  1:4917612/54‑1 (MQ=255)
                                         aatGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgc                               <  1:1656730/69‑2 (MQ=255)
                                         aGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGGCGGTCGt                               <  1:523969/71‑1 (MQ=255)
                                          atGCATCTTCTGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgca                              >  1:3327909/2‑69 (MQ=255)
                                           tGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTAt                             <  1:401398/71‑1 (MQ=255)
                                            gCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAcc                                           <  1:1254123/56‑1 (MQ=255)
                                              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATcct                          <  1:5139350/71‑2 (MQ=255)
                                              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCg                          <  1:1523005/71‑1 (MQ=255)
                                              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCg                          <  1:4443711/71‑1 (MQ=255)
                                              aTCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAACTGTGTCGGTCGTATCCg                          <  1:566046/71‑1 (MQ=255)
                                               tcttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGtgtg                                                     <  1:2499695/43‑1 (MQ=255)
                                                cttctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGa                                               <  1:4473373/48‑1 (MQ=255)
                                                 ttctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa                       >  1:5082527/1‑71 (MQ=255)
                                                  tctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTAt                      >  1:3191624/1‑71 (MQ=255)
                                                  tctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTAt                      >  1:58117/1‑71 (MQ=255)
                                                   ctGTTACCCACGAATTGAGTCCGGTCAACCGACCGTGTGCTCAGAAACCTGCgtgggtcgca                              <  1:2475869/62‑3 (MQ=255)
                                                   ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt                        <  1:675981/68‑1 (MQ=255)
                                                   ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGtt                        <  1:5263314/68‑1 (MQ=255)
                                                   ctGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTa                       <  1:5161340/69‑1 (MQ=255)
                                                      cTATCCGCGTATTTAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCt                    <  1:4448824/69‑1 (MQ=255)
                                                         tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt                   <  1:2319512/67‑1 (MQ=255)
                                                         tCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCtt                   <  1:2269130/67‑1 (MQ=255)
                                                           cgcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGt              <  1:1879847/70‑1 (MQ=255)
                                                            gcgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGc            >  1:2062906/1‑71 (MQ=255)
                                                             cgTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt           >  1:208666/1‑71 (MQ=255)
                                                                aTTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCt           >  1:2992793/1‑68 (MQ=255)
                                                                aTTGAAGCGGGTCAGACGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgtt        >  1:2292953/1‑71 (MQ=255)
                                                                 ttGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg       >  1:5677282/1‑71 (MQ=255)
                                                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg       >  1:4447281/1‑70 (MQ=255)
                                                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg       <  1:2321273/70‑1 (MQ=255)
                                                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg       >  1:3607389/1‑70 (MQ=255)
                                                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg       >  1:3838078/1‑70 (MQ=255)
                                                                  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg       >  1:866778/1‑70 (MQ=255)
                                                                   gAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGAGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa     >  1:269440/1‑71 (MQ=255)
                                                                    aaGCGTGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa     <  1:3936445/70‑1 (MQ=255)
                                                                     aGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACg   <  1:4309061/71‑1 (MQ=255)
                                                                      gCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa  >  1:4667558/1‑71 (MQ=255)
                                                                      |                                                                      
CAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGA  >  minE/835969‑836109

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: