Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,249,944 C→A 11.2% P89T (CCG→ACG)  yqjG → predicted S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,249,9440CA11.2% 157.3 / 4.4 53P89T (CCG→ACG) yqjGpredicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (28/19);  new base A (2/4);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTAT  >  W3110S.gb/3249878‑3250012
                                                                  |                                                                    
caccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                  <  1:204875/71‑1 (MQ=255)
  ccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGa                                                                <  1:2303857/71‑1 (MQ=255)
  ccgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGa                                                                <  1:2985331/71‑1 (MQ=255)
   cgcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGAt                                                               >  1:584162/1‑71 (MQ=255)
    gcacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCt                                                                  <  1:338351/67‑1 (MQ=255)
     cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAcc                                                                      >  1:2288950/1‑62 (MQ=255)
     cacGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGc                                                             >  1:3937875/1‑71 (MQ=255)
       cGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTg                                                           >  1:2380850/1‑71 (MQ=255)
        gCTTATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTg                                                                 >  1:2262564/1‑64 (MQ=255)
        gCTGATCATGTGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTg                                                                 >  1:4111460/1‑64 (MQ=255)
          tGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAAccc                                                                     <  1:939182/58‑1 (MQ=255)
             tCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTg                                                                 >  1:3897614/1‑59 (MQ=255)
               aTGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACg                                                    <  1:2265292/70‑1 (MQ=255)
                 gcgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGtt                                                 >  1:1953738/1‑71 (MQ=255)
                  cgcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTg                                                >  1:3058861/1‑71 (MQ=255)
                   gcAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTgg                                               <  1:2234052/71‑1 (MQ=255)
                     aaaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTgg                                               >  1:3841645/1‑69 (MQ=255)
                      aaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                      >  1:2596192/1‑61 (MQ=255)
                      aaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                      >  1:239156/1‑61 (MQ=255)
                      aaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                      >  1:1441694/1‑61 (MQ=255)
                      aaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACTGTTGGAc                                             >  1:410261/1‑70 (MQ=255)
                      aaCTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTgg                                               <  1:2832862/68‑1 (MQ=255)
                           aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                         <  1:2636795/69‑1 (MQ=255)
                           aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                       >  1:2904029/1‑71 (MQ=255)
                           aaaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCttt                                         <  1:4234561/69‑1 (MQ=255)
                            aaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCtt                                          <  1:29112/67‑1 (MQ=255)
                            aaGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtgat                                      <  1:927012/71‑1 (MQ=255)
                              ggACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTtgatga                                    <  1:3594449/71‑1 (MQ=255)
                                 cTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTg                                                                 >  1:1162835/1‑39 (MQ=255)
                                  tGGAACCGTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCttt                                         <  1:900911/62‑1 (MQ=255)
                                    gAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGttt                              >  1:2637311/1‑71 (MQ=255)
                                     aaCCGTTTATTTCCGTTTCCGTAGTGACCCCGCTGATGCt                                                            >  1:960087/1‑40 (MQ=255)
                                      aCCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTc                            <  1:1848751/71‑1 (MQ=255)
                                       ccGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGaaaa                                                      >  1:3042598/1‑44 (MQ=255)
                                         gTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCttt                                         <  1:2857634/55‑1 (MQ=255)
                                         gTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                          >  1:3337269/1‑70 (MQ=255)
                                         gTTTATTTCCGTTTCCGTAGTGAACACGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCgg                         <  1:1311937/71‑1 (MQ=255)
                                          tttATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACgg                                                   <  1:2809233/44‑1 (MQ=255)
                                            tATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAg                      >  1:1589721/1‑71 (MQ=255)
                                             aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACa                                  <  1:501148/58‑1 (MQ=255)
                                             aTTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACa                                  <  1:2289747/58‑1 (MQ=255)
                                                  cgtttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTtga                                       >  1:1619652/1‑48 (MQ=255)
                                                    tttccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga              <  1:2894599/71‑1 (MQ=255)
                                                      tccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCg                >  1:1291767/1‑67 (MQ=255)
                                                      tccgtAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGgaga            >  1:4103912/1‑71 (MQ=255)
                                                            gTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGga              >  1:3303520/1‑63 (MQ=255)
                                                            gTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGc       <  1:2241171/70‑1 (MQ=255)
                                                             tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCg                          <  1:1699026/50‑1 (MQ=255)
                                                             tGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGctc     <  1:3787682/71‑1 (MQ=255)
                                                              gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACg        >  1:1495185/1‑67 (MQ=255)
                                                              gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGctc     >  1:3081304/1‑70 (MQ=255)
                                                              gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGctc     >  1:638074/1‑70 (MQ=255)
                                                              gAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGctc     >  1:676709/1‑70 (MQ=255)
                                                               aaCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTa   >  1:1221881/1‑71 (MQ=255)
                                                                aCCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTAt  >  1:1872823/1‑71 (MQ=255)
                                                                  ccGCTGATGCTGGAAAACGGTTGGACCTTTGATGACa                                  >  1:3444243/1‑37 (MQ=255)
                                                                  |                                                                    
CACCGCACGCTGATCATGCGCAAACTCAAAGGACTGGAACCGTTTATTTCCGTTTCCGTAGTGAACCCGCTGATGCTGGAAAACGGTTGGACCTTTGATGACAGTTTTCCGGGAGCAACCGGAGACACGCTCTAT  >  W3110S.gb/3249878‑3250012

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: ‑