Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,111,0950C→G100.0% 8.4 17G94A (GGC→GCC) insF‑1putative transposase for IS3
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (7/0):  total (17/0)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.80e-02

CTGCTCCAACAGATTTTCTGACACAGGCAGGCCGTGTGCGCGGTAG  >  REL606/1111064‑1111109
                               |              
cTGCTCCAACATATTTTCTGACACAGGCAGGGCgg             >  1:3758104/1‑34 (MQ=1)
cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg             >  1:3272260/1‑34 (MQ=0)
cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg             >  1:2759660/1‑34 (MQ=0)
cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg             >  1:15628429/1‑34 (MQ=0)
cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg             >  1:6226455/1‑34 (MQ=0)
cTGCTCCAACAGATTTTCTGACACAGGCAGGGCgg             >  1:15216934/1‑34 (MQ=0)
cTGCTCCAACAGATTTTCTGACACAGGCAGGACgg             >  1:7660244/1‑34 (MQ=0)
cTGCTCCAACAGATGTTCTGACACAGGCAGGGCgg             >  1:1746972/1‑34 (MQ=1)
 tGCTCCAACAGATTTTCTGACACAGGCAGGGCggg            >  1:16962353/1‑33 (MQ=1)
 tGCTCCAACAGATTTTCTGACACAGGCAGGGCggg            >  1:3559215/1‑33 (MQ=1)
 tGCTCCAACAGATTTTCGGACACAGGCAGGGCggg            >  1:2862191/1‑33 (MQ=1)
  gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg           >  1:14867958/1‑32 (MQ=1)
  gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg           >  1:13680243/1‑32 (MQ=1)
  gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg           >  1:15733221/1‑32 (MQ=1)
  gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg           >  1:231342/1‑32 (MQ=1)
  gCTCCAACAGATTTTCTGACACAGGCAGGGCgggg           >  1:12393622/1‑32 (MQ=1)
   cTCCAACAGATTTTCTGACACAGGCAGGGCgggta          >  1:13067323/1‑34 (MQ=1)
     ccAACAGATTTTCTGACACAGGCAGGCCgggtttg        >  1:7478450/1‑32 (MQ=1)
     ccAACAGATTTTCTGACACAGGCAGGCCGGGGgcg        >  1:5149351/1‑35 (MQ=1)
     ccAACAGATTTTCTGACACAGGCAGGCCGGGGgcg        >  1:14935944/1‑35 (MQ=1)
     ccAACAGATTTTCTGACACAGGCAGGCCGGGGgcg        >  1:4661932/1‑35 (MQ=1)
        aCAGATTTTCTGACACAGGCAGGGCGGGGGGGCgg     >  1:2483216/1‑35 (MQ=0)
        aCAGATTTTCTGACACAGGCAGGGCGGGGGGGCgg     >  1:11222377/1‑35 (MQ=0)
        aCAGATTTTCTGACACAGGCAGGCCGGGGGGGCgg     >  1:5833624/1‑35 (MQ=1)
        aCAGATTTTCTGACACAGGCAGGCCGGGGGGGCgg     >  1:9143553/1‑35 (MQ=1)
         cAGATTTTCTGACACAGGCAGGGCGGGGGGGCggg    >  1:602470/1‑34 (MQ=0)
         cAGATTTTCTGACACAGGCAGGCCGCGGGGGCggg    >  1:13327073/1‑34 (MQ=1)
          aGATTTTCTGACACAGGCAGGCCGGGTGCGGgggg   >  1:13730080/1‑33 (MQ=1)
           gATTTTCTGACACAGGCAGGGCGGGGGCGGggggg  >  1:18010159/1‑32 (MQ=0)
           gATTTTCTGACACAGGCAGGCCGTGGGCGGgggaa  >  1:1273468/1‑32 (MQ=1)
                               |              
CTGCTCCAACAGATTTTCTGACACAGGCAGGCCGTGTGCGCGGTAG  >  REL606/1111064‑1111109

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.