Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	3,894,997	Δ6,138 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	3893587–3894993	3901134	6142–7548	33 [30]	[0] 152	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	= 3894996	NA (NA)	148 (1.900)	69/98	0.0	100%	noncoding (1443/1443 nt)	IS150	repeat region
	?	REL606	3901135 =	0 (0.000)					coding (1285/1428 nt)	yieO	predicted multidrug or homocysteine efflux system

CGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATA  >  REL606/3893536‑3893601                                                   |                tgAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                 <  1:4644043/50‑1 (MQ=255) cGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                 <  1:1657444/51‑1 (MQ=255) cGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATg                 <  1:5284079/51‑1 (MQ=255)  gAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATgg                <  1:352128/51‑1 (MQ=255)   aaaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGa               >  1:307292/1‑51 (MQ=255)    aaTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGaa              <  1:2194596/51‑1 (MQ=255)     aTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt             >  1:5283302/1‑51 (MQ=255)     aTCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt             >  1:4350688/1‑51 (MQ=255)      tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:5499956/51‑1 (MQ=255)      tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2203008/51‑1 (MQ=255)      tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2270214/51‑1 (MQ=255)      tCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2063118/51‑1 (MQ=255)       cTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAg           >  1:4635512/1‑51 (MQ=255)         aaCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATa            <  1:2925144/48‑1 (MQ=255)          aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        >  1:4461796/1‑51 (MQ=255)          aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        <  1:7191698/51‑1 (MQ=255)          aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        <  1:7007078/51‑1 (MQ=255)          aCGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc        >  1:4292986/1‑51 (MQ=255)           cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc       >  1:6529167/1‑51 (MQ=255)           cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc       >  1:6687293/1‑51 (MQ=255)           cGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc       >  1:6720716/1‑51 (MQ=255)            gCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt      <  1:5765465/51‑1 (MQ=255)             ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa     >  1:4956105/1‑51 (MQ=40)             ccAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa     >  1:4020073/1‑51 (MQ=40)              cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa    <  1:894067/51‑1 (MQ=38)              cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa    >  1:5837419/1‑51 (MQ=38)              cAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa    <  1:5529550/51‑1 (MQ=38)               aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   >  1:6644406/1‑51 (MQ=35)               aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:2958245/51‑1 (MQ=35)               aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:2664828/51‑1 (MQ=35)               aGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat   <  1:776138/51‑1 (MQ=35)                gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata  <  1:2816290/51‑1 (MQ=33)                gACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata  <  1:5119266/51‑1 (MQ=33)                                                   |                CGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATA  >  REL606/3893536‑3893601

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