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	34 [0]	[0] 33	[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)	33 (1.220)	25/96	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

TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAA  >  NC_000913/167435‑167527                                                 |                                             taatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                              <  1:401689/49‑1 (MQ=255)   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                            <  1:406868‑M1/49‑3 (MQ=255)    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                           <  1:1186255‑M1/49‑4 (MQ=255)    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                           >  1:331022‑M1/1‑46 (MQ=255)      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                         >  1:404836‑M1/1‑44 (MQ=255)      aTTCTCGTTTACGTTAACATTCACTTTACATCAGAGATATACCAtacc                                          <  1:1354244‑M1/48‑5 (MQ=255)         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                       <  1:1205441‑M1/48‑8 (MQ=255)         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                       <  1:1658203‑M1/48‑8 (MQ=255)            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                   >  1:826444‑M1/1‑38 (MQ=255)             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                  <  1:598017‑M1/49‑13 (MQ=255)             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                  >  1:445685‑M1/1‑37 (MQ=255)             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                  <  1:188798‑M1/49‑13 (MQ=255)                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                             <  1:1257853‑M1/49‑18 (MQ=255)                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                             <  1:647258‑M1/49‑18 (MQ=255)                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                            >  1:995966‑M1/1‑30 (MQ=255)                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                         <  1:183268‑M1/49‑22 (MQ=255)                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                          <  1:402410‑M1/48‑21 (MQ=255)                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                      <  1:1587674‑M1/49‑25 (MQ=255)                            aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                    >  1:296447‑M1/1‑22 (MQ=255)                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                  >  1:1678083‑M1/1‑21 (MQ=255)                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga              >  1:1944758‑M1/1‑17 (MQ=255)                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga              >  1:1287392‑M1/1‑17 (MQ=255)                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat             <  1:72292‑M1/49‑34 (MQ=255)                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag         >  1:537066‑M1/1‑12 (MQ=255)                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag         >  1:1384172‑M1/1‑12 (MQ=255)                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag         >  1:745666‑M1/1‑12 (MQ=255)                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag         >  1:1041533‑M1/1‑12 (MQ=255)                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagc        >  1:1743411‑M1/1‑11 (MQ=255)                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcg       >  1:334918‑M1/1‑10 (MQ=255)                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgt      >  1:1854702‑M1/1‑9 (MQ=255)                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt     >  1:284888‑M1/1‑8 (MQ=255)                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt     <  1:1515884‑M1/49‑42 (MQ=255)                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac  >  1:930974‑M1/1‑5 (MQ=255)                                             taCCAttccatataacggccctgtacgcctggacggataaagcgttta   <  1:1469368‑M1/48‑44 (MQ=255)                                                 |                                             TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAA  >  NC_000913/167435‑167527

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: ‑
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.