Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,306,738 C→G 19.6% R401P (CGC→CCC)  yfaW ← predicted enolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,306,7380C→G19.6% 5.8 ‑3.6 0.0341202yfaWpredicted enolase
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (13/14):  total (20/14)
Fisher's exact test strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACAGCAACGCCCGGAGTGATTAGTGGCTGTAGGGGCG.TTTCAAATTGCAGTCACGATTGAGTTCGACGCC  >  REL606/2306703‑2306772
                                   |                                   
acaGCAACGCCCGGAGTGATTAGTGGCTGTAGGGGc                                     <  1:616398/36‑1 (MQ=255)
  aGCAACGCCCGGAGTGATTAGTGGCTGTGGGGGCgg                                   >  1:5328411/1‑35 (MQ=255)
    cAACGCCCGGAGTGATTAGTGGCTGAGGGGGGG.ttt                                >  1:7355666/1‑36 (MQ=37)
      aCGCCCGGAGTGATTAGTGGCTGTAGGGGCG.TTTCa                              <  1:5770934/36‑1 (MQ=255)
       cGCCCGGAGTGATTAGTGGCTGTAGGGGCG.TTTCaa                             <  1:4646896/36‑1 (MQ=255)
       cGCCCGGAGTGATTAGTGGCTGTAGGGGCG.TTTCaa                             <  1:4480618/36‑1 (MQ=255)
         cccGGAGTGATTAGTGGCTGTAGGGGGG.TTTTAAAt                           >  1:2240668/1‑36 (MQ=38)
         cccGGAGTGATTAGTGGCTGTAGGGGCG.TTTCAAAt                           <  1:5577298/36‑1 (MQ=255)
          ccGGAGTGATTAGTGGCTGTGGGGGCG.TTTCAAAtt                          >  1:2892032/1‑36 (MQ=255)
           cGGAGTGATTAGTGGCTGTAGGGGCG.TTTCAAATTg                         <  1:3177474/36‑1 (MQ=255)
           cGGAGTGATTAGTGGCTGGCGGGGGGCTTGCAAAtt                          >  1:3538383/1‑36 (MQ=9)
            ggAGTGATTAGTGGCTGTAGGGGGG.GTTTAAAtttt                        >  1:2212722/1‑34 (MQ=25)
            ggAGTGATTAGTGGCTGTAGGGGCG.TTTCAAATTGc                        <  1:4382636/36‑1 (MQ=255)
             gAGTGATTAGTGGCTGTATGGGGG.TTTTAAAtttta                       >  1:2414669/1‑33 (MQ=19)
              aGTGATTAGTGGCTGTAGGGGGG.TTTTTAATTTCAg                      >  1:3295511/1‑36 (MQ=16)
              aGTGATTAGTGGCTGTAGGGGCG.TTTCAAATTGCAg                      <  1:186152/36‑1 (MQ=255)
               gTGATTAGTGGCTGTAGGGGGG.TTTTAAATTTCAGt                     >  1:3870844/1‑36 (MQ=37)
               gTGATTAGTGGCTGTAGGGGCG.TTTCAAATTGCAGt                     <  1:4220269/36‑1 (MQ=255)
                 gATTAGTGGCTGTAGGGGCG.TTTCAAATTTCAGTcc                   >  1:2371492/1‑35 (MQ=255)
                 gATTAGTGGCTGTAGGGGCG.TTTCAAATTGCAGTCa                   <  1:2747792/36‑1 (MQ=255)
                  aTTAGTGGCTGTANGGGCG.TTTTAAATTTCAGTCAc                  >  1:4318239/1‑36 (MQ=37)
                      gTGGCTGTAGGGGCG.TTTCAAATTTCAGTCACGAtt              >  1:5279462/1‑36 (MQ=255)
                      gTGGCTGTAGGGGCG.TTTCAAATTGCAGTCACGAtt              >  1:7248424/1‑36 (MQ=255)
                        ggCTGTAGGTGCG.TTTCAAATTGCAGTCACGATTGa            <  1:6645187/36‑1 (MQ=255)
                        ggCTGTAGGGTCG.ATTCAAATTGCAGTCACGATTGa            >  1:4394770/1‑36 (MQ=38)
                        ggCTGTAGGGGCG.TTTCAAATTGCAGTCACGATTGa            >  1:720087/1‑36 (MQ=255)
                         gCTGTAGGGGCG.TTTCTAATTGCAGTCACGATTGAg           >  1:4722214/1‑36 (MQ=255)
                         gCTGTAGGGGCG.TTTCAAATTGCAGTCACGATTGAg           >  1:3897428/1‑36 (MQ=255)
                         gCTGTAGGGGCG.TTTCAAATTGCAGTCACGATTGAg           <  1:1968067/36‑1 (MQ=255)
                            gTAGGGGCG.TTTCAAATTGCAGTCACGATTGAGTTc        >  1:3352697/1‑36 (MQ=255)
                            gTAGGGGCG.TTTCAAATTGCAGTCACGATTGAGTTc        >  1:1135799/1‑36 (MQ=255)
                              aGGGGCG.TTTCAAATTGCAGTCACGATTGAGTTCGa      <  1:1742916/36‑1 (MQ=255)
                                gggCG.TTTCAAATTGCAGTCGCGATTGAGTTCGACg    <  1:3941333/36‑1 (MQ=255)
                                 ggCG.TTTCAAATTGCAGTCACGATTGAGTTCGACGc   >  1:714656/1‑36 (MQ=255)
                                  gCG.TTTCAAATTGCAGTCACGATTGAGTTCGACGcc  <  1:4402834/36‑1 (MQ=255)
                                   |                                   
ACAGCAACGCCCGGAGTGATTAGTGGCTGTAGGGGCG.TTTCAAATTGCAGTCACGATTGAGTTCGACGCC  >  REL606/2306703‑2306772

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG