Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,617,694 A→C 31.7% G93G (GGT→GGG yhhQ ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,617,6940AC31.7% 84.5 / 74.7 101G93G (GGT→GGGyhhQconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (23/46);  new base C (17/15);  total (40/61)
Fisher's exact test for biased strand distribution p-value = 8.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

Only 100 of 103 total aligned reads displayed.

TGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGATAA  >  minE/2617630‑2617761
                                                                |                                                                   
tGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCa                                                                 <  1:1069378/69‑1 (MQ=255)
tGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCa                                                                 <  1:3295755/69‑1 (MQ=255)
tGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCa                                                                 <  1:3241249/69‑1 (MQ=255)
    ggTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCata                                                               <  1:232170/67‑1 (MQ=255)
     gTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAg                                                            <  1:1020359/69‑1 (MQ=255)
     gTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAg                                                            <  1:1676501/69‑1 (MQ=255)
     gTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAg                                                            <  1:3219115/69‑1 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:2758470/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:2455131/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:2950630/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:2223850/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:376334/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:2060895/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:1953119/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:3208974/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:585023/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:2981795/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:1506173/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:135644/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:1353402/62‑2 (MQ=255)
            tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                            <  1:1209263/62‑2 (MQ=255)
                 gctaCAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:300215/66‑1 (MQ=255)
                 gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:642899/69‑1 (MQ=255)
                 gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:2905312/69‑1 (MQ=255)
                 gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:33337/69‑1 (MQ=255)
                 gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:2367590/69‑1 (MQ=255)
                 gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:35063/69‑1 (MQ=255)
                  cGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:1229325/68‑1 (MQ=255)
                  cGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                                <  1:536454/68‑1 (MQ=255)
                            gTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:1143512/69‑1 (MQ=255)
                            gTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:1664173/69‑1 (MQ=255)
                            gTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:2249286/69‑1 (MQ=255)
                               gAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGc                                                      <  1:703543/49‑1 (MQ=255)
                               gAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGc                                                      <  1:684833/49‑1 (MQ=255)
                               gAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGc                                                      <  1:58399/49‑1 (MQ=255)
                               gAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGc                                                      <  1:2360554/49‑1 (MQ=255)
                                   tGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:2590432/62‑1 (MQ=255)
                                     ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGCATAGCGACGAGATGACGTAGga                                     <  1:3142348/60‑1 (MQ=255)
                                     ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:171409/60‑1 (MQ=255)
                                     ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:3376493/60‑1 (MQ=255)
                                     ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:450812/60‑1 (MQ=255)
                                     ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:3090544/60‑1 (MQ=255)
                                     ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:2974658/60‑1 (MQ=255)
                                     ggCGACTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:2453329/60‑1 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                     <  1:3153321/54‑1 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:82201/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:2594126/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:617806/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:3072595/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:612177/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:2231138/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:1474799/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                        >  1:1992714/1‑67 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAgg                      >  1:2576192/1‑69 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAgg                      >  1:1253656/1‑69 (MQ=255)
                                           tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAgg                      >  1:1979635/1‑69 (MQ=255)
                                               ccGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGAt                   <  1:2111666/68‑1 (MQ=255)
                                               ccGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGAt                   <  1:2829793/68‑1 (MQ=255)
                                               ccGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGAt                   <  1:1598137/68‑1 (MQ=255)
                                               ccGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGAt                   <  1:1236117/68‑1 (MQ=255)
                                                  aaTCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAtt               <  1:1255979/69‑1 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGGAATGCGGGGATCa                 >  1:1429148/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATc                  >  1:310434/1‑64 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:325766/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:3372276/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:10266/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:678474/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:858748/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:3206630/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:1072707/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:1110254/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:1128437/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:2705635/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:1146064/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:2520323/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:1190698/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:1891244/1‑65 (MQ=255)
                                                    tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa                 >  1:190258/1‑65 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACgaga                                               >  1:1457972/1‑35 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTa              >  1:2442406/1‑68 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc             >  1:2840864/1‑69 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc             >  1:1629162/1‑69 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc             >  1:165292/1‑69 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc             >  1:3204965/1‑69 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc             >  1:2169768/1‑69 (MQ=255)
                                                    tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc             >  1:3316298/1‑69 (MQ=255)
                                                      ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt                >  1:533646/1‑64 (MQ=255)
                                                      ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt                >  1:702083/1‑64 (MQ=255)
                                                      ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt                >  1:2421347/1‑64 (MQ=255)
                                                      ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt                >  1:356670/1‑64 (MQ=255)
                                                            aGGAACCCATATAGAATATCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:614695/69‑1 (MQ=255)
                                                            aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATCACGCAGGGATCATTACCGCGAAGa     <  1:3235735/69‑1 (MQ=255)
                                                            aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:1562420/69‑1 (MQ=255)
                                                            aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:628811/69‑1 (MQ=255)
                                                            aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:2431780/69‑1 (MQ=255)
                                                            aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:2576512/69‑1 (MQ=255)
                                                            aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:2789802/69‑1 (MQ=255)
                                                             ggAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa     <  1:2178731/68‑1 (MQ=255)
                                                               aaCCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGATaa  <  1:845478/69‑1 (MQ=255)
                                                               aaCCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGATaa  <  1:2025757/69‑1 (MQ=255)
                                                                |                                                                   
TGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGATAA  >  minE/2617630‑2617761

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: