Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 76,644 G→A 47.2% E164E (GAG→GAA ftsI → transpeptidase involved in septal peptidoglycan synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE76,6440GA47.2% 9.1 / 78.1 91E164E (GAG→GAAftsItranspeptidase involved in septal peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base A (22/21);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTC  >  minE/76585‑76711
                                                           |                                                                   
ccTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGt                                                            <  1:2109683/69‑1 (MQ=255)
      aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTa                                                       <  1:765170/68‑1 (MQ=255)
      aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTa                                                       <  1:1622731/68‑1 (MQ=255)
      aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTa                                                       <  1:2025542/68‑1 (MQ=255)
      aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAt                                                      >  1:1255953/1‑69 (MQ=255)
      aTGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAt                                                      >  1:2101407/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGt                                                                    >  1:210042/1‑53 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACt                                                        >  1:2390101/1‑65 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:1485901/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:1041248/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:1075223/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:3125999/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:2263910/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:219542/1‑69 (MQ=255)
        ggCGGACTACATCAAAAAACTGAAACTGACGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTAtcc                                                    >  1:2799384/1‑69 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:535651/1‑68 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:1631488/1‑68 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:3362911/1‑68 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:3095836/1‑68 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:1002303/1‑68 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:553399/1‑68 (MQ=255)
               tACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCgg                                              >  1:1283810/1‑68 (MQ=255)
                    cAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGt                                        >  1:2673180/1‑69 (MQ=255)
                       aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                     <  1:1015712/69‑1 (MQ=255)
                       aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                     <  1:934593/69‑1 (MQ=255)
                       aaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                     <  1:1196489/69‑1 (MQ=255)
                        aaactgaaactgCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGAc                                     <  1:1680330/68‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTTAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:780613/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:2013892/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:918616/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:2641038/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:1308214/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:1831536/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:1819949/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:428820/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:1389752/69‑1 (MQ=255)
                                actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:1711784/69‑1 (MQ=255)
                                 ctgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:3015986/68‑1 (MQ=255)
                                 ctgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                            <  1:1107494/68‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:2868887/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:384056/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:2381278/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:2825082/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:2608426/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:1170104/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:2120674/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:2937030/69‑1 (MQ=255)
                                        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                    <  1:844886/69‑1 (MQ=255)
                                         gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                   <  1:2435747/69‑1 (MQ=255)
                                         gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                   <  1:3393720/69‑1 (MQ=255)
                                         gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                   <  1:448353/69‑1 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:994081/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:286462/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:56160/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:870030/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:3376729/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:374175/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:377119/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:793880/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:695189/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:1580895/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:1875093/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:1961626/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2126888/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2321556/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2357503/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:139449/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2433302/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:1586463/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2593161/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATACATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2445032/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTATTATACATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:2077632/1‑68 (MQ=255)
                                           ttCATCTGCGTGAAGAATCTCGCCGTTAATATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                  >  1:865401/1‑68 (MQ=255)
                                              aTCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACt               <  1:2953729/68‑1 (MQ=255)
                                               tCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTa              <  1:1998388/68‑1 (MQ=255)
                                               tCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTa              <  1:1113879/68‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTACCg           <  1:1466386/65‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg           <  1:1494070/65‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg           <  1:1282612/65‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg           <  1:249047/65‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg           <  1:2701702/65‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc         >  1:1107666/1‑67 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc         >  1:1004820/1‑67 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc         <  1:2675384/67‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc         <  1:169959/67‑1 (MQ=255)
                                                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc         <  1:402147/67‑1 (MQ=255)
                                                       aaGAGTCTTGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa     >  1:817909/1‑69 (MQ=255)
                                                       aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa     >  1:1774297/1‑69 (MQ=255)
                                                       aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa     >  1:1130336/1‑69 (MQ=255)
                                                       aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa     >  1:3111406/1‑69 (MQ=255)
                                                       aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa     >  1:976760/1‑69 (MQ=255)
                                                          agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc  <  1:1600629/69‑1 (MQ=255)
                                                          agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc  <  1:171839/69‑1 (MQ=255)
                                                          agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc  <  1:3104636/69‑1 (MQ=255)
                                                           |                                                                   
CCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTC  >  minE/76585‑76711

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: