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	35 [0]	[0] 34	[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)	34 (1.120)	32/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

TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167438‑167483 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173448‑173496                                                                                                 TCATTCTCGTTTACGTTATCATTCGCTTTACATCAGAGATATACCATTCC                                               <  1:771122/50‑1  CATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                              >  1:2374602/1‑50   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                              >  1:2223860/1‑49   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                              <  1:1143513/49‑1    TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                            <  1:836883/50‑1      CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                           >  1:751828/1‑49      CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                           <  1:2458619/49‑1        CGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                        >  1:284845/1‑50         GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAA                                         <  1:1946778/48‑1            TACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGC                                     <  1:1107042/49‑1               GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                 <  1:936049/50‑1                  ATCATTCACTTGACATCAGAGATATACCATTCCATATAACGGCCCTGTA                               <  1:2594135/49‑1                   TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                              <  1:1658837/49‑1                      TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCT                          >  1:250433/1‑50                      TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCT                          <  1:2315193/50‑1                            TTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACG                     >  1:1500062/1‑49                             TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                   >  1:763135/1‑50                              ACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                   >  1:1027362/1‑49                               CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGAT                  >  1:2210179/1‑49                                ATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                 >  1:1124643/1‑49                                ATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                 >  1:1035576/1‑49                                 TCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAA               <  1:1910860/50‑1                                 TCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                 >  1:2693368/1‑48                                  CAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG              >  1:2458812/1‑50                                   AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC             >  1:2293292/1‑50                                     AGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGT           <  1:1603671/50‑1                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTT         <  1:1194531/50‑1                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT          >  1:2125934/1‑49                                        TATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTA        <  1:2184264/50‑1                                          TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACG      >  1:1229444/1‑50                                            CCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGC     <  1:1653313/49‑1                                              ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  <  1:1225413/50‑1                                              ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  <  1:1223613/50‑1                                              ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA   >  1:567954/1‑49                                                                                                 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167438‑167483 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173448‑173496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ 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.