Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,103,846 T→C 30.3% 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,8460TC30.3% 103.1 / 64.6 76Y5H (TAC→CAC) ompRDNA‑binding response regulator in two‑component regulatory system with EnvZ
Reads supporting (aligned to +/- strand):  ref base T (31/22);  new base C (15/8);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACA  >  W3110S.gb/4103778‑4103913
                                                                    |                                                                   
aGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACt                                                                     <  1:163228/69‑1 (MQ=255)
       ttAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGAt                                                              >  1:2676024/1‑69 (MQ=255)
       ttAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGAt                                                              >  1:3165773/1‑69 (MQ=255)
        tAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGAt                                                              <  1:617886/68‑1 (MQ=255)
            aaTACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGa                                                               <  1:2386615/63‑1 (MQ=255)
            aaTACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCtg                                                          <  1:846081/68‑1 (MQ=255)
              tACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCtggtg                                                       >  1:710344/1‑69 (MQ=255)
              tACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCtggtg                                                       >  1:2317006/1‑69 (MQ=255)
              tACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCtggtg                                                       >  1:271940/1‑69 (MQ=255)
              tACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCtggtg                                                       >  1:598699/1‑69 (MQ=255)
              tACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCtggtg                                                       >  1:403129/1‑69 (MQ=255)
              tACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCtggt                                                        >  1:1054750/1‑68 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTc                                                            >  1:3646640/1‑62 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCtggtg                                                       >  1:414461/1‑67 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCtggtg                                                       >  1:3550139/1‑67 (MQ=255)
                cacGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCtggtgg                                                      >  1:2345046/1‑68 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                    <  1:192652/68‑1 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                    <  1:1586654/68‑1 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                    <  1:619445/68‑1 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTc                                                    <  1:2645028/68‑1 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTc                                                    <  1:3097517/68‑1 (MQ=255)
                  cGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTc                                                    <  1:603595/68‑1 (MQ=255)
                   gCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTc                                                    <  1:2534332/67‑1 (MQ=255)
                        cAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgac                                              >  1:1764235/1‑68 (MQ=255)
                         aaaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgac                                              <  1:3691467/67‑1 (MQ=255)
                         aaaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgac                                              <  1:3036390/67‑1 (MQ=255)
                          aaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            <  1:438446/68‑1 (MQ=255)
                          aaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            <  1:1572644/68‑1 (MQ=255)
                          aaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATgacga                                            <  1:2908401/68‑1 (MQ=255)
                          aaTTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATgacga                                            <  1:103249/68‑1 (MQ=255)
                            ttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCtgg                                                         >  1:287606/1‑53 (MQ=255)
                            ttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            >  1:1061953/1‑66 (MQ=255)
                            ttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            >  1:663922/1‑66 (MQ=255)
                            ttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            >  1:1100117/1‑66 (MQ=255)
                            ttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            >  1:2552281/1‑66 (MQ=255)
                            ttgttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATgacga                                            >  1:1709023/1‑66 (MQ=255)
                              gttgCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATgc                                       >  1:735194/1‑69 (MQ=255)
                                  cGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGa                                                               <  1:740758/41‑1 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCg                                                   >  1:3092997/1‑50 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCg                                                   >  1:630631/1‑50 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCg                                                   >  1:1349399/1‑50 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCg                                                   >  1:2963555/1‑50 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCg                                                   >  1:2924411/1‑50 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcg                                <  1:3690433/69‑1 (MQ=255)
                                     aCCTTTGGGAGTACAAACAATGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcg                                <  1:912826/69‑1 (MQ=255)
                                        tttGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcgtgc                             >  1:2583578/1‑69 (MQ=255)
                                        tttGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCtgcgtgc                             >  1:2466581/1‑69 (MQ=255)
                                                 acaaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCg                                                   <  1:2702991/38‑1 (MQ=255)
                                                 acaaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCg                                                   <  1:2020308/38‑1 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTc                                                    >  1:1196887/1‑36 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTc                                                    >  1:2414050/1‑36 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTc                                                    >  1:2376274/1‑36 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                           <  1:12075/45‑1 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                           <  1:1248125/45‑1 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATgacgac                                           <  1:1559358/45‑1 (MQ=255)
                                                  caaacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGAGGACGACATGCGCCTGCGTGCGCTGCTGGaa                   >  1:2325689/1‑69 (MQ=255)
                                                  caaacaaTGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGaa                   >  1:1202672/1‑69 (MQ=255)
                                                  caaacaaTGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGaa                   >  1:3324104/1‑69 (MQ=255)
                                                  caaacaaTGCAAGAGAACCACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGaa                   >  1:857086/1‑69 (MQ=255)
                                                    aacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCgct                          >  1:2915868/1‑60 (MQ=255)
                                                    aacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCgct                          >  1:2126381/1‑60 (MQ=255)
                                                    aacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCgct                          >  1:3151839/1‑60 (MQ=255)
                                                    aacaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAAc                  >  1:3465860/1‑68 (MQ=255)
                                                     acaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTgcgc                           >  1:3709745/1‑58 (MQ=255)
                                                     acaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTgcgc                           >  1:1180300/1‑58 (MQ=255)
                                                     acaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGt                <  1:401592/69‑1 (MQ=255)
                                                     acaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGt                >  1:2455705/1‑69 (MQ=255)
                                                     acaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGt                <  1:3390938/69‑1 (MQ=255)
                                                     acaaTGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGt                >  1:3336717/1‑69 (MQ=255)
                                                           cAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTAtctc          >  1:1495251/1‑69 (MQ=255)
                                                           cAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTAtctc          <  1:600070/69‑1 (MQ=255)
                                                           cAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTAtctc          <  1:120799/69‑1 (MQ=255)
                                                              gagaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTgg                     <  1:1236527/55‑1 (MQ=255)
                                                              gagaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTgg                     <  1:1232007/55‑1 (MQ=255)
                                                                gaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCg      >  1:2709279/1‑68 (MQ=255)
                                                                gaACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCg      >  1:3491993/1‑68 (MQ=255)
                                                                   cTACAAGAATCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACa  >  1:1075370/1‑69 (MQ=255)
                                                                    tACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACa  <  1:1339800/68‑1 (MQ=255)
                                                                    |                                                                   
AGCTTTTTTAAGAATACACGCTTACAAATTGTTGCGAACCTTTGGGAGTACAAACAATGCAAGAGAACTACAAGATTCTGGTGGTCGATGACGACATGCGCCTGCGTGCGCTGCTGGAACGTTATCTCACCGAACA  >  W3110S.gb/4103778‑4103913

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: