Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,208,832 A→G 13.6% H126R (CAC→CGC)  yrdD → predicted DNA topoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,208,8320AG13.6% 106.9 / 10.5 44H126R (CAC→CGC) yrdDpredicted DNA topoisomerase
Reads supporting (aligned to +/- strand):  ref base A (18/20);  new base G (3/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATTACATGCCCCCAATGTCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGGAGAATGCCCTG  >  W3110S.gb/4208765‑4208900
                                                                   |                                                                    
aTTACATGCCCCCAATGTCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCActc                                                                   >  1:1196568/1‑71 (MQ=255)
  tACATGCCCCCAATGTCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCtt                                                                 >  1:872611/1‑71 (MQ=255)
      tGCCCCCAATGTCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGa                                                             >  1:966117/1‑71 (MQ=255)
        cccccAATGTCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATc                                                           >  1:1946295/1‑71 (MQ=255)
                 tCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCGCTCTTGTGATCGCTACCCgg                                                  <  1:464851/71‑1 (MQ=255)
                 tCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGAt                                                            >  1:560593/1‑61 (MQ=255)
                 tCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCgg                                                  <  1:289264/71‑1 (MQ=255)
                  cGGACGGGCCATCTGGTCCAGCGCCTCTCCCGTTATGGCAAAACATTTCGCTCTTGTGATCg                                                          >  1:2423516/1‑62 (MQ=255)
                    gACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTtgt                                                               <  1:1442887/55‑1 (MQ=255)
                    gACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAgt                                               <  1:652923/71‑1 (MQ=255)
                      cGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGa                                                             >  1:1054458/1‑55 (MQ=255)
                        ggCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCgg                                                  <  1:2321401/64‑1 (MQ=255)
                             tCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCGCTCTtgtg                                                              <  1:1016875/47‑1 (MQ=255)
                             tCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCGCTCTTGTGATCGCTa                                                       <  1:969397/54‑1 (MQ=255)
                               tGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCGCTCTTGTGATCGCTACCCGGAGTGTCAATTTGc                                     >  1:475675/1‑70 (MQ=255)
                                 gTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCgg                                                  <  1:2265023/55‑1 (MQ=255)
                                 gTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGc                                     <  1:1391599/68‑1 (MQ=255)
                                  tCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCAtt                                 <  1:1153467/71‑1 (MQ=255)
                                  tCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCAtt                                 <  1:543820/71‑1 (MQ=255)
                                  tCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCAtt                                 <  1:2444926/71‑1 (MQ=255)
                                    cAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTa                                <  1:1661609/70‑1 (MQ=255)
                                      gcgcCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGc                                     <  1:2797904/63‑1 (MQ=255)
                                        gccgcTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACtt                            >  1:1922077/1‑70 (MQ=255)
                                              cccGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAAcc                      >  1:449525/1‑70 (MQ=255)
                                               ccGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCa                    <  1:2389376/71‑1 (MQ=255)
                                                  ttATGGCAAAACATTTCACTCTTGTGATCGCTACCCg                                                   <  1:2280521/37‑1 (MQ=255)
                                                   tATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCaa                         >  1:499813/1‑62 (MQ=255)
                                                    aTGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAg                                                >  1:2085249/1‑38 (MQ=255)
                                                    aTGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTg                                      >  1:870215/1‑48 (MQ=255)
                                                     tGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAg                                                >  1:1823471/1‑37 (MQ=255)
                                                       gCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAAcc                      <  1:1426117/61‑1 (MQ=255)
                                                       gCAAAACATTTCACTCTTGTGATCGCTAACCGGAGTGTCAATTTGCCATTAACTCCAAAcc                      <  1:574596/61‑1 (MQ=255)
                                                        cAAAACATTTCGCTCTTGTGATCGCTACCCGGAGTGTCAAtt                                        >  1:2692951/1‑42 (MQ=255)
                                                        cAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGgag           <  1:1172840/71‑1 (MQ=255)
                                                        cAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGgag           <  1:1503266/71‑1 (MQ=255)
                                                        cAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGgag           <  1:2339788/71‑1 (MQ=255)
                                                            aCATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTaa                               >  1:2630917/1‑47 (MQ=255)
                                                            aCATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAAc                       >  1:654793/1‑55 (MQ=255)
                                                             cATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAAc                       >  1:1927414/1‑54 (MQ=255)
                                                               tttCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGgag           <  1:1211470/64‑1 (MQ=255)
                                                                 tCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGGAGAATGCCCt   >  1:1806457/1‑70 (MQ=255)
                                                                 tCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGGAGAATGCCCt   >  1:2062590/1‑70 (MQ=255)
                                                                  cACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGGAGAATGCCCt   >  1:590155/1‑69 (MQ=255)
                                                                  cACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGGAGAATGCCCTg  <  1:1224801/70‑1 (MQ=255)
                                                                   |                                                                    
ATTACATGCCCCCAATGTCGGACGGGCCATCTGGTCCAGCGCCGCTCCCGTTATGGCAAAACATTTCACTCTTGTGATCGCTACCCGGAGTGTCAATTTGCCATTAACTTCAAACCCATAGCTGGAGAATGCCCTG  >  W3110S.gb/4208765‑4208900

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: