Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 76,677 T→C 50.0% T175T (ACT→ACC ftsI → transpeptidase involved in septal peptidoglycan synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE76,6770TC50.0% ‑3.3 / 78.8 86T175T (ACT→ACCftsItranspeptidase involved in septal peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base T (7/36);  new base C (22/21);  total (29/57)
Fisher's exact test for biased strand distribution p-value = 1.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

ACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAG  >  minE/76617‑76734
                                                            |                                                         
actgCCGGGGATTCATCTGCGTTAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:780613/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:2641038/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:918616/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:1831536/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:1819949/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:1711784/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:2013892/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:428820/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:1308214/69‑1 (MQ=255)
actgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:1389752/69‑1 (MQ=255)
 ctgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:1107494/68‑1 (MQ=255)
 ctgCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCt                                                   <  1:3015986/68‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:2120674/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:384056/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:1170104/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:2381278/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:2937030/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:2868887/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:2608426/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:2825082/69‑1 (MQ=255)
        ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCt                                           <  1:844886/69‑1 (MQ=255)
         gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                                          <  1:448353/69‑1 (MQ=255)
         gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                                          <  1:2435747/69‑1 (MQ=255)
         gATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCtt                                          <  1:3393720/69‑1 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:2433302/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:2593161/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:994081/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:2357503/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:286462/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:56160/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:695189/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:793880/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:870030/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:2321556/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:2126888/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:139449/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:1961626/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:1875093/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:1580895/1‑68 (MQ=255)
           ttCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACCGCTCACCTCATCGGCttt                                         >  1:1586463/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:2701702/65‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                                  <  1:1494070/65‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                                  <  1:249047/65‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                                  <  1:1282612/65‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                                <  1:402147/67‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                                >  1:1004820/1‑67 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                                <  1:169959/67‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                                <  1:2675384/67‑1 (MQ=255)
                     tGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTc                                >  1:1107666/1‑67 (MQ=255)
                       aaGAGTCTTGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                            >  1:817909/1‑69 (MQ=255)
                       aaGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATa                            >  1:3111406/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)
                          agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc                         <  1:3104636/69‑1 (MQ=255)
                          agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc                         <  1:171839/69‑1 (MQ=255)
                          agTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTc                         <  1:1600629/69‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:477988/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:1108104/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:3269866/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:3258919/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:3061179/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:1789328/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:2496160/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:2130973/66‑1 (MQ=255)
                                  cGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAggg                    <  1:2287579/66‑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:1727681/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)
                                                  gCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGagaag  <  1:404446/68‑1 (MQ=255)
                                                   cGAAGTGACTGCTCACCTCATCGGCTTTACTAACg                                  <  1:964976/35‑1 (MQ=255)
                                                            |                                                         
ACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAG  >  minE/76617‑76734

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: