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 3893614–3894994 3901930 6937–8317 15 [13] [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)61 (1.380) 43/146 0.1 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 0 (0.000)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

AAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893538‑3893626
                                                                           |             
aaaTCTAACGCCAGACTGTACTNCACCCATTTTATTGGNCGATGAAACGGAACAGCCCCTCATAGGTCGANACCaa               >  1:28125/1‑76 (MQ=255)
aaaTCTAACGCCAGACTGTACTGCACCCATATTGTTGGACGATGAAATGGAATAGCCCCTAANATGTCAAAGCCaa               >  2:1321125/1‑76 (MQ=17)
 aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa              <  1:271113/76‑1 (MQ=255)
 aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa              <  2:967866/76‑1 (MQ=255)
         gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt      >  2:837594/1‑76 (MQ=39)
          ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCttt     >  1:1211756/1‑76 (MQ=38)
           cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  1:425728/1‑76 (MQ=35)
           cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  1:437429/1‑76 (MQ=35)
            aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   >  2:349593/1‑76 (MQ=34)
            aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   <  2:974760/76‑1 (MQ=34)
             gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  1:392373/1‑76 (MQ=32)
             gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  1:720796/1‑76 (MQ=32)
             gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  1:818686/1‑76 (MQ=32)
             gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  2:572340/1‑76 (MQ=32)
             gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  2:683128/1‑76 (MQ=32)
                                                                           |             
AAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893538‑3893626

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