Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,090,347	Δ708 bp	coding (1‑708/708 nt)	endA →	DNA‑specific endonuclease I

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3090347	3091054	708	34 [0]	[0] 34	endA	DNA‑specific endonuclease I

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3090346	0 (0.000)	34 (1.130)	26/98	0.0	100%	intergenic (+94/‑1)	yggI/endA	conserved protein YggI/DNA‑specific endonuclease I
	?	NC_000913	3091055 =	0 (0.000)					intergenic (+1/‑79)	endA/rsmE	DNA‑specific endonuclease I/16S rRNA m(3)U1498 methyltransferase

CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3090298‑3090346 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTA  >  NC_000913/3091055‑3091099                                                                                                CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCC                                             <  1:1183629/50‑1    AACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTA                                           >  1:1913277/1‑49      CAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACAC                                        >  1:1865159/1‑50       AGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACAC                                        >  1:1336511/1‑49           TTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCG                                   >  1:1124583/1‑50           TTCGCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCG                                   >  1:202875/1‑50              GCTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGA                                >  1:1511848/1‑50               CTACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGAT                               <  1:20500/50‑1                TACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATT                              >  1:1260981/1‑50                 ACGTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTC                             >  1:1480530/1‑50                   GTTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTT                           <  1:1042813/50‑1                    TTGCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTT                          >  1:1544527/1‑50                      GCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTT                        >  1:679274/1‑50                      GCTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTT                         <  1:974142/49‑1                       CTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTT                        <  1:69237/49‑1                       CTGGCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTT                        <  1:1992066/49‑1                          GCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTA                    <  1:1588398/50‑1                          GCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTA                    <  1:641802/50‑1                          GCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTA                    <  1:2109256/50‑1                          GCTCGTTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTT                     <  1:385195/49‑1                               TTTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCC               >  1:1076157/1‑50                                TTTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCC               >  1:1787023/1‑49                                 TTAACACGGAG‑AAGTGCCTACGCTAGCGGGATTCTTTTTGTTAACCCCT              <  1:2119301/49‑1                                 TTAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCT              <  1:2208914/49‑1                                  TAACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTA             <  1:2031307/49‑1                                   AACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACC           >  1:408566/1‑50                                    ACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACC           <  1:1325966/49‑1                                    ACACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACC           <  1:1328159/49‑1                                     CACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCC         <  1:1275652/50‑1                                     CACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCC         <  1:2202855/50‑1                                      ACGGAG‑AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCA        <  1:24276/50‑1                                           agAAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCG    >  1:2164781/3‑50                                             AAGTGCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTA  >  1:1898039/1‑50                                              AGTGCCTACACTAGCGGGATTCTTTTTGTTAAACCCTACCCCACGCGTA  >  1:1207978/1‑49                                                                                                CAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3090298‑3090346 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTA  >  NC_000913/3091055‑3091099

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