Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,375,828	IS186 (+) +4 bp	coding (132‑135/963 nt)	menC ←	o‑succinylbenzoate synthase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2375828 =	0 (0.000)	25 (1.010)	24/484	0.2	100%	coding (135/963 nt)	menC	o‑succinylbenzoate synthase
	?	NC_000913	= 2515617	NA (NA)					noncoding (1345/1345 nt)	IS186	repeat region
*	?	NC_000913	= 2375831	1 (0.040)	18 (0.730)	18/484	0.5	94.7%	coding (132/963 nt)	menC	o‑succinylbenzoate synthase
	?	NC_000913	2514273 =	NA (NA)					noncoding (1/1345 nt)	IS186	repeat region

TTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCA  >  NC_000913/2515372‑2515847                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       ttAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAtt                                                                                                                                                                                                                                   >  1:176385/1‑251 (MQ=31)            cATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGATTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                                                                                                                                                                                                        <  2:177761/251‑1 (MQ=25)                  ccAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACg                                                                                                                                                                                                                   <  1:54804/250‑1 (MQ=38)                     gCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                                                                                                                                                                                                               >  1:236929/1‑251 (MQ=40)                                                      aTCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaa                                                                                                                                                                              <  2:76011/251‑1 (MQ=255)                                                        ccGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa                                                                                                                                                                             <  2:50081/250‑1 (MQ=255)                                                                                                                           ttCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTc                                                                                                         >  1:157573/1‑251 (MQ=255)                                                                                                                                              gccAGCCTACAAAACCCAATTGTACACCGCGAGGCATCGTCTTAATGAGGAGCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTGAAGTGACCGCTGTTGGGATAATTCCCGGGTTTTTGTACGCCTGGTAATCAGCTAATGGTGGCGGGCGACCAATCAAATCCCCTCGTATATATGTCCCCCCGAGCTTATGCATCAGCGCCTTCTCCTTTTGGGTCTTTGCTAACTCCGCGACGACCCTCAATg                                                                                        <  2:103382/247‑1 (MQ=12)                                                                                                                                                                                        tAACGAGGGACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACTCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGc                                             <  2:133212/250‑1 (MQ=255)                                                                                                                                                                                        tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCg                                              <  1:34834/249‑1 (MQ=255)                                                                                                                                                                                                   cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTATGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGtt                                 >  1:59985/1‑251 (MQ=255)                                                                                                                                                                                                    gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTc                                <  2:49329/251‑1 (MQ=255)                                                                                                                                                                                                     aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCg                               <  1:123393/251‑1 (MQ=255)                                                                                                                                                                                                                                   ttAAGTTAGCGCTTATGGGATAATTCCCAGGTTTTTTTACGCCTGTTAATCAGCTACTGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCa  >  2:12689/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       TTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCTTGCAGATATTGCACCCCGAGCTTATGCAATAGCGCCTGCTGCTGTTGCGTCTCGACAAACTCCGCGACCACACTCAATGACTTCGCTTTCGCCAGATCGGTAATTGATCTCACAATCATCGCATCCAGCGTGTTCGTGACAATATCTTTCACAAAGACGCCATCA  >  NC_000913/2515372‑2515847

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