Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,389,357 G→A 13.5% S217N (AGT→AAT)  yjfJ → predicted transcriptional regulator effector protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,389,3570G→A13.5% 4.7 ‑3.0 ‑0.037650yjfJpredicted transcriptional regulator effector protein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (15/17):  total (20/17)
Fisher's exact test strand distribution p-value = 4.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

GCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGCAACGCC  >  REL606/4389322‑4389392
                                   |                                   
gCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAg                                     <  1:1736444/36‑1 (MQ=255)
 cTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGt                                    <  1:3128673/36‑1 (MQ=255)
    gCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAgc                                 <  1:6761730/36‑1 (MQ=255)
    gCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAgc                                 <  1:4086195/36‑1 (MQ=255)
      ggAAGCCGGAATTGGCGGTGGCAATAAAAATAacgg                               >  1:1087971/1‑35 (MQ=21)
       gAAGCCGGAATTGGCGGTAGCAATAAAAGAAGCGcc                              >  1:107391/1‑36 (MQ=255)
          gCCGGAATTGGCGGTAGCAATAAAAAAAGCGccccg                           >  1:7758067/1‑34 (MQ=25)
            cGGAATTGGCGGTAGCAATAAAAATAGGGCCCAggg                         >  1:1718878/1‑35 (MQ=25)
            cGGAATTGGCGGTAGCAATAAAAAAAGGGCCCCGGa                         >  1:1485756/1‑36 (MQ=19)
             ggAATTGGCGGTAGTAATAAAAATAGCGCACAgggt                        >  1:4134196/1‑34 (MQ=25)
                aTTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTa                     <  1:4709761/36‑1 (MQ=255)
                 ttGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTAt                    <  1:7914272/36‑1 (MQ=255)
                  tGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTAtt                   <  1:7246789/36‑1 (MQ=255)
                   ggCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTa                  <  1:3062089/36‑1 (MQ=255)
                   ggCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTa                  >  1:3822491/1‑36 (MQ=255)
                   ggCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTa                  <  1:1466099/36‑1 (MQ=255)
                   ggCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTa                  <  1:1451113/36‑1 (MQ=255)
                   ggCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTa                  >  1:447869/1‑36 (MQ=255)
                     cGGTAGCAATAAAAGTAGCGCCCCGGATGTATTAGc                <  1:2807521/36‑1 (MQ=255)
                     cGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGc                >  1:1314260/1‑36 (MQ=255)
                     cGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGc                >  1:6523341/1‑36 (MQ=255)
                      ggTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCa               >  1:5429678/1‑36 (MQ=255)
                       gTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCaa              >  1:7814991/1‑36 (MQ=255)
                         aGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGa            >  1:4000883/1‑36 (MQ=255)
                          gCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGAc           >  1:5451044/1‑36 (MQ=255)
                          gCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGAc           <  1:7844621/36‑1 (MQ=255)
                           cAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACt          >  1:7958410/1‑36 (MQ=255)
                           cAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACt          >  1:3144233/1‑36 (MQ=255)
                           cAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACt          >  1:3019513/1‑36 (MQ=255)
                           cAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACt          >  1:8293655/1‑36 (MQ=255)
                            aaTAAAAGTAGCGCCCAGGATGTATTAGCAAGACTg         >  1:6387803/1‑36 (MQ=255)
                            aaTAAAAGTAGCGCCCAGGATGTATTAGCAAGACTg         <  1:1795798/36‑1 (MQ=255)
                            aaTAAAAGTAGCGCCCAGGATGTATTAGCAACACTg         >  1:4444059/1‑36 (MQ=255)
                             aTAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGc        <  1:2514376/36‑1 (MQ=255)
                              tAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGCa       <  1:8054710/36‑1 (MQ=255)
                                 aaGTAGCGCCCAGGATGTATTAGCAAGACTGCAACg    <  1:4272525/36‑1 (MQ=255)
                                 aaGTAGCGCCCAGGATGTATTAGCAAGACTGCAACg    <  1:1135688/36‑1 (MQ=255)
                                  aGTAGCGCCCAGGATGTATTAGCAAGACTGCAACGc   <  1:1953610/36‑1 (MQ=255)
                                   gTAGCGCCCAGGATGTATTAGCAAGACTGCAACGcc  >  1:5562540/1‑36 (MQ=255)
                                   gTAGCGCCCAGGATGTATTAGCAAGACTGCAACGcc  >  1:849303/1‑36 (MQ=255)
                                   |                                   
GCTGGCGGAAGCCGGAATTGGCGGTAGCAATAAAAGTAGCGCCCAGGATGTATTAGCAAGACTGCAACGCC  >  REL606/4389322‑4389392

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG