Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,912 C→G 5.1% V210V (GTC→GTG 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,9120CG 5.1% 246.9 / 2.6 78V210V (GTC→GTGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (31/43);  new base G (2/2);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAC  >  minE/835843‑835981
                                                                     |                                                                     
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)
    aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                     >  1:757825/1‑68 (MQ=255)
     aCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                 >  1:3624603/1‑71 (MQ=255)
      cccGGCATGTGTGGCGACCTGCGCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                 >  1:4018316/1‑70 (MQ=255)
      cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                 >  1:3070801/1‑70 (MQ=255)
        cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                              <  1:4018509/71‑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:1035948/58‑1 (MQ=255)
              gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                     <  1:3808084/58‑1 (MQ=255)
               tgtgGCGACCTGCCCGAGCGGTTCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                       <  1:3843982/71‑1 (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)
                     gACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                                                <  1:2757191/56‑1 (MQ=255)
                       ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                                               >  1:1070637/1‑55 (MQ=255)
                       ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                              >  1:2224999/1‑56 (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)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                              <  1:2819306/48‑1 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                              <  1:1439225/48‑1 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                     >  1:697619/1‑57 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                     >  1:2295813/1‑57 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctg                                          >  1:727784/1‑68 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                       >  1:1832239/1‑71 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                       >  1:1429452/1‑71 (MQ=255)
                               aGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                       <  1:3693529/71‑1 (MQ=255)
                                 cGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                     <  1:1031246/71‑1 (MQ=255)
                                 cGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGtt                                     <  1:3105730/71‑2 (MQ=255)
                                    tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                       <  1:2986992/50‑1 (MQ=255)
                                    tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                       <  1:3485686/66‑1 (MQ=255)
                                     gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                 >  1:3744958/1‑71 (MQ=255)
                                      cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                                                  >  1:578630/1‑37 (MQ=255)
                                      cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                                                  >  1:3754100/1‑37 (MQ=255)
                                      cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                 <  1:3394863/70‑1 (MQ=255)
                                         tttACAAGCGTGAAGAAGATGGCATCGTCCTGATGGACCAGGATAAATGCCGt                                               <  1:4032851/53‑1 (MQ=255)
                                         tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                              <  1:2983292/54‑1 (MQ=255)
                                         tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                             <  1:1977388/71‑1 (MQ=255)
                                           tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGa                                                               >  1:4123587/1‑35 (MQ=255)
                                           tACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                           >  1:1857657/1‑71 (MQ=255)
                                           tACAAGAGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa                            >  1:539242/1‑70 (MQ=255)
                                           tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                                            >  1:1667466/1‑54 (MQ=255)
                                           tACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggc                                            >  1:2352836/1‑54 (MQ=255)
                                            aCAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt                          >  1:2751332/1‑71 (MQ=255)
                                             cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                         <  1:2122629/71‑1 (MQ=255)
                                             cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                         <  1:1445361/71‑1 (MQ=255)
                                             cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                         <  1:2085022/71‑1 (MQ=255)
                                             cAAGCGCGAGGAGGACGGCATCGTCCTGATCGACCAGGACAAGTGCCGCGGCTGGCGGATGTGCa                               <  1:1434330/65‑1 (MQ=255)
                                              aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                       <  1:1203018/40‑1 (MQ=255)
                                              aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                        <  1:770933/71‑1 (MQ=255)
                                              aaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                        <  1:3258100/71‑1 (MQ=255)
                                                gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                     >  1:3891022/1‑40 (MQ=255)
                                                gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                     >  1:3903465/1‑40 (MQ=255)
                                                gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                       <  1:3126887/54‑1 (MQ=255)
                                                 cGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGtggc                                            <  1:1464355/48‑1 (MQ=255)
                                                 cGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                       <  1:3506649/53‑1 (MQ=255)
                                                   tGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGt                                      <  1:917889/52‑1 (MQ=255)
                                                     aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGc                                <  1:1409670/56‑1 (MQ=255)
                                                     aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc                  >  1:836903/1‑70 (MQ=255)
                                                     aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                 >  1:2578843/1‑71 (MQ=255)
                                                     aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                 >  1:1113775/1‑71 (MQ=255)
                                                       gaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                                         >  1:328541/1‑45 (MQ=255)
                                                          gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                        >  1:166444/1‑59 (MQ=255)
                                                          gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACa             >  1:3400660/1‑70 (MQ=255)
                                                           aTGGGATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                         <  1:1931453/57‑1 (MQ=255)
                                                             ggCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                         <  1:2478135/55‑1 (MQ=255)
                                                              gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                           >  1:1487657/1‑52 (MQ=255)
                                                               cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaaa       >  1:3841752/1‑71 (MQ=255)
                                                                 tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc     <  1:1790183/71‑1 (MQ=255)
                                                                   gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTa   <  1:1253443/71‑1 (MQ=255)
                                                                   gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTa   <  1:608344/71‑1 (MQ=255)
                                                                    tCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATc                             <  1:342487/44‑1 (MQ=255)
                                                                    tCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAc  <  1:949667/71‑1 (MQ=255)
                                                                     ccTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                     >  1:909908/1‑35 (MQ=255)
                                                                     |                                                                     
CCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAC  >  minE/835843‑835981

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: