Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 285,523 T→C 56.2% F126F (TTT→TTC clpP → proteolytic subunit of ClpA‑ClpP and ClpX‑ClpP ATP‑dependent serine proteases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE285,5230TC56.2% 19.1 / 54.4 73F126F (TTT→TTCclpPproteolytic subunit of ClpA‑ClpP and ClpX‑ClpP ATP‑dependent serine proteases
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base C (28/13);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 4.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCTCGATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGAT  >  minE/285471‑285589
                                                    |                                                                  
cggcCTCGATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTcgc                                                    >  1:2564828/1‑69 (MQ=255)
cggcCTCGATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTcgc                                                    >  1:2063067/1‑69 (MQ=255)
      cGATGGGTGCTTTCTTGCTGACTGCAGGGGCAAAAGGTAAACGTTTCtgcc                                                                >  1:794538/1‑51 (MQ=255)
      cGATGGGTGCTTTCTTGCTGACTGCAGGGGCAAAAGGTAAACGTTTCtgcc                                                                >  1:1919444/1‑51 (MQ=255)
      cGATGGGTGCTTTCTTGCTGACTGCAGGGGCAAAAGGTAAACGTTTCtgcc                                                                >  1:1334192/1‑51 (MQ=255)
      cGATGGGTGCTTTCTTGCTGACTGCAGGGGCAAAAGGTAAACGTTTCtgcc                                                                >  1:1787675/1‑51 (MQ=255)
      cGATGGGTGCTTTCTTGCTGACTGCAGGGGCAAAAGGTAAACGTTTCtgcc                                                                >  1:1460889/1‑51 (MQ=255)
      cGATGGGTGCTTTCTTGCTGACTGCAGGGGCAAAAGGTAAACGTTTCtgcc                                                                >  1:3303365/1‑51 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         <  1:2734328/69‑1 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         <  1:2646528/69‑1 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         <  1:1839263/69‑1 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         <  1:1685265/69‑1 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         >  1:923736/1‑69 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         <  1:1106756/69‑1 (MQ=255)
           ggCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGAtt                                         <  1:1282487/69‑1 (MQ=255)
                            cGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGggcg                         >  1:2149242/1‑68 (MQ=255)
                            cGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGggcg                         >  1:120605/1‑68 (MQ=255)
                            cGCAGGGGCAAAAGGTAAACGTTTTTGCCTGACGAATTCGCGGGTGATGATTCACCAACCGTTGggcg                         >  1:1861600/1‑68 (MQ=255)
                             gCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGggcggc                       <  1:228052/69‑1 (MQ=255)
                              cAGGGTCAAAAGGTAAACGTTTCTGCCTGCCGAAtt                                                       >  1:236979/1‑36 (MQ=37)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2933870/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:1676276/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:842085/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2410691/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2825842/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2824300/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2814415/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2666778/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2679605/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2742283/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:3435069/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:3373894/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:1737025/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:3122085/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:3063064/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:1388422/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:1808890/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2001051/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2051032/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGt                                                >  1:2055122/1‑43 (MQ=255)
                              cAGGGGCAAAAGGTAAACGTTTCTGCCTGACGAAATc                                                      >  1:488502/1‑37 (MQ=25)
                                ggggCAAAAGGTAAACGTTTCTGCCTGGCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:2321910/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTTATGATTCACCAGCCGTTGGGCGGCt                      <  1:3018764/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:2866279/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:2888374/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:3094659/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:479302/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:1039705/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:2575272/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:1113582/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:2075028/67‑1 (MQ=255)
                                ggggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:1782464/67‑1 (MQ=255)
                                 gggCAAAAGGTAAACGTTTCTGCCTGCCGAATTCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:1767711/66‑1 (MQ=255)
                                   gCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:2523617/66‑1 (MQ=255)
                                   gCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:2399699/66‑1 (MQ=255)
                                   gCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:2063521/66‑1 (MQ=255)
                                   gCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:1825511/66‑1 (MQ=255)
                                   gCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:1757764/66‑1 (MQ=255)
                                   gCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:1594534/66‑1 (MQ=255)
                                    cAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTAc                    <  1:2688866/65‑1 (MQ=255)
                                         ggTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:711613/69‑1 (MQ=255)
                                         ggTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:1175955/69‑1 (MQ=255)
                                         ggTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:1532528/69‑1 (MQ=255)
                                         ggTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:714777/69‑1 (MQ=255)
                                         ggTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:836937/69‑1 (MQ=255)
                                           tAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAgg          >  1:2596654/1‑68 (MQ=255)
                                           tAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGc         >  1:1097105/1‑69 (MQ=255)
                                           tAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGc         >  1:528034/1‑69 (MQ=255)
                                           tAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGc         >  1:3373435/1‑69 (MQ=255)
                                           tAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGc         >  1:1561348/1‑69 (MQ=255)
                                           tAAACGTTTCTGCCTGCCGAAATCGCGCGTGATGATTCACCAGCCGTTGGGCGGCt                      <  1:950006/56‑1 (MQ=255)
                                            aaaCGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:2951801/66‑1 (MQ=255)
                                               cGTTTTTGCCTGCCGATTTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAg           <  1:2319577/63‑1 (MQ=255)
                                                  ttttGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGat  <  1:2232957/69‑1 (MQ=255)
                                                    |                                                                  
CGGCCTCGATGGGCGCTTTCTTGCTGACCGCAGGGGCAAAAGGTAAACGTTTTTGCCTGCCGAATTCGCGCGTGATGATTCACCAACCGTTGGGCGGCTACCAGGGCCAGGCGACCGAT  >  minE/285471‑285589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: