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 3893623–3894996 3901930 6935–8308 19 [17] [0] 64 [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)63 (1.020) 43/146 0.3 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 0 (0.000)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

GCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893547‑3893626
                                                                           |    
gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt      >  1:407030/1‑76 (MQ=39)
gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt      <  2:415877/76‑1 (MQ=39)
 ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCttt     <  1:886167/76‑1 (MQ=38)
 ccAGACTGTACTGCACCCATTTGGTTGGACGATGAAATGGAATAGCCCCTAATATTTCAAAGCCAAAATACCCttt     <  2:168988/76‑1 (MQ=255)
  cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCACAATACCCtttt    >  1:803317/1‑76 (MQ=25)
  cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  2:898131/1‑76 (MQ=35)
  cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  1:803308/1‑76 (MQ=35)
  cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  2:1686104/1‑76 (MQ=35)
  cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  2:1876836/1‑76 (MQ=35)
  cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  2:200246/1‑76 (MQ=35)
   aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   <  2:559071/76‑1 (MQ=34)
   aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   >  1:1585617/1‑76 (MQ=34)
   aGACTGTACTGCACCCATTTTGTTGGACGATGAAATCGACTAGCCCCTNATATGTCAAAGCCAAAATACCCTTTTg   <  1:1783116/76‑1 (MQ=255)
    gACTGTACTGCACCNATTTTGTTGGGCGATGAAATGGAATCGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  1:1979507/1‑76 (MQ=255)
    gACTGTACTGCACCNATTTTGTTGGACGATGAAATGGAATAGCCCCTAATCTGTCAAAGCCAAAATACCCTTTTGa  >  1:4329/1‑76 (MQ=12)
    gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  1:872249/1‑76 (MQ=32)
    gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  <  2:1896974/76‑1 (MQ=32)
    gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  <  1:723214/76‑1 (MQ=32)
    gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  <  1:1649465/76‑1 (MQ=32)
                                                                           |    
GCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893547‑3893626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: