Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,493,301 +A coding (851/1080 nt) frwC ← predicted enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,493,2981.A98.9% 275.2 / ‑3.4 92coding (854/1080 nt)frwCpredicted enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (49/42);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 103 total aligned reads displayed.

GGCAGCACAATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGCC  >  W3110S.gb/3493229‑3493363
                                                                      |                                                                  
ggCAGCACAATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGcaccac                                                                     >  1:2908520/1‑70 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaact                                                              >  1:1356472/1‑67 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaact                                                              >  1:4494431/1‑67 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaact                                                              >  1:1569096/1‑67 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaact                                                              >  1:3673124/1‑67 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAAcc                                                         <  1:2186587/71‑1 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAAcc                                                         >  1:2930281/1‑71 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAAcc                                                         <  1:1848577/71‑1 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAAc                                                          >  1:3144808/1‑70 (MQ=255)
           cAGACAACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAAcc                                                         >  1:2974530/1‑70 (MQ=255)
             gACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGAc                                                      <  1:1135207/71‑1 (MQ=255)
             gACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGAc                                                      <  1:4581590/71‑1 (MQ=255)
               ccaccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCa                                                    >  1:1081914/1‑71 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCTACCat                                                   <  1:1367852/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACTat                                                   <  1:413529/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCat                                                   <  1:3172858/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCat                                                   <  1:854461/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCat                                                   <  1:2137931/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAACCTGAACCGACCat                                                   <  1:3337242/71‑1 (MQ=255)
                    ccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCatgatg                                               >  1:530754/1‑71 (MQ=255)
                    ccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCatgatg                                               >  1:1334673/1‑71 (MQ=255)
                    ccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCatgatg                                               >  1:943172/1‑71 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCTACCATGATggaagg                                          <  1:2187113/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:3728017/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:1089566/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:4431773/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:2982784/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:754701/71‑1 (MQ=255)
                           gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaaggaa                                        >  1:2315479/1‑71 (MQ=255)
                           gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaaggaa                                        >  1:2580879/1‑71 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagga                                         <  1:1535159/69‑1 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:2846081/68‑1 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAAt                                       <  1:654904/71‑1 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGAAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAAt                                       <  1:2672746/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGa                                     <  1:3617619/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGa                                     <  1:1835821/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGa                                     <  1:3265386/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGa                                     <  1:2807371/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGa                                     <  1:808666/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGa                                     >  1:2326174/1‑71 (MQ=255)
                                gCACTGCGCACCGCCCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGaca                                   >  1:3776659/1‑71 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGaa                                                           >  1:2690187/1‑47 (MQ=38)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGaa                                                           >  1:4584204/1‑47 (MQ=38)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGa                                                       >  1:2276868/1‑51 (MQ=38)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCatg                                                  >  1:1959357/1‑56 (MQ=39)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaaggaa                                        >  1:1689045/1‑66 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAAt                                       >  1:2803396/1‑67 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGaca                                   >  1:2788639/1‑71 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGaca                                   >  1:2434097/1‑71 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGaca                                   >  1:303765/1‑71 (MQ=255)
                                   cTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGc                                <  1:3789032/71‑1 (MQ=255)
                                    tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:3898809/60‑1 (MQ=255)
                                    tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:1940941/60‑1 (MQ=255)
                                     gcgcACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAg                              >  1:1275068/1‑71 (MQ=255)
                                     gcgcACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                                >  1:4445190/1‑36 (MQ=255)
                                      cgcACGGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGc                             >  1:2493045/1‑71 (MQ=255)
                                      cgcACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGc                             >  1:566173/1‑71 (MQ=255)
                                       gcACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:144932/57‑1 (MQ=39)
                                           cGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATc                        <  1:4206556/71‑1 (MQ=255)
                                           cGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATc                        <  1:753602/71‑1 (MQ=255)
                                           cGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATc                        <  1:1181284/71‑1 (MQ=255)
                                                 gcgcCGCAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgc                  <  1:1895033/71‑1 (MQ=255)
                                                  cgccgcAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcg                 >  1:3426649/1‑71 (MQ=255)
                                                  cgccgcAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcg                 >  1:2097316/1‑71 (MQ=255)
                                                  cgccgcAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcg                 >  1:3165307/1‑71 (MQ=255)
                                                   gccgcAGTTACTGCACCACAAAACTGAACCGACCATGATggaaggaa                                        >  1:1454200/1‑47 (MQ=39)
                                                   gccgcAGTTACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcgg                >  1:490150/1‑71 (MQ=255)
                                                       cAGTTACTGCACCACAAAACTGAACCGACCATGATggaagga                                         <  1:1367407/42‑1 (MQ=255)
                                                       cAGTTACTGCACCACAAAACTGAACCGACCATGATggaagga                                         <  1:3712414/42‑1 (MQ=255)
                                                       cAGTTACTGCACCACAAAACTGAACCGACCATGATggaagga                                         <  1:3150647/42‑1 (MQ=255)
                                                          ttACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:2321156/38‑1 (MQ=38)
                                                          ttACTGCACCACAAAACTGAACCGACCATGATggaagg                                          <  1:4010313/38‑1 (MQ=38)
                                                           tACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:2215352/1‑70 (MQ=255)
                                                           tACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:2520118/1‑70 (MQ=255)
                                                           tACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:713655/1‑70 (MQ=255)
                                                           tACTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:4101341/1‑70 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGaaa           >  1:3883251/1‑67 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAg          >  1:3107740/1‑68 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:642337/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3491298/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:2603950/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:2385138/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3100736/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3425972/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3135558/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3407848/1‑71 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:4350867/1‑70 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:3551975/70‑1 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:883358/70‑1 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:1701069/70‑1 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:415576/70‑1 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:2060833/1‑70 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:1969076/1‑70 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:4545293/1‑70 (MQ=255)
                                                            aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:4439089/70‑1 (MQ=255)
                                                                caccacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGc   <  1:4118535/71‑1 (MQ=255)
                                                                caccacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGc   <  1:3584047/71‑1 (MQ=255)
                                                                 accacaAAACTGAACCGACCATTATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGcc  <  1:1392130/71‑1 (MQ=255)
                                                                 accacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGcc  <  1:3864295/71‑1 (MQ=255)
                                                                 accacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGcc  <  1:1807212/71‑1 (MQ=255)
                                                                      |                                                                  
GGCAGCACAATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGCC  >  W3110S.gb/3493229‑3493363

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: