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 3893609–3894994 3901930 6937–8322 4 [3] [0] 23 [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)23 (1.070) 20/146 0.2 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 0 (0.000)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

GGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/3652457‑3652533
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGG  >  REL606/3901931‑3901993
                                                                                                                                            
GGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTAC                                                                  <  2:592794/76‑1
GGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTAC                                                                  <  2:311489/76‑1
GGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTAC                                                                  <  2:276993/76‑1
 GTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                                                 <  2:94743/76‑1
 GTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                                                 <  1:83613/76‑1
     GACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                                             <  2:636323/76‑1
        TCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCC                                                          >  2:376269/1‑76
        TCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCC                                                          >  1:106627/1‑76
           AATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAG                                                       >  1:582872/1‑76
                AATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGcc                                                  >  2:325246/1‑74
                  TATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCA                                                <  1:683802/76‑1
                  TATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCA                                                >  2:669821/1‑76
                   ATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAG                                               >  1:423398/1‑76
                   ATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAG                                               <  1:277201/76‑1
                          TCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCG                                        <  2:699102/76‑1
                          TCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTG                                                    <  1:709744/64‑1
                          TCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTG                                                    >  2:709744/1‑64
                              ACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAA                                    <  2:202792/76‑1
                              ACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAA                                    <  2:375011/76‑1
                                      GCCTCGTGTTTAACTGTCCACCTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGG                            <  2:553099/76‑1
                                       CCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGC                           <  2:312002/76‑1
                                       CCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGC                           <  2:60745/76‑1
                                             GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAA                     >  1:174690/1‑76
                                                  ACTGTCCAACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAATTTGGC                >  2:702613/1‑76
                                                         AACTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCA         >  1:461326/1‑76
                                                           CTTTTTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGA       >  2:470682/1‑76
                                                               TTGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCG   <  2:148658/76‑1
                                                                TGGGGTCAGTACAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGG  <  2:452230/76‑1
                                                                                                                                            
GGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/3652457‑3652533
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGG  >  REL606/3901931‑3901993

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