Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,679,130 C→T 40.9% R294R (CGG→CGA pbpC ← fused transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,679,1300CT40.9% 32.9 / 53.1 78R294R (CGG→CGApbpCfused transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (8/38);  new base T (15/17);  total (23/55)
Fisher's exact test for biased strand distribution p-value = 6.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

ACGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCC  >  minE/1679076‑1679184
                                                      |                                                      
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:3015094/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:560807/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:529958/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:2042763/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:575061/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:1822558/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:1717184/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:616519/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:1645727/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:1536038/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:821813/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:27470/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:2960034/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:884073/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:1480869/69‑1 (MQ=255)
acGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAATCGTCCTTTCCAGtt                                          <  1:858003/69‑1 (MQ=255)
       tCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGcc                                   >  1:2117483/1‑69 (MQ=255)
               aTCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGttct                             >  1:2459389/1‑67 (MQ=255)
               aTCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGttct                             >  1:597508/1‑67 (MQ=255)
               aTCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGttct                             >  1:2359707/1‑67 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:44806/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:493832/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:1431260/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:2377504/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:3254899/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:2875580/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:109695/69‑1 (MQ=255)
                      aCGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCACTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:2681446/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgt                    <  1:2740007/64‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:220941/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:2041467/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:1049718/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:994524/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:1246398/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:1277543/69‑1 (MQ=255)
                           catcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:744356/69‑1 (MQ=255)
                            atcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:2119599/68‑1 (MQ=255)
                            atcGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTg               <  1:2192660/68‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:623025/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:2414910/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:1600445/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:389425/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:3354140/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:1597643/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:802687/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:1483979/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:3341254/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:583771/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:2833197/69‑1 (MQ=255)
                              cGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:2273450/69‑1 (MQ=255)
                               gCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAg            <  1:1067652/68‑1 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGAcgtc                   >  1:1519527/1‑56 (MQ=38)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCgg       >  1:778881/1‑68 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:2252758/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:1089744/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:1491581/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:1608054/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:1704944/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:1741820/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:491290/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:2603456/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:393056/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:2785853/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:2665815/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:3407129/1‑69 (MQ=255)
                                    tGAGCTGCGCGGTGGCAATCGTCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGc      >  1:3407133/1‑69 (MQ=255)
                                     gAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCa     >  1:161942/1‑69 (MQ=255)
                                        cTGCGCGGTGGCAATCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:1688747/69‑1 (MQ=255)
                                        cTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:3228117/69‑1 (MQ=255)
                                        cTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:1356539/69‑1 (MQ=255)
                                        cTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:260397/69‑1 (MQ=255)
                                        cTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:203053/69‑1 (MQ=255)
                                        cTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:2200783/69‑1 (MQ=255)
                                        cTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:475440/69‑1 (MQ=255)
                                                 ggCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATcc  <  1:2505453/60‑1 (MQ=255)
                                                  gCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCa                         >  1:2764068/1‑36 (MQ=255)
                                                  gCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCa                         >  1:3097974/1‑36 (MQ=255)
                                                  gCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCa                         >  1:1483406/1‑36 (MQ=255)
                                                      |                                                      
ACGCATATCGGTATGATCAACCACGATCATCGCCAGTGAGCTGCGCGGTGGCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCC  >  minE/1679076‑1679184

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: