Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,617,691 G→T 32.0% S94S (TCC→TCA yhhQ ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,617,6910GT32.0% 121.8 / 63.4 100S94S (TCC→TCAyhhQconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/45);  new base T (17/15);  total (40/60)
Fisher's exact test for biased strand distribution p-value = 8.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 103 total aligned reads displayed.

TGAAACTGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGA  >  minE/2617624‑2617758
                                                                   |                                                                   
tGAAACTGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAgg                                                                    <  1:2488903/69‑1 (MQ=255)
      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:2950630/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:2981795/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:2455131/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:376334/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:2223850/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:2060895/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:1953119/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:585023/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:3208974/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:1506173/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:1209263/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:1353402/62‑2 (MQ=255)
                  tACGGGCGACAAACAGGTTGAAGTGGGCGAGCGCGCCGAATCCTTGCCATGACCCCAtataa                                                         <  1:135644/62‑2 (MQ=255)
                       gctaCAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:300215/66‑1 (MQ=255)
                       gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:642899/69‑1 (MQ=255)
                       gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:35063/69‑1 (MQ=255)
                       gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:2367590/69‑1 (MQ=255)
                       gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:33337/69‑1 (MQ=255)
                       gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:2905312/69‑1 (MQ=255)
                       gCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:910197/69‑1 (MQ=255)
                        cGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:536454/68‑1 (MQ=255)
                        cGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACgag                                             <  1:1229325/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:58399/49‑1 (MQ=255)
                                     gAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGc                                                   <  1:703543/49‑1 (MQ=255)
                                     gAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGc                                                   <  1:684833/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:2974658/60‑1 (MQ=255)
                                           ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                  <  1:3090544/60‑1 (MQ=255)
                                           ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                  <  1:3376493/60‑1 (MQ=255)
                                           ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                  <  1:171409/60‑1 (MQ=255)
                                           ggCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGga                                  <  1:450812/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:3072595/1‑67 (MQ=255)
                                                 tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                     >  1:612177/1‑67 (MQ=255)
                                                 tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                     >  1:1992714/1‑67 (MQ=255)
                                                 tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                     >  1:617806/1‑67 (MQ=255)
                                                 tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                     >  1:2594126/1‑67 (MQ=255)
                                                 tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                     >  1:2231138/1‑67 (MQ=255)
                                                 tGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCa                     >  1:1474799/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:1598137/68‑1 (MQ=255)
                                                     ccGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGAt                <  1:2829793/68‑1 (MQ=255)
                                                     ccGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGAt                <  1:1236117/68‑1 (MQ=255)
                                                        aaTCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAtt            <  1:2083190/69‑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:858748/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:678474/1‑65 (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:1072707/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:1110254/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:1128437/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:1146064/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:2705635/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:1190698/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:2520323/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:1891244/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:190258/1‑65 (MQ=255)
                                                          tCCTTGCCATGACCCCATATAAAACAGCGATGAGATGACGTAGGAGATCAGCAATGCGGGGATCa              >  1:3206630/1‑65 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACgaga                                            >  1:1457972/1‑35 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTa           >  1:2442406/1‑68 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc          >  1:2169768/1‑69 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc          >  1:1629162/1‑69 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc          >  1:165292/1‑69 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc          >  1:3316298/1‑69 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc          >  1:2840864/1‑69 (MQ=255)
                                                          tCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTAc          >  1:3204965/1‑69 (MQ=255)
                                                            ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt             >  1:533646/1‑64 (MQ=255)
                                                            ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt             >  1:2421347/1‑64 (MQ=255)
                                                            ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt             >  1:356670/1‑64 (MQ=255)
                                                            ccTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCAt             >  1:702083/1‑64 (MQ=255)
                                                                  aGGAACCCATATAGAATATCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa  <  1:614695/69‑1 (MQ=255)
                                                                  aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATCACGCAGGGATCATTACCGCGAAGa  <  1:3235735/69‑1 (MQ=255)
                                                                  aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa  <  1:2431780/69‑1 (MQ=255)
                                                                  aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa  <  1:628811/69‑1 (MQ=255)
                                                                  aGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGa  <  1:1562420/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)
                                                                   |                                                                   
TGAAACTGGCGGTGGCGATACGGGCGACAAACAGGTTGAAGTGGGCGAGTGCGCCGAATCCCTGCCAGGAACCCATATAGAATAGCGACGAGATGACGTAGGAGATTAATAACGCAGGGATCATTACCGCGAAGA  >  minE/2617624‑2617758

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: