Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,612,205 G→A 84.1% G316G (GGC→GGT 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,2050GA84.1% 90.4 / 7.4 56G316G (GGC→GGTeutAreactivating factor for ethanolamine ammonia lyase
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base A (31/16);  total (38/18)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGGCGCG  >  minE/1612148‑1612272
                                                         |                                                                   
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:102936/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:998095/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:866840/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:3334848/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:3152731/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:313294/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:1878216/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:1849652/1‑69 (MQ=255)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAAt                                                          >  1:1774492/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)
cGGCAGTTGTACGCCCTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGTATGAGCACCCGCGCCAATc                                                         >  1:1583815/1‑69 (MQ=255)
   cAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCAc                                                       <  1:3023099/69‑1 (MQ=255)
    aGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCAc                                                       <  1:2871924/68‑1 (MQ=255)
        gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg                                                  >  1:3017057/1‑69 (MQ=255)
        gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg                                                  >  1:764178/1‑69 (MQ=255)
        gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg                                                  >  1:1602923/1‑69 (MQ=255)
        gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg                                                  >  1:2911726/1‑69 (MQ=255)
        gTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTgg                                                  >  1:1398582/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:784915/1‑65 (MQ=255)
             cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc                                                 >  1:1758738/1‑65 (MQ=255)
             cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc                                                 >  1:307398/1‑65 (MQ=255)
             cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc                                                 >  1:403387/1‑65 (MQ=255)
             cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc                                                 >  1:2656726/1‑65 (MQ=255)
             cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc                                                 >  1:2733493/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:404915/1‑65 (MQ=255)
             cccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGgc                                                 >  1:3176119/1‑65 (MQ=255)
              ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa                                            <  1:1206676/69‑1 (MQ=255)
              ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa                                            <  1:2138397/69‑1 (MQ=255)
              ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa                                            <  1:1368181/69‑1 (MQ=255)
              ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa                                            <  1:2479775/69‑1 (MQ=255)
              ccTCCAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCa                                            <  1:863950/69‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:579800/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:346073/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:1666298/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:3182869/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:1109757/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:1399003/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:1281533/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:1887369/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:3182068/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:2232249/67‑1 (MQ=255)
                 ccaGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGAGCACCCGCGCCAATCACCGTGGCGCGCACgg                                         <  1:2405097/67‑1 (MQ=255)
                                                 gTATGCGCGCCCGCGCCAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGc         >  1:1975564/1‑69 (MQ=255)
                                                 gTATGCGCGCCCGCGCCAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGc         >  1:88357/1‑69 (MQ=255)
                                                      cgcgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGgcg    >  1:2506186/1‑69 (MQ=255)
                                                      cgcgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGgcg    >  1:1045077/1‑69 (MQ=255)
                                                      cgcgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGgcg    >  1:77917/1‑69 (MQ=255)
                                                      cgcgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGgcg    >  1:2498585/1‑69 (MQ=255)
                                                      cgcgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGgcg    >  1:3382777/1‑69 (MQ=255)
                                                      cgcgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGACGGAAACTGTACGTTCATCTCACGCAGgcg    >  1:1103679/1‑69 (MQ=255)
                                                        cgcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGGcgcg  <  1:1081036/69‑1 (MQ=255)
                                                         gcccgcgccAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGGcgcg  <  1:2980047/68‑1 (MQ=255)
                                                         |                                                                   
CGGCAGTTGTACGCCCTCAAGCCAGATTGTGCTGCCAGAGAGCGAAAGGGTATGCGCGCCCGCGCCAATCACCGTGGCGCGCACGGTTTGCGCCGGAAACTGTACGTTCATCTCACGCAGGCGCG  >  minE/1612148‑1612272

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: