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	49 [0]	[0] 50	[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)	49 (1.310)	37/98	0.0	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

CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173399‑173496                                                  |                                                 taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                  <  1:314773‑M2/1‑1 (MQ=255)   atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                <  1:2395423‑M2/3‑1 (MQ=255)   atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                >  1:1275352‑M2/48‑50 (MQ=255)    tcattctcgtttacgttatcattcactttacatcagagatataccaTTcc                                               <  1:1929200‑M2/4‑1 (MQ=255)    tcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                >  1:1692686‑M2/47‑49 (MQ=255)       ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            <  1:745008‑M2/7‑1 (MQ=255)       ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            >  1:693083‑M2/44‑50 (MQ=255)       ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            <  1:710117‑M2/7‑1 (MQ=255)        tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            >  1:1190388‑M2/43‑49 (MQ=255)         ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                           <  1:560261‑M2/8‑1 (MQ=255)         ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                           <  1:1472260‑M2/8‑1 (MQ=255)         ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                           <  1:734138‑M2/8‑1 (MQ=255)            gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                        <  1:1155770‑M2/11‑1 (MQ=255)             tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                      <  1:2324366‑M2/13‑1 (MQ=255)             tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                      >  1:2178806‑M2/38‑50 (MQ=255)             tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                      <  1:326160‑M2/13‑1 (MQ=255)             tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                      <  1:1434060‑M2/13‑1 (MQ=255)                  gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                 >  1:1424584‑M2/33‑50 (MQ=255)                  gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                 >  1:79612‑M2/33‑50 (MQ=255)                  gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                 <  1:2529079‑M2/18‑1 (MQ=255)                   ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                 >  1:208619‑M2/32‑49 (MQ=255)                   ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                <  1:386472‑M2/19‑1 (MQ=255)                   ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                <  1:1359119‑M2/19‑1 (MQ=255)                     atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                              >  1:1688284‑M2/30‑50 (MQ=255)                     atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                              >  1:822311‑M2/30‑50 (MQ=255)                      tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                              <  1:854417‑M2/21‑1 (MQ=255)                      tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                              <  1:2082984‑M2/21‑1 (MQ=255)                      tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                             <  1:1943769‑M2/22‑1 (MQ=255)                       cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                            >  1:410145‑M2/28‑50 (MQ=255)                        attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGcc                           >  1:360876‑M2/27‑50 (MQ=255)                         ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                          >  1:1372035‑M2/26‑50 (MQ=255)                          tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                         >  1:1828355‑M2/25‑50 (MQ=255)                           cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                         >  1:283749‑M2/24‑49 (MQ=255)                            actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                       <  1:495500‑M2/28‑1 (MQ=255)                              tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                     >  1:2319112‑M2/21‑50 (MQ=255)                                tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                   >  1:36191‑M2/19‑50 (MQ=255)                                  catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                  <  1:2218473‑M2/33‑1 (MQ=255)                                  catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                  <  1:2195786‑M2/33‑1 (MQ=255)                                  catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                 <  1:159698‑M2/34‑1 (MQ=255)                                  catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                 <  1:606347‑M2/34‑1 (MQ=255)                                   atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                 >  1:598383‑M2/16‑49 (MQ=255)                                     cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg              >  1:2017384‑M2/14‑50 (MQ=255)                                     cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg              >  1:797677‑M2/14‑50 (MQ=255)                                          atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt          >  1:422850‑M2/9‑49 (MQ=255)                                            ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc       >  1:1091676‑M2/7‑50 (MQ=255)                                             taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcg      >  1:1933297‑M2/6‑50 (MQ=255)                                             taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc       <  1:1167749‑M2/44‑1 (MQ=255)                                              accaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc     >  1:769645‑M2/5‑50 (MQ=255)                                                caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa   >  1:1132325‑M2/3‑50 (MQ=255)                                                  ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc  >  1:806656/1‑49 (MQ=255)                                                  |                                                 CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173399‑173496

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