Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,744,613 T→C 19.4% D272G (GAT→GGT)  waaG ← glucosyltransferase I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,744,6130T→C19.4% 3.1 ‑1.9 ‑0.330815waaGglucosyltransferase I
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (8/16):  total (8/22)
Fisher's exact test strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.49e-01

GCTGCTTCCTGGTAGGCAGGATGCAGTAATAAATCCGCCGCCGCCATTAATTCCGAGACATCGTTGC  >  REL606/3744580‑3744646
                                 |                                 
gctgctTCCTGGTAGGCAGGATGCAGTAATAAATcc                                 <  1:7158505/36‑1 (MQ=255)
  tgctTCCTGGTAGGCAGGATGCAGTAATAAATccgc                               <  1:5671827/36‑1 (MQ=255)
   gctTCCTGGTAGGCAGGATGCAGTAATAAATccgcc                              <  1:356032/36‑1 (MQ=255)
   gctTCCTGGTAGGCAGGATGCAGTAATAAATccgcc                              <  1:3830040/36‑1 (MQ=255)
   gctTCCAGGTAGGCAGGATTCAGTAATAAATccgcc                              <  1:6409206/36‑1 (MQ=37)
    ctTCCTGGTAGGCAGGATGCAGTAATAAATccgccg                             >  1:5320153/1‑36 (MQ=255)
       ccTGGTAGGCAGGATGCAGTAATAAATccgccgccg                          <  1:3385972/36‑1 (MQ=255)
       ccTGGTAGGCAGGATGCAGTAATAAATccgccgccg                          <  1:4347912/36‑1 (MQ=255)
        cTGGTAGGCAGGATGCAGTAATAAATCcgccgccgc                         >  1:6356912/1‑36 (MQ=255)
        cTGGTAGGCAGGATGCAGTAATAAATCcgccgccgc                         <  1:3486462/36‑1 (MQ=255)
         tGGTAGGCAGGATGCAGTAATAAATCCgccgccgcc                        <  1:2593722/36‑1 (MQ=255)
          ggTAGGCAGGATGCAGTAATAAATCCGCCGCCGCCa                       >  1:5667389/1‑36 (MQ=255)
           gTAGGCAGGATGCAGTAATAAATCCGCCGCCGCCAt                      <  1:5970015/36‑1 (MQ=255)
           gTAGGCAGGATGCAGTAATAAATCCGCCGCCGCCAt                      >  1:5288455/1‑36 (MQ=255)
             aGGCAGGATGAAGTAATAAATCCGCCGCCGCCATTa                    <  1:672568/36‑1 (MQ=255)
             aGGCAGAATGCAGAAAAAAACCCGCCGCCGCCATTa                    <  1:5634835/36‑1 (MQ=16)
              ggCAGGATGCAGTAATAAACCGGCCGCCGCCATTaa                   <  1:3040215/36‑1 (MQ=37)
              ggCAGGATGCAGAAATAAACCCGCCGCCGCCATTaa                   <  1:4973318/36‑1 (MQ=37)
               gCAGGATGCAGTAATAAATCCGCCGCCGCCATTAAt                  <  1:5978741/36‑1 (MQ=255)
               gCAGGATGCAGTAATAAATCCGCCGCCGCCATTAAt                  <  1:1558568/36‑1 (MQ=255)
                 agaaTGCAGTAATAAATCCGCCGCCGCCATTAATTc                <  1:291129/33‑1 (MQ=255)
                 aGGATGCATTACTAAATCCGCGGCCGCCATTAATTc                <  1:5700382/36‑1 (MQ=25)
                  ggATGCAGTAATAAATCCGCCGCCGCCATTAATTcc               >  1:5703100/1‑36 (MQ=255)
                     tGCAGTAATAAATCCGCCGCCGCCATTAATTCCgag            >  1:3460191/1‑36 (MQ=255)
                        aGTAATAAATCCGCCGCCGCCATTAATTCCGAGACa         >  1:4824989/1‑36 (MQ=255)
                          taataaACCCGCCGCCGCCATTAATTCCGAGACATc       <  1:1869096/36‑1 (MQ=255)
                             taaATCCGCCGCCGCCATTAATTCCGAGACATCGtt    >  1:953946/1‑36 (MQ=255)
                             taaACCCGCCGCCGCCATTAATTCCGAGACATCGtt    <  1:3903655/36‑1 (MQ=255)
                             aaaaCCCGCCGCCGCCATTAATTCCGAGACATCGtt    <  1:6861692/35‑1 (MQ=255)
                               aaTCCGCCGCCGCCATTAATTCCGAGACATCGTTGc  <  1:994562/36‑1 (MQ=255)
                                 |                                 
GCTGCTTCCTGGTAGGCAGGATGCAGTAATAAATCCGCCGCCGCCATTAATTCCGAGACATCGTTGC  >  REL606/3744580‑3744646

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG