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–3894995 3901930 6936–8322 11 [10] [0] 34 [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)32 (0.890) 27/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

TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAG  >  REL606/3652461‑3652607
                                                                        |                                                                          
tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                                                         >  1:521308/1‑76 (MQ=32)
tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                                                         <  2:267264/76‑1 (MQ=32)
 gACTCCACTTGAATATCGGAATCAGACCTATATGNCTCGTGTTTAACTGTCNAACTTTTTGGGGACAGNACAAACt                                                                        >  2:1056032/1‑76 (MQ=11)
  aCTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                                                                       >  1:49264/1‑76 (MQ=35)
   cTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACctt                                                                      >  1:561512/1‑76 (MQ=17)
     ccAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGa                                                                    <  2:444706/76‑1 (MQ=40)
         ttGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTta                                                                <  2:458568/76‑1 (MQ=255)
           gAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTtata                                                              <  2:803161/76‑1 (MQ=255)
             atatCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                                                            >  1:785106/1‑76 (MQ=255)
                   gAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTg                                                      <  1:221650/76‑1 (MQ=255)
                       cAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGc                                                  >  2:440048/1‑76 (MQ=255)
                        aGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt                                                 >  2:23567/1‑76 (MQ=255)
                         gACCTATATGCCTCGTGTTTAACTGTCCGACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCtt                                                >  2:1155856/1‑76 (MQ=255)
                            cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc                                             >  1:415238/1‑76 (MQ=255)
                            cTATATGCCTCGTGTTTAACTGACCNACTTTTAGGGATCAGTATAAACTTTGNTTTATAGTCNGGTGGGGCntttn                                             >  1:1022984/1‑75 (MQ=21)
                               tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACCAACTTTGATTTATAGTCAGGTGGGGCTTTTctct                                          >  2:553527/1‑76 (MQ=255)
                               tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTNTctct                                          >  2:957936/1‑76 (MQ=255)
                                atGCCTCGTGTTTAACTGTCCAGCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtctcgg                                         <  1:860935/76‑3 (MQ=255)
                                   ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtct                                      <  1:485059/76‑1 (MQ=255)
                                   ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtct                                      <  2:419408/76‑1 (MQ=255)
                                      cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGcc                                   >  2:525054/1‑76 (MQ=255)
                                         gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt                                >  1:1601/1‑76 (MQ=255)
                                           ttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg                              <  1:528997/76‑1 (MQ=255)
                                           ttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg                              <  1:528987/76‑1 (MQ=255)
                                            tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg                             >  1:71925/1‑76 (MQ=255)
                                              aCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTg                           <  1:991579/76‑1 (MQ=255)
                                              aCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTg                           >  1:649515/1‑76 (MQ=255)
                                                tGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGaa                         <  2:334001/76‑1 (MQ=255)
                                                 gTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAAt                        >  1:761578/1‑76 (MQ=255)
                                                 gTCCAACTTATTGGGGTCAGTACAAACTTTGATTTATATTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAAt                        >  1:833320/1‑76 (MQ=255)
                                                  tCCAACTTTTTGGGGTCANTACAAACTTTGATTAATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATa                       <  1:935405/76‑1 (MQ=255)
                                                    cAACTTTTTGGGTATAGGACAAACNTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCNGCCTTTCGGTGAATAcc                     <  2:990127/76‑1 (MQ=255)
                                                        tttttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTgag                 >  1:622627/1‑76 (MQ=255)
                                                        tttttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGATTTTCTCTGTCTGCCTTTCGGTGAATACCTgag                 >  1:772205/1‑76 (MQ=255)
                                                            tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACaa             <  2:9661/76‑1 (MQ=255)
                                                              gggTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAAc           <  2:1186426/76‑1 (MQ=255)
                                                               ggTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACa          <  1:273117/76‑1 (MQ=255)
                                                                 tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGt        <  2:71925/76‑1 (MQ=255)
                                                                   aGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATAGCTGAGACAAACAGtct      >  2:767448/1‑76 (MQ=255)
                                                                      aCAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCaa   >  2:273876/1‑76 (MQ=255)
                                                                      aCAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCaa   <  2:746064/76‑1 (MQ=255)
                                                                      aCAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCaa   >  2:224965/1‑76 (MQ=255)
                                                                       cAAACTTTGACTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCGCAAg  <  1:163832/76‑1 (MQ=255)
                                                                        |                                                                          
TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAG  >  REL606/3652461‑3652607

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