Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,184,782	Δ6,916 bp		ubiK–yqiI	ubiK, ygiL, yqiG, insC‑5, insD‑5, yqiG, yqiH, yqiI

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3184782	3191697	6916	36 [0]	[0] 36	ubiK–yqiI	ubiK,ygiL,yqiG,insC‑5,insD‑5,yqiG,yqiH,yqiI

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3184781	0 (0.000)	35 (1.300)	28/94	0.0	100%	intergenic (‑315/‑59)	ribB/ubiK	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/ubiquinone biosynthesis accessory factor UbiK
	?	NC_000913	3191698 =	0 (0.000)					intergenic (+2/+41)	yqiI/glgS	putative fimbrial protein YqiI/surface composition regulator

GACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTG  >  NC_000913/3191653‑3191744                                              |                                               atgattagcggcagcaacgcatagcttcacataattctggtttatATCt                                             <  1:1455265‑M2/4‑1 (MQ=255)  tgattagcggcagcaacgcatagcttcacataattctggtttatATCt                                             >  1:1107642‑M2/45‑48 (MQ=255)    attagcggcagcaacgcatagcttcacataattctggtttatATCTAAt                                          <  1:3663‑M2/7‑1 (MQ=255)     ttagcggcagcaacgcatagcttcacataattctggtttatATCTAATa                                         <  1:1784341‑M2/8‑1 (MQ=255)      tagcggcagcaacgcatagcttcacataattctggtttatATCTAATa                                         >  1:1142778‑M2/41‑48 (MQ=255)      tagcggcagcaacgcatagcttcacataattctggtttatATCTAATa                                         >  1:125194‑M2/41‑48 (MQ=255)        gcggcagcaacgcatagcttcacataattctggtttatATCTAATACa                                       <  1:1902354‑M2/10‑1 (MQ=255)         cggcagcaacgcatagcttcacataattctggtttatATCTAATACaaa                                     <  1:1269833‑M2/12‑1 (MQ=255)          ggcagcaacgcatagcttcacataattctggtttatATCTAATACaaa                                     >  1:551287‑M2/37‑48 (MQ=255)          ggcagcaacgcatagcttcacataattctggtttatATCTAATACAAac                                    <  1:1802611‑M2/13‑1 (MQ=255)          ggcagcaacgcatagcttcacataattctggtttatATCTAATACAAac                                    <  1:486193‑M2/13‑1 (MQ=255)           gcagcaacgcatagcttcacataattctggtttatATCTAATACAAaca                                   <  1:854038‑M2/14‑1 (MQ=255)            cagcaacgcatagcttcacataattctggtttatATCTAATACAAacac                                  >  1:1980559‑M2/35‑49 (MQ=255)            cagcaacgcatagcttcacataattctggtttatATCTAATACAAacac                                  >  1:930976‑M2/35‑49 (MQ=255)              gcaacgcatagcttcacataattctggtttatATCTAATACAAACACTa                                <  1:155225‑M2/17‑1 (MQ=255)                aacgcatagcttcacataattctggtttatATCTAATACAAACACTaaa                              <  1:356582‑M2/19‑1 (MQ=255)                aacgcatagcttcacataattctggtttatATCTAATACAAACACTaaa                              <  1:952673‑M2/19‑1 (MQ=255)                    catagcttcacataattctggtttatATCTAATACAAACACTAAAACgg                          <  1:1967735‑M2/23‑1 (MQ=255)                    catagcttcacataattctggtttatATCTAATACAAACACTAAAACgg                          >  1:1145395‑M2/27‑49 (MQ=255)                    catagcttcacataattctggtttatATCTAATACAAACACTAAAACgg                          >  1:301642‑M2/27‑49 (MQ=255)                     atagcttcacataattctggtttatATCTAATACAAACACTAAAACggg                         >  1:604284‑M2/26‑49 (MQ=255)                      tagcttcacataattctggtttatATCTAATACAAACACTAAAACGGGc                        <  1:1220741‑M2/25‑1 (MQ=255)                      tagcttcacataattctggtttatATCTAATACAAACACTAAAACGGGc                        >  1:79209‑M2/25‑49 (MQ=255)                       agcttcacataattctggtttatATCTAATACAAACACTAAAACGGGc                        >  1:1269895‑M2/24‑48 (MQ=255)                        gcttcacataattctggtttatATCTAATACAAACACTAAAACGGGcc                       <  1:733732‑M2/26‑1 (MQ=255)                           tcacataattctggtttatATCTAATACAAACACTAAAACGGGCCAt                     >  1:226418‑M2/20‑47 (MQ=255)                            cacataattctggtttatATCTAATACAAACACTAAAACGGGCCATCa                   <  1:797194‑M2/30‑1 (MQ=255)                               ataattctggtttatATCTAATACAAACACTAAAACGGGCCATCAGGcc               <  1:850261‑M2/34‑1 (MQ=255)                                  attctggcttatATCTAATACAAACACTAAAACGGGCCATCAGGCCCGt            >  1:1774499‑M2/13‑49 (MQ=255)                                   ttctggtttatATCTAATACAAACACTAAAACGGGCCATCAGGCCCGtt           <  1:481676‑M2/38‑1 (MQ=255)                                   ttctggtttatATCTAATACAAACACTAAAACGGGCCATCAGGCCCGt            >  1:1524626‑M2/12‑48 (MQ=255)                                   ttctggtttatATCTAATACAAACACTAAAACGGGCCATCAGGCCCGt            >  1:1722272‑M2/12‑48 (MQ=255)                                     ctggtttatATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTta          >  1:1039740‑M2/10‑48 (MQ=255)                                         tttatATCTAATACAAACACTAAAACTGGCCATCAGGCCCGTTatat       >  1:1195540‑M2/6‑47 (MQ=255)                                         tttatATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTatat       >  1:1566472‑M2/6‑47 (MQ=255)                                             tATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTg  <  1:851554‑M2/47‑1 (MQ=255)                                              |                                               GACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTG  >  NC_000913/3191653‑3191744

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.