Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,391,592 A→G 30.8% E137G (GAG→GGG)  yjfC → predicted synthetase/amidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,391,5920A→G30.8% 3.2 ‑2.0 0.018410yjfCpredicted synthetase/amidase
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (4/8):  total (9/9)
Fisher's exact test strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGAT  >  REL606/4391558‑4391626
                                  |                                  
tACGAGTCGGCTTATTTCCAGTGGCTGTGGCTggag                                   <  1:2683378/36‑1 (MQ=255)
 aCGCGTCGGCTTATTTCCAGTGGCTGTGGCTggagg                                  <  1:3212666/36‑1 (MQ=255)
 aCGAGTCGGCTTATTTCCAGTGGCTGGGGCTggggg                                  >  1:1331249/1‑36 (MQ=37)
  cGAGTCGTCTTATTTCCAGTGGCTGTGGCTggagga                                 <  1:5193419/36‑1 (MQ=255)
  cGAGTCGGCTTATTTCCAGTGGCTGTGGCTggagga                                 <  1:2474697/36‑1 (MQ=255)
   gAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGAt                                <  1:4538921/36‑1 (MQ=255)
     gTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGc                              <  1:4693818/36‑1 (MQ=255)
       cGGCTTATTTCCAGTGGCTGTGGGTGGGGGGTGgcc                            >  1:1117667/1‑36 (MQ=19)
         gCTTATTTCCAGTGGCTGTGGGTTGGGGGTGTCCgg                          >  1:3381354/1‑36 (MQ=16)
          cTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGgc                         <  1:491341/36‑1 (MQ=255)
             aTTTCCAGTGGCTGTGGGTGGGGGGTGGCCGgggcg                      >  1:1634270/1‑35 (MQ=9)
              tttCCAGTGGCTGTGGGTTGGGGGTTCCCGgggcgg                     >  1:5779592/1‑34 (MQ=9)
              tttCCAGTGGCTGTGGCTGGGGGATGCCCCgcgccg                     >  1:3323607/1‑34 (MQ=25)
                  cAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGc                 <  1:1368698/36‑1 (MQ=255)
                     tggctgtggctgGAGGATGCCCGGCGCAGCGGCatt              <  1:874716/36‑1 (MQ=255)
                        ctgtggctgGAGGATGCCCGGCGCAGCGGCattatt           <  1:4003770/36‑1 (MQ=255)
                        ctgtggataGAGGATGCCCGGCGCAGCGGCattatt           <  1:2249043/36‑1 (MQ=37)
                         tgtggctgGAGGATGCCCGGCGCAGCGGGATTTTTc          >  1:3145843/1‑36 (MQ=25)
                          gtggctgGAGGATGCCCGGGGCAGCGGCATTATTcc         <  1:852610/36‑1 (MQ=255)
                           tggctgGAGGATGCCCGGAGCAGCGGCATTATTCcg        >  1:4734695/1‑36 (MQ=255)
                            ggctgGGNNNTGCCCGGCGCCGCGGCATTATTCcgc       <  1:1987197/36‑1 (MQ=25)
                             gctgGAGGATGCCAGGCGCAGCGGCATTATTCcgcg      >  1:3853969/1‑36 (MQ=255)
                               tgGAGGATGCCCGGCGCAGAGGCATTATTCCGCGTg    >  1:1363793/1‑36 (MQ=255)
                                ggaggaTGCCCGGCGCAGCGGCATTATTCCGCGTGa   >  1:774518/1‑36 (MQ=255)
                                ggaggaTGCCCGGCGCAGCGGCATTATTCCGCGTGa   >  1:83479/1‑36 (MQ=255)
                                 gaggaTGCCCGGCGCAGGGGCATTATTCCGCGTGAt  >  1:564206/1‑36 (MQ=255)
                                 gaggaTGCCCGGCGCAGCGGCATTATTCCGCGTGAt  <  1:4994325/36‑1 (MQ=255)
                                  |                                  
TACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGAT  >  REL606/4391558‑4391626

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