Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,354 T→G 13.7% V24V (GTT→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,3540TG13.7% 149.4 / 28.2 73V24V (GTT→GTGnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (34/29);  new base G (2/8);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAC  >  minE/835284‑835419
                                                                      |                                                                 
tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                                   <  1:2563524/71‑1 (MQ=255)
  aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGTCTGCCACACCTGTTCAGTTAc                                                                 >  1:1714717/1‑71 (MQ=255)
  aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc                                                                 >  1:525186/1‑71 (MQ=255)
         gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                          <  1:3949171/71‑1 (MQ=255)
         gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                          <  1:1209399/71‑1 (MQ=255)
         gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                          <  1:584646/71‑1 (MQ=255)
         gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                          <  1:2609513/71‑1 (MQ=255)
          ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa                                                          <  1:2759346/70‑1 (MQ=255)
            cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc                                                       >  1:218643/1‑71 (MQ=255)
             acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                                                       <  1:1398341/70‑1 (MQ=255)
             acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc                                                       <  1:1864425/70‑1 (MQ=255)
              caAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt                                                     <  1:2737055/71‑1 (MQ=255)
               aaGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCt                                                               >  1:3289879/1‑60 (MQ=255)
                 gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                                  <  1:142995/71‑1 (MQ=255)
                 gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                                  <  1:3414327/71‑1 (MQ=255)
                 gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg                                                  >  1:2017075/1‑71 (MQ=255)
                  tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg                                                 <  1:1611548/71‑1 (MQ=255)
                   cGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg                                                 >  1:435241/1‑70 (MQ=255)
                    ggCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAc                                               <  1:3251052/71‑1 (MQ=255)
                          gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGcc                                          <  1:377194/70‑1 (MQ=255)
                          gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg                                         >  1:829858/1‑71 (MQ=255)
                            gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa                                                          <  1:2487693/52‑1 (MQ=255)
                             cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                       >  1:696405/1‑70 (MQ=255)
                             cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                       <  1:2688763/70‑1 (MQ=255)
                             cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                       >  1:701179/1‑70 (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:3667292/1‑35 (MQ=255)
                                    tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                                   >  1:739812/1‑35 (MQ=255)
                                    tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt                                                                   <  1:2048618/35‑1 (MQ=255)
                                    tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTa                                                                  >  1:535111/1‑36 (MQ=255)
                                    tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaaaa                                                        >  1:1745641/1‑46 (MQ=255)
                                      gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                               <  1:612248/69‑1 (MQ=255)
                                      gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                               <  1:1161257/69‑1 (MQ=255)
                                      gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg                               <  1:687478/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:3887534/1‑71 (MQ=255)
                                                  ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt                 >  1:3936105/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)
                                                          tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg         >  1:244932/2‑71 (MQ=255)
                                                          cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg         >  1:2469016/1‑71 (MQ=255)
                                                            ccTGTTCGGTGACCTGTAAAAACGTCTGTACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc       >  1:2059072/1‑70 (MQ=255)
                                                            ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                          >  1:1167491/1‑52 (MQ=255)
                                                            ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg       <  1:3461557/71‑1 (MQ=255)
                                                             cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa                            <  1:3328891/49‑1 (MQ=255)
                                                              tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt                  <  1:2760596/58‑1 (MQ=255)
                                                              tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa           >  1:2640198/1‑65 (MQ=255)
                                                               gTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa           >  1:3949634/1‑64 (MQ=255)
                                                                ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATg                               <  1:940882/43‑1 (MQ=16)
                                                                ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt                 <  1:3319294/57‑1 (MQ=255)
                                                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg                                       <  1:512008/35‑1 (MQ=255)
                                                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                          >  1:2868468/1‑48 (MQ=255)
                                                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                >  1:2777811/1‑58 (MQ=255)
                                                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                >  1:197828/1‑58 (MQ=255)
                                                                ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa    <  1:1224112/70‑1 (MQ=255)
                                                                 tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac          >  1:2321343/1‑63 (MQ=255)
                                                                  cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg       >  1:2559532/1‑65 (MQ=255)
                                                                  cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAc  >  1:1769383/1‑70 (MQ=255)
                                                                    gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg                              <  1:1442120/40‑1 (MQ=255)
                                                                    gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata                          >  1:2612771/1‑44 (MQ=255)
                                                                      tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAat                           >  1:3478885/1‑41 (MQ=255)
                                                                      tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTgg                   >  1:236573/1‑49 (MQ=255)
                                                                      |                                                                 
TGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAC  >  minE/835284‑835419

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: