Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,334,786 Δ1 bp 5.8% coding (23/279 nt) sfsB → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,334,7860G. 5.8% 217.5 / 4.4 69coding (23/279 nt)sfsBDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base . (2/2);  total (37/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.16e-01

GAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAAT  >  W3110S.gb/3334717‑3334856
                                                                     |                                                                      
gaagaaGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTgg                                                                       <  1:1946325/71‑1 (MQ=255)
     agTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATcc                                                                  >  1:100943/1‑71 (MQ=255)
     agTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATcc                                                                  >  1:902793/1‑71 (MQ=255)
         ttCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGc                                                               >  1:1304623/1‑70 (MQ=255)
         ttCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCt                                                              <  1:3192993/71‑1 (MQ=255)
         ttCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCt                                                              <  1:3249326/71‑1 (MQ=255)
             gaatgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAca                                                          >  1:1544893/1‑71 (MQ=255)
             gaatgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAca                                                          >  1:3619621/1‑71 (MQ=255)
              aatgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcat                                                         >  1:2589230/1‑71 (MQ=255)
              aatgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAc                                                           <  1:2589077/69‑1 (MQ=255)
               atgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatc                                                        <  1:387174/71‑1 (MQ=255)
               atgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatc                                                        <  1:2121511/71‑1 (MQ=255)
                tgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatca                                                       >  1:2670545/1‑71 (MQ=255)
                tgaaAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatca                                                       >  1:690223/1‑71 (MQ=255)
                  aaaGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCAtt                                                     <  1:580289/71‑1 (MQ=255)
                     gCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCg                                                  <  1:2893693/71‑1 (MQ=255)
                        aCTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGc                                                   <  1:570807/67‑1 (MQ=255)
                        aCTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGc                                                   <  1:589303/67‑1 (MQ=255)
                        aCTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGt                                               >  1:2150784/1‑71 (MQ=255)
                         cTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCa                                                                     >  1:1242694/1‑48 (MQ=255)
                            aCTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTgc                                           >  1:3668007/1‑71 (MQ=255)
                            aCTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCAGCt                                                              >  1:1169562/1‑52 (MQ=255)
                            aCTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGCATCCCGCTGACATCATTGCGGGttt                                             <  1:3603679/68‑1 (MQ=255)
                            aCTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGCATCCCGCTGACATCATTGCGGGttt                                             <  1:3358779/68‑1 (MQ=255)
                            aCTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGCATCCCGCTGACATCATTGCGGGTTTgc                                           >  1:1590824/1‑70 (MQ=255)
                            aCTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGCATCCCGCTGACATCATTGCGGGTTTgc                                           >  1:3677638/1‑70 (MQ=255)
                                aaGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCga                                                               >  1:1480275/1‑46 (MQ=255)
                                aaGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTgcgc                                         >  1:3567277/1‑69 (MQ=255)
                                aaGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTgcgc                                         >  1:2513810/1‑69 (MQ=255)
                                aaGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTgcgc                                         >  1:2590504/1‑69 (MQ=255)
                                 aGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTgcgc                                         >  1:3419341/1‑68 (MQ=255)
                                 aGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaag                                      <  1:3787783/71‑1 (MQ=255)
                                 aGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaag                                      <  1:3783241/71‑1 (MQ=255)
                                  ggAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaaga                                     <  1:3677035/71‑1 (MQ=255)
                                   gAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaagaa                                    <  1:985296/71‑1 (MQ=255)
                                   gAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaag                                      >  1:1501300/1‑69 (MQ=255)
                                    aGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaagaag                                   <  1:3589097/71‑1 (MQ=255)
                                    aGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaagaag                                   <  1:3387611/71‑1 (MQ=255)
                                       gAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAggg                                 <  1:453118/70‑1 (MQ=255)
                                       gAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAggg                                 >  1:5292/1‑70 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCTCTGACATCATTGCGGGTTTGCGCAAGAAGGGaa                               >  1:2462483/1‑71 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCTCTGACATCATTGCGGGTTTGCGCAAGAAGGGaa                               >  1:3753808/1‑71 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatca                                                       >  1:303149/1‑47 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatca                                                       >  1:688801/1‑47 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGAcatca                                                       >  1:315130/1‑47 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGa                                <  1:620448/70‑1 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCAGCTGAcatca                                                       >  1:790499/1‑47 (MQ=255)
                                        aGGAAAAATGGAAAGTAATTTCACTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGaa                               >  1:3032352/1‑71 (MQ=255)
                                           aaaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaa                                       >  1:2316716/1‑60 (MQ=255)
                                           aaaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaa                                       >  1:474747/1‑60 (MQ=255)
                                           aaaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaa                                       >  1:2184004/1‑60 (MQ=255)
                                            aaaaTTGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:1712915/71‑1 (MQ=255)
                                            aaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:315634/71‑1 (MQ=255)
                                            aaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:1544629/71‑1 (MQ=255)
                                            aaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:2304915/71‑1 (MQ=255)
                                            aaaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:3068802/71‑1 (MQ=255)
                                             aaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCa                          <  1:1081902/71‑1 (MQ=255)
                                             aaaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCa                          <  1:3693990/71‑1 (MQ=255)
                                              aaTGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCa                          <  1:715183/70‑1 (MQ=255)
                                                      gTAATTTCATTGACTGGCATCCCGCTGACATCATTg                                                    >  1:955212/1‑36 (MQ=255)
                                                       tAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:3133860/60‑1 (MQ=255)
                                                       tAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:2160677/60‑1 (MQ=255)
                                                       tAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTc                           <  1:3732092/60‑1 (MQ=255)
                                                        aaTTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCaagaag                                   >  1:1980720/1‑51 (MQ=255)
                                                        aaTTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGcggcgg                >  1:2900285/1‑70 (MQ=255)
                                                            tCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAAtct           >  1:3726175/1‑71 (MQ=255)
                                                            tCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAAtc            >  1:3687291/1‑70 (MQ=255)
                                                                 gACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAgg                                  >  1:1754046/1‑43 (MQ=255)
                                                                    tGGCATCCCGCTGACATCATTGCGGGTTTGCTCAAGAAgg                                  <  1:1836438/40‑1 (MQ=255)
                                                                    tGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAgg                                  <  1:3316977/40‑1 (MQ=255)
                                                                     ggCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAAt  >  1:2541301/1‑71 (MQ=255)
                                                                     ggCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAAt  >  1:1325601/1‑71 (MQ=255)
                                                                     |                                                                      
GAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAAGGGAACGTCAATGGCGGCGGAATCTCGCAGAAAT  >  W3110S.gb/3334717‑3334856

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: