Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 67,373 T→G 36.6% I299I (ATA→ATC leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE67,3730TG36.6% 26.0 / 43.8 52I299I (ATA→ATCleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base T (6/27);  new base G (8/11);  total (14/38)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

ACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCA  >  minE/67310‑67434
                                                               |                                                             
aCCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATAc                                                           <  1:2104952/68‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATATCggagg                                                     <  1:441469/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2845526/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:444822/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2957173/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:1494236/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2572708/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:390579/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:3244746/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2475742/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:1769173/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:1938951/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2360249/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2255160/69‑1 (MQ=255)
     tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:2077608/69‑1 (MQ=255)
     aataTTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCggagg                                                     <  1:981394/66‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:2002410/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:3208649/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:3177830/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:1988453/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:2407626/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:1727194/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:2595626/68‑1 (MQ=255)
                            tAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGcc                               <  1:1333375/68‑1 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:3189514/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:40713/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:2920317/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:113701/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:1146017/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:1151966/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:1714183/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            >  1:2282342/1‑69 (MQ=255)
                              gTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            <  1:1084642/69‑1 (MQ=255)
                              gTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaac                            <  1:3207759/69‑1 (MQ=255)
                                 ttCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaacaa                          >  1:1658292/1‑68 (MQ=255)
                                 ttCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAAt                         >  1:44899/1‑69 (MQ=255)
                                 ttCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAAt                         >  1:2942100/1‑69 (MQ=255)
                                 ttCGCGGTTTTTCAGCACGCAATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCcaacaa                          >  1:3133332/1‑68 (MQ=255)
                                         ttttCAGTACACTATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:2807979/44‑1 (MQ=25)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:2852562/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:791427/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:2795113/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:3428923/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:1993936/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:2340772/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:3147708/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:3120243/44‑1 (MQ=37)
                                         ttttCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:2666574/44‑1 (MQ=37)
                                         ttttCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCtttgt                 >  1:2769140/1‑69 (MQ=255)
                                         ttttCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCtttgt                 >  1:2848113/1‑69 (MQ=255)
                                            tCAGTACACCATCCTGGTGGATACCGGAGGAGTGTGCGAAt                                          <  1:2901649/41‑1 (MQ=37)
                                                        ccTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCa  <  1:1707229/69‑1 (MQ=255)
                                                               |                                                             
ACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCA  >  minE/67310‑67434

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: