Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,597,394 C→T 9.0% R352R (CGC→CGT tsr → methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,597,3940CT 9.0% 155.7 / 4.0 55R352R (CGC→CGTtsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base T (2/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGTTAAGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCA  >  W3110S.gb/4597324‑4597464
                                                                      |                                                                      
tGGCGTTAAGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATgcgc                                                                        >  1:3105448/1‑71 (MQ=255)
tGGCGTTAAGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATgcgc                                                                        >  1:1710268/1‑71 (MQ=255)
  gCGTTAAGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGa                                                                      <  1:2946532/71‑1 (MQ=255)
       aaGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATAtct                                                                 >  1:733763/1‑71 (MQ=255)
       aaGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATAtct                                                                 >  1:3706043/1‑71 (MQ=255)
            cTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTccac                                                             >  1:620620/1‑70 (MQ=255)
               cTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGt                                                         >  1:1378973/1‑71 (MQ=255)
                  aaaCGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCg                                                      >  1:446876/1‑71 (MQ=255)
                      ggCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGat                                                                     <  1:3201397/52‑1 (MQ=255)
                      ggCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGat                                                                     <  1:2618154/52‑1 (MQ=255)
                       gcgcAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATAtct                                                                 >  1:458722/1‑55 (MQ=255)
                       gcgcAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGaa                                                 <  1:2802602/71‑1 (MQ=255)
                        cgcAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTccacc                                                            >  1:862681/1‑59 (MQ=255)
                        cgcAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTccacc                                                            >  1:3417120/1‑59 (MQ=255)
                        cgcAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGaaa                                                <  1:2046764/71‑1 (MQ=255)
                            aCGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAAt                                              <  1:3766773/69‑1 (MQ=255)
                             cgcgGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGaaaa                                               >  1:2115862/1‑67 (MQ=255)
                              gcggcgGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                                          <  1:3377699/71‑1 (MQ=255)
                                ggcggTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGa                                        <  1:1252446/71‑1 (MQ=255)
                                 gcggTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGa                                        >  1:1259417/1‑70 (MQ=255)
                                 gcggTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGa                                        >  1:1466849/1‑70 (MQ=255)
                                 gcggTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGa                                        <  1:1812781/70‑1 (MQ=255)
                                     tAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGat                                                                     <  1:1140984/37‑1 (MQ=255)
                                       aaGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTc                                                       >  1:3770276/1‑49 (MQ=255)
                                       aaGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTc                                                       >  1:3533349/1‑49 (MQ=255)
                                       aaGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGAtattat                                  >  1:583220/1‑70 (MQ=255)
                                       aaGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGAtattat                                  >  1:1877425/1‑70 (MQ=255)
                                        aGTGGTAGATAACGTGGTGCAGACTATGCGTGATATCTCCACCAGTTCGCAGAAAATCGCCGatat                                     >  1:429082/1‑66 (MQ=255)
                                        aGTGGTAGATAACGTGGTGCAGACTATGCGTGATATCTCCACCAGTTCGCAGAAAATCGCCGatat                                     >  1:1785290/1‑66 (MQ=255)
                                            gTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCg                             >  1:2224756/1‑70 (MQ=255)
                                            gTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGt                            <  1:78546/71‑1 (MQ=255)
                                              aGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGAtattat                                  >  1:1872734/1‑63 (MQ=255)
                                              aGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGAtattat                                  >  1:821497/1‑63 (MQ=255)
                                               gATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGaa                                                 <  1:1327442/47‑1 (MQ=255)
                                               gATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGc                              >  1:3364210/1‑66 (MQ=255)
                                                 tAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGt                            <  1:1611491/66‑1 (MQ=255)
                                                 tAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGt                            <  1:3502493/66‑1 (MQ=255)
                                                 tAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCAGATATTATCAGCGTAAtt                        >  1:3451745/1‑70 (MQ=255)
                                                  aaCGTGGTGCAGACTATGGGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGa                      >  1:312540/1‑71 (MQ=255)
                                                  aaCGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGa                      >  1:3245621/1‑71 (MQ=255)
                                                   aCGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAAtt                        <  1:2602787/68‑1 (MQ=255)
                                                   aCGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGAc                     >  1:3354751/1‑71 (MQ=255)
                                                      tggtgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCg                             <  1:2916975/60‑1 (MQ=255)
                                                      tggtgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCg                             <  1:2606687/60‑1 (MQ=255)
                                                       ggtgCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGcc                                          <  1:2119906/46‑1 (MQ=255)
                                                        gtgCAGACTATGCGTGATATTTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGt                            <  1:567889/59‑1 (MQ=38)
                                                        gtgCAGACTATGCGTGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGt                            <  1:3707658/59‑1 (MQ=255)
                                                        gtgCAGACTATGCGTGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGt                            <  1:790464/59‑1 (MQ=255)
                                                           cAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGc             >  1:2024536/1‑71 (MQ=255)
                                                                  tGCGCTATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAAt                         >  1:2104720/1‑52 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAAt                         >  1:1201884/1‑52 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGcc            >  1:393095/1‑65 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGcc            >  1:2414788/1‑65 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTcc        <  1:3679471/69‑1 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTc         >  1:3179336/1‑68 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAg      <  1:3608430/71‑1 (MQ=255)
                                                                  tGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAg      <  1:2646814/71‑1 (MQ=255)
                                                                    cgcgATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGc                  <  1:2000465/57‑1 (MQ=255)
                                                                    cgcgATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGc                  <  1:427130/57‑1 (MQ=255)
                                                                    cgcgATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGc                  <  1:1431113/57‑1 (MQ=255)
                                                                      cgATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCa  <  1:1080973/71‑1 (MQ=255)
                                                                      |                                                                      
TGGCGTTAAGTGCTTCTGAAACGGCGCAACGCGGCGGTAAAGTGGTAGATAACGTGGTGCAGACTATGCGCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCA  >  W3110S.gb/4597324‑4597464

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: