Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,852 G→A 10.7% P190P (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,8520GA10.7% 190.8 / 11.2 67P190P (CCG→CCAnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base A (4/3);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

GATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAC  >  minE/835783‑835921
                                                                      |                                                                     
gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGACCACTGCCTGAACCCg                                                                       >  1:3526373/1‑70 (MQ=255)
gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg                                                                       >  1:2089863/1‑70 (MQ=255)
 aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                                                      >  1:2079504/1‑70 (MQ=255)
    tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg                                                                  <  1:1048690/71‑1 (MQ=255)
     ataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg                                                                  <  1:3191144/70‑1 (MQ=255)
      taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                                                                >  1:1669244/1‑71 (MQ=255)
        gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt                                                               >  1:2683894/1‑70 (MQ=255)
          cAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                                                            >  1:761217/1‑71 (MQ=255)
             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                                                                <  1:1640584/64‑1 (MQ=255)
             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                                                         <  1:904302/71‑1 (MQ=255)
             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                                                         <  1:3718816/71‑1 (MQ=255)
               tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCt                                                       >  1:1518126/2‑71 (MQ=255)
               cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                                                       >  1:2410564/1‑71 (MQ=255)
                gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                                                       <  1:751524/70‑1 (MQ=255)
                   aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGCGACCTGcc                                                    >  1:2527627/1‑69 (MQ=255)
                    acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg                                                  <  1:3555968/71‑1 (MQ=255)
                       ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGc                                               <  1:3346729/71‑1 (MQ=255)
                        cTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGc                                               <  1:117428/70‑1 (MQ=255)
                           cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc                                                         >  1:47455/1‑57 (MQ=255)
                           cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCgg                                             <  1:3516571/69‑1 (MQ=255)
                           cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCgg                                             <  1:1394333/69‑1 (MQ=255)
                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                           >  1:3527305/1‑69 (MQ=255)
                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                           >  1:3600629/1‑69 (MQ=255)
                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                           >  1:117734/1‑69 (MQ=255)
                            atgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGAGGTGc                                          >  1:3149782/1‑71 (MQ=255)
                                  tATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGccc                                                   >  1:1951244/1‑56 (MQ=255)
                                  tATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAttt                                     <  1:3044142/70‑1 (MQ=255)
                                   aTTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                                                                >  1:3106394/1‑42 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                                                                >  1:1395579/1‑41 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAtt                                      <  1:1673720/67‑1 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                   >  1:3463761/1‑70 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                   >  1:1111159/1‑70 (MQ=255)
                                     ctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTActa                                 <  1:1585897/70‑3 (MQ=255)
                                      tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGa                                                          >  1:1564231/1‑45 (MQ=255)
                                      tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                >  1:1829469/1‑71 (MQ=255)
                                      tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                >  1:4045314/1‑71 (MQ=255)
                                          gcgcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                       <  1:624428/60‑1 (MQ=255)
                                          gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                            <  1:373081/71‑3 (MQ=255)
                                          gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                            <  1:2385307/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:2590296/69‑1 (MQ=255)
                                                 tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                       <  1:669598/69‑1 (MQ=255)
                                                 tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                      <  1:2782525/70‑1 (MQ=255)
                                                   cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                   >  1:3143601/1‑71 (MQ=255)
                                                   cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                   >  1:1010523/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:3533152/1‑51 (MQ=255)
                                                                      ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                     >  1:1090911/1‑51 (MQ=255)
                                                                      |                                                                     
GATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAC  >  minE/835783‑835921

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: