Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 119,269 G→A 39.1% A48T (GCA→ACA) ‡ yacL → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE119,2690GA39.1% 38.0 / 47.2 74A48T (GCA→ACA) ‡yacLconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/14);  new base A (15/14);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAA  >  minE/119203‑119337
                                                                  |                                                                    
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:377636/69‑1 (MQ=255)
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:28425/69‑1 (MQ=255)
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:2338726/69‑1 (MQ=255)
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:1113669/69‑1 (MQ=255)
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:128945/69‑1 (MQ=255)
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:1285103/69‑1 (MQ=255)
cACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgca                                                                    <  1:1201085/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:311925/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:1703337/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:2991846/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:2241304/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:3092621/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:504103/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:745314/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:883358/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:2843058/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAAGACGCAGCGCGCACGCTCAAAgg                                                            <  1:1991038/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAgg                                                            <  1:967701/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAgg                                                            <  1:3363958/69‑1 (MQ=255)
        tAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGAGCACGCACTGAAAgg                                                            <  1:1836801/69‑1 (MQ=255)
          aTGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGt                                                           <  1:1882498/68‑1 (MQ=255)
          aTGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGt                                                           <  1:2655601/68‑1 (MQ=255)
          aTGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGt                                                           <  1:1550179/68‑1 (MQ=255)
          aTGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGt                                                           <  1:2092501/68‑1 (MQ=255)
          aTGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGt                                                           <  1:2426418/68‑1 (MQ=255)
          aTGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGt                                                           <  1:2650255/68‑1 (MQ=255)
             aaGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGcgcacgc                                                                     >  1:287461/1‑55 (MQ=255)
             aaGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCg                                                       >  1:982881/1‑69 (MQ=255)
             aaGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCg                                                       <  1:1423017/69‑1 (MQ=255)
             aaGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCg                                                       >  1:2169544/1‑69 (MQ=255)
                   tGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGt                                                 <  1:2207585/69‑1 (MQ=255)
                   tGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGt                                                 <  1:1318576/69‑1 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:1008155/1‑69 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:1980752/1‑69 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:1102441/1‑69 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:2920728/1‑69 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:1294101/1‑69 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:1765534/1‑69 (MQ=255)
                    gaaagaaaACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTc                                                >  1:1440057/1‑69 (MQ=255)
                                         tgatgaAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATg                            <  1:2819121/68‑1 (MQ=255)
                                         tgatgaAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATg                            <  1:2254637/68‑1 (MQ=255)
                                             gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa                        <  1:2544201/68‑1 (MQ=255)
                                             gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa                        <  1:1900362/68‑1 (MQ=255)
                                             gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa                        <  1:309013/68‑1 (MQ=255)
                                             gaAGTGGAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATa                        <  1:1570894/68‑1 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATacactc                   >  1:980139/1‑66 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:2696449/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:2742559/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:2976015/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:3065903/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:492682/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:947099/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:1013215/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:1465438/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:2074315/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:2017031/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:1913342/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:1131717/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:1673265/1‑69 (MQ=255)
                                                  ggAAGACGCAGCGCGCACGCTCAAAGGCAGCGAACGATCCTGGCAGCGGGCAGGGCATGAATACActct                  >  1:1087898/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:2619783/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:126952/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:363438/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:1482998/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:3155806/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:1636572/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:1668587/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:1789710/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:104995/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:1891720/1‑69 (MQ=255)
                                                     aCAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCtgtg               >  1:2180959/1‑69 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:3240947/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:761141/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:542744/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:1139403/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:463874/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:2660614/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:2813884/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:1060095/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:1902514/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:1872224/1‑57 (MQ=255)
                                                        agcagcGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATa                        >  1:1759761/1‑57 (MQ=255)
                                                                acgcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTg    <  1:3229804/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:1041267/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:1134053/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:724980/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:3184376/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:1691232/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:3051024/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:1881490/69‑1 (MQ=255)
                                                                  gcacTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGaa  <  1:2299148/69‑1 (MQ=255)
                                                                  |                                                                    
CACTGGTTTAATGAAGAGGTGAAAGAAAACCTGGCCTTGCTTGATGAAGTGGAACAAGCAGCGCACGCACTGAAAGGTAGCGAACGGTCCTGGCAACGGGCAGGGCATGAATACACCCTGTGGATGGACGGTGAA  >  minE/119203‑119337

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: