Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,291,258 C→T 22.1% 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,2580CT22.1% 179.4 / 45.8 86G326D (GGC→GAC) yejOpredicted autotransporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base T (14/5);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 3.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

ATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGA  >  W3110S.gb/2291191‑2291327
                                                                   |                                                                     
atttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAg                                                                    >  1:1510954/1‑71 (MQ=255)
atttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                    <  1:3836804/71‑1 (MQ=255)
atttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                    <  1:688525/71‑1 (MQ=255)
atttatttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                    >  1:1406722/1‑71 (MQ=255)
atttatttGCGGCTTCATGACCTCAAGAATGACATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                    >  1:390152/1‑71 (MQ=255)
          ggCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGtt                                                          <  1:873679/71‑1 (MQ=255)
             ttCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGtt                                                          <  1:1961444/68‑1 (MQ=255)
             ttCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGtt                                                          <  1:3647569/68‑1 (MQ=255)
              tCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAgc                                                       >  1:1166176/1‑70 (MQ=255)
                   aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCgg                                                   >  1:425819/1‑69 (MQ=255)
                   aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCgg                                                   >  1:240296/1‑69 (MQ=255)
                   aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGc                                                  <  1:639565/70‑1 (MQ=255)
                    ccTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGcc                                                 <  1:2692357/70‑1 (MQ=255)
                       cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGcc                                                                      >  1:1368258/1‑46 (MQ=255)
                       cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttg                                             >  1:597159/1‑71 (MQ=255)
                       cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttg                                             >  1:3845873/1‑71 (MQ=255)
                       cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttg                                             >  1:2723746/1‑71 (MQ=255)
                       cAAGAATGCCATCATTCCCTCTGACAGAAACGTTAACCATTGTGcc                                                                      >  1:2163773/1‑46 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:3577867/1‑71 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:3556835/1‑71 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:4054018/1‑71 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttg                                             <  1:866564/70‑1 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:1914550/1‑71 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:2063970/1‑71 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:3119246/1‑71 (MQ=255)
                        aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                              <  1:906101/69‑1 (MQ=255)
                         aGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgtt                                           <  1:3249729/71‑1 (MQ=255)
                         aGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgtt                                           <  1:978114/71‑1 (MQ=255)
                         aGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttgtt                                           <  1:4002852/71‑1 (MQ=255)
                          gAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          <  1:3211482/71‑1 (MQ=255)
                          gAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          <  1:3405355/71‑1 (MQ=255)
                          gAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          <  1:466798/71‑1 (MQ=255)
                          gAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgttg                                          <  1:2175019/71‑1 (MQ=255)
                            aTGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGa                                         >  1:469088/1‑70 (MQ=255)
                                  tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGAt                                        <  1:3506122/65‑1 (MQ=255)
                                  tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGAt                                        <  1:3930963/65‑1 (MQ=255)
                                  tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAAc                                     <  1:2578963/68‑1 (MQ=255)
                                  tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCat                                  <  1:2841277/71‑1 (MQ=255)
                                   catTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAt                                               >  1:2387504/1‑57 (MQ=255)
                                   catTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAt                                               >  1:3735726/1‑57 (MQ=255)
                                   catTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGAt                                        <  1:2851752/64‑1 (MQ=255)
                                   catTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCata                                 >  1:305417/1‑71 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:3344535/1‑59 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCatat                                >  1:1798640/1‑71 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCatat                                >  1:62889/1‑71 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCata                                 >  1:1233237/1‑70 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCa                                   >  1:4006200/1‑68 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCa                                   >  1:3161638/1‑68 (MQ=255)
                                    atTCCCTCTGACTGAAACGTTAACAATTGTGCCAGCCCAGACGTTAGCGCGGCCAttgt                                            >  1:107768/1‑59 (MQ=255)
                                       cccTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtttt                             >  1:1623475/1‑71 (MQ=255)
                                       cccTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtttt                             >  1:1795460/1‑71 (MQ=255)
                                       cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtttt                             >  1:494950/1‑71 (MQ=255)
                                       cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAtttt                             >  1:2807425/1‑71 (MQ=255)
                                       cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCAATTGTTGATAACCATAtttt                             >  1:706179/1‑71 (MQ=255)
                                        ccTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATaa                                      >  1:2508411/1‑61 (MQ=255)
                                        ccTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATaa                                      >  1:3306353/1‑61 (MQ=255)
                                         ctctGACTGAAACGTTAACCATTGTGCCAGCCCAGACGtt                                                          <  1:1560102/40‑1 (MQ=255)
                                           ctgactgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:1871888/1‑69 (MQ=255)
                                           ctgactgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:2161234/1‑69 (MQ=255)
                                           ctgactgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:2283477/1‑69 (MQ=255)
                                           ctgactgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:1375700/1‑69 (MQ=255)
                                           ctgactgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:4026478/1‑69 (MQ=255)
                                            tgactgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                           >  1:772565/1‑68 (MQ=255)
                                              actgaAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:3031240/68‑1 (MQ=255)
                                              actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:3027385/68‑1 (MQ=255)
                                              actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:1781568/68‑1 (MQ=255)
                                              actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:1309238/68‑1 (MQ=255)
                                              actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGAt                      <  1:836628/71‑1 (MQ=255)
                                               ctgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATg                     <  1:669344/71‑1 (MQ=255)
                                               ctgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATg                     <  1:190365/71‑1 (MQ=255)
                                                tgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                         <  1:1704810/66‑1 (MQ=255)
                                                    aCGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCttt                >  1:1442150/1‑71 (MQ=255)
                                                     cGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAgaga                       <  1:2485672/63‑1 (MQ=255)
                                                     cGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAgaga                       <  1:3036438/63‑1 (MQ=255)
                                                     cGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGc                   >  1:799919/1‑67 (MQ=255)
                                                     cGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCtttt               <  1:224552/71‑1 (MQ=255)
                                                      gTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTg              <  1:2239529/71‑1 (MQ=255)
                                                      gTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTg              <  1:93666/71‑1 (MQ=255)
                                                      gTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTg              <  1:93096/71‑1 (MQ=255)
                                                      gTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTg              <  1:847583/71‑1 (MQ=255)
                                                       ttAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                              <  1:783461/38‑1 (MQ=255)
                                                       ttAACCATTGTGCCAGACCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCtttt               >  1:1383695/1‑69 (MQ=255)
                                                          aCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAgg           <  1:3031230/70‑1 (MQ=255)
                                                           ccATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAAc                                     <  1:152542/43‑1 (MQ=255)
                                                                gtgtCAGCCCAGACGTTAGCGCGGCCATTGTTGATaa                                      >  1:178682/1‑37 (MQ=255)
                                                                gtgtCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGt    >  1:3011863/1‑71 (MQ=255)
                                                                gtgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGt    >  1:2729456/1‑71 (MQ=255)
                                                                gtgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGACGAATACGt    >  1:3922380/1‑71 (MQ=255)
                                                                 tgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGt    <  1:3316919/70‑1 (MQ=255)
                                                                 tgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTAGGAATACGTGa  <  1:4084165/71‑1 (MQ=255)
                                                                   ccAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGAt                      <  1:3857783/50‑1 (MQ=255)
                                                                   ccAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGa  <  1:1874439/70‑1 (MQ=255)
                                                                   ccAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGa  <  1:2797666/70‑1 (MQ=255)
                                                                   |                                                                     
ATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGTGA  >  W3110S.gb/2291191‑2291327

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: