Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,516,285 T→C 6.9% intergenic (+20/‑47) hslU → / → menA molecular chaperone and ATPase component of HslUV protease/1,4‑dihydroxy‑2‑naphthoate octaprenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,516,2850TC 6.9% 156.6 / 4.1 58intergenic (+20/‑47)hslU/menAmolecular chaperone and ATPase component of HslUV protease/1,4‑dihydroxy‑2‑naphthoate octaprenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/30);  new base C (2/2);  total (26/32)
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.97e-01

GATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTAGCCGA  >  W3110S.gb/3516221‑3516355
                                                                |                                                                      
gATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTtcatca                                                                      <  1:1437602/67‑1 (MQ=255)
gATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTtcacca                                                                      <  1:2978020/67‑1 (MQ=255)
  tGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGttt                                                                <  1:215914/71‑1 (MQ=255)
    cGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                              >  1:438653/1‑71 (MQ=255)
       tggtggCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATg                                                            >  1:2890528/1‑70 (MQ=255)
            gCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCt                                                       <  1:3146136/70‑1 (MQ=255)
            gCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCt                                                       >  1:3898787/1‑70 (MQ=255)
            gCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCt                                                       >  1:79438/1‑70 (MQ=255)
             cAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGa                                                     <  1:2951828/71‑1 (MQ=255)
             cAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTGa                                                     <  1:3737196/71‑1 (MQ=255)
              aGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGttt                                                                <  1:3637369/59‑1 (MQ=255)
                aTGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                              <  1:2281288/59‑1 (MQ=255)
                aTGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                              <  1:2412586/59‑1 (MQ=255)
                aTGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGa                                                     <  1:2201308/68‑1 (MQ=255)
                aTGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAg                                                  <  1:4146224/71‑1 (MQ=255)
                 tGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTtcat                                                                        >  1:4095946/1‑48 (MQ=255)
                  gAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATg                                                            >  1:1908225/1‑59 (MQ=255)
                  gAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTg                                                               >  1:2367583/1‑56 (MQ=255)
                   aaGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGcc                                               >  1:4074496/1‑71 (MQ=255)
                    aGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGccc                                              >  1:416518/1‑71 (MQ=255)
                    aGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGccc                                              <  1:1675723/71‑1 (MQ=255)
                    aGATCTGAGCCGTTTTATCCTATAATAGCGTTCAATCATTTTCATCATTGtt                                                                 >  1:434557/1‑52 (MQ=255)
                     gATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGt                                                                  <  1:3501227/50‑1 (MQ=255)
                     gATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGccc                                              >  1:880921/1‑70 (MQ=255)
                        cTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTcatcat                                                                     <  1:2900485/44‑1 (MQ=255)
                         tGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCt                                                       <  1:2308331/57‑1 (MQ=255)
                           aGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGa                                                     >  1:2777919/1‑57 (MQ=255)
                           aGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTGAAAGGCCCCAtttt                                        <  1:2007189/70‑1 (MQ=255)
                           aGCCGTTTTATCCTATAATCGCGTTCAATCAATTTCATCATTGTTTGATGGGGCTGAAAGGCCCCAttt                                         <  1:188169/69‑1 (MQ=255)
                             ccGTTTTATCCTATAATCGCGTTCAATCATTTtcatc                                                                       <  1:164558/37‑1 (MQ=255)
                                ttttATCCTATAATCGCGTTCAATCATTTtcatca                                                                      <  1:2059367/35‑1 (MQ=255)
                                ttttATCCTATAATCGCGTTCAATCATTTTCATCATTGttt                                                                <  1:1086917/41‑1 (MQ=255)
                                ttttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAgg                                                 <  1:572518/56‑1 (MQ=255)
                                  ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                              >  1:2350633/1‑41 (MQ=255)
                                  ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                              >  1:1331317/1‑41 (MQ=255)
                                  ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTg                                   >  1:4172576/1‑68 (MQ=255)
                                  ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgcg                                >  1:4052078/1‑71 (MQ=255)
                                  ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgcg                                >  1:3514465/1‑71 (MQ=255)
                                  ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgcg                                >  1:897329/1‑71 (MQ=255)
                                  ttATCCGATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgcg                                >  1:554997/1‑71 (MQ=255)
                                       cTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGa                                                     <  1:2177963/45‑1 (MQ=255)
                                        tataATCGCGTTCAATCATTTTCATCATTGTTTGa                                                              >  1:2806529/1‑35 (MQ=255)
                                          taATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTg                                                      <  1:3255290/41‑1 (MQ=255)
                                          taATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTg                                                      <  1:3680459/41‑1 (MQ=255)
                                          taATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGtattat                        >  1:3572236/1‑71 (MQ=255)
                                             tCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGAc                     >  1:1518040/1‑71 (MQ=255)
                                               gcgTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGAc                     >  1:2990999/1‑69 (MQ=255)
                                                cgTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGa                  <  1:3484980/71‑1 (MQ=255)
                                                  ttCAATCATTTTCATCATTGTTTGATGGGGCTGaa                                                    <  1:3861214/35‑1 (MQ=255)
                                                    cAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCAttttt                                       <  1:680663/46‑1 (MQ=255)
                                                    cAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa              >  1:1390418/1‑71 (MQ=255)
                                                    cAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa              >  1:58832/1‑71 (MQ=255)
                                                     aaTCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa              >  1:2036275/1‑70 (MQ=255)
                                                     aaTCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa              <  1:2965603/70‑1 (MQ=255)
                                                       tCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGtatta                         >  1:3105436/1‑57 (MQ=255)
                                                       tCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGa                      <  1:4175133/60‑1 (MQ=255)
                                                       tCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAAcaacaa           <  1:2501052/71‑1 (MQ=255)
                                                         aTTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACaaa          <  1:2345647/70‑1 (MQ=255)
                                                         aTTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACaaa          <  1:1806381/70‑1 (MQ=255)
                                                         aTTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAAt         <  1:2253187/71‑1 (MQ=255)
                                                         aTTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAAt         <  1:2747931/71‑1 (MQ=255)
                                                         aTTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAAt         <  1:664353/71‑1 (MQ=255)
                                                          ttttCATCATTGTTTGATGGGGCGGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAAtt        <  1:1538468/71‑1 (MQ=255)
                                                            ttCATCATTGTTTGATGGGGCTGAAAGGCCCCAtttt                                        <  1:3810744/37‑1 (MQ=255)
                                                               atcatTGTTTGATGGGGCTGAAAGGCCTCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTAGcc    >  1:2018330/1‑70 (MQ=255)
                                                               atcatTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTAGCCg   <  1:1890408/71‑1 (MQ=255)
                                                                tcatTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTAGCCGa  <  1:381101/71‑1 (MQ=255)
                                                                |                                                                      
GATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTAGCCGA  >  W3110S.gb/3516221‑3516355

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: