Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,097,921 C→T 67.3% V242I (GTC→ATC)  garL ← alpha‑dehydro‑beta‑deoxy‑D‑glucarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,097,9210CT67.3% 42.6 / 29.0 58V242I (GTC→ATC) garLalpha‑dehydro‑beta‑deoxy‑D‑glucarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (14/5);  new base T (17/22);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTGCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGATAACGACGCGCATCGGCTTCGA  >  minE/2097859‑2097989
                                                              |                                                                    
ttgttgCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACg                                                                     >  1:2562841/1‑64 (MQ=255)
ttgttgCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACg                                                                     >  1:623677/1‑64 (MQ=255)
ttgttgCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACg                                                                     >  1:2230860/1‑64 (MQ=255)
  gttgCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCgag                                                                <  1:3306153/67‑1 (MQ=255)
  gttgCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCgag                                                                <  1:976501/67‑1 (MQ=255)
          ggtggtTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGcc                                                       >  1:3124601/1‑68 (MQ=255)
          ggtggtTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGcc                                                       >  1:1979741/1‑68 (MQ=255)
           gtggtTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGa                                                     <  1:1559194/69‑1 (MQ=255)
           gtggtTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGa                                                     <  1:1823144/69‑1 (MQ=255)
              gtTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACa                                                   >  1:858181/1‑68 (MQ=255)
                    ttttAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCcaca                                             >  1:319387/1‑68 (MQ=255)
                    ttttAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACaa                                            >  1:396805/1‑69 (MQ=255)
                    ttttAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACaa                                            >  1:1193086/1‑69 (MQ=255)
                    ttttAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACaa                                            >  1:2119134/1‑69 (MQ=255)
                                 aGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCa                               >  1:3118493/1‑69 (MQ=255)
                                 aGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCa                               >  1:277335/1‑69 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:2156562/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:2482155/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:2203022/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:542563/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:277093/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:2805958/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:3017684/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:3169717/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1122660/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1941076/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:193402/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1870005/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:799086/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1749663/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:841261/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1641866/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1478178/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCAAAAACGTcgc                                     <  1:1751369/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTCCCGACAGCCACAAACGTcgc                                     <  1:2238422/60‑1 (MQ=255)
                                    cAGTTTCTGAGTGGCAGAGCGGAAGATGCCGAGATCGCTACCGACAGCCACAAACGTcgc                                     <  1:662172/60‑1 (MQ=255)
                                         tCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGAt                       <  1:475778/69‑1 (MQ=255)
                                           tGAGTGGCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:1737885/53‑1 (MQ=255)
                                                 gCAGAGCGGATGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:97994/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACCGCCACAAACGTcgca                                    >  1:1196846/1‑47 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:493725/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:879925/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:2692846/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:382977/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:664944/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:3225814/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:3143848/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:3039775/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:1824136/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:1399795/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:1388502/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:1226909/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:1099951/1‑68 (MQ=255)
                                                 gCAGAGCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAac                    >  1:130264/1‑64 (MQ=255)
                                                 gCAGAGAGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTCGCACCCCATTCCAGATAacgacg                >  1:2197797/1‑68 (MQ=255)
                                                   agagCGGAAGATGCCGAGATCGCTGCCGACAGCCACAAACGTcgc                                     <  1:2834954/45‑1 (MQ=255)
                                                            gACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGATAACGACGCGCATCGGCTTc    >  1:2763067/1‑69 (MQ=255)
                                                            gACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGATAACGACGCGCATCGGCTTc    >  1:2442894/1‑69 (MQ=255)
                                                              cGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGATAACGACGCGCATCGGCTTCGa  >  1:491191/1‑69 (MQ=255)
                                                              cGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGATAACGACGCGCATCGGCTTCGa  >  1:1055495/1‑69 (MQ=255)
                                                              |                                                                    
TTGTTGCGGTGGTGGTTATTTTTTAAAGGTATCAGCCAGTTTCTGAGTGGCAGAGCGGAAGACGCCGAGATCGCTGCCGACAGCCACAAACGTCGCGCCCCATTCCAGATAACGACGCGCATCGGCTTCGA  >  minE/2097859‑2097989

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: