Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,892,019:1 +GC 100% intergenic (‑437/+54) gudP ← / ← ygdL predicted D‑glucarate transporter/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,892,0191.G100.0% 209.3 / NA 60intergenic (‑437/+54)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (33/27);  total (33/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*minE1,892,0192.C100.0% 208.5 / NA 60intergenic (‑437/+54)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (33/27);  total (33/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAA  >  minE/1891951‑1892076
                                                                     ||                                                         
cTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCag                                                            <  1:789688/70‑2 (MQ=255)
   gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                          <  1:1225496/69‑4 (MQ=255)
          ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                                  >  1:414005/1‑70 (MQ=255)
            tGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                               >  1:2487044/1‑71 (MQ=255)
            tGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                               >  1:2364049/1‑71 (MQ=255)
             gAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                              >  1:913381/1‑71 (MQ=255)
             gAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                              >  1:3827292/1‑71 (MQ=255)
              aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                              <  1:3916246/70‑1 (MQ=255)
               aaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                            >  1:335679/1‑71 (MQ=255)
                 aGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGt                                          <  1:285135/71‑1 (MQ=255)
                  gattcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTc                                         >  1:1044333/1‑71 (MQ=255)
                  gattcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTc                                         >  1:122424/1‑71 (MQ=255)
                  gattcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTc                                         >  1:2337461/1‑71 (MQ=255)
                    ttcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCgg                                       <  1:1180183/71‑1 (MQ=255)
                    ttcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCgg                                       <  1:3999776/71‑1 (MQ=255)
                    ttcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCg                                        >  1:42657/1‑70 (MQ=255)
                     tcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                      >  1:220326/1‑71 (MQ=255)
                       gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                     <  1:2595183/70‑1 (MQ=255)
                       gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                     >  1:1632627/1‑70 (MQ=255)
                       gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                    <  1:1289629/71‑1 (MQ=255)
                       gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                    <  1:4113041/71‑1 (MQ=255)
                        attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGcaca                                                   >  1:3111225/1‑55 (MQ=255)
                        attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGcaca                                                   >  1:1596669/1‑55 (MQ=255)
                        attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                   <  1:64744/71‑1 (MQ=255)
                         ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                     >  1:308400/1‑68 (MQ=255)
                         ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATAAGCTAAGCACAGAACCGGGTCGGATaaa                                  >  1:3507146/1‑70 (MQ=255)
                          tcAGCACCTCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                 >  1:836157/1‑71 (MQ=255)
                          tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAg                                  >  1:2323771/1‑70 (MQ=255)
                           cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                 <  1:3947748/70‑1 (MQ=255)
                           cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc                                >  1:907608/1‑71 (MQ=255)
                           cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc                                >  1:1707263/1‑71 (MQ=255)
                              cACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                             >  1:640098/1‑71 (MQ=255)
                                 cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                          <  1:1336674/71‑1 (MQ=255)
                                  gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                   >  1:2442540/1‑61 (MQ=255)
                                   cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTc                                         <  1:460337/54‑1 (MQ=255)
                                   cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCg                                        <  1:2883854/55‑1 (MQ=255)
                                    tAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                     <  1:2452528/57‑1 (MQ=255)
                                     aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                      >  1:2008661/1‑71 (MQ=255)
                                      aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                                  <  1:4094262/42‑1 (MQ=255)
                                      aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                                  <  1:150945/42‑1 (MQ=255)
                                      aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                     <  1:4092920/71‑1 (MQ=255)
                                      aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                      >  1:1143293/1‑70 (MQ=255)
                                       aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                             >  1:486456/1‑46 (MQ=255)
                                        aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                     <  1:1200719/53‑1 (MQ=255)
                                        aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                     <  1:487517/53‑1 (MQ=255)
                                        aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                          >  1:2580096/1‑64 (MQ=255)
                                        aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCAt                   <  1:1076995/71‑1 (MQ=255)
                                        aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCAt                   <  1:1657670/71‑1 (MQ=255)
                                         cGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc                  >  1:1708637/1‑71 (MQ=255)
                                            cATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                           <  1:1326638/43‑1 (MQ=255)
                                             aTTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg                <  1:1061905/69‑1 (MQ=255)
                                             aTTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGAc              <  1:745883/71‑1 (MQ=255)
                                              tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGGTa                                    >  1:3461465/1‑48 (MQ=255)
                                                tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGCTa                                    >  1:2983141/1‑46 (MQ=255)
                                                tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt            <  1:1835633/70‑1 (MQ=255)
                                                tACCGCTCGCTGAACATATCAGCTAAGCACAGAACAGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt            >  1:833959/1‑70 (MQ=255)
                                                   cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGa          >  1:3084193/1‑69 (MQ=255)
                                                    gctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTg       <  1:2490622/71‑1 (MQ=39)
                                                      tcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGtt     >  1:2917436/1‑71 (MQ=255)
                                                        gctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGttt    <  1:3200323/70‑1 (MQ=38)
                                                        gctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGtt     <  1:1473672/69‑1 (MQ=255)
                                                         ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTaa  >  1:2653013/1‑71 (MQ=255)
                                                                     ||                                                         
CTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAA  >  minE/1891951‑1892076

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: