Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,349,295 C→T 45.2% 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,2950CT45.2% 6.4 / 38.4 31A124T (GCC→ACC) mtgAbiosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (10/7);  new base T (8/6);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

AATCCGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCT  >  W3110S.gb/3349234‑3349362
                                                             |                                                                   
aaTCCGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGttg                                                                   <  1:1537543/64‑1 (MQ=255)
aaTCCGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGttg                                                                   <  1:2067125/64‑1 (MQ=255)
aaTCCGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctgtctg                                                               >  1:3275878/1‑68 (MQ=255)
aaTCCGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctg                                                                   >  1:2170724/1‑64 (MQ=255)
    cGGCTTCCAGCCCTTTTCGCACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctg                                                                   >  1:61313/1‑60 (MQ=255)
      gCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGaa                                                        >  1:121987/1‑69 (MQ=255)
      gCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGaa                                                        >  1:651516/1‑69 (MQ=255)
      gCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGaa                                                        >  1:1773995/1‑69 (MQ=255)
                   ttCGAACCCATCTACGCCCATCCCATAAAAAGAGATTTTTGGttg                                                                   >  1:1816014/1‑45 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctg                                                                   <  1:2689203/45‑1 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctg                                                                   <  1:2742569/45‑1 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctg                                                                   <  1:1821583/45‑1 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGctg                                                                   <  1:1803178/45‑1 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGc                                             >  1:228541/1‑67 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTgaga                                                         <  1:1801445/55‑1 (MQ=255)
                   ttCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAg                                              >  1:2194807/1‑66 (MQ=255)
                         cccAGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGaa                                     <  1:90415/69‑1 (MQ=255)
                            aGCTACGCCCATCCCATAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATAc                                  <  1:3042723/69‑1 (MQ=255)
                               tACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGt                               >  1:2578965/1‑69 (MQ=255)
                                aCGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGtt                              >  1:2231621/1‑69 (MQ=255)
                                aCGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGtt                              >  1:2961329/1‑69 (MQ=255)
                                   cccatcccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGtttt                            <  1:1166052/68‑1 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCg                                                    >  1:1246847/1‑38 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTgcg                     >  1:3423828/1‑69 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAAt                                                      <  1:2031164/36‑1 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAAt                                                      <  1:295820/36‑1 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTg                                                 >  1:199453/1‑41 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGt                               >  1:1057728/1‑59 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTgcg                     >  1:3287547/1‑69 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTgcg                     >  1:2541098/1‑69 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTgcg                     >  1:2597875/1‑69 (MQ=255)
                                         ccatAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTgcg                     >  1:760180/1‑69 (MQ=255)
                                            tAAAAAGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTc                  <  1:1943906/69‑1 (MQ=255)
                                              aaaaGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCa                          >  1:2423802/1‑59 (MQ=255)
                                              aaaaGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCa                          >  1:2703492/1‑59 (MQ=255)
                                              aaaaGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGt                <  1:2886276/69‑1 (MQ=255)
                                              aaaaGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGt                <  1:2994667/69‑1 (MQ=255)
                                              aaaaGAGATTTTTGGTTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGt                <  1:3127298/69‑1 (MQ=255)
                                                   agaTTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTgc           <  1:2992710/69‑1 (MQ=255)
                                                     aTTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTgcgc         <  1:2372243/69‑1 (MQ=255)
                                                     aTTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTgcgc         <  1:2570096/69‑1 (MQ=255)
                                                             ctgtctgtTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCt  <  1:3156291/68‑1 (MQ=255)
                                                             |                                                                   
AATCCGGCTTCCAGCCCTTTTCGAACCCAGCTACGCCCATCCCATAAAAAGAGATTTTTGGCTGTCTGTTGAGAAATCGTTGAAGCACCGCGAATACGGTTTTCATTGCGCTCGTTGTGCGCCAGGGCT  >  W3110S.gb/3349234‑3349362

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: