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

GCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAC  >  minE/1891954‑1892077
                                                                  ||                                                          
gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                           >  1:1538546/1‑66 (MQ=255)
gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                           <  1:2080918/69‑4 (MQ=255)
gCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCagct                                                           >  1:1792924/1‑66 (MQ=255)
       ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                  <  1:2582492/71‑1 (MQ=255)
       ccTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                  <  1:2226297/71‑1 (MQ=255)
        cTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGaa                                                 <  1:2093584/71‑1 (MQ=255)
          gAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAcc                                               >  1:1766394/1‑71 (MQ=255)
           ataaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                <  1:48800/67‑1 (MQ=255)
            aaaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCgg                                             >  1:2586235/1‑71 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCa                                                              >  1:2058927/1‑53 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                            >  1:2232077/1‑71 (MQ=255)
             aaGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCggg                                            >  1:431714/1‑71 (MQ=255)
               gattcgattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTc                                          >  1:1038423/1‑71 (MQ=255)
                    gattcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGAt                                      <  1:243470/70‑1 (MQ=255)
                     attcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                    <  1:2201003/71‑1 (MQ=255)
                       tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                  >  1:1770963/1‑71 (MQ=255)
                       tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                  >  1:1986905/1‑71 (MQ=255)
                        cAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                  <  1:1319074/70‑1 (MQ=255)
                         aGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACg                                <  1:3030300/71‑1 (MQ=255)
                            aCCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTc                             <  1:695310/71‑1 (MQ=255)
                              cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATaa                                    <  1:90008/62‑1 (MQ=255)
                              cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                           <  1:1842036/71‑1 (MQ=255)
                              cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgc                           <  1:2255568/71‑1 (MQ=255)
                               gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGt                                           >  1:272461/1‑54 (MQ=255)
                               gCTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                          <  1:2068352/71‑1 (MQ=255)
                                cTAAAACGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGAc                                 >  1:1516061/1‑63 (MQ=255)
                                  aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                >  1:315656/1‑46 (MQ=255)
                                  aaaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                >  1:2881680/1‑46 (MQ=255)
                                  aaaaCGACATTTACCGCTCGCTGAACATAACAGCTAAGCACAGAAc                                                >  1:1059935/1‑46 (MQ=25)
                                   aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGtctc                      <  1:1190695/71‑1 (MQ=255)
                                   aaaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGt                         >  1:2150965/1‑68 (MQ=255)
                                    aaCGACATTTACCGCTCGCTGAACATATCAGCTAAGCACAGAAc                                                >  1:151205/1‑44 (MQ=255)
                                          aTTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGa                                  >  1:1444769/1‑52 (MQ=255)
                                          aTTTACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGAc               <  1:2778303/71‑1 (MQ=255)
                                           tttACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                  <  1:1543297/35‑1 (MQ=25)
                                           tttACCGCTCGCTGAACATATCAGCTAAGCACAGa                                                  <  1:65477/35‑1 (MQ=25)
                                             tACCGCTCGCTGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt             >  1:2122003/1‑70 (MQ=255)
                                                cgctcgctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTcgcg                          <  1:2915207/54‑1 (MQ=255)
                                                     gctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGttt     <  1:1745506/70‑1 (MQ=38)
                                                     gctGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGttt     <  1:715490/70‑1 (MQ=38)
                                                       tGAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAc  >  1:1068912/1‑71 (MQ=39)
                                                        gAACATATCAGCTAAGCACAGAACCGGGTCGGATAAGACGt                               >  1:2416937/1‑41 (MQ=37)
                                                                  ||                                                          
GCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCA‑‑TAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAC  >  minE/1891954‑1892077

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