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	3893611–3894994	3901930	6937–8320	19 [18]	[0] 75	[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	= 3652533	NA (NA)	72 (1.250)	53/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

TCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893535‑3893626                                                                            |                 tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGc                  >  2:1554282/1‑76 (MQ=255)    aaaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa               >  1:892084/1‑76 (MQ=255)     aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa              <  2:614820/76‑1 (MQ=255)       tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAAt            >  1:422424/1‑76 (MQ=255)        cTAACGCCAGACTGTGCTACACCCATTTTGTTGGACGGTGAAATGGAATAGCCCCTAATATGTCAAAGCCAACATa           >  2:1028534/1‑76 (MQ=14)        cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa           >  2:673905/1‑76 (MQ=255)        cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa           >  2:617124/1‑76 (MQ=255)        cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa           >  2:1028508/1‑76 (MQ=255)        cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa           >  2:1028524/1‑76 (MQ=255)         tAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc          >  1:129147/1‑76 (MQ=255)          aaCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAACATAcc         >  2:343529/1‑76 (MQ=255)            cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt       >  2:143380/1‑76 (MQ=40)            cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt       >  2:194282/1‑76 (MQ=40)             gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtt      >  2:20123/1‑76 (MQ=39)               cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    >  2:1038002/1‑76 (MQ=35)               cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCtttt    <  2:1000700/76‑1 (MQ=35)                aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   >  2:235112/1‑76 (MQ=34)                aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg   >  2:1705437/1‑76 (MQ=34)                 gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa  >  1:40589/1‑76 (MQ=32)                                                                            |                 TCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGA  >  REL606/3893535‑3893626

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