Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 76,662 G→A 45.7% P170P (CCG→CCA ftsI → transpeptidase involved in septal peptidoglycan synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE76,6620GA45.7% 17.6 / 76.8 94P170P (CCG→CCAftsItranspeptidase involved in septal peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base G (15/36);  new base A (22/21);  total (37/57)
Fisher's exact test for biased strand distribution p-value = 3.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGG  >  minE/76600‑76724
                                                              |                                                              
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:3095836/1‑68 (MQ=255)
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:1631488/1‑68 (MQ=255)
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:553399/1‑68 (MQ=255)
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:535651/1‑68 (MQ=255)
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:1283810/1‑68 (MQ=255)
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:1002303/1‑68 (MQ=255)
tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                                           >  1:3362911/1‑68 (MQ=255)
     cAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGt                                                     >  1:2673180/1‑69 (MQ=255)
        aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                                  <  1:1015712/69‑1 (MQ=255)
        aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                                  <  1:1196489/69‑1 (MQ=255)
        aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                                  <  1:934593/69‑1 (MQ=255)
         aaactgaaactgCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                                  <  1:1680330/68‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTTAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:780613/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:2641038/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:1711784/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:428820/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:1389752/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:1308214/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:1819949/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:1831536/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:2013892/69‑1 (MQ=255)
                 actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:918616/69‑1 (MQ=255)
                  ctgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:3015986/68‑1 (MQ=255)
                  ctgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                         <  1:1107494/68‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:2937030/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:2381278/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:2120674/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:1170104/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:2608426/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:844886/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:384056/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:2825082/69‑1 (MQ=255)
                         ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                 <  1:2868887/69‑1 (MQ=255)
                          gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                                <  1:448353/69‑1 (MQ=255)
                          gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                                <  1:3393720/69‑1 (MQ=255)
                          gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                                <  1:2435747/69‑1 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:994081/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:3376729/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:870030/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:286462/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2593161/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2433302/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:695189/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2321556/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:1586463/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:1580895/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:56160/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:139449/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:793880/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:1875093/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2357503/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:1961626/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2126888/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:377119/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:374175/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATACATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2077632/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTATTATACATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:2445032/1‑68 (MQ=255)
                            ttCATCTGCGTGAAGAATCTCGCCGTTAATATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                               >  1:865401/1‑68 (MQ=255)
                               aTCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACt                            <  1:2953729/68‑1 (MQ=255)
                                tCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTa                           <  1:1113879/68‑1 (MQ=255)
                                tCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTa                           <  1:1998388/68‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTACCg                        <  1:1466386/65‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                        <  1:249047/65‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                        <  1:1282612/65‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                        <  1:1494070/65‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                        <  1:2701702/65‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                      >  1:1004820/1‑67 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                      <  1:169959/67‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                      >  1:1107666/1‑67 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                      <  1:2675384/67‑1 (MQ=255)
                                      tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                      <  1:402147/67‑1 (MQ=255)
                                        aaGAGTCTTGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                  >  1:817909/1‑69 (MQ=255)
                                        aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                  >  1:1130336/1‑69 (MQ=255)
                                        aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                  >  1:1774297/1‑69 (MQ=255)
                                        aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                  >  1:976760/1‑69 (MQ=255)
                                        aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                  >  1:3111406/1‑69 (MQ=255)
                                           agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc               <  1:171839/69‑1 (MQ=255)
                                           agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc               <  1:1600629/69‑1 (MQ=255)
                                           agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc               <  1:3104636/69‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:2287579/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:2130973/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:1789328/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:2496160/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:3269866/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:3061179/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:477988/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:3258919/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg          <  1:1108104/66‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGAtt       <  1:1727681/69‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGAtt       <  1:1644161/69‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGAtt       <  1:3280127/69‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGAtt       <  1:3180286/69‑1 (MQ=255)
                                                   cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGAtt       <  1:2736009/69‑1 (MQ=255)
                                                         tATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAggg  <  1:1150894/68‑1 (MQ=255)
                                                           tccgtccgGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc               <  1:1225050/53‑1 (MQ=255)
                                                           tccgtccgGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc               <  1:95700/53‑1 (MQ=255)
                                                              |                                                              
TACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGG  >  minE/76600‑76724

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: