Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	167,484	Δ5,964 bp		[fhuA]–fhuB	[fhuA], fhuC, fhuD, fhuB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	167484	173447	5964	29 [0]	[0] 29	[fhuA]–fhuB	[fhuA],fhuC,fhuD,fhuB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 167483	0 (0.000)	27 (0.900)	23/98	0.1	100%	intergenic (+219/‑1)	mrcB/fhuA	peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
	?	NC_000913	173448 =	0 (0.000)					intergenic (+4/+154)	fhuB/hemL	iron(III) hydroxamate ABC transporter membrane subunit/glutamate‑1‑semialdehyde aminotransferase

CGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173402‑173497                                               |                                                  tcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 >  1:602813‑M2/47‑49 (MQ=255)  cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:1965786‑M2/46‑50 (MQ=255)    ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             >  1:1920772‑M2/44‑50 (MQ=255)         gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:676686‑M2/9‑1 (MQ=255)         gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:1230760‑M2/9‑1 (MQ=255)         gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:1830658‑M2/11‑1 (MQ=255)         gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                        >  1:2034387‑M2/39‑50 (MQ=255)          tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       >  1:1456448‑M2/38‑50 (MQ=255)          tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       >  1:1959950‑M2/38‑50 (MQ=255)            tacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      <  1:625345‑M2/14‑1 (MQ=255)             acgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                    >  1:2201642‑M2/35‑50 (MQ=255)                ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:540361‑M2/32‑49 (MQ=255)                 tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                >  1:897479‑M2/31‑50 (MQ=255)                   tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               <  1:1831087‑M2/21‑1 (MQ=255)                   tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               <  1:1362911‑M2/21‑1 (MQ=255)                      ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           >  1:1705883‑M2/26‑50 (MQ=255)                      ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           <  1:1500503‑M2/25‑1 (MQ=255)                             tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:704119‑M2/19‑50 (MQ=255)                              acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:288724‑M2/18‑49 (MQ=255)                                atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                 <  1:2094072‑M2/35‑1 (MQ=255)                                   agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg               <  1:2121491‑M2/37‑1 (MQ=255)                                   agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:1812992‑M2/38‑1 (MQ=255)                                   agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              >  1:1149498‑M2/13‑50 (MQ=255)                                    gagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCg             >  1:1449230‑M2/12‑50 (MQ=255)                                    gagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCg             >  1:1130813‑M2/12‑50 (MQ=255)                                       atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:668100‑M2/9‑49 (MQ=255)                                          taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcg       >  1:395912‑M2/6‑50 (MQ=255)                                               ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:1772005/50‑1 (MQ=255)                                               ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:1017196/50‑1 (MQ=255)                                               |                                                  CGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173402‑173497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.