Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,966,285 T→C 21.2% E269G (GAG→GGG)  ptsP ← fused PEP‑protein phosphotransferase (enzyme I) of PTS system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,966,2850TC21.2% 63.2 / 14.4 33E269G (GAG→GGG) ptsPfused PEP‑protein phosphotransferase (enzyme I) of PTS system
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base C (3/4);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAG  >  W3110S.gb/2966217‑2966346
                                                                    |                                                             
tGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCActct                                                             >  1:2114498/1‑71 (MQ=255)
tGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCAct                                                               <  1:4470276/69‑1 (MQ=255)
     tAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAc                                                        >  1:4491620/1‑71 (MQ=255)
      aGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAcc                                                       >  1:4015731/1‑71 (MQ=255)
       gTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCAct                                                               >  1:868419/1‑62 (MQ=255)
       gTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCAcc                                                       <  1:4002295/70‑1 (MQ=255)
         tATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGc                                                           <  1:4581940/64‑1 (MQ=255)
          aTCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAg                                                    >  1:2119916/1‑70 (MQ=255)
          aTCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGc                                                   >  1:1449984/1‑71 (MQ=255)
           tCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCActc                                                              >  1:301444/1‑59 (MQ=255)
             aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCtt                                                <  1:1264046/71‑1 (MQ=255)
             aCTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCtt                                                >  1:3595462/1‑71 (MQ=255)
                tAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACCCTGCCACCGAGCCTTTAt                                             >  1:346319/1‑71 (MQ=255)
                 aGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTatat                                             >  1:2469741/1‑70 (MQ=255)
                  gcgcGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATc                                            <  1:3015051/70‑1 (MQ=255)
                                tCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACaa                             >  1:819216/1‑71 (MQ=255)
                                 cGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAAt                            >  1:1853701/1‑71 (MQ=255)
                                  ggCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAAtt                           <  1:4051661/71‑1 (MQ=255)
                                   gCAAATTTTTCAATGACCGTTTTTACCGCCCACCCTGCCACCGAGCCTTTATc                                            <  1:4594297/53‑1 (MQ=255)
                                   gCAAATTTTTCAATGACAGTTTTTACCGCCCACCCTGCCACCGAGCCTTTATCAACCTCGGCAAACAAtt                           <  1:2819224/70‑1 (MQ=255)
                                    cAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCtt                                                <  1:222976/48‑1 (MQ=255)
                                     aaaTTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTcgc                        <  1:2873604/71‑1 (MQ=255)
                                     aaaTTTTTCAATGACCGCTTTTACCGCCCACCCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTcgc                        <  1:2178333/71‑1 (MQ=255)
                                       aTTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTAt                                             <  1:3676792/48‑1 (MQ=255)
                                        tttttCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGAc                     <  1:76397/71‑1 (MQ=255)
                                        tttttCAATGACCGTTTTTACCGCCCACCCTGCCACCGAGCCTTTATCaa                                          >  1:4328975/1‑50 (MQ=255)
                                        tttttCAATGACCGTTTTTACCGCCCACCCTGCCACCGAGCCTTTATCa                                           <  1:2824688/49‑1 (MQ=255)
                                         ttttCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTa                                              >  1:1420772/1‑45 (MQ=255)
                                              aaTGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGcc               <  1:4412355/71‑1 (MQ=255)
                                               aTGACCGTTTTTACCGCCCACTCTGCCACCGAGCCtt                                                <  1:3331263/37‑1 (MQ=255)
                                                tGACCGTTTTTACCGCCCACTCTGCCACCGAGCCtt                                                <  1:2463189/36‑1 (MQ=255)
                                                tGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATc                                            <  1:517060/40‑1 (MQ=255)
                                                  aCCGTTTTTACCGCCCACCCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGc                >  1:118011/1‑66 (MQ=255)
                                                    cGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATc                                            >  1:426775/1‑36 (MQ=255)
                                                    cGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCaa                                          >  1:2072823/1‑38 (MQ=255)
                                                      tttttACCGCCCACTCTGCCACCGAGCCTTTATCAACCt                                       <  1:2947527/39‑1 (MQ=255)
                                                          tACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGaa    <  1:1185424/70‑1 (MQ=255)
                                                                   ctctGCCACCGAGCCTTTATCAACCTCGGCAAACAATTcgc                        <  1:1012789/41‑1 (MQ=255)
                                                                   ctctGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAg  <  1:2543656/63‑1 (MQ=255)
                                                                    |                                                             
TGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAG  >  W3110S.gb/2966217‑2966346

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: