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% 183.1 / NA 52intergenic (‑437/+54)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (20/32);  total (20/32)
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% 181.6 / NA 52intergenic (‑437/+54)gudP/ygdLpredicted D‑glucarate transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (20/32);  total (20/32)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAC  >  minE/1891950‑1892077
                                                                      ||                                                          
aCTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCa                                                              >  1:740926/1‑70 (MQ=255)
 cTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagc                                                            >  1:1629793/1‑69 (MQ=255)
    gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                           <  1:1745533/69‑4 (MQ=255)
           ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                                   <  1:2957963/70‑1 (MQ=255)
               aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                               <  1:143175/70‑1 (MQ=255)
               aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                              >  1:336790/1‑71 (MQ=255)
                aaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCag                                                             >  1:536265/1‑54 (MQ=255)
                 aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                            >  1:3010603/1‑71 (MQ=255)
                 aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                            <  1:393377/71‑1 (MQ=255)
                   gattcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTc                                          <  1:312230/71‑1 (MQ=255)
                     ttcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCgg                                        <  1:557825/71‑1 (MQ=255)
                     ttcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCg                                         >  1:3286337/1‑70 (MQ=255)
                      tcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCgg                                        <  1:3942488/70‑1 (MQ=255)
                       cgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                       <  1:1175395/70‑1 (MQ=255)
                          ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                    <  1:2652314/70‑1 (MQ=255)
                           tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                             >  1:3517014/1‑60 (MQ=255)
                           tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGGTa                                     >  1:3352307/1‑68 (MQ=255)
                                aCCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGcaca                                                    >  1:829440/1‑48 (MQ=255)
                                 ccGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                >  1:1228281/1‑51 (MQ=255)
                                  cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                           <  1:3899473/71‑1 (MQ=255)
                                   gCTAAAACGACATTTACCGCTCGCTGAACATATCag                                                             >  1:837498/1‑35 (MQ=255)
                                   gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                               >  1:370766/1‑50 (MQ=255)
                                   gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                               >  1:3669917/1‑50 (MQ=255)
                                   gCTAAAACGACATTTAACGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                          >  1:2022162/1‑71 (MQ=255)
                                   gCTAAAACGACATTTAACGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                          >  1:2026961/1‑71 (MQ=255)
                                    cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                >  1:793316/1‑48 (MQ=255)
                                    cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGt                         <  1:3546789/71‑1 (MQ=255)
                                    cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGt                         >  1:702564/1‑71 (MQ=255)
                                      aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                                   >  1:2819454/1‑43 (MQ=255)
                                       aaaCGACATTTACCTCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                      <  1:2006219/71‑1 (MQ=255)
                                       aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                       <  1:539328/70‑1 (MQ=255)
                                       aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                       >  1:3718435/1‑70 (MQ=255)
                                       aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                       <  1:1878400/70‑1 (MQ=255)
                                       aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                       <  1:1545806/70‑1 (MQ=255)
                                        aaCGACATTTACCGCTCGCTTAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                      <  1:1148882/54‑1 (MQ=255)
                                        aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCg                                         <  1:2246974/51‑1 (MQ=255)
                                        aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCg                                         <  1:1898920/51‑1 (MQ=255)
                                        aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                       <  1:572696/69‑1 (MQ=255)
                                         aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                  <  1:2584143/41‑1 (MQ=255)
                                         aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCa                     >  1:1584394/1‑70 (MQ=255)
                                          cGACATTTTCCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGGCGGATAAGACGTTCGCGTCTCATc                   <  1:378932/71‑1 (MQ=255)
                                           gACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATcc                  <  1:2899970/71‑1 (MQ=255)
                                            aCATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc                                 <  1:3135646/55‑1 (MQ=255)
                                             cATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg                 <  1:2989768/70‑1 (MQ=255)
                                             cATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGa                <  1:841773/71‑1 (MQ=255)
                                               tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                      <  1:2889541/63‑1 (MQ=255)
                                               tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg                 <  1:686545/68‑1 (MQ=255)
                                                ttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt             >  1:3285646/1‑71 (MQ=255)
                                                 tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                     >  1:1318459/1‑46 (MQ=255)
                                                  aCCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                          >  1:3978994/1‑56 (MQ=255)
                                                   ccGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTg            <  1:2408022/69‑1 (MQ=255)
                                                    cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc                                 >  1:4141772/1‑47 (MQ=255)
                                                    cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGAtt         <  1:3892839/71‑1 (MQ=255)
                                                      ctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt       <  1:1032648/71‑1 (MQ=255)
                                                      ctcgctGAACATATCAGCTAAGCAC‑‑AACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGttt     <  1:1684202/71‑1 (MQ=255)
                                                       tcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                       <  1:1604156/38‑1 (MQ=25)
                                                             aaCATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAc  <  1:881496/69‑1 (MQ=38)
                                                                      ||                                                          
ACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAC  >  minE/1891950‑1892077

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