Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 192,898 T→C 25.9% A9A (GCT→GCC frr → ribosome recycling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb192,8980TC25.9% 88.7 / 34.1 54A9A (GCT→GCCfrrribosome recycling factor
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (4/10);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GACTGTCTGGTCTGACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTG  >  W3110S.gb/192829‑192965
                                                                     |                                                                   
gACTGTCTGGTCTGACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCt                                                                     >  1:2491511/1‑70 (MQ=255)
 aCTGTCTGGTCTGACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTg                                                                    >  1:971488/1‑70 (MQ=255)
 aCTGTCTGGTCTGACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTg                                                                    <  1:83701/70‑1 (MQ=255)
 aCTGTCTGGTCTGACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTg                                                                    <  1:1198933/70‑1 (MQ=255)
    gtctggtctgACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAg                                                                 <  1:3541528/70‑1 (MQ=255)
     tctggtctgACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAg                                                                 <  1:526258/69‑1 (MQ=255)
       tggtctgACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACg                                                             >  1:141028/1‑71 (MQ=255)
         gtctgACTGAGACAAGTTTTCAAGGATACGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAGTACg                                                             <  1:3410076/69‑1 (MQ=255)
            tgactgaGACAAGTTTTCAATGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAg                                                                 <  1:1779154/62‑1 (MQ=255)
            tgactgaGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAg                                                                 >  1:3644489/1‑62 (MQ=255)
            tgactgaGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACg                                                             <  1:4414413/66‑1 (MQ=255)
            tgactgaAACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATg                                                         >  1:1832831/1‑70 (MQ=255)
               ctgaGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATg                                                         <  1:2146117/67‑1 (MQ=255)
                  agaCAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAAt                                                  >  1:2752484/1‑71 (MQ=255)
                    aCAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAGTACGCATGGACAAATGc                                                <  1:1476151/71‑1 (MQ=255)
                     cAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCg                                               <  1:2010665/71‑1 (MQ=255)
                     cAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCg                                               <  1:382479/71‑1 (MQ=255)
                      aaGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGa                                                                   <  1:4080533/50‑1 (MQ=255)
                      aaGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACg                                                             <  1:3807868/56‑1 (MQ=255)
                      aaGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAg                                                                 >  1:1062750/1‑52 (MQ=255)
                      aaGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAg                                                                 >  1:894089/1‑52 (MQ=255)
                       aGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATgag                                               >  1:1769555/1‑67 (MQ=255)
                       aGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAGTACGCATGGACAAATGCGTa                                             <  1:2552439/71‑1 (MQ=255)
                         ttttCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAg                                            >  1:101443/1‑70 (MQ=255)
                           ttCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAg                                                                 >  1:4170230/1‑47 (MQ=255)
                           ttCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGaa                                          >  1:1750202/1‑70 (MQ=255)
                           ttCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAg                                         >  1:2328034/1‑71 (MQ=255)
                                 gATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACa                                                     <  1:648651/53‑1 (MQ=255)
                                 gATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCa                                   <  1:3309649/71‑1 (MQ=255)
                                  aTTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTa                                             <  1:889526/60‑1 (MQ=255)
                                   ttCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCaa                                  >  1:4163289/1‑70 (MQ=255)
                                     cGTAACGTGATTAGCGATATCAGAAAAGATGCCGAAGTACGCATGGACAAATGCGTAGAAGCGTTCaaaa                                <  1:1629616/70‑1 (MQ=255)
                                       tAACGTGATTAGCGATATCAGAAAAGATGCCGAAGTACGCATGGACa                                                     <  1:3331243/47‑1 (MQ=255)
                                       tAACGTGATTAGCGATATCAGAAAAGATGCCGAAGTACGCATGGACa                                                     <  1:4159213/47‑1 (MQ=255)
                                           gTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATg                                                 >  1:557803/1‑47 (MQ=255)
                                           gTGATTAGCGATATCAGAAAAGATGCCGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAAt                         <  1:3293195/71‑1 (MQ=255)
                                            tGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATc                        <  1:2744160/71‑1 (MQ=255)
                                               ttAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAAt                         <  1:110168/67‑1 (MQ=255)
                                                tAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCaaa                          >  1:686269/1‑65 (MQ=255)
                                                tAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAg                      >  1:1962367/1‑69 (MQ=255)
                                                 aGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCaa                                  >  1:3491868/1‑56 (MQ=255)
                                                   cGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCaaaa                 <  1:3434962/71‑1 (MQ=255)
                                                   cGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCaaaa                 <  1:4193651/71‑1 (MQ=255)
                                                       atCAGAAAAGATGCCGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACg             >  1:2368869/1‑71 (MQ=255)
                                                        tCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTa                                             <  1:3449361/38‑1 (MQ=255)
                                                        tCAGAAAAGATGCCGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACg             >  1:3657474/1‑70 (MQ=255)
                                                        tCAGAAAAGATGCCGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACg             <  1:3509761/70‑1 (MQ=255)
                                                         cAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTa                                             >  1:3576827/1‑37 (MQ=255)
                                                         cAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATc                        <  1:2246943/58‑1 (MQ=255)
                                                          aGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAAt                         >  1:1662604/1‑56 (MQ=255)
                                                              aaGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAg                      >  1:4394623/1‑55 (MQ=255)
                                                              aaGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGt      >  1:886711/1‑71 (MQ=255)
                                                              aaGATGCCGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGt      <  1:2471707/71‑1 (MQ=255)
                                                                 aTGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGt   <  1:316217/71‑1 (MQ=255)
                                                                   gCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTg  <  1:4388973/70‑1 (MQ=255)
                                                                     |                                                                   
GACTGTCTGGTCTGACTGAGACAAGTTTTCAAGGATTCGTAACGTGATTAGCGATATCAGAAAAGATGCTGAAGTACGCATGGACAAATGCGTAGAAGCGTTCAAAACCCAAATCAGCAAAATACGCACGGGTCGTG  >  W3110S.gb/192829‑192965

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: