Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,890,961 G→A 69.4% 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,9610GA69.4% 65.2 / 39.9 49P208S (CCG→TCG) gudPpredicted D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base G (5/10);  new base A (18/16);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGAC  >  minE/1890893‑1891020
                                                                    |                                                           
ttGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGa                                                            <  1:3305279/70‑1 (MQ=255)
ttGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGa                                                            <  1:2310563/70‑1 (MQ=255)
       tCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGAt                                                     <  1:3390391/70‑1 (MQ=255)
           tatTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGc                                                >  1:766399/1‑71 (MQ=255)
           tatTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGc                                                >  1:1252820/1‑71 (MQ=255)
           tatTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGc                                                >  1:2031621/1‑71 (MQ=255)
           tatTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGc                                                >  1:2060957/1‑71 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTc                                              >  1:738350/1‑71 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTc                                              >  1:4244216/1‑71 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCt                                               <  1:2809255/70‑1 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTc                                              >  1:1139546/1‑71 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTc                                              >  1:4188039/1‑71 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTc                                              >  1:3651435/1‑71 (MQ=255)
             ttGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTc                                              >  1:190255/1‑71 (MQ=255)
                aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTg                                           <  1:987301/71‑1 (MQ=255)
                aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTg                                           <  1:2711960/71‑1 (MQ=255)
                aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTg                                           <  1:2179619/71‑1 (MQ=255)
                aTCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAAGGTTGATTTGGCTCGTg                                           <  1:3732988/71‑1 (MQ=255)
                    gcgcACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                       <  1:1898702/71‑1 (MQ=255)
                    gcgcACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                       <  1:180146/71‑1 (MQ=255)
                     cgcACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                       <  1:167243/70‑1 (MQ=255)
                          caccCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTc                                 <  1:416127/71‑1 (MQ=255)
                           accCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCa                                <  1:3782310/71‑1 (MQ=255)
                            cccGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTg                                           <  1:3374103/59‑1 (MQ=255)
                            cccGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCaa                               >  1:4020611/1‑71 (MQ=255)
                            cccGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGtt                                           <  1:3219993/59‑2 (MQ=255)
                                        tACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCt                   >  1:2783679/1‑71 (MQ=255)
                                        tACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCt                   <  1:469009/71‑1 (MQ=255)
                                         aCTCCAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACtt                                   >  1:1868096/1‑54 (MQ=255)
                                            ccAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGc                    >  1:1210315/1‑66 (MQ=255)
                                             cAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatga              >  1:1027094/1‑71 (MQ=255)
                                             cAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatg               <  1:3558120/70‑1 (MQ=255)
                                             cAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatg               <  1:1280107/70‑1 (MQ=255)
                                             cAGCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCtgatg               <  1:3514757/70‑1 (MQ=255)
                                               gCTCTTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACa            <  1:3231462/71‑1 (MQ=255)
                                                ctctTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGAc                                     >  1:3596497/1‑45 (MQ=255)
                                                ctctTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGAc                                     >  1:370118/1‑45 (MQ=255)
                                                ctctTTCTTATTTACCCCCGAATGTTGATTTGGCTCGTGGATGAc                                     >  1:4206692/1‑45 (MQ=255)
                                                   tttcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                       <  1:1764310/40‑1 (MQ=255)
                                                   tttcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATg                                       <  1:2976259/40‑1 (MQ=255)
                                                     tcttATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATccc       >  1:1777746/1‑70 (MQ=255)
                                                     tcttATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATc         <  1:3256431/68‑1 (MQ=255)
                                                     tcttATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATc         <  1:3079609/68‑1 (MQ=255)
                                                      cttATTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATccc       >  1:3856114/1‑69 (MQ=255)
                                                         aTTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCa     <  1:1298944/68‑1 (MQ=255)
                                                         aTTTACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGAc  >  1:3625218/1‑71 (MQ=255)
                                                          tttACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCa     <  1:3512657/67‑1 (MQ=255)
                                                           ttACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACaa           >  1:2822508/1‑60 (MQ=255)
                                                           ttACCCCCGAATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACaa           >  1:1806485/1‑60 (MQ=255)
                                                                    |                                                           
TTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGAC  >  minE/1890893‑1891020

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: