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 3893600–3894994 3901930 6937–8331 15 [14] [1] 61 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)59 (1.350) 45/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

GGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTG  >  REL606/3893524‑3893625
                                                                           |                          
gggTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa                            <  1:1394504/76‑1 (MQ=255)
 ggTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat                           <  1:552737/76‑1 (MQ=255)
   tACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                         >  2:17406/1‑76 (MQ=255)
   tACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatat                         <  2:414624/76‑1 (MQ=255)
     cTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGt                       <  2:1159826/76‑1 (MQ=255)
       cAGGTCTAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGGCGATGAAATGGAATAGCCCCTAATATGTCa                     <  1:1305809/76‑1 (MQ=255)
       cAGGTCGAAATCTAACGCCAGAATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCa                     >  2:859696/1‑76 (MQ=255)
           tCGAAATCTAACGCCAGACTGTACTGCACTCATTTTCTTGGACGATGAAATGGAATANCCCCTAATATGTCAAAGc                 <  1:1361972/76‑1 (MQ=18)
              aaaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa              >  2:60109/1‑76 (MQ=255)
               aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa             <  1:270366/76‑1 (MQ=255)
                aTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaaa            >  2:676936/1‑76 (MQ=255)
                 tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAAt           <  1:1086262/76‑1 (MQ=255)
                  cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa          >  1:942267/1‑76 (MQ=255)
                         cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt   >  2:521843/1‑76 (MQ=35)
                          aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg  <  2:581610/76‑1 (MQ=34)
                                                                           |                          
GGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTG  >  REL606/3893524‑3893625

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