Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,909 C→T 14.4% I209I (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,9090CT14.4% 191.1 / 20.0 76I209I (ATC→ATTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base T (6/5);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATC  >  minE/835841‑835978
                                                                    |                                                                     
tGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGt                                                                     >  1:499022/1‑71 (MQ=255)
  ccTGAACCCGGCATGTGTGGCGAGCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTc                                                                    <  1:3077938/70‑1 (MQ=255)
   cTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt                                                                  <  1:3603997/71‑1 (MQ=255)
      aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGa                                                                <  1:1373010/70‑1 (MQ=255)
        cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg                                                             <  1:982906/71‑1 (MQ=255)
          cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc                                                           >  1:3930504/1‑71 (MQ=255)
            gCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTc                                                                    >  1:3153570/1‑60 (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:1143689/71‑1 (MQ=255)
                   tgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                  <  1:3425337/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)
                            gCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGa                                                                >  1:3035260/1‑48 (MQ=255)
                            gCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGAt                                                               >  1:951557/1‑49 (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)
                               cGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc                                                                       <  1:2830272/38‑1 (MQ=255)
                                  gcggcgCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGt                                   <  1:423219/71‑1 (MQ=255)
                                  gcggcgCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGt                                   <  1:2485156/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:2288057/1‑68 (MQ=255)
                                        cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg                                >  1:2190091/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:2449836/71‑1 (MQ=255)
                                                  gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                   <  1:2291124/71‑1 (MQ=255)
                                                  gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                   <  1:364089/71‑1 (MQ=255)
                                                  gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                   <  1:2115226/71‑1 (MQ=255)
                                                     tGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCg                                    <  1:1785532/51‑1 (MQ=255)
                                                      gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:2415818/64‑1 (MQ=255)
                                                      gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg                      <  1:1989575/64‑1 (MQ=255)
                                                      gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc               >  1:1734274/1‑71 (MQ=255)
                                                      gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc               >  1:3641783/1‑71 (MQ=255)
                                                       aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg              <  1:986271/71‑1 (MQ=255)
                                                       aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg              <  1:3582159/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:635731/71‑1 (MQ=255)
                                                                gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa     <  1:3279261/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)
                                                                    |                                                                     
TGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATC  >  minE/835841‑835978

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