Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,918 C→T 14.8% I212I (ATC→ATT 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,9180CT14.8% 183.2 / 20.7 74I212I (ATC→ATTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base T (6/5);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTT  >  minE/835849‑835983
                                                                     |                                                                 
cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                                                  <  1:982906/71‑1 (MQ=255)
  cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                                >  1:3930504/1‑71 (MQ=255)
     cATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAg                                                             >  1:1687624/1‑71 (MQ=255)
         tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAgaga                                                          >  1:1414510/1‑67 (MQ=255)
         tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                          <  1:3621278/70‑1 (MQ=255)
         cgtgGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATa                                                         <  1:2255555/70‑1 (MQ=255)
          gtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATc                                                                   >  1:804926/1‑60 (MQ=255)
           tgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                       <  1:3425337/71‑1 (MQ=255)
           tgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                       <  1:1143689/71‑1 (MQ=255)
            ggCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAAt                                                      >  1:1969010/1‑71 (MQ=255)
            ggCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAAt                                                      >  1:820397/1‑71 (MQ=255)
             gCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATg                                                     >  1:2882119/1‑71 (MQ=255)
              cGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGc                                                    <  1:3429905/71‑1 (MQ=255)
                aCCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCg                                                  >  1:728049/1‑71 (MQ=255)
                aCCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCg                                                  >  1:1984304/1‑71 (MQ=255)
                 ccTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGt                                                 <  1:529350/71‑1 (MQ=255)
                     cccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGcc                                                   >  1:2063342/1‑65 (MQ=255)
                      ccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                       >  1:354336/1‑60 (MQ=255)
                       caaGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggctgg                                           <  1:3365225/69‑1 (MQ=255)
                       caaGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggctg                                            >  1:780153/3‑70 (MQ=255)
                       caaGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggctg                                            >  1:1106065/3‑70 (MQ=255)
                          gcggcgCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGt                                        <  1:2485156/71‑1 (MQ=255)
                          gcggcgCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGt                                        <  1:423219/71‑1 (MQ=255)
                           cGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGt                                        <  1:3323047/70‑1 (MQ=255)
                              tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                       <  1:1792878/68‑1 (MQ=255)
                              tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                       <  1:2532702/68‑1 (MQ=255)
                               gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                   >  1:3399541/1‑71 (MQ=255)
                                cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                   <  1:2331074/70‑1 (MQ=255)
                                cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg                                   >  1:3524678/1‑70 (MQ=255)
                                cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg                                     >  1:2190091/1‑68 (MQ=255)
                                cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg                                     >  1:2288057/1‑68 (MQ=255)
                                cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAt                                      >  1:943696/1‑67 (MQ=255)
                                  aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt                                                          <  1:1692177/45‑1 (MQ=255)
                                  aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt                                >  1:163782/1‑71 (MQ=255)
                                   tttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATa                                                         >  1:3559595/1‑45 (MQ=255)
                                   tttAAAAGCGTGAAGAAGATGGCATCGTCCTGATCGCCCAGGATa                                                         >  1:1291911/1‑45 (MQ=255)
                                    ttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaa                                                        >  1:3580415/1‑45 (MQ=255)
                                    ttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa                              >  1:2083247/1‑71 (MQ=255)
                                      actaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc                             <  1:1885904/67‑1 (MQ=255)
                                      aCAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGt                                                 <  1:2037880/50‑1 (MQ=255)
                                       cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt                            <  1:293862/70‑1 (MQ=255)
                                          gCGTGAAGAAGATGGCATCGTCCTGATCGACCAgg                                                            <  1:545146/35‑1 (MQ=255)
                                          gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                        <  1:364089/71‑1 (MQ=255)
                                          gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                        <  1:2115226/71‑1 (MQ=255)
                                          gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                        <  1:2449836/71‑1 (MQ=255)
                                          gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                        <  1:2291124/71‑1 (MQ=255)
                                             tGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                         <  1:1785532/51‑1 (MQ=255)
                                              gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                           <  1:1989575/64‑1 (MQ=255)
                                              gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                           <  1:2415818/64‑1 (MQ=255)
                                              gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc                    >  1:3641783/1‑71 (MQ=255)
                                              gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc                    >  1:1734274/1‑71 (MQ=255)
                                               aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                   <  1:3582159/71‑1 (MQ=255)
                                               aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                   <  1:986271/71‑1 (MQ=255)
                                                agaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg                   >  1:3266908/1‑70 (MQ=255)
                                                  aagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGtggctgg                                           <  1:1841925/44‑1 (MQ=255)
                                                  aagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg                          <  1:1780459/61‑1 (MQ=255)
                                                   agaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGa                         >  1:933712/1‑61 (MQ=255)
                                                    gacggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACa               >  1:2196522/4‑70 (MQ=255)
                                                    gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaa              <  1:781931/71‑1 (MQ=255)
                                                      tGGCATCGTCCTGATCGACCAGGATAAATGCCGtg                                                >  1:1358608/1‑35 (MQ=255)
                                                      tGGCATCGTCCTGAACGACCAGGATAAATGCCGTGGCTGGCGTAtgt                                    >  1:2608690/1‑47 (MQ=255)
                                                      cggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACaaa             >  1:155852/2‑70 (MQ=255)
                                                       ggCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaa            >  1:1867789/1‑70 (MQ=255)
                                                        gCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACaaaaaa          <  1:71648/71‑1 (MQ=255)
                                                        gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa          <  1:3279261/71‑1 (MQ=255)
                                                        gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa          <  1:635731/71‑1 (MQ=255)
                                                         cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAt                                      >  1:3091753/1‑42 (MQ=255)
                                                         cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaaa         <  1:177164/71‑1 (MQ=255)
                                                           tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc       >  1:3691146/1‑71 (MQ=255)
                                                             gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCt      <  1:2782806/70‑1 (MQ=255)
                                                              tCCTGATCGACCAGGATAAATGCCGTGGCTGGCTTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAc    <  1:925842/71‑1 (MQ=255)
                                                                cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt   <  1:219070/70‑1 (MQ=255)
                                                                cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt   <  1:100622/70‑1 (MQ=255)
                                                                 tGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACtt  >  1:2159777/1‑70 (MQ=255)
                                                                  gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc                    >  1:256887/1‑51 (MQ=255)
                                                                     |                                                                 
CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTT  >  minE/835849‑835983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: