Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,103,846 T→C 22.7% Y5H (TAC→CAC)  ompR → DNA‑binding response regulator in two‑component regulatory system with EnvZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,103,8460TC22.7% 145.5 / 43.3 75Y5H (TAC→CAC) ompRDNA‑binding response regulator in two‑component regulatory system with EnvZ
Reads supporting (aligned to +/- strand):  ref base T (36/22);  new base C (9/8);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACAAG  >  W3110S.gb/4103778‑4103915
                                                                    |                                                                     
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTa                                                                      >  1:1204097/1‑70 (MQ=255)
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTa                                                                      >  1:3493733/1‑70 (MQ=255)
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTa                                                                      >  1:355330/1‑70 (MQ=255)
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTa                                                                      >  1:897473/1‑70 (MQ=255)
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTa                                                                      >  1:3782704/1‑70 (MQ=255)
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTa                                                                      >  1:3884705/1‑70 (MQ=255)
     ttttAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGa                                                                 <  1:1430120/70‑1 (MQ=255)
     ttttAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGAt                                                                >  1:921706/1‑71 (MQ=255)
        tAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAAAGCAAGAGAACTACAAGATTc                                                              <  1:1787134/70‑1 (MQ=255)
                cacGCTTACAAATTGTTGGGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:4077879/1‑70 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCtg                                                            >  1:3552604/1‑64 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTggtggt                                                       >  1:123544/1‑69 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:2353049/1‑70 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      <  1:3548081/70‑1 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:82094/1‑70 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:1013224/1‑70 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:1858048/1‑70 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      <  1:142689/70‑1 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:1612851/1‑70 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                      >  1:907735/1‑70 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGAt                                                   <  1:2162238/71‑1 (MQ=255)
                    cTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTc                                                      <  1:1281255/66‑1 (MQ=255)
                      tACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgac                                                >  1:2270468/1‑70 (MQ=255)
                      tACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGa                                                    >  1:3698/1‑66 (MQ=255)
                      tACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGa                                                    >  1:3054623/1‑66 (MQ=255)
                      tACAAATTGTTCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGa                                                    >  1:3026582/1‑65 (MQ=255)
                       aCAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgac                                                <  1:72562/69‑1 (MQ=255)
                          tcttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                             <  1:3088741/67‑1 (MQ=255)
                          aaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACAt                                           >  1:3470809/1‑71 (MQ=255)
                          aaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACAt                                           >  1:1485019/1‑71 (MQ=255)
                              gttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                             <  1:457413/65‑1 (MQ=255)
                                  cGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtg                                    <  1:451899/70‑1 (MQ=255)
                                  cGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCtg                                    <  1:1341848/70‑1 (MQ=255)
                                   gAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcg                                  <  1:3783376/71‑1 (MQ=255)
                                   gAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcg                                  <  1:167626/71‑1 (MQ=255)
                                   gAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcg                                  <  1:2424997/71‑1 (MQ=255)
                                      ccTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTc                                                              <  1:645225/40‑1 (MQ=255)
                                      ccTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCg                                                     <  1:2883406/49‑1 (MQ=255)
                                      ccTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCg                                                     <  1:467261/49‑1 (MQ=255)
                                         ttGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACAt                                           >  1:2238913/1‑56 (MQ=255)
                                         ttGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACAt                                           >  1:2989936/1‑56 (MQ=255)
                                         ttGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcg                                  >  1:3892082/1‑65 (MQ=255)
                                              aGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTgg                       <  1:3623795/71‑1 (MQ=255)
                                               gTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgaa                                             <  1:697942/48‑2 (MQ=255)
                                               gTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTgg                       >  1:3379137/1‑70 (MQ=255)
                                                tACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                             <  1:3749935/47‑1 (MQ=255)
                                                tACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                             <  1:754345/47‑1 (MQ=255)
                                                tACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACAt                                           <  1:3741231/49‑1 (MQ=255)
                                                tACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACAt                                           <  1:1926093/49‑1 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTggtggt                                                       >  1:3156640/1‑35 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATgac                                                >  1:3510133/1‑42 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACg                   >  1:1071540/1‑71 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACg                   >  1:2069205/1‑71 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACg                   >  1:2495389/1‑71 (MQ=255)
                                                  caaacaaTGCAAGACAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACg                   >  1:2101537/1‑71 (MQ=255)
                                                   aaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCgctg                           <  1:3532998/62‑1 (MQ=255)
                                                   aaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCgctg                           <  1:3044147/62‑1 (MQ=255)
                                                      caaTGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTa                >  1:3277977/1‑70 (MQ=255)
                                                      caaTGCAAGAGAACCACAAGATTCTGGTGGTCCATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTa                >  1:2292174/1‑70 (MQ=255)
                                                       aaTGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTAtc              >  1:124136/1‑71 (MQ=255)
                                                       aaTGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTAtc              >  1:2014590/1‑71 (MQ=255)
                                                           cAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCAc          >  1:2657112/1‑71 (MQ=255)
                                                           cAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCAc          >  1:3927170/1‑71 (MQ=255)
                                                           cAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCAc          >  1:2684501/1‑71 (MQ=255)
                                                              gagaACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCg        >  1:3481481/1‑70 (MQ=255)
                                                               agaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGaa      >  1:1319950/1‑71 (MQ=255)
                                                                gaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAAc     <  1:2078177/71‑1 (MQ=255)
                                                                gaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAAc     <  1:1368596/71‑1 (MQ=255)
                                                                 aaCCACAAGATTCTGGTGGTCTATGACGACATGCGCCTGCTTGCGCTGCTGGAACGTTATCTCACCGAACa    <  1:2395806/71‑1 (MQ=255)
                                                                 aaCCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACa    <  1:2726200/71‑1 (MQ=255)
                                                                 aaCCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACa    <  1:3426862/71‑1 (MQ=255)
                                                                   cTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACAAg  >  1:2574206/1‑71 (MQ=255)
                                                                   cTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACAAg  >  1:2426192/1‑71 (MQ=255)
                                                                   cTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACAAg  >  1:3047949/1‑71 (MQ=255)
                                                                   cTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACAAg  >  1:3507511/1‑71 (MQ=255)
                                                                    tACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTAt               <  1:2049023/57‑1 (MQ=255)
                                                                    |                                                                     
AGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACAAG  >  W3110S.gb/4103778‑4103915

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: