Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,291,258 C→T 23.4% G326D (GGC→GAC)  yejO ← predicted autotransporter outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,291,2580CT23.4% 152.6 / 43.0 77G326D (GGC→GAC) yejOpredicted autotransporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (6/12);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGA  >  W3110S.gb/2291188‑2291327
                                                                      |                                                                     
aTAATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGc                                                                       >  1:4045432/1‑71 (MQ=255)
  aaTTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGc                                                                       >  1:4328899/1‑69 (MQ=255)
   atttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                    >  1:3256631/1‑71 (MQ=255)
   atttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                    >  1:2587983/1‑71 (MQ=255)
     ttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGcc                                                                  >  1:2045718/1‑71 (MQ=255)
        tttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCa                                                                <  1:3576583/70‑1 (MQ=255)
        tttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCa                                                                <  1:4601644/70‑1 (MQ=255)
              gCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGAc                                                             <  1:1453312/67‑1 (MQ=255)
               cTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAg                                                        <  1:709/71‑1 (MQ=255)
                 tCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAgcg                                                      >  1:408303/1‑71 (MQ=255)
                 tCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAgcg                                                      >  1:2160692/1‑71 (MQ=255)
                  cATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAgcgc                                                     >  1:1201610/1‑71 (MQ=255)
                    tGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCa                                                                <  1:1346626/58‑1 (MQ=255)
                      aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGcc                                                 <  1:4313123/71‑1 (MQ=255)
                      aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGc                                                  >  1:4081959/1‑70 (MQ=255)
                      aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGc                                                  <  1:3446055/70‑1 (MQ=255)
                      aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGc                                                  <  1:1879775/70‑1 (MQ=255)
                        cTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGtt                                                          <  1:1434345/60‑1 (MQ=255)
                        cTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAt                                               <  1:4151299/71‑1 (MQ=255)
                         tCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAg                                                                    <  1:1668066/49‑1 (MQ=255)
                         tCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAg                                                                    <  1:556394/49‑1 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGc                                                                       >  1:2627038/1‑44 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:2593419/1‑71 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:2942508/1‑71 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGa                                                              >  1:2510234/1‑53 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            <  1:3861234/71‑1 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:2765253/1‑71 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                              <  1:1272775/69‑1 (MQ=255)
                           aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                              <  1:528165/69‑1 (MQ=255)
                               aTGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGa                                         <  1:4218956/70‑1 (MQ=255)
                                  ccATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGa                                                              >  1:1801246/1‑46 (MQ=255)
                                  ccATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAAc                                     >  1:1495918/1‑71 (MQ=255)
                                   catcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          <  1:2670103/65‑1 (MQ=255)
                                    atcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGcc                                                                  <  1:2793040/40‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGtt                                                          <  1:3025314/47‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGGGCGGCCATTGTTGATAACCat                                  <  1:1645828/71‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGcc                                                                  <  1:4237263/39‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTa                                                         <  1:2956005/48‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAgc                                                       <  1:4296432/50‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttg                                             <  1:4332204/60‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttg                                             <  1:98545/60‑1 (MQ=255)
                                     tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCat                                  <  1:1834119/71‑1 (MQ=255)
                                      catTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATaa                                      >  1:4010711/1‑66 (MQ=255)
                                      catTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATaa                                      >  1:3624197/1‑66 (MQ=255)
                                       atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCa                                                                >  1:2911856/1‑39 (MQ=255)
                                       atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCata                                 >  1:895429/1‑70 (MQ=255)
                                        ttCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAg                                                        >  1:2846185/1‑46 (MQ=255)
                                        ttCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtt                               <  1:3472151/71‑1 (MQ=255)
                                          cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAg                                                        >  1:1730185/1‑44 (MQ=255)
                                          cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAg                                                        >  1:4213168/1‑44 (MQ=255)
                                          cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtttt                             >  1:3507751/1‑71 (MQ=255)
                                          cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtttt                             >  1:2316861/1‑71 (MQ=255)
                                           ccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTc                            <  1:1861526/71‑1 (MQ=255)
                                             tctGACTGAAACGTTAACCATTGTGCCAGCCCAGACg                                                            <  1:4331634/37‑1 (MQ=255)
                                              ctgactgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:2705237/1‑69 (MQ=255)
                                               tgactgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:1816351/1‑68 (MQ=255)
                                               tgactgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:2396941/1‑68 (MQ=255)
                                                 actgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:4490949/68‑1 (MQ=255)
                                                 actgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:236632/68‑1 (MQ=255)
                                                 actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:2640029/68‑1 (MQ=255)
                                                 actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:3394212/68‑1 (MQ=255)
                                                  ctgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGcgcg                                                    <  1:1041205/40‑1 (MQ=255)
                                                  ctgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATg                     <  1:3579630/71‑1 (MQ=255)
                                                   tgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          <  1:602510/49‑1 (MQ=255)
                                                    gaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGc                                                  >  1:3054729/1‑40 (MQ=255)
                                                    gaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATg                     >  1:1668584/1‑69 (MQ=255)
                                                    gaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATg                     >  1:3500197/1‑69 (MQ=255)
                                                        cGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCtttt               <  1:3223245/71‑1 (MQ=255)
                                                          ttAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCtttt               <  1:4180700/69‑1 (MQ=255)
                                                          ttAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCtttt               >  1:2509442/1‑69 (MQ=255)
                                                            aaCCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATa                                       >  1:4234835/1‑43 (MQ=255)
                                                            aaCCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATa                                       >  1:1595324/1‑43 (MQ=255)
                                                            aaCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           <  1:4607376/55‑1 (MQ=255)
                                                             aCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAgg           <  1:4478488/70‑1 (MQ=255)
                                                             aCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAgg           >  1:148214/1‑70 (MQ=255)
                                                             aCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGa          <  1:3257122/71‑1 (MQ=255)
                                                               cATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          >  1:3485328/1‑37 (MQ=255)
                                                                aTTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATa       <  1:3807033/71‑1 (MQ=255)
                                                                   gtgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         >  1:1556226/1‑50 (MQ=255)
                                                                    tgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGt    <  1:1006008/70‑1 (MQ=255)
                                                                      ccAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGa  <  1:217784/70‑1 (MQ=255)
                                                                      |                                                                     
ATAATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGA  >  W3110S.gb/2291188‑2291327

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: