Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,921 C→T 15.0% D213D (GAC→GAT 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,9210CT15.0% 180.9 / 20.3 73D213D (GAC→GATnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base T (6/5);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTG  >  minE/835851‑835989
                                                                      |                                                                    
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)
         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)
                   cccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGcc                                                         >  1:2063342/1‑65 (MQ=255)
                    ccGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATaaa                                                             >  1:354336/1‑60 (MQ=255)
                     caaGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggctgg                                                 <  1:3365225/69‑1 (MQ=255)
                     caaGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggctg                                                  >  1:1106065/3‑70 (MQ=255)
                     caaGCGGCGCTATCTACAAACGCGAAGAAGACGGCATTGTGCTGATTGATCAGGATAAATGCCGtggctg                                                  >  1:780153/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:2532702/68‑1 (MQ=255)
                            tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa                                             <  1:1792878/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:2449836/71‑1 (MQ=255)
                                        gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                              <  1:364089/71‑1 (MQ=255)
                                        gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt                              <  1:2291124/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: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)
                                                                    gACCAGGATAAATGCCGTGGGTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTg  <  1:2925096/71‑1 (MQ=255)
                                                                      |                                                                    
CGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTG  >  minE/835851‑835989

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: