Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,612,166 A→C 57.0% L329L (CTT→CTG eutA ← reactivating factor for ethanolamine ammonia lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,612,1660AC57.0% 34.0 / 55.9 65L329L (CTT→CTGeutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base A (26/2);  new base C (32/5);  total (58/7)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTGGCGCGCA  >  minE/1612108‑1612230
                                                          |                                                                
ccGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtg                                                        >  1:591625/1‑69 (MQ=255)
ccGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtg                                                        >  1:531204/1‑69 (MQ=255)
ccGTTTCATAAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATtg                                                        >  1:2625591/1‑69 (MQ=255)
         caATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGtgctgc                                                  >  1:1489821/1‑66 (MQ=255)
         caATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGtgctgc                                                  >  1:2241363/1‑66 (MQ=255)
         caATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGtgctgc                                                  >  1:188545/1‑66 (MQ=255)
         caATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGtgctgc                                                  >  1:1920795/1‑66 (MQ=255)
                                       gCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:152056/1‑51 (MQ=255)
                                       gCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:922416/1‑51 (MQ=255)
                                       gCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:2488898/1‑51 (MQ=255)
                                       gCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:2434346/1‑51 (MQ=255)
                                       gCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:2023468/1‑51 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGATAGGGTATGAGCACCCGCGCCaa                 >  1:3105535/1‑68 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCaa                 >  1:327672/1‑68 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCaa                 >  1:218769/1‑68 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCaa                 >  1:1852602/1‑68 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:313294/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:102936/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:3152731/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:3334848/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:866840/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:998095/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1703795/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1000425/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1186498/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1300854/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1774492/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1878216/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                >  1:1849652/1‑69 (MQ=255)
                                        cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGTATGAGCACCCGCGCCAATc               >  1:1583815/1‑69 (MQ=255)
                                           cAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCAc             <  1:3023099/69‑1 (MQ=255)
                                            aGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCAc             <  1:2871924/68‑1 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:3160942/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:2864271/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:45151/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:1948440/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:430853/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:3344259/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:30788/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:3102812/1‑43 (MQ=255)
                                               tGTACGCCCTCAAGCCAAATTGTGCTGCCAGAGAGCGAAAggg                                   >  1:1247307/1‑43 (MQ=255)
                                                gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg        >  1:1398582/1‑69 (MQ=255)
                                                gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg        >  1:1602923/1‑69 (MQ=255)
                                                gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg        >  1:3017057/1‑69 (MQ=255)
                                                gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg        >  1:2911726/1‑69 (MQ=255)
                                                gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg        >  1:764178/1‑69 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGagcacc                         >  1:2556952/1‑47 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATc               >  1:1813478/1‑57 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTgg        >  1:610329/1‑64 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:3176119/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:307398/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:403387/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:404915/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:425035/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:3028572/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:293919/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:2733493/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:2656726/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:1758738/1‑65 (MQ=255)
                                                     cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc       >  1:784915/1‑65 (MQ=255)
                                                      ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa  <  1:2138397/69‑1 (MQ=255)
                                                      ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa  <  1:2479775/69‑1 (MQ=255)
                                                      ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa  <  1:1206676/69‑1 (MQ=255)
                                                      ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa  <  1:863950/69‑1 (MQ=255)
                                                      ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa  <  1:1368181/69‑1 (MQ=255)
                                                          |                                                                
CCGTTTCATCAATCGGGATCGCCACCGGCAAATTGCGCAGCGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTGGCGCGCA  >  minE/1612108‑1612230

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: