Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,994,992	Δ3,877 bp		[ygeL]–ygeQ	[ygeL], ygeN, ygeO, insD‑4, insC‑4, ygeP, ygeQ

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2994992	2998868	3877	31 [0]	[0] 32	[ygeL]–ygeQ	[ygeL],ygeN,ygeO,insD‑4,insC‑4,ygeP,ygeQ

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2994991	0 (0.000)	30 (1.140)	28/94	0.0	100%	pseudogene (102/156 nt)	ygeL	protein YgeL
	?	NC_000913	2998869 =	0 (0.000)					intergenic (‑41/+115)	ygeQ/glyU	protein YgeQ/tRNA‑Gly

CCCTTTCATAAAATATTCCATCAATCGTTTGTTAAAAGCGTAGATAACCACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAA  >  NC_000913/2998820‑2998918                                                  |                                                  tgtagatttgttacactccgctatagctttcttcagatctttatagccca                                                   >  1:2302770‑M2/50‑50 (MQ=255)  gtagatttgttacactccgctatagctttcttcagatctttatagcccAc                                                  >  1:699471‑M2/49‑50 (MQ=255)   tagatttgttacactccgctatagctttcttcagatctttatagccca                                                   <  1:1893077‑M2/1‑1 (MQ=255)   tagatttgttacactccgctatagctttcttcagatctttatagcccACt                                                 >  1:140452‑M2/48‑50 (MQ=255)    agatttgttacactccgctatagctttcttcagatctttatagcccACt                                                 >  1:433928‑M2/47‑49 (MQ=255)     gatttgttacactccgctatagctttcttcagatctttatagcccACTTa                                               <  1:31956‑M2/5‑1 (MQ=255)      atttgttacactccgctatagctttcttcagatctttatagcccACTTAt                                              >  1:541379‑M2/45‑50 (MQ=255)      atttgttacactccgctatagctttcttcagatctttatagcccACTTAt                                              <  1:2205511‑M2/6‑1 (MQ=255)          gttacactccgctatagctttcttcagatctttatagcccACTTATTTa                                           <  1:811536‑M2/9‑1 (MQ=255)            tacactccgctatagctttcttcagatctttatagcccACTTATTTaac                                         >  1:60328‑M2/39‑49 (MQ=255)             acactccgctatagctttcttcagatctttatagcccACTTATTTaaca                                        <  1:1392208‑M2/12‑1 (MQ=255)              cactccgctatagctttcttcagatctttatagcccACTTATTTaacaac                                      <  1:1188460‑M2/14‑1 (MQ=255)              cactccgctatagctttcttcagatctttatagcccACTTATTTaacaac                                      <  1:656033‑M2/14‑1 (MQ=255)              cactccgctatagctttcttcagatctttatagcccACTTATTTaacaa                                       <  1:2189622‑M2/13‑1 (MQ=255)               actccgctatagctttcttcagatctttatagcccACTTATTTAACAACt                                     >  1:983801‑M2/36‑50 (MQ=255)                ctccgctatagctttcttcagatctttatagcccACTTATTTAACAACtt                                    <  1:840133‑M2/16‑1 (MQ=255)                   cgctatagctttcttcagatctttatagcccACTTATTTAACAACTTGaa                                 <  1:1357440‑M2/19‑1 (MQ=255)                    gctatagctttcttcagatctttatagcccACTTATTTAACAACTTGaa                                 >  1:1219797‑M2/31‑49 (MQ=255)                     ctatagctttcttcagatctttatagcccACTTATTTAACAACTTGaaaa                               <  1:1409820‑M2/21‑1 (MQ=255)                         agctttcttcagatctttatagcccACTTATTTAACAACTTGAAAATAAg                           >  1:901107‑M2/26‑50 (MQ=255)                            tttcttcagatctttatagcccACTTATTTAACAACTTGAAAATAAGgtg                        >  1:1571173‑M2/23‑50 (MQ=255)                               cttcagatctttatagcccACTTATTTAACAACTTGAAAATAAGGTGtt                      <  1:2322872‑M2/30‑1 (MQ=255)                                  cagatctttatagcccACTTATTTAACAACTTGAAAATAAGGTGTTTTa                   <  1:610394‑M2/33‑1 (MQ=255)                                   agatctttatagcccACTTATTTAACAACTTGAAAATAAGGTGTTTTAcc                 <  1:2317584‑M2/35‑1 (MQ=255)                                       ctttatagcccACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTggg             >  1:978882‑M2/12‑50 (MQ=255)                                          tatagcccACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGgtt           <  1:1994203‑M2/41‑1 (MQ=255)                                           atagcccACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGgttgt         >  1:1026641‑M2/8‑50 (MQ=255)                                           atagcccACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGgttgt         >  1:67000‑M2/8‑50 (MQ=255)                                             agcccACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGGTTGTTa       >  1:815421‑M2/6‑50 (MQ=255)                                              gcccACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGGTTGTTa       <  1:939733‑M2/45‑1 (MQ=255)                                                 cACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACaa    >  1:2014706‑M2/2‑49 (MQ=255)                                                  aCTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACaaaa  <  1:1193023/50‑1 (MQ=255)                                                  |                                                  CCCTTTCATAAAATATTCCATCAATCGTTTGTTAAAAGCGTAGATAACCACTTATTTAACAACTTGAAAATAAGGTGTTTTACCTGGGTTGTTACAAAA  >  NC_000913/2998820‑2998918

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.