Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,291,258 C→T 9.2% 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,2580CT 9.2% 159.1 / 5.5 54G326D (GGC→GAC) yejOpredicted autotransporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (2/3);  total (25/29)
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.85e-01

TAATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGT  >  W3110S.gb/2291189‑2291325
                                                                     |                                                                   
tAATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGc                                                                     >  1:725707/1‑70 (MQ=255)
       tttGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCa                                                              <  1:648162/70‑1 (MQ=255)
           cGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGt                                                         <  1:2219975/71‑1 (MQ=255)
                tCATGACCTCAAGGATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAgcg                                                    >  1:41466/1‑71 (MQ=255)
                tCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGcc                                                                >  1:1805831/1‑59 (MQ=255)
                tCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGc                                                                 >  1:2432877/1‑58 (MQ=255)
                   tGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGc                                                                     >  1:2168924/1‑51 (MQ=255)
                     aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGc                                                >  1:2032711/1‑70 (MQ=255)
                     aCCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGc                                                <  1:2311375/70‑1 (MQ=255)
                       cTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAt                                             <  1:2718050/71‑1 (MQ=255)
                       cTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAt                                             <  1:442024/71‑1 (MQ=255)
                         cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGcc                                                                >  1:1612794/1‑50 (MQ=255)
                         cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGcc                                                                >  1:313546/1‑50 (MQ=255)
                         cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttg                                           >  1:2561256/1‑71 (MQ=255)
                         cAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAttg                                           >  1:951454/1‑71 (MQ=255)
                          aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAg                                                                  >  1:412094/1‑47 (MQ=255)
                          aaGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                            <  1:982168/69‑1 (MQ=255)
                               tGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAg                                                      <  1:2222739/54‑1 (MQ=255)
                                gCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGa                                                            >  1:2115879/1‑47 (MQ=255)
                                 ccATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGt                                                         <  1:41568/49‑1 (MQ=255)
                                 ccATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGt                                                         <  1:1799191/49‑1 (MQ=255)
                                  catcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCa                                                              <  1:1092334/43‑1 (MQ=255)
                                    tcatTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGtt                                                        <  1:2114305/47‑1 (MQ=255)
                                    tcatTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGtt                                                        <  1:1715278/47‑1 (MQ=255)
                                    tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAAcc                                  <  1:2520321/69‑1 (MQ=255)
                                    tcatTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCat                                <  1:2046574/71‑1 (MQ=255)
                                    tcatTACCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCa                                 >  1:2027812/1‑70 (MQ=255)
                                      atTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCatat                              >  1:940439/1‑71 (MQ=255)
                                      atTCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCatat                              >  1:273648/1‑71 (MQ=255)
                                      atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGa                                                            >  1:870360/1‑41 (MQ=255)
                                      atTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGa                                                            >  1:2683160/1‑41 (MQ=255)
                                       ttCCCTCTGACTGAAACGTTAACCATTGTGTCAGCCCAGACGTTAGCGCGGCCAttg                                           <  1:1473950/57‑1 (MQ=255)
                                       ttCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                            <  1:681638/56‑1 (MQ=255)
                                       ttCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCAtt                                            <  1:1232020/56‑1 (MQ=255)
                                         cccTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATAttt                            >  1:464107/1‑70 (MQ=255)
                                                actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                       <  1:358649/68‑1 (MQ=255)
                                                actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                       <  1:1261987/68‑1 (MQ=255)
                                                actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                       <  1:1610278/68‑1 (MQ=255)
                                                actgaAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGAt                    >  1:1155654/1‑71 (MQ=255)
                                                       cGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCtttt             >  1:358324/1‑71 (MQ=255)
                                                        gTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTg            <  1:367615/71‑1 (MQ=255)
                                                        gTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTg            <  1:612519/71‑1 (MQ=255)
                                                         ttAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATa                                     >  1:390301/1‑45 (MQ=255)
                                                         ttAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATa                                     >  1:2128822/1‑45 (MQ=255)
                                                         ttAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCa                         <  1:838718/57‑1 (MQ=255)
                                                           aaCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                       >  1:606246/1‑57 (MQ=255)
                                                           aaCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                       >  1:1243433/1‑57 (MQ=255)
                                                           aaCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCaga                       >  1:412515/1‑57 (MQ=255)
                                                            aCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAgg         <  1:1595116/70‑1 (MQ=255)
                                                            aCCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAgg         <  1:224436/70‑1 (MQ=255)
                                                             ccATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGaa       <  1:2131132/71‑1 (MQ=255)
                                                               aTTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATa     <  1:1539336/71‑1 (MQ=255)
                                                                ttGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCa                                 <  1:1047314/42‑1 (MQ=255)
                                                                  gtgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGt  >  1:218795/1‑71 (MQ=255)
                                                                  gtgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATATTCAGAGATGGCATTTGAGGAATACg   >  1:1070686/1‑70 (MQ=255)
                                                                   tgCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGt  <  1:839154/70‑1 (MQ=255)
                                                                     |                                                                   
TAATTTATTTGCGGCTTCATGACCTCAAGAATGCCATCATTCCCTCTGACTGAAACGTTAACCATTGTGCCAGCCCAGACGTTAGCGCGGCCATTGTTGATAACCATATTTTCAGAGATGGCTTTTGAGGAATACGT  >  W3110S.gb/2291189‑2291325

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