Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 4,415,710 IS150 (–) +3 bp coding (848‑850/1413 nt) cycA → D‑alanine/D‑serine/glycine transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3893554 =NA (NA)71 (1.100) 57/200 0.1 100% noncoding (1/1443 nt) IS150 repeat region
?REL606 4415710 = 0 (0.000)coding (848/1413 nt) cycA D‑alanine/D‑serine/glycine transporter
* ? REL606 = 3894996NA (NA)80 (1.260) 64/198 0.0 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 4415712 0 (0.000)coding (850/1413 nt) cycA D‑alanine/D‑serine/glycine transporter

AGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCC  >  REL606/3894899‑3895005
                                                                                                 |         
agaaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc        <  2:1677981/101‑1 (MQ=32)
agaaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc        <  2:340227/101‑1 (MQ=32)
      aGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAggctgtact  <  2:115273/101‑9 (MQ=2)
                                                                                                 |         
AGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCC  >  REL606/3894899‑3895005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ 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.