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

GGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGT  >  minE/1891953‑1892072
                                                                   ||                                                     
ggCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                      >  1:3243778/1‑67 (MQ=255)
ggCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                      >  1:2319323/1‑67 (MQ=255)
ggCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                      >  1:2284766/1‑67 (MQ=255)
ggCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagc                                                       >  1:1021943/1‑67 (MQ=255)
 gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                      >  1:3479105/1‑66 (MQ=255)
 gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                      <  1:4616267/69‑4 (MQ=255)
        ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                             <  1:4203993/71‑1 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                          <  1:596780/70‑1 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                          <  1:3405244/70‑1 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                          <  1:2517624/70‑1 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                          <  1:2674983/70‑1 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                         >  1:3434423/1‑71 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                         >  1:2256463/1‑71 (MQ=255)
            aaaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCg                                         >  1:3115388/1‑71 (MQ=255)
             aaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                        <  1:518472/71‑1 (MQ=255)
                 attcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCg                                    <  1:4189847/71‑1 (MQ=255)
                    cgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                  <  1:1824193/70‑1 (MQ=255)
                    cgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                  <  1:1920085/70‑1 (MQ=255)
                    cgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                  <  1:2282658/70‑1 (MQ=255)
                    cgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGa                                  <  1:1399967/70‑1 (MQ=255)
                       ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                 >  1:4846794/1‑68 (MQ=255)
                       ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                 >  1:942911/1‑68 (MQ=255)
                         cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                             <  1:2575501/70‑1 (MQ=255)
                                   aaaaCGGCATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                >  1:670825/1‑57 (MQ=255)
                                   aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                        >  1:2688104/1‑49 (MQ=255)
                                   aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcg                       >  1:2599723/1‑66 (MQ=255)
                                   aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                  >  1:4972076/1‑71 (MQ=255)
                                   aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtc                   >  1:2932908/1‑70 (MQ=255)
                                    aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                 <  1:1704781/71‑1 (MQ=255)
                                    aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                 <  1:4618371/71‑1 (MQ=255)
                                    aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                 <  1:2011673/71‑1 (MQ=255)
                                     aaCGACATTTACCGCTCGCTGAACATATCAGCTACGCACAGAACCg                                         >  1:4619551/1‑46 (MQ=255)
                                     aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCa                >  1:52402/1‑71 (MQ=255)
                                       cGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc              >  1:568446/1‑71 (MQ=255)
                                       cGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc              <  1:4239018/71‑1 (MQ=255)
                                       cGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc              >  1:3480084/1‑71 (MQ=255)
                                       cGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc              <  1:3356810/71‑1 (MQ=255)
                                       cGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATc              <  1:2538763/71‑1 (MQ=255)
                                         aCATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg            <  1:1804166/71‑1 (MQ=255)
                                         aCATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCg            <  1:2513307/71‑1 (MQ=255)
                                            tttCCCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:1816485/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                           <  1:2441151/37‑1 (MQ=25)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                           <  1:2076340/37‑1 (MQ=25)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                           <  1:874808/37‑1 (MQ=25)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                           <  1:4680449/37‑1 (MQ=25)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:4529556/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:2849025/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:304020/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:881375/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:2683892/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:3688164/55‑1 (MQ=255)
                                            tttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:2114002/55‑1 (MQ=255)
                                             ttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCa                >  1:5158971/1‑63 (MQ=255)
                                             ttACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt        >  1:2060506/1‑71 (MQ=39)
                                              tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                             >  1:3274483/1‑49 (MQ=255)
                                              tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                  >  1:3471539/1‑60 (MQ=255)
                                              tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtct                  >  1:2918064/1‑60 (MQ=255)
                                                   ctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt  <  1:3101642/71‑1 (MQ=255)
                                                    tcgctGAACATATCAGCTAAGCACAGACCCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt  <  1:2782315/70‑1 (MQ=255)
                                                    tcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGt  <  1:3936276/70‑1 (MQ=255)
                                                     cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                <  1:840843/39‑1 (MQ=37)
                                                     cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                <  1:1059162/39‑1 (MQ=37)
                                                     cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                <  1:1468492/39‑1 (MQ=37)
                                                     cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                <  1:4824030/39‑1 (MQ=37)
                                                     cgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATa                                <  1:4068402/39‑1 (MQ=37)
                                                      gctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGtt                         <  1:2766192/45‑1 (MQ=255)
                                                                   ||                                                     
GGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGT  >  minE/1891953‑1892072

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: