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	33 [0]	[0] 32	[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)	31 (1.130)	28/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

AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167436‑167483 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495                                                                                                  AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATT                                                >  1:2270713/1‑50  ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               <  1:978209/50‑1  ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATT                                                <  1:236210/49‑1   TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:947332/1‑49    CATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                             >  1:1841605/1‑50     ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                             >  1:1077550/1‑49        CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          >  1:1071664/1‑49        CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:1644934/49‑1        CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:1134290/49‑1           GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                      <  1:1773715/50‑1           GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                      >  1:1896621/1‑50           GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                      >  1:1118892/1‑50              TACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGC                                    <  1:1323561/49‑1               ACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCC                                   <  1:267948/49‑1                CGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCT                                 <  1:2014195/50‑1                    ATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTA                              <  1:63024/49‑1                     TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACG                            >  1:1550883/1‑50                          CACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                        <  1:781300/49‑1                           ACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                      <  1:108570/50‑1                               TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                  >  1:529952/1‑50                                ACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                  <  1:1226676/49‑1                                     AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG             <  1:1397737/49‑1                                       AGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGT          >  1:677985/1‑50                                       AGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGT          >  1:1491939/1‑50                                        GATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT         <  1:605795/50‑1                                         ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTT        <  1:1496175/50‑1                                         ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT         >  1:405844/1‑49                                          TATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTA       >  1:2178415/1‑50                                             ACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGC    >  1:2245114/1‑50                                              CCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCG   <  1:2315053/50‑1                                               CATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  1:2145504/1‑50                                                                                                  AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167436‑167483 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495

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.