Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,513,100 T→C 11.4% N81D (AAT→GAT)  ypdG ← predicted enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,513,1000TC11.4% 88.2 / 4.9 35N81D (AAT→GAT) ypdGpredicted enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (11/20);  new base C (2/2);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGA  >  W3110S.gb/2513030‑2513165
                                                                      |                                                                 
ccgccgAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTAtt                                                                   <  1:3634062/71‑1 (MQ=255)
ccgccgAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTAtt                                                                   <  1:306361/71‑1 (MQ=255)
ccgccgAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTAtt                                                                   <  1:1883819/71‑1 (MQ=255)
  gccgAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATCaa                                                                 <  1:3684093/71‑1 (MQ=255)
     gAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAagag                                                              >  1:4610575/1‑71 (MQ=255)
          aCTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATcc                                                         <  1:3912924/71‑1 (MQ=255)
            tGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATCAAGAGCATCCat                                                       >  1:4311385/1‑71 (MQ=255)
              cGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCat                                                       <  1:5028565/69‑1 (MQ=255)
              cGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCat                                                       <  1:3030998/69‑1 (MQ=255)
                       ggCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGAtct                                            <  1:2378342/71‑1 (MQ=255)
                        gCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTaa                                                                 >  1:2264338/1‑49 (MQ=255)
                        gCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGAtct                                            >  1:3764245/1‑70 (MQ=255)
                        gCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGAtct                                            >  1:4197692/1‑70 (MQ=255)
                          aaaTTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATCAAGAGCATCCATGATGCCGAtctc                                           >  1:3358618/1‑69 (MQ=255)
                           aaTTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTAtt                                                                   <  1:3959838/44‑1 (MQ=255)
                             tttAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCa                                                            >  1:1873178/1‑49 (MQ=255)
                                     agGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTa                              <  1:2298813/71‑1 (MQ=255)
                                     agGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTa                              <  1:306149/71‑1 (MQ=255)
                                      ggAGAGGGCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTaa                             <  1:4323338/71‑1 (MQ=255)
                                         gagGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGcc                          <  1:4778600/71‑1 (MQ=255)
                                         gagGTCGAAACCGGAGAATTTCCCGCTATCAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGcc                          <  1:401520/71‑1 (MQ=255)
                                          agGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCa                         <  1:4339443/71‑1 (MQ=255)
                                            gTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAgg                                     <  1:4396309/57‑1 (MQ=255)
                                            gTCGAAACCGGAGAATTTCCAGCTATTAAGAGCATCCATGATGCCGATCTCAGCAgg                                     <  1:3921021/57‑1 (MQ=255)
                                             tCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAAt                                  <  1:2575614/59‑1 (MQ=255)
                                             tCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTaa                             <  1:2826729/64‑1 (MQ=255)
                                              cGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCa                                                        >  1:2566928/1‑36 (MQ=255)
                                                aaaCCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCtat                    <  1:1502484/70‑1 (MQ=255)
                                                     ggAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCaa              >  1:2550299/1‑71 (MQ=255)
                                                      gagaATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAAtttt                               <  1:2331974/53‑1 (MQ=255)
                                                        gaATTTCCCGCTATTAAGAGCATCAATGATGCCGATCTcagc                                        >  1:2297429/1‑42 (MQ=255)
                                                         aaTTTCCCGCTATTAAGAGCATCCATGATGCCGATCTcag                                         <  1:4169918/40‑1 (MQ=255)
                                                         aaTTTCCCGCTATTAAGAGCATCCATGATGCCGATCTcag                                         <  1:168436/40‑1 (MQ=255)
                                                         aaTTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAg                            >  1:2726670/1‑53 (MQ=255)
                                                          aTTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAg          <  1:3702978/70‑1 (MQ=255)
                                                                 gCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGa  <  1:748929/71‑1 (MQ=255)
                                                                   tATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTa                              >  1:1430888/1‑41 (MQ=255)
                                                                   tATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAg                 <  1:1458141/54‑1 (MQ=255)
                                                                   tATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGa  >  1:4599775/1‑69 (MQ=255)
                                                                     ttAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAAt             <  1:3634658/56‑1 (MQ=255)
                                                                      |                                                                 
CCGCCGAAGGACTGCGACAGCCAGGCAAATTTAAGTAAGGAGAGGTCGAAACCGGAGAATTTCCCGCTATTAAGAGCATCCATGATGCCGATCTCAGCAGGAATTTTAAGCCAGCTATAAGCAATCAGCTGGCTGA  >  W3110S.gb/2513030‑2513165

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: