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→CCA) narHnitrate 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

GATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAC  >  minE/835783‑835921
                                                                      |                                                                     
gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGACCACTGCCTGAA‑CCCg                                                                       >  1:3526373/1‑70 (MQ=255)
gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCg                                                                       >  1:2089863/1‑70 (MQ=255)
 aTGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCgg                                                                      >  1:2079504/1‑70 (MQ=255)
    tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtg                                                                  <  1:1048690/71‑1 (MQ=255)
     ataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtg                                                                  <  1:3191144/70‑1 (MQ=255)
      taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                                                                >  1:1669244/1‑71 (MQ=255)
        gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtgt                                                               >  1:2683894/1‑70 (MQ=255)
          cAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc                                                            >  1:761217/1‑71 (MQ=255)
             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                                                                <  1:1640584/64‑1 (MQ=255)
             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                                                         <  1:904302/71‑1 (MQ=255)
             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                                                         <  1:3718816/71‑1 (MQ=255)
               tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt                                                       >  1:1518126/2‑71 (MQ=255)
               cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                                                       >  1:2410564/1‑71 (MQ=255)
                gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt                                                       <  1:751524/70‑1 (MQ=255)
                   aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGT‑GCGACCTGcc                                                    >  1:2527627/1‑69 (MQ=255)
                    acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCg                                                  <  1:3555968/71‑1 (MQ=255)
                       ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGc                                               <  1:3346729/71‑1 (MQ=255)
                        cTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGc                                               <  1:117428/70‑1 (MQ=255)
                           cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGAc                                                         >  1:47455/1‑57 (MQ=255)
                           cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCgg                                             <  1:3516571/69‑1 (MQ=255)
                           cATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCgg                                             <  1:1394333/69‑1 (MQ=255)
                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg                                           >  1:3527305/1‑69 (MQ=255)
                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg                                           >  1:3600629/1‑69 (MQ=255)
                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg                                           >  1:117734/1‑69 (MQ=255)
                            atgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGAGGTGc                                          >  1:3149782/1‑71 (MQ=255)
                                  tATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGccc                                                   >  1:1951244/1‑56 (MQ=255)
                                  tATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAttt                                     <  1:3044142/70‑1 (MQ=255)
                                   aTTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                                                                >  1:3106394/1‑42 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg                                                                >  1:1395579/1‑41 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAtt                                      <  1:1673720/67‑1 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                   >  1:3463761/1‑70 (MQ=255)
                                    tttGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                   >  1:1111159/1‑70 (MQ=255)
                                     ctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTActa                                 <  1:1585897/70‑3 (MQ=255)
                                      tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGa                                                          >  1:1564231/1‑45 (MQ=255)
                                      tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                >  1:1829469/1‑71 (MQ=255)
                                      tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                                >  1:4045314/1‑71 (MQ=255)
                                          gcgcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                       <  1:624428/60‑1 (MQ=255)
                                          gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                            <  1:373081/71‑3 (MQ=255)
                                          gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACgcg                            <  1:2385307/71‑3 (MQ=255)
                                           cgcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                            >  1:1196786/1‑70 (MQ=255)
                                             ccTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                         <  1:2951842/71‑1 (MQ=255)
                                               tgtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                         >  1:3680443/1‑69 (MQ=255)
                                               tgtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                         >  1:1926426/1‑69 (MQ=255)
                                                 tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                       <  1:2590296/69‑1 (MQ=255)
                                                 tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                       <  1:669598/69‑1 (MQ=255)
                                                 tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                      <  1:2782525/70‑1 (MQ=255)
                                                   cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                   >  1:3143601/1‑71 (MQ=255)
                                                   cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATg                   >  1:1010523/1‑71 (MQ=255)
                                                     aaCACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                  >  1:1651483/1‑70 (MQ=255)
                                                     aaCACTGCCTCAA‑TCCCAGATGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACgg                  >  1:2350301/1‑70 (MQ=255)
                                                      acacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                       >  1:2775066/1‑64 (MQ=255)
                                                      acacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                          >  1:3369353/1‑61 (MQ=255)
                                                      acacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCa                <  1:71421/71‑1 (MQ=255)
                                                       cacTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGa                                        <  1:3839322/46‑1 (MQ=255)
                                                        acTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                         <  1:3717996/60‑1 (MQ=255)
                                                          tGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                       <  1:294023/60‑1 (MQ=255)
                                                          tGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt            >  1:3699780/1‑71 (MQ=255)
                                                          tGCCTGAA‑‑CCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                       <  1:39978/59‑1 (MQ=255)
                                                            ccTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc          >  1:4021566/1‑71 (MQ=255)
                                                            ccTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTcc          >  1:911336/1‑71 (MQ=255)
                                                             cTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt         <  1:1308139/71‑1 (MQ=255)
                                                             cTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt         <  1:2695292/71‑1 (MQ=255)
                                                                aa‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt         >  1:757825/1‑68 (MQ=255)
                                                                 a‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                     >  1:1386404/1‑55 (MQ=255)
                                                                 a‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc     >  1:3624603/1‑71 (MQ=255)
                                                                 a‑CCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                     >  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)
                                                                      |                                                                     
GATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAC  >  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: ‑