Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,221,439 C→T 37.4% N422N (AAC→AAT
T5M (ACG→ATG) 
yebS →
yebT →
conserved inner membrane protein
conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,221,4390CT37.4% 43.0 / 42.5 72N422N (AAC→AAT
T5M (ACG→ATG) 
yebS
yebT
conserved inner membrane protein
conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base T (4/23);  total (24/48)
Fisher's exact test for biased strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

ATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCTTT  >  minE/1221371‑1221493
                                                                    |                                                      
aTTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAAc                                                        >  1:3394272/1‑69 (MQ=255)
aTTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAAc                                                        >  1:3108986/1‑69 (MQ=255)
aTTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAAc                                                        >  1:300017/1‑69 (MQ=255)
   ttGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGcc                                                     >  1:3376413/1‑69 (MQ=255)
   ttGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGcc                                                     >  1:2610646/1‑69 (MQ=255)
                       aTGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAgcg                                  <  1:1007822/68‑1 (MQ=255)
                       aTGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAgcg                                  <  1:124066/68‑1 (MQ=255)
                       aTGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAgcg                                  <  1:1396912/68‑1 (MQ=255)
                       aTGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGGCAGGAAACGCCCGCTTCGACGACTGAAgcg                                  <  1:2881540/68‑1 (MQ=255)
                        tGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCa                                <  1:778710/69‑1 (MQ=255)
                        tGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCa                                <  1:3223806/69‑1 (MQ=255)
                        tGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCa                                <  1:2150746/69‑1 (MQ=255)
                        tGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCa                                <  1:2958387/69‑1 (MQ=255)
                        tGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCa                                <  1:1685404/69‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:2689536/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:2252267/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:3239290/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:1756180/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:183389/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:1410448/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:598245/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:150412/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:515479/68‑1 (MQ=255)
                                cAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTaaa                         <  1:151445/68‑1 (MQ=255)
                                   ccGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATa                     <  1:1199363/69‑1 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:1717181/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:803362/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:2752282/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:252400/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:3129555/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:3166115/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:2263602/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:400211/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATaa                    >  1:396295/1‑69 (MQ=255)
                                    cGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATa                     <  1:2673349/68‑1 (MQ=255)
                                    cGCTTACTTAGGGATGCACATGAGTCAGCAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATa                     >  1:1653123/1‑68 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGATATGACCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2273487/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:3209293/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:3125429/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:3256524/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:3262128/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2539387/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:392604/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:800982/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:37622/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2944655/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:989177/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2871874/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1520924/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1560359/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1605075/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1198839/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1798636/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1909823/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2209810/61‑1 (MQ=255)
                                        tACTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2375064/61‑1 (MQ=255)
                                        tACTATGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:307113/61‑1 (MQ=255)
                                         aCTTTGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:1391614/60‑1 (MQ=255)
                                             tGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACTACTGAAGCGCAGATTAAAAATAAACGCCGTAt           >  1:1583560/1‑69 (MQ=255)
                                             tGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTAt           >  1:1961068/1‑69 (MQ=255)
                                             tGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTAt           >  1:1141445/1‑69 (MQ=255)
                                             tGGGATGCACATGAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTAt           >  1:2467538/1‑69 (MQ=255)
                                                    cacaTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCAcc     >  1:3379367/1‑68 (MQ=255)
                                                    cacaTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCAcc     >  1:2602093/1‑68 (MQ=255)
                                                    cacaTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCt    >  1:2316056/1‑69 (MQ=255)
                                                    cacaTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCt    >  1:2060204/1‑69 (MQ=255)
                                                    cacaTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCt    >  1:1457815/1‑69 (MQ=255)
                                                      caTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCttt  >  1:2343638/1‑69 (MQ=255)
                                                      caTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCttt  >  1:1701221/1‑69 (MQ=255)
                                                      caTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCttt  >  1:814744/1‑69 (MQ=255)
                                                      caTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCtt   <  1:2764428/68‑1 (MQ=255)
                                                      caTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCtt   <  1:1676902/68‑1 (MQ=255)
                                                       aTGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCtt   <  1:3376103/67‑1 (MQ=255)
                                                         gAGTCAGGAAATGCCCGCTTCGACGACTGAAGCGCAGATTaaaa                        <  1:2570841/44‑1 (MQ=255)
                                                           gTCAGCAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCtt   <  1:1106761/63‑1 (MQ=255)
                                                                    |                                                      
ATTTTGACTATTCTTGCTGTGGAATGGCTGGACAGCCGCTTACTTTGGGATGCACATGAGTCAGGAAACGCCCGCTTCGACGACTGAAGCGCAGATTAAAAATAAACGCCGTATCTCACCTTT  >  minE/1221371‑1221493

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: