Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,318,439 T→A 18.4% G47G (GGA→GGT cld ← regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,318,4390TA18.4% 222.9 / 36.4 113G47G (GGA→GGTcldregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base T (59/33);  new base A (13/8);  total (72/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTG  >  minE/1318374‑1318504
                                                                 |                                                                 
tGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1095078/1‑68 (MQ=255)
tGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:2042219/1‑68 (MQ=255)
tGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:637739/1‑68 (MQ=255)
tGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1421625/1‑68 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAAt                                                              >  1:2124638/1‑66 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:501967/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:516915/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:2928556/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:3284642/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:2420501/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:2930475/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:1425524/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:2720911/1‑67 (MQ=255)
     cACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATa                                                             >  1:1537715/1‑67 (MQ=255)
           gggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGc                                                     <  1:2311544/69‑1 (MQ=255)
           gggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGc                                                     <  1:2021327/69‑1 (MQ=255)
           gggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATAGCCAGGGc                                                     <  1:216774/69‑1 (MQ=255)
           gggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATAGCCAGGGc                                                     <  1:1667328/69‑1 (MQ=255)
           gggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATAGCCAGGGc                                                     <  1:750310/69‑1 (MQ=255)
           gggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATAGCCAGGGc                                                     <  1:2921220/69‑1 (MQ=255)
            ggCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATAGCCAGGGc                                                     <  1:3171003/68‑1 (MQ=255)
              cTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCaa                                                   >  1:1143262/1‑68 (MQ=255)
              cTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAAt                                                  >  1:3319307/1‑69 (MQ=255)
                    gATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCcac                                            >  1:706369/1‑69 (MQ=255)
                    gATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCcac                                            >  1:2196885/1‑69 (MQ=255)
                    gATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCcac                                            >  1:136167/1‑69 (MQ=255)
                        aTTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1371197/1‑44 (MQ=255)
                        aTTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1742315/1‑44 (MQ=255)
                        aTTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:3163530/1‑44 (MQ=255)
                        aTTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:3136980/1‑44 (MQ=255)
                         ttGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATAACCAATAGCCAGGGc                                                     <  1:3242806/55‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1423162/1‑35 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:2287737/1‑35 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1622616/1‑35 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:1483251/1‑35 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATcc                                                                 >  1:452288/1‑35 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:1882637/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:2507851/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:2790729/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:3156994/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:1353378/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:2027178/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:141058/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:1323383/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:1810686/68‑1 (MQ=255)
                                 gACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAat                                <  1:1687239/68‑1 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGc                                                           <  1:3235433/35‑1 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:2361361/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:2795290/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:1600974/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:3179774/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:1634001/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:372978/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   <  1:3251524/59‑1 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   <  1:1709039/59‑1 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   <  1:1184234/59‑1 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:862/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:332288/1‑59 (MQ=255)
                                       cATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGa                                   >  1:3008909/1‑59 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGAtt                         >  1:1874532/1‑67 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:815073/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:667002/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:891189/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:79931/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:3250242/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:624083/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:1132932/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:434206/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:3376343/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:2835672/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:1159970/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:1439938/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:1489190/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:1639720/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:1811447/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:2215033/1‑69 (MQ=255)
                                         tttCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGt                       >  1:288889/1‑69 (MQ=255)
                                         tttCACCTTCGCTACTGCCAAATATCCAATAGCCAGGGCTATAGCCACAATGACGGAAATGAAGATTGt                       >  1:1161249/1‑69 (MQ=255)
                                           tctcCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCa                     >  1:795531/1‑69 (MQ=255)
                                           tctcCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCa                     >  1:2758137/1‑69 (MQ=255)
                                           tctcCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCa                     >  1:3417228/1‑69 (MQ=255)
                                           tctcCTTCGCTAATGCCAAATATCCAATAGCCAGGGCAATAGc                                               >  1:1879474/1‑43 (MQ=255)
                                                 tCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgc               <  1:2500146/69‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:1716973/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:2336848/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:2328247/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:1305057/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:3314890/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:1549091/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:261363/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:2105449/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:2555767/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:1972856/68‑1 (MQ=255)
                                                    cTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgcca             <  1:2101489/68‑1 (MQ=255)
                                                         gCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTTTCATCTTgccacgccac       <  1:1996949/69‑1 (MQ=255)
                                                         gCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgccac       <  1:2484937/69‑1 (MQ=255)
                                                         gCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgccac       <  1:2132644/69‑1 (MQ=255)
                                                              atatCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTg  <  1:1758329/69‑1 (MQ=255)
                                                                 tCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTg  >  1:2814922/1‑66 (MQ=255)
                                                                 tCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTg  >  1:349761/1‑66 (MQ=255)
                                                                 |                                                                 
TGCCCCACGTCGGGCTGAGTGATAATTGCTGTTGACGTCCATTTCTCCTTCGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTG  >  minE/1318374‑1318504

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: