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.

GGCAGCACAATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑AAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGCC  >  W3110S.gb/3493229‑3493363
                                                                      |                                                                  
ggCAGCACAATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGcaccac                                                                     >  1:2908520/1‑70 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑aaaact                                                              >  1:1356472/1‑67 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑aaaact                                                              >  1:4494431/1‑67 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑aaaact                                                              >  1:1569096/1‑67 (MQ=255)
      acaATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑aaaact                                                              >  1:3673124/1‑67 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAAcc                                                         <  1:2186587/71‑1 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAAcc                                                         >  1:2930281/1‑71 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAAcc                                                         <  1:1848577/71‑1 (MQ=255)
          tCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAAc                                                          >  1:3144808/1‑70 (MQ=255)
           cAGACAACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAAcc                                                         >  1:2974530/1‑70 (MQ=255)
             gACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGAc                                                      <  1:1135207/71‑1 (MQ=255)
             gACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGAc                                                      <  1:4581590/71‑1 (MQ=255)
               ccaccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCa                                                    >  1:1081914/1‑71 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCTACCat                                                   <  1:1367852/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACTat                                                   <  1:413529/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCat                                                   <  1:3172858/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCat                                                   <  1:854461/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCat                                                   <  1:2137931/71‑1 (MQ=255)
                caccCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑AAACCTGAACCGACCat                                                   <  1:3337242/71‑1 (MQ=255)
                    ccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCatgatg                                               >  1:530754/1‑71 (MQ=255)
                    ccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCatgatg                                               >  1:1334673/1‑71 (MQ=255)
                    ccAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCatgatg                                               >  1:943172/1‑71 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCTACCATGATggaagg                                          <  1:2187113/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:3728017/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:1089566/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:4431773/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:2982784/71‑1 (MQ=255)
                         cTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:754701/71‑1 (MQ=255)
                           gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaaggaa                                        >  1:2315479/1‑71 (MQ=255)
                           gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaaggaa                                        >  1:2580879/1‑71 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagga                                         <  1:1535159/69‑1 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:2846081/68‑1 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAAt                                       <  1:654904/71‑1 (MQ=255)
                            cATAGCACTGCGCACCGACCAGCGCCGAAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAAt                                       <  1:2672746/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGa                                     <  1:3617619/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGa                                     <  1:1835821/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGa                                     <  1:3265386/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGa                                     <  1:2807371/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGa                                     <  1:808666/71‑1 (MQ=255)
                              tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGa                                     >  1:2326174/1‑71 (MQ=255)
                                gCACTGCGCACCGCCCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGaca                                   >  1:3776659/1‑71 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGaa                                                           >  1:2690187/1‑47 (MQ=38)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGaa                                                           >  1:4584204/1‑47 (MQ=38)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGa                                                       >  1:2276868/1‑51 (MQ=38)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCatg                                                  >  1:1959357/1‑56 (MQ=39)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaaggaa                                        >  1:1689045/1‑66 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAAt                                       >  1:2803396/1‑67 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGaca                                   >  1:2788639/1‑71 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGaca                                   >  1:2434097/1‑71 (MQ=255)
                                gCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGaca                                   >  1:303765/1‑71 (MQ=255)
                                   cTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGc                                <  1:3789032/71‑1 (MQ=255)
                                    tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:3898809/60‑1 (MQ=255)
                                    tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:1940941/60‑1 (MQ=255)
                                     gcgcACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAg                              >  1:1275068/1‑71 (MQ=255)
                                     gcgcACCGACCAGCGCCGCAGTTACTGCACCAC‑aaaa                                                                >  1:4445190/1‑36 (MQ=255)
                                      cgcACGGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGc                             >  1:2493045/1‑71 (MQ=255)
                                      cgcACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGc                             >  1:566173/1‑71 (MQ=255)
                                       gcACCGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:144932/57‑1 (MQ=39)
                                           cGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATc                        <  1:4206556/71‑1 (MQ=255)
                                           cGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATc                        <  1:753602/71‑1 (MQ=255)
                                           cGACCAGCGCCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATc                        <  1:1181284/71‑1 (MQ=255)
                                                 gcgcCGCAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgc                  <  1:1895033/71‑1 (MQ=255)
                                                  cgccgcAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcg                 >  1:3426649/1‑71 (MQ=255)
                                                  cgccgcAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcg                 >  1:2097316/1‑71 (MQ=255)
                                                  cgccgcAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcg                 >  1:3165307/1‑71 (MQ=255)
                                                   gccgcAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaaggaa                                        >  1:1454200/1‑47 (MQ=39)
                                                   gccgcAGTTACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAgcgg                >  1:490150/1‑71 (MQ=255)
                                                       cAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagga                                         <  1:1367407/42‑1 (MQ=255)
                                                       cAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagga                                         <  1:3712414/42‑1 (MQ=255)
                                                       cAGTTACTGCACCACAAAA‑CTGAACCGACCATGATggaagga                                         <  1:3150647/42‑1 (MQ=255)
                                                          ttACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:2321156/38‑1 (MQ=38)
                                                          ttACTGCACCACAAAA‑CTGAACCGACCATGATggaagg                                          <  1:4010313/38‑1 (MQ=38)
                                                           tACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:2215352/1‑70 (MQ=255)
                                                           tACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:2520118/1‑70 (MQ=255)
                                                           tACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:713655/1‑70 (MQ=255)
                                                           tACTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg         >  1:4101341/1‑70 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGaaa           >  1:3883251/1‑67 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAg          >  1:3107740/1‑68 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:642337/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3491298/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:2603950/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:2385138/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3100736/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3425972/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3135558/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGaa       >  1:3407848/1‑71 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:4350867/1‑70 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:3551975/70‑1 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:883358/70‑1 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:1701069/70‑1 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:415576/70‑1 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:2060833/1‑70 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:1969076/1‑70 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        >  1:4545293/1‑70 (MQ=255)
                                                            aCTGCACCACAAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGa        <  1:4439089/70‑1 (MQ=255)
                                                                caccacaAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGc   <  1:4118535/71‑1 (MQ=255)
                                                                caccacaAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGc   <  1:3584047/71‑1 (MQ=255)
                                                                 accacaAAA‑CTGAACCGACCATTATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGcc  <  1:1392130/71‑1 (MQ=255)
                                                                 accacaAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGcc  <  1:3864295/71‑1 (MQ=255)
                                                                 accacaAAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGcc  <  1:1807212/71‑1 (MQ=255)
                                                                      |                                                                  
GGCAGCACAATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCAC‑AAA‑CTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGGAATCGCC  >  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: ‑