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	3893610–3894994	3901930	6937–8321	18 [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	= 3652533	NA (NA)	62 (1.140)	45/146	0.2	100%	noncoding (1443/1443 nt)	IS150	repeat region
	?	REL606	3901931 =	0 (0.000)					coding (489/1428 nt)	yieO	predicted multidrug or homocysteine efflux system

GTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAA  >  REL606/3893534‑3893627                                                                            |                   gTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg                    <  1:304452/76‑1 (MQ=255)   cGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGcc                  <  1:1037657/76‑1 (MQ=255)     aaaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaa                >  1:241934/1‑76 (MQ=255)      aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa               <  1:1209716/76‑1 (MQ=255)         cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa            >  1:1484514/1‑76 (MQ=255)           aaCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAcc          >  2:1524894/1‑76 (MQ=255)            aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAccc         >  1:308674/1‑76 (MQ=255)             cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCt        <  1:1423404/76‑1 (MQ=40)               ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCttt      <  2:873049/76‑1 (MQ=38)                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg    >  1:253206/1‑76 (MQ=34)                 aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg    <  2:1469255/76‑1 (MQ=34)                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  1:1490261/1‑76 (MQ=32)                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  1:1472925/1‑76 (MQ=32)                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  2:1504654/1‑76 (MQ=32)                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  1:1331093/1‑76 (MQ=32)                  gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa   >  2:830216/1‑76 (MQ=32)                   aCTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaa  >  2:1232328/1‑76 (MQ=31)                   aCTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaa  <  1:1167155/76‑1 (MQ=31)                                                                            |                   GTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAA  >  REL606/3893534‑3893627

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