Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,508,088 C→A 23.5% G615V (GGC→GTC)  ypdD ← fused predicted PTS enzymes Hpr component, enzyme I component, and enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,508,0880CA23.5% 139.3 / 21.1 55G615V (GGC→GTC) ypdDfused predicted PTS enzymes Hpr component, enzyme I component, and enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (27/15);  new base A (7/6);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCG  >  W3110S.gb/2508023‑2508158
                                                                 |                                                                      
ggCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGt                                                                   >  1:3500630/1‑71 (MQ=255)
       catactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCg                                                             >  1:4498549/1‑70 (MQ=255)
       catactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGACCAGttacc                                                               >  1:3244432/1‑65 (MQ=255)
        atactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCa                                                                >  1:2012922/1‑66 (MQ=255)
        atactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCa                                                                >  1:2018249/1‑66 (MQ=255)
        atactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCAcc                                                              >  1:2831667/1‑68 (MQ=255)
        atactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGc                                                            <  1:417601/70‑1 (MQ=255)
              aGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGcc                                                     >  1:4137061/1‑71 (MQ=255)
              aGCTCATCAAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCaa                                                          >  1:781101/1‑66 (MQ=255)
                 tCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCaaa                                                         >  1:4503226/1‑64 (MQ=255)
                  cATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTAcc                                                 >  1:987755/1‑71 (MQ=255)
                  cATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGACCAg                                                                    >  1:4213323/1‑52 (MQ=255)
                    tCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGACCAGTTCACCGCa                                                           <  1:4525653/59‑1 (MQ=255)
                     ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCaaa                                                         >  1:937576/1‑60 (MQ=255)
                     ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATg                                                       >  1:162272/1‑62 (MQ=255)
                     ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTAc                                                  >  1:4505914/1‑67 (MQ=255)
                     ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCAt                                              >  1:1201818/1‑71 (MQ=255)
                      cAGGCCCAGCCCAAGCAGTAGCGGCATATAACGGCTTTCACCGACCAGTTCACCGCa                                                           <  1:3995013/57‑1 (MQ=255)
                           ccagcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGcc                                                     >  1:3877366/1‑58 (MQ=255)
                           ccagcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGACCAGTTCACCGCaaa                                                         >  1:323351/1‑54 (MQ=255)
                              gcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGcc                                                                       <  1:1132791/37‑1 (MQ=255)
                               cccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCt                                    <  1:2555127/71‑1 (MQ=255)
                                ccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTg                                   <  1:3062116/71‑1 (MQ=255)
                                 caAGCAGTAGCGGCAGATAACGGCTTTCACCGACCAGTTCACCGCaaa                                                         >  1:3728082/1‑48 (MQ=255)
                                    gCAGTAGCGGCAGATAACGGCTTTCACCGACCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCt                                    <  1:2814217/66‑1 (MQ=255)
                                      aGTAGCGGCAGATAACGGCTTTCACCGACCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATgcgc                             <  1:4523455/71‑1 (MQ=255)
                                       gTAGCGGCAGATAAGGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCAGCTGATGCGCa                            <  1:3178810/71‑1 (MQ=255)
                                        tAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAg                           >  1:841079/1‑71 (MQ=255)
                                           cGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTAccc                                                <  1:2852069/47‑1 (MQ=255)
                                           cGGCAGATAACGGCTTTCACCGACCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTgg                        <  1:89602/71‑1 (MQ=255)
                                             gCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGt                       >  1:3765581/1‑70 (MQ=255)
                                                 aTAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTa                  >  1:992983/1‑71 (MQ=255)
                                                 aTAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTa                  >  1:1464028/1‑71 (MQ=255)
                                                  tAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAt                 <  1:2861335/71‑1 (MQ=255)
                                                  tAACGGCTTTCACCGACCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAt                 >  1:2329399/1‑71 (MQ=255)
                                                   aaCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCAtt                                             >  1:1519614/1‑42 (MQ=255)
                                                   aaCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATgcg                              >  1:598859/1‑57 (MQ=255)
                                                    aCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCa                            >  1:1766849/1‑58 (MQ=255)
                                                    aCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAtt                <  1:3730444/70‑1 (MQ=255)
                                                    aCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAtt                <  1:748707/70‑1 (MQ=255)
                                                    aCGGCTTTCACCGACCAGTTCACCGCAAATGCCTTCCCATTTGCCCCGCt                                    <  1:3313116/50‑1 (MQ=255)
                                                     cGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtg              <  1:2655117/71‑1 (MQ=255)
                                                      ggCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCa                            >  1:3988327/1‑56 (MQ=255)
                                                      ggCTTTCACCGACCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGAt                                 >  1:2321157/1‑51 (MQ=255)
                                                        cTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGa                                  >  1:3858910/1‑48 (MQ=255)
                                                        cTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGa                                  >  1:3853674/1‑48 (MQ=255)
                                                         tttCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCt                                    <  1:2007578/45‑1 (MQ=255)
                                                          ttCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTtgctgc          >  1:623726/1‑70 (MQ=255)
                                                           tCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAtt                >  1:707730/1‑63 (MQ=255)
                                                           tCACCGACCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACt                   >  1:4168716/1‑60 (MQ=255)
                                                            cACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACt                   >  1:53677/1‑59 (MQ=255)
                                                            cACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCaa        <  1:3191622/70‑1 (MQ=255)
                                                            cACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCaa        <  1:3085012/70‑1 (MQ=255)
                                                             aCCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTAt                 <  1:507000/60‑1 (MQ=255)
                                                             aCCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAAc       <  1:2813402/70‑1 (MQ=255)
                                                                gCCCAGTTCACCCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATgcg  <  1:2173346/71‑1 (MQ=255)
                                                                 cccAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATgcg  <  1:3666101/71‑1 (MQ=255)
                                                                 |                                                                      
GGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCG  >  W3110S.gb/2508023‑2508158

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: