Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,533,467 C→T 24.0% G47S (GGT→AGT)  ypeB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,533,4670CT24.0% 97.0 / 27.8 50G47S (GGT→AGT) ypeBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/20);  new base T (7/5);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTA  >  W3110S.gb/2533399‑2533535
                                                                    |                                                                    
cTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                                                                    <  1:1907231/71‑1 (MQ=255)
  tCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGc                                                                  >  1:2661795/1‑71 (MQ=255)
              ccAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCggg                                                      <  1:1768952/71‑1 (MQ=255)
                    gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGc                                                                  <  1:2837296/53‑1 (MQ=255)
                    gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGc                                                                  <  1:1926656/53‑1 (MQ=255)
                    gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTcc                                                              >  1:65453/1‑57 (MQ=255)
                    gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTcc                                                              >  1:237454/1‑57 (MQ=255)
                    gCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTc                                                 >  1:2069709/1‑70 (MQ=255)
                       tCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTc                                                               <  1:3402901/53‑1 (MQ=255)
                        cAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCtt                                            >  1:4385573/1‑71 (MQ=255)
                              cTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCggg                                      <  1:2411033/71‑1 (MQ=255)
                              cTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCg                                        >  1:3189625/1‑69 (MQ=255)
                               tCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCggg                                      >  1:890692/1‑70 (MQ=255)
                                cGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTccc                                                             >  1:1808397/1‑46 (MQ=255)
                                cGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCg                                        >  1:390321/1‑67 (MQ=255)
                                 ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                    >  1:5050538/1‑70 (MQ=255)
                                 ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCa                                    <  1:1103081/70‑1 (MQ=255)
                                 ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACTGAGATTCCCTGCCGGGAACTCATCTTTGCGGGCa                                    >  1:3692009/1‑70 (MQ=255)
                                 ggTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGACGGGAACTCATCTTTGCGGGCaa                                   <  1:3316390/71‑1 (MQ=255)
                                       aaTCAGCAGCGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACcac                              <  1:4986763/70‑1 (MQ=255)
                                        aTCAGCAGCTTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                        >  1:2979655/1‑43 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                        >  1:1564464/1‑43 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                        >  1:3718882/1‑43 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                        >  1:2831761/1‑43 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                        >  1:1827705/1‑43 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCg                                                        >  1:2541293/1‑43 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAAc                                 <  1:3252389/66‑1 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAAc                                 <  1:3258189/66‑1 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAAc                                 <  1:1914050/66‑1 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcag                            <  1:2916341/71‑1 (MQ=255)
                                        aTCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcag                            <  1:2544133/71‑1 (MQ=255)
                                           agcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcag                            >  1:2218603/1‑68 (MQ=255)
                                            gcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                        <  1:2979501/71‑1 (MQ=255)
                                            gcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                        <  1:886930/71‑1 (MQ=255)
                                            gcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                        <  1:2005145/71‑1 (MQ=255)
                                            gcagcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                        <  1:3785808/71‑1 (MQ=255)
                                             cagcGTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                        >  1:5017505/1‑70 (MQ=255)
                                               gcGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCaa                                   <  1:2223844/57‑1 (MQ=255)
                                                 gTGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCt                                             >  1:922462/1‑45 (MQ=255)
                                                 gTGATTTGCATTAGCTCACCGAGCTTCCCTGCCggg                                                      <  1:3669585/36‑1 (MQ=255)
                                                  tGATTTGCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCggg                                      >  1:327517/1‑51 (MQ=255)
                                                     tttGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCAt                                               <  1:1969232/39‑1 (MQ=255)
                                                       tGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGc                                         <  1:3117388/43‑1 (MQ=255)
                                                       tGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCgg             >  1:4923951/1‑71 (MQ=255)
                                                        gCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcagcag                         >  1:2821339/1‑58 (MQ=255)
                                                        gCATTAGCTCACTGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCAcagcag                         >  1:3159397/1‑58 (MQ=255)
                                                                 cACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGAt    >  1:866461/1‑70 (MQ=255)
                                                                  aCCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGAtt   <  1:2387472/70‑1 (MQ=255)
                                                                  aCCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGAtt   <  1:3285422/70‑1 (MQ=255)
                                                                  aCCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTa  <  1:3460087/71‑1 (MQ=255)
                                                                   ccGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGa                        <  1:2539775/48‑1 (MQ=255)
                                                                    |                                                                    
CTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCATCTTTGCGGGCAAACCACAGCAGATACTCCTCCGGCAGGTCGATTA  >  W3110S.gb/2533399‑2533535

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: