Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,698,757 C→T 39.4% L213L (CTA→TTA)  csiE → stationary phase inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,698,7570CT39.4% 38.6 / 48.9 76L213L (CTA→TTA) csiEstationary phase inducible protein
Reads supporting (aligned to +/- strand):  ref base C (33/13);  new base T (19/11);  total (52/24)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGGACAGACAA  >  minE/1698690‑1698819
                                                                   |                                                              
aaTTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                                                               >  1:3011443/1‑69 (MQ=255)
aaTTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                                                               >  1:1237520/1‑69 (MQ=255)
aaTTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                                                               >  1:733238/1‑69 (MQ=255)
aaTTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                                                               >  1:2322499/1‑69 (MQ=255)
aaTTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                                                               >  1:660090/1‑69 (MQ=255)
aaTTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCact                                                               >  1:2483451/1‑69 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:2432042/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:183187/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:2857065/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:284687/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:207466/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:3289979/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:3378408/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:854848/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:933724/68‑1 (MQ=255)
       cGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCGCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:2299322/68‑1 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:2499701/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:1651592/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:2183761/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:2317392/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:2904236/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:1747179/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:3124902/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:521133/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:654374/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:681409/1‑69 (MQ=255)
           aCTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCa                                                    >  1:1005066/1‑69 (MQ=255)
            atGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCAt                                                   >  1:3331989/2‑69 (MQ=255)
                                  ccTGTACCTCCTGATGAGCCACTGTTTATGGCATTACtttt                                                         <  1:1704320/41‑1 (MQ=38)
                                        ccaccTGATGAGCCACTGTTTATGGGACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGAt                       >  1:508826/1‑69 (MQ=255)
                                        ccaccTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCg                           >  1:1087641/1‑65 (MQ=255)
                                        ccaccTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGa                        >  1:103936/1‑68 (MQ=255)
                                        ccaccTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGAt                       >  1:1047604/1‑69 (MQ=255)
                                        ccaccTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGAt                       >  1:1338768/1‑69 (MQ=255)
                                        ccaccTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGAt                       >  1:2120799/1‑69 (MQ=255)
                                        ccaccTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGAt                       >  1:3054973/1‑69 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:1579385/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:228960/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:3247973/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:3158507/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:3431004/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:3075347/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:447560/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:1551540/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:1937776/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:570393/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:2008897/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:2030510/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:586194/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:606187/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:140758/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:2502898/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:2079224/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:875405/1‑39 (MQ=255)
                                           ccTGATGAGCCACTGTTTATGGCATTACTTTTTTCCATg                                                  >  1:2223944/1‑39 (MQ=255)
                                                 gaGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGAtg                                  >  1:3143104/1‑48 (MQ=255)
                                                 gaGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATc                                  >  1:70019/1‑49 (MQ=255)
                                                 gaGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATc                                  >  1:1555274/1‑49 (MQ=255)
                                                 gaGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATc                                  >  1:1580574/1‑49 (MQ=255)
                                                    ccACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:2288537/69‑1 (MQ=255)
                                                    ccACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:3427098/69‑1 (MQ=255)
                                                    ccACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:1414934/69‑1 (MQ=255)
                                                    ccACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:3075073/69‑1 (MQ=255)
                                                    ccACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:641823/69‑1 (MQ=255)
                                                    ccACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:2867983/69‑1 (MQ=255)
                                                     cACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCgg           <  1:1407817/68‑1 (MQ=255)
                                                          tttATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacag      >  1:2706473/1‑68 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCTGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:1500030/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCAt                 >  1:2513667/1‑55 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   <  1:526421/69‑1 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:3195018/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:405417/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:1640567/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:694784/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   <  1:2142161/69‑1 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:2407271/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   <  1:777441/69‑1 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   <  1:112458/69‑1 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   >  1:2344131/1‑69 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGgacagaca   <  1:3402707/69‑1 (MQ=255)
                                                            tATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGACCCATTGCGGGATGCGCATCAGCGGgacagaca   <  1:165600/69‑1 (MQ=255)
                                                            tATGGCACTACTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGGACAGACaa  >  1:2508102/1‑69 (MQ=255)
                                                                   |                                                              
AATTGGACGCCACTGGCAGCGTCGGGCGCTCCAGCCTGTACCACCTGATGAGCCACTGTTTATGGCACTACTTTTTTCCATGTTGCGGGTTCCCGATCCATTGCGGGATGCGCATCAGCGGGACAGACAA  >  minE/1698690‑1698819

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: