Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 3,894,997 Δ6,934 bp IS150‑mediated rbsD[yieO] rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3893616–3894994 3901930 6937–8315 8 [5] [0] 20 [insJ‑5]–[yieO] [insJ‑5],insK‑5,rbsD,rbsA,rbsC,rbsB,rbsK,rbsR,[yieO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3652533NA (NA)18 (0.750) 18/146 0.4 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 0 (0.000)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

TGGGCATATGTTTGCGCGCGTAAAGAAGGCCTGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGACAAAA  >  REL606/3901858‑3902006
                                                                         |                                                                           
tgactccaattgaatatcggaatcagacctatatgcctcgtgtttaattgtccaactttttggggtcagtacaGGt                                                                           <  2:763780‑M2/3‑1 (MQ=255)
tgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGt                                                                           >  1:360449‑M2/74‑76 (MQ=255)
      caattgaatatcggaatcagacctntatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCa                                                                     >  2:131797‑M2/68‑76 (MQ=255)
         ttgaatatcggaatcagccctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGt                                                                  >  2:173969‑M2/65‑76 (MQ=255)
          tgaatatctgaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTg                                                                 <  2:419729‑M2/13‑1 (MQ=255)
              tatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCa                                                             >  2:599145‑M2/60‑76 (MQ=255)
                            ctatatgcctagtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGa                                               >  1:111532‑M2/46‑76 (MQ=255)
                                 tgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACggg                                          >  1:780229‑M2/41‑76 (MQ=255)
                                          tttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAAt                                 >  2:195744‑M2/32‑76 (MQ=255)
                                                   ccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCccac                        >  1:393132‑M2/23‑76 (MQ=255)
                                                     aactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACAGGGCCTAAAATTGGCCccacca                      <  2:485619‑M2/56‑1 (MQ=255)
                                                          tttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGAccc                 >  2:372846‑M2/16‑76 (MQ=255)
                                                              gggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCa             <  1:697581‑M2/65‑1 (MQ=255)
                                                                gtcagtacaGGTTGCCCAGGTCACCAGTACGCCGGCAAGAACGGGGCCTAAAATTGGCCCCACCAGACACGGCATg           <  2:704975‑M2/67‑1 (MQ=255)
                                                                 tcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCANACCCGGTATgg          <  1:118289‑M2/68‑1 (MQ=255)
                                                                   agtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCg        >  2:433243‑M2/7‑76 (MQ=255)
                                                                     tacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGAc      >  1:386970‑M2/5‑76 (MQ=255)
                                                                      gcaGGTTGCCTAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCGGAGCCGGCATGGCGACa     <  1:114937‑M2/73‑1 (MQ=255)
                                                                      acaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGACa     <  2:340387‑M2/73‑1 (MQ=255)
                                                                         ggTTGCCCAGGTGACCAGCATGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGACaaaa  <  2:663486/76‑1 (MQ=255)
                                                                         |                                                                           
TGGGCATATGTTTGCGCGCGTAAAGAAGGCCTGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGGCATGGCGACAAAA  >  REL606/3901858‑3902006

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