Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,679,847 G→A 29.9% R55R (CGC→CGT pbpC ← fused transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,679,8470GA29.9% 47.5 / 19.8 50R55R (CGC→CGTpbpCfused transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (9/26);  new base A (12/3);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 5.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

ACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGG  >  minE/1679787‑1679909
                                                            |                                                              
aCGTGGAGAAACATCTTCGATTGTTACTGGATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCagag                                                        <  1:1013791/69‑1 (MQ=255)
aCGTGGAGAAACATCTTCGATTGTTACTGGATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCagag                                                        <  1:1183915/69‑1 (MQ=255)
aCGTGGAGAAACATCTTCGATTGTTACTGGATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCagag                                                        <  1:585360/69‑1 (MQ=255)
                         aCCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGAc                                <  1:1314839/68‑1 (MQ=255)
                         aCCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGAc                                <  1:910781/68‑1 (MQ=255)
                         aCCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGAc                                <  1:69460/68‑1 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCAt                                             >  1:1169010/3‑54 (MQ=38)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCAt                                             >  1:2762369/3‑54 (MQ=38)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGccaccacca                              >  1:65101/3‑69 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGccaccacca                              >  1:1790444/3‑69 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGccaccacca                              >  1:1876962/3‑69 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGccaccacca                              >  1:2915755/3‑69 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGccaccacca                              >  1:945285/3‑69 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGccaccacca                              >  1:2467055/3‑69 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGcc                                     >  1:1343055/3‑62 (MQ=39)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGCCaccacc                               >  1:3148372/3‑68 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTGCGCCaccacc                               >  1:2615319/3‑68 (MQ=255)
                          ctggATAACGCCAGATGCCGTCAGCATCGGCGAAACGCCAGAGTGGCGTACCATCCTACGccaccacca                              >  1:2077950/3‑69 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTg                        >  1:2660592/1‑63 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCggg                    >  1:391637/1‑67 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:1363845/1‑68 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:3307787/1‑68 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:3433260/1‑68 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:1915854/1‑68 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:751581/1‑68 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:893258/1‑68 (MQ=255)
                                      aGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGa                   >  1:2365515/1‑68 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:2209589/69‑1 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:2663726/69‑1 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:2752940/69‑1 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:1409724/69‑1 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:1534050/69‑1 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:3159159/69‑1 (MQ=255)
                                           ccGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACt             <  1:2431300/69‑1 (MQ=255)
                                             gTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:972870/68‑1 (MQ=255)
                                             gTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:367347/68‑1 (MQ=255)
                                             gTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:3224782/68‑1 (MQ=255)
                                             gTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:2672010/68‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:3106685/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:571439/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:1913514/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:2612499/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:86819/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:2482072/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:1329388/62‑1 (MQ=255)
                                                   aTCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:2297204/62‑1 (MQ=255)
                                                    tCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACtt            <  1:1883650/61‑1 (MQ=255)
                                                      ggCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCgg  <  1:1751540/69‑1 (MQ=255)
                                                      ggCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCgg  <  1:685693/69‑1 (MQ=255)
                                                      ggCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCgg  <  1:1340985/69‑1 (MQ=255)
                                                            |                                                              
ACGTGGAGAAACATCTTCGATTGTTACCGGATAACGCCAGATGCCGTCAGCATCGGCGAAGCGCCAGAGCGGCGTACCATCCTGCGCCACGACCACTCGTGCGGGATTGACTTCATGCAGCGG  >  minE/1679787‑1679909

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: