Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 923,181 C→T 12.5% intergenic (‑169/‑154) cspD ← / → clpS cold shock protein homolog/regulatory protein for ClpA substrate specificity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb923,1810CT12.5% 218.3 / 19.2 80intergenic (‑169/‑154)cspD/clpScold shock protein homolog/regulatory protein for ClpA substrate specificity
Reads supporting (aligned to +/- strand):  ref base C (29/39);  new base T (7/3);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTA  >  W3110S.gb/923113‑923242
                                                                    |                                                             
gagaATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAg                                                             <  1:1620880/71‑1 (MQ=255)
  gaATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGtt                                                           <  1:1219831/71‑1 (MQ=255)
  gaATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGtt                                                           <  1:865346/71‑1 (MQ=255)
     ttAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAAc                                                        <  1:3569960/71‑1 (MQ=255)
             aGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGa                                                >  1:3431847/1‑71 (MQ=255)
             aGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGa                                                <  1:2182985/71‑1 (MQ=255)
               ccGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGTAGTTAACGCTATTGa                                                >  1:4139358/1‑69 (MQ=255)
                 gTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGTAGTTAACGCTATTGACAg                                             >  1:4281601/1‑70 (MQ=255)
                  tCAAGCTTTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAgg                                            <  1:1925908/70‑1 (MQ=255)
                  tCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGc                                                      <  1:2397905/60‑1 (MQ=255)
                  tCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGa                                           <  1:1671367/71‑1 (MQ=255)
                   cAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTa                                                    >  1:3353623/1‑61 (MQ=255)
                    aaGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGGAGTTAACGCTATTGACAGGAAt                                         >  1:3436666/1‑71 (MQ=255)
                      tcATTTTACAGGGGACAAGTGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGa                                           <  1:2724188/66‑1 (MQ=255)
                      gCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGTAGTTAACGCTATTGACAGGAAtgt                                       <  1:1798161/71‑1 (MQ=255)
                      gCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGtt                                                           <  1:4322343/51‑1 (MQ=255)
                      gCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAtgt                                       <  1:4019185/71‑1 (MQ=255)
                         tttGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGTAGtt                                                           <  1:355784/48‑1 (MQ=255)
                         tttGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGtt                                                           <  1:739801/48‑1 (MQ=255)
                          ttGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGtt                                                           <  1:4148766/47‑1 (MQ=255)
                           tGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTaa                                                         >  1:1492474/1‑48 (MQ=255)
                           tGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGa                                     >  1:3342933/1‑68 (MQ=255)
                           tGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACa                                   <  1:2023527/70‑1 (MQ=255)
                                 gggACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCg                            <  1:4628889/71‑1 (MQ=255)
                                  ggACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTaa                                                         >  1:2926475/1‑41 (MQ=255)
                                  ggACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGc                           <  1:4219735/71‑1 (MQ=255)
                                   gACAAGGGGCAGGTATTAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCt                          <  1:2790959/71‑1 (MQ=255)
                                    aCAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGt                                                            <  1:105475/36‑1 (MQ=255)
                                    aCAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCt                          >  1:1267678/1‑70 (MQ=255)
                                    aCAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTg                         <  1:723903/71‑1 (MQ=255)
                                      aaGGGGCAGGTATGAATCAAAAATTTGAAGTAGTTaa                                                         >  1:1101697/1‑37 (MQ=255)
                                      aaGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGAc                                               >  1:1323067/1‑47 (MQ=255)
                                        ggggCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAt                                         >  1:1823563/1‑51 (MQ=255)
                                         gggCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGcc                    <  1:2369995/71‑1 (MQ=255)
                                         gggCAGATATGAATCAAAAATTTGAAGTAGTTAACGCTATTGACAGGAATGTGACAGa                                 >  1:5904/1‑58 (MQ=255)
                                           gCAGGTATGAATCAAAAATTTGAAGCAGTTAACGc                                                      <  1:401144/35‑1 (MQ=255)
                                           gCAGGTATGAATCAAAAATTTGAAGCAGTTAACGc                                                      <  1:3272225/35‑1 (MQ=255)
                                           gCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTa                                                    <  1:3641838/37‑1 (MQ=255)
                                           gCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCaa                  <  1:1844448/71‑1 (MQ=255)
                                            cAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGcc                    <  1:992118/68‑1 (MQ=255)
                                              ggTATGAATCAAAAATTTGAAGTAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               <  1:4146950/71‑1 (MQ=255)
                                              ggTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               >  1:2500861/1‑71 (MQ=255)
                                               gTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGa                                                <  1:3222781/37‑1 (MQ=255)
                                               gTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGAt                                <  1:2607814/53‑1 (MQ=255)
                                               gTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGt                              <  1:4378642/55‑1 (MQ=255)
                                                 attaaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               >  1:4152487/4‑68 (MQ=255)
                                                 aTGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAg                                             >  1:4498034/1‑38 (MQ=255)
                                                 aTGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGa                                           >  1:12698/1‑40 (MQ=255)
                                                 aTGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGa                                 >  1:1764577/1‑50 (MQ=255)
                                                 aTGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               >  1:3861127/1‑68 (MQ=255)
                                                 aTGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               >  1:908796/1‑68 (MQ=255)
                                                  tGAATCAAAAATTTGAAGTAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               >  1:1590167/1‑67 (MQ=255)
                                                  tGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAg                                             <  1:2242117/37‑1 (MQ=255)
                                                  tGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAg                                             <  1:1485268/37‑1 (MQ=255)
                                                  tGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAtgt                                       <  1:2737845/43‑1 (MQ=255)
                                                  tGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCt                          <  1:394450/56‑1 (MQ=255)
                                                  tGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGa               >  1:711894/1‑67 (MQ=255)
                                                   gAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGa                                           <  1:2876190/38‑1 (MQ=255)
                                                   gAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgat          >  1:183152/1‑71 (MQ=255)
                                                   gAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgat          >  1:1702948/1‑71 (MQ=255)
                                                    aaTCAAAAATTTGAAGTAGTTAACGCTATTGACAGGAATGTGACAGa                                 >  1:1718934/1‑47 (MQ=255)
                                                    aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGaa                                          >  1:906680/1‑38 (MQ=255)
                                                    aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGa                                 >  1:3743927/1‑47 (MQ=255)
                                                    aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgat          >  1:3547038/1‑70 (MQ=255)
                                                    aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgat          >  1:3015422/1‑70 (MQ=255)
                                                    aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgat          <  1:2783506/70‑1 (MQ=255)
                                                     aTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatga        <  1:2672409/71‑1 (MQ=255)
                                                     aTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatg         <  1:1635311/70‑1 (MQ=255)
                                                      tCAAAAATTTGAAGTAGTTAAAGCTATTGACAGGAAtg                                        >  1:1813017/1‑38 (MQ=38)
                                                      tCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAgatgat       <  1:996394/71‑1 (MQ=255)
                                                        aaaaaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACa                                   <  1:951336/41‑1 (MQ=255)
                                                        aaaaaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATa      >  1:82944/1‑70 (MQ=255)
                                                        aaaaaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAg     >  1:2624707/1‑71 (MQ=255)
                                                        aaaaaTTTGAAGAAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCaa                  >  1:255139/1‑58 (MQ=255)
                                                           aaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCa                   >  1:4487717/1‑54 (MQ=255)
                                                           aaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCa                   >  1:3866684/1‑54 (MQ=255)
                                                           aaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCa                   >  1:3512818/1‑54 (MQ=255)
                                                           aaTTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTta  >  1:695802/1‑71 (MQ=255)
                                                               tGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGa                        >  1:4438046/1‑45 (MQ=255)
                                                                gAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCt                          <  1:753257/42‑1 (MQ=255)
                                                                    |                                                             
GAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTA  >  W3110S.gb/923113‑923242

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: