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 3893618–3894994 3901930 6937–8313 15 [14] [0] 62 [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)58 (1.150) 48/146 0.0 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 0 (0.000)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

CTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893542‑3893626
                                                                           |         
cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa           >  1:1546876/1‑76 (MQ=255)
 tAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc          <  1:251548/76‑1 (MQ=255)
 tAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc          <  1:280346/76‑1 (MQ=255)
 tAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc          <  1:2820/76‑1 (MQ=255)
 tAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc          >  1:493123/1‑76 (MQ=255)
  aaCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAcc         >  1:1123196/1‑76 (MQ=255)
    tgCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt       >  2:766909/2‑76 (MQ=39)
    cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt       <  1:1621780/76‑1 (MQ=40)
     gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt      >  2:1118039/1‑76 (MQ=39)
       cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    <  2:149696/76‑1 (MQ=35)
       cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  2:984689/1‑76 (MQ=35)
        aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   >  1:1608259/1‑76 (MQ=34)
        aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   >  1:28357/1‑76 (MQ=34)
        aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACACTTTTg   <  1:444020/76‑1 (MQ=21)
         gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  <  1:527513/76‑1 (MQ=32)
                                                                           |         
CTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893542‑3893626

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