Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,119,901 A→T 6.5% T63S (ACA→TCA) ‡ hofQ → predicted fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,119,9010AT 6.5% 433.7 / 17.9 153T63S (ACA→TCA) ‡hofQpredicted fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base A (69/74);  new base T (4/6);  total (73/80)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 163 total aligned reads displayed.

TGGCTGAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGAAGG  >  W3110S.gb/4119833‑4119971
                                                                    |                                                                      
  gctgAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAg                                                                     <  1:1161132/70‑1 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:5784882/70‑1 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATg                                                                  <  1:3941294/69‑1 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                <  1:6715431/71‑1 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:7273056/1‑71 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:5614435/1‑71 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:1874807/1‑71 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:1864862/1‑71 (MQ=255)
      aaCAGGAGAAGTTGAACCTGGGCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:4221674/70‑1 (MQ=255)
        cAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTTTCGTTACATCTAACAGATGTTcc                                                              >  1:3838553/1‑71 (MQ=255)
        cAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTc                                                               >  1:6606830/1‑70 (MQ=255)
          ggAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            >  1:1028221/1‑71 (MQ=255)
          ggAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:4253403/71‑1 (MQ=255)
            agaAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:2891174/1‑65 (MQ=255)
               aGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGtaa                                                               >  1:6110285/1‑61 (MQ=38)
               aGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGtaa                                                               >  1:397929/1‑61 (MQ=38)
               aGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGtaa                                                               >  1:2260818/1‑61 (MQ=38)
                gTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGt                                                                 >  1:5680683/1‑60 (MQ=255)
                 ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:6230181/64‑1 (MQ=255)
                 ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:884058/64‑1 (MQ=255)
                 ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCa                                                     <  1:4918012/71‑1 (MQ=255)
                 ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCa                                                     <  1:2603603/71‑1 (MQ=255)
                  tGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAg                                                    >  1:7923536/1‑71 (MQ=255)
                   gAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTcc                                                              <  1:833801/60‑1 (MQ=255)
                    aaCCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAgg                                                   <  1:5215682/70‑1 (MQ=255)
                      ccTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAg                                                    <  1:7860940/67‑1 (MQ=255)
                       cTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGa                                                                    <  1:3955327/50‑1 (MQ=255)
                       cTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGa                                                                    <  1:1276064/50‑1 (MQ=255)
                        tGGTCGTGTCGCCAGACGTCAGCGGTACTGTGTCGTTACATCTAACAGa                                                                    <  1:119330/49‑1 (MQ=255)
                        tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:7813468/52‑1 (MQ=37)
                        tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:4596050/52‑1 (MQ=37)
                        tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAgg                                                   >  1:1157202/1‑66 (MQ=255)
                        tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:1923321/1‑68 (MQ=255)
                           tCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:5075195/1‑70 (MQ=255)
                           tCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:7491781/1‑70 (MQ=255)
                              tgtCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:2598720/51‑1 (MQ=255)
                              tgtCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACt                                        <  1:7566058/71‑1 (MQ=255)
                                tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGc                                                      <  1:6348585/55‑1 (MQ=255)
                                tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:2759606/1‑65 (MQ=255)
                                tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:4606978/1‑65 (MQ=255)
                                  gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                <  1:5299571/59‑1 (MQ=255)
                                  gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                <  1:4827381/59‑1 (MQ=255)
                                  gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                <  1:6789828/59‑1 (MQ=255)
                                  gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACaaa                                          <  1:3730596/65‑1 (MQ=255)
                                  gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTa                                     <  1:1614352/70‑1 (MQ=255)
                                    cAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTg                                  <  1:5089195/71‑1 (MQ=255)
                                     aGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTg                                  >  1:6196928/1‑70 (MQ=255)
                                      gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:830633/71‑1 (MQ=255)
                                      gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:5783512/71‑1 (MQ=255)
                                      gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:5854511/71‑1 (MQ=255)
                                      gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:4401420/71‑1 (MQ=255)
                                      gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:4063159/71‑1 (MQ=255)
                                      gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:2530827/71‑1 (MQ=255)
                                       aCGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:5904929/1‑53 (MQ=255)
                                       aCGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaa                               <  1:6460168/71‑1 (MQ=255)
                                        cGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              <  1:2362923/71‑1 (MQ=255)
                                        cGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              <  1:26902/71‑1 (MQ=255)
                                            aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTa                                              <  1:2894625/51‑1 (MQ=255)
                                            aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACaaa                                          >  1:1434582/1‑55 (MQ=255)
                                            aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACaaa                                          >  1:6526558/1‑55 (MQ=255)
                                            aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGt                                      <  1:4065147/59‑1 (MQ=255)
                                            aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGt                                   >  1:2811964/1‑62 (MQ=255)
                                            aGCGCTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACaaa                                          >  1:4167514/1‑55 (MQ=255)
                                                ggaCGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCCCTACAGTTTGTAGTgc                                 >  1:2759731/3‑59 (MQ=255)
                                                gTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAGTTTGTAGTgc                                 >  1:2045013/1‑59 (MQ=38)
                                                gTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGa                      >  1:5477323/1‑71 (MQ=255)
                                                 tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:6546734/1‑43 (MQ=255)
                                                 tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:3522496/1‑43 (MQ=255)
                                                 tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:7884690/1‑43 (MQ=255)
                                                 tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGAc                     >  1:1677194/1‑71 (MQ=255)
                                                  aCGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACt                    >  1:5560762/1‑71 (MQ=255)
                                                  aCGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACt                    >  1:2204757/1‑71 (MQ=255)
                                                   cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                >  1:5125792/1‑42 (MQ=255)
                                                   cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:2058481/58‑1 (MQ=255)
                                                    ggTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGa                  >  1:1720993/1‑71 (MQ=255)
                                                      tgtCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:5326709/48‑1 (MQ=255)
                                                       gtCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACt                    <  1:6823061/66‑1 (MQ=255)
                                                        tCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACt                                        <  1:5644004/45‑1 (MQ=255)
                                                         cGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGa                  >  1:6986587/1‑66 (MQ=255)
                                                         cGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGa                  >  1:7296561/1‑66 (MQ=255)
                                                          gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcgc            >  1:7293633/1‑71 (MQ=255)
                                                          gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcg             <  1:3702694/70‑1 (MQ=255)
                                                          gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcg             <  1:2598065/70‑1 (MQ=255)
                                                           ttACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:1238199/43‑1 (MQ=255)
                                                           ttACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:7683829/43‑1 (MQ=255)
                                                           ttACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:3568588/43‑1 (MQ=255)
                                                             aCATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCgg          >  1:7548555/1‑70 (MQ=255)
                                                             aCATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCgg          >  1:7801026/1‑70 (MQ=255)
                                                             aCATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCgg          >  1:1615460/1‑70 (MQ=255)
                                                              cATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              >  1:1062207/1‑49 (MQ=255)
                                                              cATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              >  1:1901965/1‑49 (MQ=255)
                                                              cATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCa        <  1:7123507/71‑1 (MQ=255)
                                                               aTCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAg       <  1:2610923/71‑1 (MQ=255)
                                                                 cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              <  1:471168/46‑1 (MQ=255)
                                                                 cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAg                             <  1:4261020/47‑1 (MQ=255)
                                                                 cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAgg      >  1:4145598/1‑70 (MQ=255)
                                                                 cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAgg      >  1:1394960/1‑70 (MQ=255)
                                                                   aaCAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGc         <  1:1381532/65‑1 (MQ=255)
                                                                    aCAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAgg      <  1:6306518/67‑1 (MQ=255)
                                                                    aCAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGaa    <  1:31473/69‑1 (MQ=255)
                                                                    |                                                                      
TGGCTGAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGAAGG  >  W3110S.gb/4119833‑4119971

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: