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

GCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTG  >  minE/1891954‑1892067
                                                                  ||                                                
gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                 <  1:1434005/69‑4 (MQ=255)
       ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                        <  1:324613/71‑1 (MQ=255)
       ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                        <  1:1261442/71‑1 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                  >  1:1002852/1‑71 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                  >  1:3647799/1‑71 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                  >  1:3110810/1‑71 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                  >  1:2372123/1‑71 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                  >  1:407874/1‑71 (MQ=255)
                     attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                          <  1:491213/71‑1 (MQ=255)
                     attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                          <  1:2372405/71‑1 (MQ=255)
                     attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                          <  1:2425032/71‑1 (MQ=255)
                      ttcAGCACCTCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:3877952/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:3803615/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:3516607/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:3482302/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:3152924/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:250773/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:2196878/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:584087/55‑1 (MQ=255)
                      ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:470800/55‑1 (MQ=255)
                              cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                 <  1:3871084/71‑1 (MQ=255)
                              cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                 <  1:1909959/71‑1 (MQ=255)
                                cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                          >  1:3907706/1‑60 (MQ=255)
                                cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                          >  1:2939470/1‑60 (MQ=255)
                                cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                          >  1:3837474/1‑60 (MQ=255)
                                cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                 >  1:1216600/1‑69 (MQ=255)
                                  aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAg                                         <  1:111972/43‑1 (MQ=255)
                                     aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                            >  1:760702/1‑53 (MQ=255)
                                     aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                            >  1:945508/1‑53 (MQ=255)
                                     aCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                            >  1:2098794/1‑53 (MQ=255)
                                            ttACCGCTCGCTGAACATATCAGCTGAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt   >  1:1814190/1‑71 (MQ=255)
                                             tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                  >  1:3665978/1‑39 (MQ=37)
                                             tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt   >  1:2241214/1‑70 (MQ=255)
                                             tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt   >  1:1348671/1‑70 (MQ=255)
                                              aCCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTg  <  1:3667952/70‑1 (MQ=255)
                                                    cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCa           <  1:3873461/55‑1 (MQ=255)
                                                    cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCa           <  1:2738251/55‑1 (MQ=255)
                                                      ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc         >  1:3719009/1‑55 (MQ=255)
                                                      ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc         >  1:2986920/1‑55 (MQ=255)
                                                      ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc         >  1:1678255/1‑55 (MQ=255)
                                                      ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc         >  1:1066412/1‑55 (MQ=255)
                                                      ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt   <  1:570368/61‑1 (MQ=38)
                                                      ctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt   <  1:1029431/61‑1 (MQ=38)
                                                                  ||                                                
GCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTG  >  minE/1891954‑1892067

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: