Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_pae_pgi	4,319,069	(CATCGC)2→1	coding (850‑855/1065 nt)	phnI ←	ribophosphonate triphosphate synthase complex putative catalytic subunit

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_pgi	4319060	4319065	6	12 [0]	[0] 12	phnI	ribophosphonate triphosphate synthase complex putative catalytic subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	4319066 =	0 (0.000)	12 (0.800)	12/234	0.4	100%	coding (858/1065 nt)	phnI	ribophosphonate triphosphate synthase complex putative catalytic subunit
	?	NC_000913_3_pae_pgi	= 4319059	0 (0.000)					coding (865/1065 nt)	phnI	ribophosphonate triphosphate synthase complex putative catalytic subunit

GGGAGTTTGAGGTGTGAGACAAAGCCTGCGGCTTCGACGTTGTCGGCATGTGCCAGCACGAACTCTTCATCCTGCGCCGGGCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCG  >  NC_000913_3_pae_pgi/4318929‑4319186                                                                                                                                   |                                                                                                                                gggAGTTTGAGGTGTGAGACAAAGCCTGCGGCTTCGACGTTGTCGGCATGTGCCAGCACGAACTCTTCATCCTGCGCCGGGCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcg                                                                                                              <  1:133990‑M2/150‑20 (MQ=255)      tttGAGGTGTGAGACAAAGCCTGCGGCTTCGACGTTGTCGGCATGTGCCAGCACGAACTCTTCATCCTGCGCCGGGCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcat                                                                                                         <  1:162176‑M2/150‑25 (MQ=255)                                       gTTGTCGGCATGTGCCAGCACGAACTCTTCATCCTGCGCCGGGCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtg                                                                        <  2:104660‑M2/150‑58 (MQ=255)                                             ggCATGTGCCAGCACGAACTCTTCATCCTGCGCCGGGCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcgg                                                                  >  2:114508‑M2/1‑87 (MQ=255)                                                                                ggCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgataaaaccgttaaccattt                               >  2:133990‑M2/1‑52 (MQ=255)                                                                                               cGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgataaaaccgttaaccatttcacactcagtcatca                >  1:104660‑M2/1‑37 (MQ=255)                                                                                                  cGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgat                                                <  3:93500‑M2/115‑82 (MQ=255)                                                                                                  cGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgat                                                >  4:93500‑M2/1‑34 (MQ=255)                                                                                                            ggCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgataaaaccgttaaccatttcacactcagtcatcagcagttcgccgac   >  1:130245‑M2/1‑24 (MQ=255)                                                                                                             gCTTGCAGAGCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgataaaaccgttaaccatttcacactcagtcatcagcagttcgccgacg  >  2:162176‑M2/1‑23 (MQ=255)                                                                                                                   agagCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgataaaaccgttaac                                    >  3:40958‑M2/1‑17 (MQ=255)                                                                                                                   agagCGCGGTCGACCAGcgccatcgctttgcgctcgctcatgccgaacaccagcccgtagccgcgcgtgaagtgcggcggctcacccggcgggtcgataaaaccgttaac                                    <  4:40958‑M2/110‑94 (MQ=255)                                                                                                                                   |                                                                                                                                GGGAGTTTGAGGTGTGAGACAAAGCCTGCGGCTTCGACGTTGTCGGCATGTGCCAGCACGAACTCTTCATCCTGCGCCGGGCCTGTCGCGTGCTCGCCGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCG  >  NC_000913_3_pae_pgi/4318929‑4319186

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