Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,090,238 A→G 8.3% S14S (AGT→AGC ycdR ← predicted enzyme associated with biofilm formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,090,2380AG 8.3% 194.5 / 8.6 72S14S (AGT→AGCycdRpredicted enzyme associated with biofilm formation
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base G (4/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

AGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCA  >  W3110S.gb/1090169‑1090304
                                                                     |                                                                  
aGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATAc                                                                   >  1:6487409/1‑71 (MQ=255)
    tCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCAc                                                               >  1:7480810/1‑71 (MQ=255)
    tCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCAc                                                               >  1:747065/1‑71 (MQ=255)
     cgcgATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCAc                                                               >  1:657723/1‑70 (MQ=255)
         aTCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCa                                                          >  1:1216895/1‑71 (MQ=255)
          tCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGCGAGCATAATTATACTCACCAGCAt                                                         <  1:7504601/71‑1 (MQ=255)
               tgGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAg                                                      <  1:7820088/69‑1 (MQ=255)
               tgGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAg                                                      >  1:1292783/1‑69 (MQ=255)
                   ggTATAAATGATGTTCTTGACTGGCTATTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGata                                                <  1:4385127/71‑1 (MQ=255)
                   ggTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGata                                                <  1:3270076/71‑1 (MQ=255)
                    gTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGga                                                    >  1:3843308/1‑66 (MQ=255)
                     tataAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCAc                                                               >  1:281828/1‑54 (MQ=255)
                       taAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGTATAATTATGCTCAcc                                                              >  1:1643600/1‑53 (MQ=37)
                       taAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGTATAATTATGCTCAcc                                                              >  1:490244/1‑53 (MQ=37)
                       taAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGTATAATTATGCTCAcc                                                              >  1:4160044/1‑53 (MQ=37)
                       taAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCa                                                       <  1:6863468/60‑1 (MQ=255)
                       taAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCa                                                       <  1:6982343/60‑1 (MQ=255)
                       taAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGAtattt                                             <  1:2571689/70‑1 (MQ=255)
                             atgTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCa                                      <  1:1071019/71‑1 (MQ=255)
                               gTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCAt                                     >  1:4115235/1‑70 (MQ=255)
                                 tCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTAc                                  >  1:1742910/1‑71 (MQ=255)
                                  cTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGgaga                                                  >  1:3490085/1‑54 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATc                                                        <  1:1738666/45‑1 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATatttattt                                         >  1:1710273/1‑60 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATatttattt                                         >  1:171668/1‑60 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATatttattt                                         >  1:6640012/1‑60 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGt                                <  1:5602353/69‑1 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTaa                              >  1:3195169/1‑71 (MQ=255)
                                     gACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTa                               <  1:6356448/70‑1 (MQ=255)
                                      aCTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAAc                             <  1:3048643/71‑1 (MQ=255)
                                          gCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACata                          <  1:2708432/70‑1 (MQ=255)
                                            tAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATAtt                        >  1:3583862/1‑70 (MQ=255)
                                            tAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATAtt                        >  1:5087123/1‑70 (MQ=255)
                                            tAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATAtt                        >  1:4111926/1‑70 (MQ=255)
                                              aTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTAt                     >  1:6120156/1‑71 (MQ=255)
                                              aTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTAt                     >  1:2265542/1‑71 (MQ=255)
                                              aTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTAt                     <  1:3963970/71‑1 (MQ=255)
                                              aTGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTAt                     <  1:3954628/71‑1 (MQ=255)
                                                gCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAAc                             >  1:6733245/1‑61 (MQ=255)
                                                gCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAAc                             >  1:5964540/1‑61 (MQ=255)
                                                  aCGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCt                  >  1:3525962/1‑70 (MQ=255)
                                                  aCGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCt                  >  1:4937359/1‑70 (MQ=255)
                                                   cgcgGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATAttt                       <  1:3706814/64‑1 (MQ=255)
                                                    gcgGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTa                      <  1:7524357/64‑1 (MQ=255)
                                                    gcgGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTa                >  1:6060111/1‑70 (MQ=255)
                                                      ggTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTaaaa             <  1:4302577/71‑1 (MQ=255)
                                                       gTGAGTATAATTATGCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTaaaa             >  1:3321985/1‑70 (MQ=255)
                                                       gTGAGTATAATTATGCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAAt            <  1:7792517/71‑1 (MQ=255)
                                                       gTGAGTATAATTATGCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAAt            <  1:7801709/71‑1 (MQ=255)
                                                       gTGAGCATAATTATACTCACCAGCATCAGGATATATTTATTTCCATTACGTAACATATTTATCCTTAAAAt            >  1:3250917/1‑71 (MQ=255)
                                                       gTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTcc                                       <  1:5331052/44‑1 (MQ=255)
                                                       gTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAAt            <  1:7138151/71‑1 (MQ=255)
                                                         gAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCAt                                     >  1:1263635/1‑44 (MQ=255)
                                                         gAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCt          >  1:5842865/1‑71 (MQ=255)
                                                          aGCATAATTATACTCACCAGCATCAGGAGAtatttat                                           >  1:787330/1‑37 (MQ=255)
                                                          aGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCt          >  1:5805344/1‑70 (MQ=255)
                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAAc                             >  1:3920633/1‑50 (MQ=255)
                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAAc                             >  1:2729702/1‑50 (MQ=255)
                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGa        <  1:6063648/71‑1 (MQ=255)
                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGa        <  1:4204643/71‑1 (MQ=255)
                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGa        <  1:3985412/71‑1 (MQ=255)
                                                           gCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGa        <  1:1443013/71‑1 (MQ=255)
                                                            cATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGaa       >  1:5771828/1‑71 (MQ=255)
                                                              tAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATg     >  1:963090/1‑71 (MQ=255)
                                                              tAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATg     >  1:5663366/1‑71 (MQ=255)
                                                              tAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATg     >  1:3227918/1‑71 (MQ=255)
                                                                 ttATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTaaa              <  1:7074805/59‑1 (MQ=255)
                                                                 ttATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTaaa              <  1:1830520/59‑1 (MQ=255)
                                                                 ttATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa  <  1:6057679/71‑1 (MQ=255)
                                                                 ttATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa  <  1:5384351/71‑1 (MQ=255)
                                                                 ttATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa  <  1:454058/71‑1 (MQ=255)
                                                                 ttATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa  <  1:802671/71‑1 (MQ=255)
                                                                  tataCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa  >  1:6419873/1‑70 (MQ=255)
                                                                  tataCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCa  >  1:1159295/1‑70 (MQ=255)
                                                                    taCTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGa        <  1:3713974/62‑1 (MQ=255)
                                                                     |                                                                  
AGATTCGCGATCCTGTGGCGGTATAAATGATGTTCTTGACTGGCTAATGCACGCGGTGAGCATAATTATACTCACCAGCATCAGGAGATATTTATTTCCATTACGTAACATATTTATCCTTAAAATCTGAATGTCA  >  W3110S.gb/1090169‑1090304

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: