Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,846 G→C 12.6% L188L (CTG→CTC)  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,8460GC12.6% 165.8 / 11.9 63L188L (CTG→CTC) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base C (5/3);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCT  >  minE/835780‑835914
                                                                  |                                                                     
ggCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑cc                                                                  <  1:3816693/71‑1 (MQ=255)
   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:3600629/1‑69 (MQ=255)
                              cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg                                    >  1:117734/1‑69 (MQ=255)
                              cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggcg                                    >  1:3527305/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:4045314/1‑71 (MQ=255)
                                         tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                         >  1:1829469/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:669598/69‑1 (MQ=255)
                                                    tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                <  1:2590296/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:2695292/71‑1 (MQ=255)
                                                                cTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt  <  1:1308139/71‑1 (MQ=255)
                                                                  |                                                                     
GGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCT  >  minE/835780‑835914

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: ‑