Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,873 G→A 10.4% P197P (CCG→CCA 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,8730GA10.4% 228.0 / 13.7 77P197P (CCG→CCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base A (4/4);  total (39/38)
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.73e-01

ACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGG  >  minE/835803‑835942
                                                                       |                                                                     
acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg                                                                       <  1:3555968/71‑1 (MQ=255)
   ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGc                                                                    <  1:3346729/71‑1 (MQ=255)
    cTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGc                                                                    <  1:117428/70‑1 (MQ=255)
       cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCgg                                                                  <  1:1394333/69‑1 (MQ=255)
       cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCgg                                                                  <  1:3516571/69‑1 (MQ=255)
       cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                                                >  1:117734/1‑69 (MQ=255)
       cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                                                >  1:3600629/1‑69 (MQ=255)
       cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                                                >  1:3527305/1‑69 (MQ=255)
        atgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGAGGTGc                                                               >  1:3149782/1‑71 (MQ=255)
              tATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAttt                                                          <  1:3044142/70‑1 (MQ=255)
                tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAtt                                                           <  1:1673720/67‑1 (MQ=255)
                tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                                        >  1:1111159/1‑70 (MQ=255)
                tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                                        >  1:3463761/1‑70 (MQ=255)
                 ctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTActa                                                      <  1:1585897/70‑3 (MQ=255)
                  tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                                     >  1:4045314/1‑71 (MQ=255)
                  tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                                     >  1:1829469/1‑71 (MQ=255)
                      gcgcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                                            <  1:624428/60‑1 (MQ=255)
                      gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                                                 <  1:2385307/71‑3 (MQ=255)
                      gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                                                 <  1:373081/71‑3 (MQ=255)
                       cgcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                                                 >  1:1196786/1‑70 (MQ=255)
                         ccTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                              <  1:2951842/71‑1 (MQ=255)
                           tgtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                              >  1:3680443/1‑69 (MQ=255)
                           tgtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                              >  1:1926426/1‑69 (MQ=255)
                             tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:669598/69‑1 (MQ=255)
                             tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:2590296/69‑1 (MQ=255)
                             tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                                           <  1:2782525/70‑1 (MQ=255)
                               cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                        >  1:1010523/1‑71 (MQ=255)
                               cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                                        >  1:3143601/1‑71 (MQ=255)
                                 aaCACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                                       >  1:1651483/1‑70 (MQ=255)
                                 aaCACTGCCTCAATCCCAGATGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACgg                                       >  1:2350301/1‑70 (MQ=255)
                                  acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            >  1:2775066/1‑64 (MQ=255)
                                  acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                                               >  1:3369353/1‑61 (MQ=255)
                                  acacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCa                                     <  1:71421/71‑1 (MQ=255)
                                   cacTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGa                                                             <  1:3839322/46‑1 (MQ=255)
                                    acTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                                              <  1:3717996/60‑1 (MQ=255)
                                      tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:294023/60‑1 (MQ=255)
                                      tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt                                 >  1:3699780/1‑71 (MQ=255)
                                      tGCCTGAACCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:39978/59‑1 (MQ=255)
                                        ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                               >  1:4021566/1‑71 (MQ=255)
                                        ccTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc                               >  1:911336/1‑71 (MQ=255)
                                         cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              <  1:1308139/71‑1 (MQ=255)
                                         cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              <  1:2695292/71‑1 (MQ=255)
                                            aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              >  1:757825/1‑68 (MQ=255)
                                             aCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:1386404/1‑55 (MQ=255)
                                             aCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                          >  1:3624603/1‑71 (MQ=255)
                                             aCCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:3363334/1‑55 (MQ=255)
                                               cccGGCATGTGTGGCGACCTGCGCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                          >  1:4018316/1‑70 (MQ=255)
                                               cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                          >  1:3070801/1‑70 (MQ=255)
                                                ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:3014250/51‑1 (MQ=255)
                                                 cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                       <  1:4018509/71‑1 (MQ=255)
                                                  ggCATGTGTGGCGACCTGCCCGAGCGGTGCGAttt                                                          >  1:2682848/1‑35 (MQ=255)
                                                  ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:1090911/1‑51 (MQ=255)
                                                  ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:3533152/1‑51 (MQ=255)
                                                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:3252888/45‑1 (MQ=255)
                                                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:1350903/45‑1 (MQ=255)
                                                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGc                                      <  1:1662162/51‑1 (MQ=255)
                                                      tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                             <  1:2435676/60‑1 (MQ=255)
                                                      tgtgtgGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGa                  <  1:1152876/71‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              <  1:35998/58‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              <  1:3808084/58‑1 (MQ=255)
                                                       gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                              <  1:1035948/58‑1 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTTCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                <  1:3843982/71‑1 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            >  1:1343375/1‑43 (MQ=255)
                                                        tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                 <  1:3039922/70‑1 (MQ=255)
                                                        cgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGa                  <  1:201995/68‑1 (MQ=255)
                                                           ggCGCCCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcgt                              >  1:2612787/1‑52 (MQ=255)
                                                            gCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAAt             <  1:1605920/70‑1 (MQ=255)
                                                              gACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:2066572/1‑39 (MQ=255)
                                                              gACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                         <  1:2757191/56‑1 (MQ=255)
                                                               aCCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                                            <  1:2988911/36‑1 (MQ=255)
                                                                ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                        >  1:1070637/1‑55 (MQ=255)
                                                                ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                       >  1:2224999/1‑56 (MQ=255)
                                                                 cTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:110676/1‑36 (MQ=255)
                                                                 cTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                                          >  1:2541091/1‑36 (MQ=255)
                                                                 cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                        >  1:1070854/1‑54 (MQ=255)
                                                                 cTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg       >  1:2985414/1‑71 (MQ=255)
                                                                     ccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg   <  1:1017280/71‑1 (MQ=255)
                                                                      cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg  >  1:853566/1‑71 (MQ=255)
                                                                       gAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTg                             <  1:325348/43‑1 (MQ=255)
                                                                       |                                                                     
ACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGG  >  minE/835803‑835942

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: