Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,029,877 C→G 22.2% G143A (GGC→GCC)  yccW ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,029,8770CG22.2% 47.4 / 10.3 27G143A (GGC→GCC) yccWpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/6);  new base G (2/4);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAAC  >  W3110S.gb/1029816‑1029940
                                                             |                                                               
cagTTTCATTCCGTGTTCTTAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAg                                                       >  1:3886628/1‑71 (MQ=255)
   tttCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCa                                                        <  1:1075648/68‑1 (MQ=255)
     tCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGGCCTGGGTCa                                                        >  1:1281987/1‑66 (MQ=255)
     tCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCt                                                     >  1:3412888/1‑69 (MQ=255)
          ccGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTcc                                              >  1:1751662/1‑71 (MQ=255)
          ccGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTcc                                              >  1:88529/1‑71 (MQ=255)
          ccGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTc                                               >  1:1578466/1‑70 (MQ=255)
               ttcttcAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCcttctt                                         >  1:2071679/1‑71 (MQ=255)
               ttcttcAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCcttctt                                         >  1:3735362/1‑71 (MQ=255)
               ttcttcAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCcttct                                          <  1:2211650/70‑1 (MQ=255)
                 cttcAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCtttt                                       >  1:970206/1‑71 (MQ=255)
                 cttcAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCtttt                                       >  1:1868187/1‑71 (MQ=255)
                     aaTCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCT‑‑GTCAGCTCCATTCCTTCTTTTTTACGt                                 <  1:123355/71‑1 (MQ=255)
                       tCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGt                                 <  1:754506/71‑1 (MQ=255)
                                   caggTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCg                     >  1:366892/1‑71 (MQ=255)
                                   caggTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCg                     >  1:371580/1‑71 (MQ=255)
                                      gTGGCAACTCGCCGGTGACGGGGGCCTGGGTCAGCTCCATTCCTTCttt                                        <  1:3398011/49‑1 (MQ=255)
                                      gTGGCAACTCGCCGGTGACGGGGGCCTGGGTCAGCTCCATTCCTTCttt                                        <  1:819016/49‑1 (MQ=255)
                                        ggCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGc                 >  1:3008612/1‑70 (MQ=255)
                                         gCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCct                                             >  1:1646283/1‑41 (MQ=255)
                                           aaCTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATc             >  1:3968296/1‑71 (MQ=255)
                                              tCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTa                                    >  1:3397089/1‑45 (MQ=255)
                                               cGCCGGTGACGGGGGCCTGGGTCAGCTCCATTCCTTCttt                                        >  1:212641/1‑40 (MQ=255)
                                                  cGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATc      <  1:2326198/71‑1 (MQ=255)
                                                     tGACGGGGGCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGaa   <  1:3338884/71‑1 (MQ=255)
                                                      gACGGGGGCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGaa   <  1:908578/70‑1 (MQ=255)
                                                       aCGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAAc  <  1:3295184/70‑1 (MQ=255)
                                                             |                                                               
CAGTTTCATTCCGTGTTCTTCAATCGGCAGCAGGGCAGGTGGCAACTCGCCGGTGACGGGGCCCTGGGTCAGCTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAAC  >  W3110S.gb/1029816‑1029940

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: