Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,516,335 A→T 15.4% intergenic (+166/+417) yjhR → / ← yjhS frameshift suppressor contains phospolipase D/nuclease; b4308_2/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,516,3350A→T15.4% 4.7 ‑3.4 0.026intergenic (+166/+417)yjhR/yjhSframeshift suppressor contains phospolipase D/nuclease; b4308_2/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (8/14):  total (12/14)
Fisher's exact test strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CATAATGCGCCTGGGATTAATGGC.GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAAT  >  REL606/4516301‑4516370
                                   |                                   
cATAATGCGCCTGGGATTAATGGC.GGGTGCTTTTAg                                    <  1:4067467/36‑1 (MQ=255)
        gcCTGGGATTAATGGC.GGGTGCTTTTAGGAATgaaa                            <  1:3940158/36‑1 (MQ=255)
           tGGGATTAATGGC.GGGTGCTTTTAGGAATgaaagaa                         <  1:991648/36‑1 (MQ=255)
            gggATTAATGGCGGGGGGGTTTTTGAAATgaaaaaa                         >  1:1255692/1‑36 (MQ=2)
             ggATTAATGGC.GGGTGCTTTTAGGAATGAAAGaaaa                       <  1:1713147/36‑1 (MQ=255)
             ggATTAATGGC.GGGGGGTTTTTTTAATGTAAAaaaa                       >  1:3989419/1‑36 (MQ=2)
             ggATTAATGGC.GGGGGCGTTTTTGTAAGAAAAaaaa                       >  1:1434811/1‑36 (MQ=2)
               aTTAATGGC.GGGCGCTTTTTGGTATTATAGAAAAta                     >  1:1226057/1‑35 (MQ=12)
                ttAATGGC.GGTTGCTTTTAGGAATGAAAGAAAATgc                    <  1:5334329/36‑1 (MQ=255)
                 tAATGGC.GGGTGCTTTTAGGAATGAAAGAAAATgca                   <  1:2469853/36‑1 (MQ=255)
                   aTGGC.GGGTGCTTTTAGGAATGAAAGAAAATgcaag                 <  1:2094936/36‑1 (MQ=255)
                    tGGC.GGGTGCTTTTAGGAATGAAAGAAAATgcaagc                >  1:4192863/1‑36 (MQ=255)
                      gC.GGGTGCTTTTAGGAATGAAAGAAAATgcaagcaa              <  1:3754210/36‑1 (MQ=255)
                      gC.GGGTGCTTTCAGGAATGAAAGAAAATgcaagcaa              <  1:5352616/36‑1 (MQ=255)
                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCaata            >  1:1675664/1‑36 (MQ=255)
                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCaata            <  1:5292756/36‑1 (MQ=255)
                          ggTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa           >  1:5897844/1‑36 (MQ=255)
                          ggTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa           >  1:4094807/1‑36 (MQ=255)
                           gTGCTTTTAGGAATGAAAGAAAATGCAAGCaataat          <  1:1353332/36‑1 (MQ=255)
                           atGCTTTAAGGAATTAAAGAAAATGCAAGCaataat          <  1:4200655/35‑1 (MQ=37)
                            tGCTTTTAGGAATGAAAGAAAATGCAAGCAataata         >  1:403912/1‑36 (MQ=255)
                               ttttAGGAATGAAAGAAAATGCAAGCAataatataa      <  1:4133224/36‑1 (MQ=255)
                               ttttAGGAATGAAAGAAAATGCAAGCAataatataa      >  1:5867290/1‑36 (MQ=255)
                                 ttAGGAATGAAAGAAAATGCAAGCAAtaatataata    >  1:5651891/1‑36 (MQ=255)
                                 ttAGGAATGAAAGAAAATGCAAGCAAtaatataata    <  1:2258455/36‑1 (MQ=255)
                                 ttAGGAATGAAAGAAAATGCAAGCAAtaatataata    >  1:1050353/1‑36 (MQ=255)
                                   aGGAATGAAAGAAAATGCAAGCAATAATATaataat  >  1:2735433/1‑36 (MQ=255)
                                   |                                   
CATAATGCGCCTGGGATTAATGGC.GGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAAT  >  REL606/4516301‑4516370

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