Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,349,295 C→T 61.8% A124T (GCC→ACC)  mtgA ← biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,349,2950CT61.8% 21.3 / 36.2 34A124T (GCC→ACC) mtgAbiosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (11/2);  new base T (7/14);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 5.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGG  >  W3110S.gb/3349238‑3349360
                                                         |                                                                 
cGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGaa                                                      <  1:3499495/71‑1 (MQ=255)
cGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGaa                                                      <  1:2521617/71‑1 (MQ=255)
         gCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTgaga                                                       <  1:729664/61‑1 (MQ=255)
         gCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTgaga                                                       <  1:354334/61‑1 (MQ=255)
          cccTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAg                                            <  1:2293995/71‑1 (MQ=255)
             ttttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGa                                              <  1:1241034/66‑1 (MQ=255)
               ttCGCACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACcg                                       >  1:2690599/1‑71 (MQ=255)
               ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACcg                                       >  1:983663/1‑71 (MQ=255)
               ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACcg                                       >  1:3204425/1‑71 (MQ=255)
               ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACcg                                       >  1:1808905/1‑71 (MQ=255)
                           tACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGttt                           >  1:3143489/1‑71 (MQ=255)
                             cGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTc                         <  1:1576162/71‑1 (MQ=255)
                                ccatcccatAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGc                                           <  1:3110412/50‑1 (MQ=255)
                                  atcccatAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGaa                                   <  1:752657/56‑1 (MQ=255)
                                   tcccatAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGtttt                          <  1:3681351/64‑1 (MQ=255)
                                   tcccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACcgcg                                     <  1:298640/53‑1 (MQ=255)
                                    cccatAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGtttt                          <  1:814979/63‑1 (MQ=255)
                                     ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAg                                            >  1:2748186/1‑44 (MQ=255)
                                     ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAg                                            >  1:3525901/1‑44 (MQ=255)
                                     ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTgcgc                  >  1:906226/1‑70 (MQ=255)
                                     ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCt                 >  1:1653694/1‑71 (MQ=255)
                                        tAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGt              <  1:708091/71‑1 (MQ=255)
                                          aaaaGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAg                                            >  1:557312/1‑39 (MQ=255)
                                          aaaaGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCa                                          >  1:2430797/1‑41 (MQ=255)
                                          aaaaGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTtg            >  1:3196649/1‑71 (MQ=255)
                                          aaaaGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTtg            >  1:1580610/1‑71 (MQ=255)
                                          aaaaGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTtg            >  1:2109279/1‑71 (MQ=255)
                                                gaTTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTgcgc       >  1:347542/1‑70 (MQ=255)
                                                gaTTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTgcgc       >  1:1014475/1‑70 (MQ=255)
                                                gaTTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGcc      <  1:2079535/71‑1 (MQ=255)
                                                    tttGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCa                        >  1:2163561/1‑49 (MQ=255)
                                                    tttGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCa                        >  1:2295991/1‑49 (MQ=255)
                                                     ttGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAggg  <  1:2738028/70‑1 (MQ=255)
                                                     ttGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGAGCGCCAggg  <  1:3709364/70‑1 (MQ=255)
                                                         |                                                                 
CGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGG  >  W3110S.gb/3349238‑3349360

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: