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.A96.5% 91.4 / ‑2.8 29coding (854/1080 nt)frwCpredicted enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (7/21);  total (7/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

ATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGG  >  W3110S.gb/3493238‑3493356
                                                             |                                                           
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:669066/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:2205349/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:1974138/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:748025/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:1018464/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:1879189/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:2323646/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:1676493/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:1665669/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:475993/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:870644/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:1370544/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:2905672/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:3015652/65‑2 (MQ=255)
aTCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:3071385/65‑2 (MQ=255)
               ccTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCatga                                          >  1:1607427/1‑65 (MQ=255)
               ccTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCatga                                          >  1:1733483/1‑65 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:404254/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:2891641/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:2495679/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:448361/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:1893749/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:1556570/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:1337969/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATgg                                       <  1:1049367/65‑1 (MQ=255)
                  gCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAATCTGAACCGACCATGATgg                                       <  1:208341/65‑1 (MQ=255)
                     tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaag                                    <  1:2263609/65‑1 (MQ=255)
                     tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaag                                    <  1:1539747/65‑1 (MQ=255)
                     tAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaag                                    <  1:967619/65‑1 (MQ=255)
                      aGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaag                                    <  1:844917/64‑1 (MQ=255)
                         aCTGCGCACCGACCAGCGCCGCAGTTACTGCACCACaaaa                                                         <  1:41033/40‑2 (MQ=255)
                           tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                   >  1:1848409/1‑60 (MQ=255)
                           tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                   >  1:1787999/1‑60 (MQ=255)
                           tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                   >  1:2844985/1‑60 (MQ=255)
                           tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                   >  1:112459/1‑60 (MQ=255)
                           tGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAAACTGAACCGACCATGATggaagg                                   >  1:3011527/1‑60 (MQ=255)
                                                   aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGa      <  1:2288707/65‑1 (MQ=255)
                                                   aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGa      <  1:299313/65‑1 (MQ=255)
                                                   aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGa      <  1:1951217/65‑1 (MQ=255)
                                                   aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGa      <  1:1614943/65‑1 (MQ=255)
                                                   aCTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGa      <  1:1538170/65‑1 (MQ=255)
                                                    cTGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGa      <  1:347185/64‑1 (MQ=255)
                                                     tGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGaa     <  1:2799743/64‑1 (MQ=255)
                                                     tGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGaa     <  1:2605438/64‑1 (MQ=255)
                                                     tGCACCACAAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGaa     <  1:940120/64‑1 (MQ=255)
                                                       caccacaAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAgg  >  1:1043150/1‑65 (MQ=255)
                                                             |                                                           
ATCAGACCACCCCAGCCTGCATAGCACTGCGCACCGACCAGCGCCGCAGTTACTGCACCACAAACTGAACCGACCATGATGGAAGGAATGACACGCAGCGGATCGGCAGCGGCGAAAGG  >  W3110S.gb/3493238‑3493356

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: