Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,464,140 A→G 16.5% D521G (GAT→GGT)  trg → methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,464,1400A→G16.5% 2.3 ‑3.3 0.0421562trgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (20/14):  total (28/14)
Fisher's exact test strand distribution p-value = 3.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAGAACAGGCGGCACGATTAACTGAGGCGGTGGATGTATTCCGTCTGCACAAACATTCTGTGTCGGCAG  >  REL606/1464107‑1464175
                                 |                                   
aagaaCAGGCGGCACGATTAACTGAGGCGGGGGGTg                                   >  1:3135675/1‑36 (MQ=37)
aagaaCAGGCGGCACGATTAACTGAGGCGGGGGATg                                   >  1:66778/1‑36 (MQ=255)
 agaaCAGGCGGCACGATTAACTGAGGCGGTGGATGt                                  <  1:1843841/36‑1 (MQ=255)
 aaaaCAGTCTGCACGATTAGCTGAGGCGGTGGATGt                                  <  1:2850882/34‑1 (MQ=25)
  gaaCAGGCGGCACGATTAACTGAGGCGGTGGGTgga                                 >  1:5816518/1‑34 (MQ=38)
   aaCAGGCGGCACGATTAACTGAGGCGGTGGATGTAt                                <  1:3182863/36‑1 (MQ=255)
        gcggcACGATTAACTGAGGCGGTGGGTGTATTCCGt                           >  1:4117938/1‑36 (MQ=255)
         cggcACGATTAACTGAGGCGGCGGGTGTATTCCGTc                          >  1:5251091/1‑36 (MQ=37)
          ggcACGATTAACTGAGGCGGTGGGTGTATTTCGTCt                         >  1:3444707/1‑36 (MQ=38)
          ggcACGATTAACTGAGGCGGTGGATGTTTTTCGTCt                         >  1:1214354/1‑36 (MQ=37)
          ggcACGATTAACTGAGGCGGTGGATGTATTCCGTCt                         <  1:3492602/36‑1 (MQ=255)
           gcACGATTAACTGAGGCGGTGGGTGTATTTCGGCtt                        >  1:738458/1‑35 (MQ=25)
           gcACGATTAACTGAGGCGGGGGATGTATTCCGTCtt                        >  1:657/1‑35 (MQ=255)
              cGATTAACTTAGGCGGTGGATGTATTCCGTCTGcac                     <  1:839855/36‑1 (MQ=255)
              cGATTAACTGAGGCGGTGGATGTATTCCGTCTGcac                     <  1:3943449/36‑1 (MQ=255)
               gATTAACTGAGGCGGTGGATGTATTCCGTCTGcaca                    <  1:5238095/36‑1 (MQ=255)
               gATTAACTGAGGCGGTGGATGTATTCCGTCTGcaca                    >  1:967638/1‑36 (MQ=255)
               gATTAACTGAGGCGGGGGATGTATTTCGTCTGcaca                    >  1:59610/1‑36 (MQ=37)
                aTTAACTGAGGCGGTGGGTGTATTTCTTCTGCCCaa                   >  1:3849326/1‑36 (MQ=16)
                aTTAACTGAGGCGGTGGATGTATTTCGTCTTCACaa                   >  1:4062298/1‑36 (MQ=38)
                  tAACTGAGGCGGTGGATGTATTCCGTCTGCACAAAc                 <  1:4641553/36‑1 (MQ=255)
                  tAACTGAGGCGGGGGATGTGTTTTGTCTGCTCGAAc                 >  1:3522378/1‑36 (MQ=9)
                    aCTGAGGCGGGGGGTGTTTTCCGTCTGCACAAACag               >  1:1413019/1‑35 (MQ=25)
                     cTGAGGCGGTGGAGGTATTTCGTTTGCCCAAACAtt              >  1:1794234/1‑36 (MQ=255)
                        aGGCGGGGGATGTATTCCGTCTTCACAAACATTCtg           >  1:341672/1‑36 (MQ=38)
                         ggCGGTGGATGTATTCCGTCTGCACAAACATTCtgt          <  1:4975113/36‑1 (MQ=255)
                          gCGGTGGATGTATTCCGTCTGCACAAACATTCtgtg         <  1:1005682/36‑1 (MQ=255)
                          gCGGTGGATGTATTCCGTCTGCACAAACATTCtgtg         >  1:4574886/1‑36 (MQ=255)
                          gCGGTGGATGTATTACGTCTGCACAAACATTCtgtg         >  1:2825392/1‑36 (MQ=255)
                           cGGTGGATGTATTCCGTCTGCACAAACATTCTgtgt        <  1:5151203/36‑1 (MQ=255)
                           cGGTGGATGTATTCCGTCTGCACAAACATTCTgtgt        <  1:2088690/36‑1 (MQ=255)
                            ggcggATGTATTCCGACTGCACGAACATTCTGGGTc       >  1:2762679/4‑36 (MQ=19)
                            ggTGGATGTATTCCGTCTGCACAACCATTCTGTGTc       >  1:3038980/1‑36 (MQ=255)
                            ggTGGATGTATTCCGTCTGCACAAACATTCTGTGTc       >  1:3020343/1‑36 (MQ=255)
                            ggTGGATGTATTCCGTCTGCACAAACATTCTGTGTc       >  1:328342/1‑36 (MQ=255)
                            ggTGGATGTATTCCGTCTGCACAAACATTCTGTGTc       <  1:2431942/36‑1 (MQ=255)
                             gTGGATGTATTCCGTCTGCACAAACATTCTGTGTCg      >  1:716125/1‑36 (MQ=255)
                              tGGATGTATTCCGTCTGCACAAACATTCTGTGTCgt     >  1:1888709/1‑35 (MQ=255)
                               ggATGTATTCCGTCTGCACAAACATTCTGTGTCGGc    <  1:1229126/36‑1 (MQ=255)
                                gATGTATTCCGTCTGCACAAACATTCTGTGTCggaa   >  1:3414576/1‑34 (MQ=255)
                                gATGTATTCCGTCTGCACAAACATTCTGTGTCGGCa   >  1:163106/1‑36 (MQ=255)
                                gATGTATTCCGTCTGCACAAACATTCTGTGTCGGCa   >  1:1431452/1‑36 (MQ=255)
                                 aTGTATTCCGTCTGCACAAACATTCTGTGTCGGCAg  >  1:1638740/1‑36 (MQ=255)
                                 aTGTATTCCGTCTGCACAAACATTCTGTGTCGGCAg  <  1:1475923/36‑1 (MQ=255)
                                 |                                   
AAGAACAGGCGGCACGATTAACTGAGGCGGTGGATGTATTCCGTCTGCACAAACATTCTGTGTCGGCAG  >  REL606/1464107‑1464175

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