Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,342 C→T 12.6% H20H (CAC→CAT 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,3420CT12.6% 197.3 / 17.2 72H20H (CAC→CATnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (24/38);  new base T (2/7);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACG  >  minE/835273‑835412
                                                                     |                                                                      
agagCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCcaca                                                                       <  1:2470985/71‑1 (MQ=255)
    cGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg                                                                   <  1:456784/71‑1 (MQ=255)
      tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGtt                                                                 >  1:1315777/1‑71 (MQ=255)
      aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt                                                                 >  1:1565500/2‑71 (MQ=255)
        aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                                               <  1:3013980/71‑1 (MQ=255)
        aaaTGAAAATTCGTTCACAAATCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:2036912/70‑1 (MQ=255)
         aaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCTGCTGCCACACCTGTTCa                                                               <  1:1130553/70‑1 (MQ=255)
           tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                            <  1:2563524/71‑1 (MQ=255)
           tGAAAATCCGTACACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGt                                                             >  1:760565/1‑70 (MQ=255)
             aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGTCTGCCACACCTGTTCAGTTAc                                                          >  1:1714717/1‑71 (MQ=255)
             aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                                                          >  1:525186/1‑71 (MQ=255)
                atccGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc                                                                <  1:1220197/59‑1 (MQ=255)
                    gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                   <  1:584646/71‑1 (MQ=255)
                    gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                   <  1:3949171/71‑1 (MQ=255)
                    gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                   <  1:2609513/71‑1 (MQ=255)
                    gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                   <  1:1209399/71‑1 (MQ=255)
                     ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                   <  1:2759346/70‑1 (MQ=255)
                      tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt                                                             >  1:959615/1‑59 (MQ=255)
                       cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:2091827/55‑1 (MQ=255)
                       cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:1258608/55‑1 (MQ=255)
                       cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                                >  1:218643/1‑71 (MQ=255)
                        acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                                                <  1:1864425/70‑1 (MQ=255)
                        acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                                                <  1:1398341/70‑1 (MQ=255)
                         caAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt                                              <  1:2737055/71‑1 (MQ=255)
                          aaGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCt                                                        >  1:3289879/1‑60 (MQ=255)
                            gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                           >  1:2017075/1‑71 (MQ=255)
                            gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                           <  1:3414327/71‑1 (MQ=255)
                            gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                           <  1:142995/71‑1 (MQ=255)
                             tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:3542954/49‑1 (MQ=255)
                             tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:555095/49‑1 (MQ=255)
                             tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:2673847/49‑1 (MQ=255)
                             tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg                                          <  1:1611548/71‑1 (MQ=255)
                              cGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg                                          >  1:435241/1‑70 (MQ=255)
                               ggCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAc                                        <  1:3251052/71‑1 (MQ=255)
                                   tGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:1698773/43‑1 (MQ=255)
                                     gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGcc                                   <  1:377194/70‑1 (MQ=255)
                                     gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg                                  >  1:829858/1‑71 (MQ=255)
                                      tGCTGAATCTCGATAAGTGCATCGGCTGCCAAACCTg                                                                   >  1:620591/1‑37 (MQ=255)
                                       gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa                                                   <  1:2487693/52‑1 (MQ=255)
                                        cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                <  1:2688763/70‑1 (MQ=255)
                                        cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                >  1:696405/1‑70 (MQ=255)
                                        cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                >  1:701179/1‑70 (MQ=255)
                                         tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:2638249/37‑1 (MQ=255)
                                         tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc                                                                <  1:883458/37‑1 (MQ=255)
                                         tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                                               <  1:2111351/38‑1 (MQ=255)
                                         tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa                                                               <  1:3554844/38‑1 (MQ=255)
                                          gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg                             >  1:2275248/1‑71 (MQ=255)
                                            acctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc                           <  1:3752733/69‑1 (MQ=255)
                                             tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                                                          <  1:1352567/39‑1 (MQ=255)
                                              ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                           <  1:776056/69‑1 (MQ=255)
                                              ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc                           <  1:1376149/69‑1 (MQ=255)
                                               tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                            <  1:2048618/35‑1 (MQ=255)
                                               tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                            >  1:3667292/1‑35 (MQ=255)
                                               tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                            >  1:739812/1‑35 (MQ=255)
                                               tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTa                                                           >  1:535111/1‑36 (MQ=255)
                                               tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaaaa                                                 >  1:1745641/1‑46 (MQ=255)
                                                 gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                        <  1:1161257/69‑1 (MQ=255)
                                                 gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                        <  1:687478/69‑1 (MQ=255)
                                                 gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                        <  1:612248/69‑1 (MQ=255)
                                                      gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg                 <  1:372211/71‑1 (MQ=255)
                                                       tGCATCGGCTGCCACACCTGTTCAGTTACCTGTATAAACGTCTGGACCAGcc                                   >  1:2302690/1‑52 (MQ=255)
                                                       tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                           >  1:3829178/1‑44 (MQ=255)
                                                       tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                        >  1:908627/1‑63 (MQ=255)
                                                         cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt              >  1:947244/1‑71 (MQ=255)
                                                          aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAcc                                       >  1:654250/1‑45 (MQ=255)
                                                           tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                   >  1:2318982/1‑64 (MQ=255)
                                                             ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt          >  1:3936105/1‑71 (MQ=255)
                                                             ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt          >  1:3887534/1‑71 (MQ=255)
                                                               cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca        <  1:796317/71‑1 (MQ=255)
                                                                 gCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                <  1:771015/45‑1 (MQ=255)
                                                                  ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac      <  1:815090/70‑1 (MQ=255)
                                                                  ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac      <  1:3895725/70‑1 (MQ=255)
                                                                  ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa                     <  1:1636282/55‑1 (MQ=255)
                                                                    acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                       <  1:3749559/51‑1 (MQ=255)
                                                                    acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                   >  1:3281131/1‑55 (MQ=255)
                                                                    acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                   >  1:247609/1‑55 (MQ=255)
                                                                    acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac   >  1:2683398/1‑71 (MQ=255)
                                                                     cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg  >  1:2469016/1‑71 (MQ=255)
                                                                     |                                                                      
AGAGCGTAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACG  >  minE/835273‑835412

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: