Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,890,961 G→A 14.3% P208S (CCG→TCG)  gudP ← predicted D‑glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,890,9610GA14.3% 90.5 / 7.4 35P208S (CCG→TCG) gudPpredicted D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base G (20/10);  new base A (2/3);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAG  >  minE/1890896‑1891031
                                                                 |                                                                      
cTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGtt                                                                   >  1:670733/1‑71 (MQ=255)
cTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCAGGATGtt                                                                   <  1:603133/71‑1 (MQ=255)
 tGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGAt                                                                      >  1:1283002/1‑67 (MQ=255)
   aTCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGa                                                                 <  1:946144/70‑1 (MQ=255)
        tatTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGAt                                                                <  1:1983353/66‑1 (MQ=255)
          ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTg                                                             >  1:1930399/1‑67 (MQ=255)
          ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCt                                                          >  1:483951/1‑70 (MQ=255)
             aTCAGCGCACCACCCGCGGCGATGTACTTCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCt                                                          <  1:1802614/67‑1 (MQ=255)
             aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCt                                                          <  1:1590880/67‑1 (MQ=255)
             aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTg                                                      >  1:258152/1‑71 (MQ=255)
             aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTg                                                      <  1:2188296/71‑1 (MQ=255)
              tCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTgg                                                     <  1:51996/71‑1 (MQ=255)
               cAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGAt                                                                >  1:864411/1‑59 (MQ=255)
               cAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTc                                                         >  1:1401316/1‑66 (MQ=255)
               cAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGa                                                    >  1:257454/1‑71 (MQ=255)
               cAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACACCCGGATGTTGAt                                                                >  1:1051786/1‑59 (MQ=255)
                aGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGAt                                                   <  1:214041/71‑1 (MQ=255)
                aGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGAt                                                   <  1:923021/71‑1 (MQ=255)
                aGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTAGTGGAt                                                   <  1:4189/71‑1 (MQ=255)
                        accCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTc                                            >  1:1658827/1‑70 (MQ=255)
                            gcggcgATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAAc                                         >  1:346490/1‑69 (MQ=255)
                               gcgATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGa                                     >  1:506597/1‑70 (MQ=255)
                                     tACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCt                              <  1:1858650/71‑1 (MQ=255)
                                       cTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGAc                                                <  1:2085174/51‑1 (MQ=255)
                                          cAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatga                         >  1:1605921/1‑71 (MQ=255)
                                           aGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatga                         >  1:1020219/1‑70 (MQ=255)
                                             ctctTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGAc                                                <  1:825928/45‑1 (MQ=255)
                                             ctctTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGAGTTTCaa                                          >  1:1830379/1‑51 (MQ=255)
                                              tctttcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAAt                     >  1:1742090/1‑71 (MQ=255)
                                               ctttcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACttt                                             >  1:575220/1‑46 (MQ=255)
                                               ctttcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATc                    >  1:1101378/1‑71 (MQ=255)
                                               ctttcttATTTACCCCCGGATGTTGAATTGGCTCGTGGATGACTTTCAAACAGATAAAGCTGATGACAATc                    >  1:17344/1‑71 (MQ=255)
                                                  tcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATgg                                                 >  1:2168546/1‑38 (MQ=255)
                                                  tcttATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACaa                      >  1:302990/1‑66 (MQ=255)
                                                         tACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAAcc                                        >  1:38343/1‑41 (MQ=255)
                                                                 gATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAaag  <  1:1166216/71‑1 (MQ=255)
                                                                 gATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAaag  <  1:771309/71‑1 (MQ=255)
                                                                 |                                                                      
CTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAG  >  minE/1890896‑1891031

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: