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→CTCnarHnitrate 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

GGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCT  >  minE/835780‑835914
                                                                  |                                                                     
ggCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc                                                                  <  1:3816693/71‑1 (MQ=255)
   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:3600629/1‑69 (MQ=255)
                              cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                    >  1:117734/1‑69 (MQ=255)
                              cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggcg                                    >  1:3527305/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:4045314/1‑71 (MQ=255)
                                         tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAg                         >  1:1829469/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:669598/69‑1 (MQ=255)
                                                    tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                <  1:2590296/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:2695292/71‑1 (MQ=255)
                                                                cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt  <  1:1308139/71‑1 (MQ=255)
                                                                  |                                                                     
GGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCT  >  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: