Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,239,316	IS5 (–) +4 bp	intergenic (‑27/+31)	mglA ← / ← mglB	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit/D‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2239316 =	1 (0.070)	18 (1.160)	13/308	0.1	94.6%	intergenic (‑27/+34)	mglA/mglB	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit/D‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region
*	?	NC_000913	= 2239319	0 (0.000)	12 (0.780)	10/304	0.3	100%	intergenic (‑30/+31)	mglA/mglB	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit/D‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein
	?	NC_000913	2288919 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

GCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATA  >  NC_000913/2289865‑2290355                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   ncATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCATTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                                                                                                                                                                                                  >  1:91349/2‑251 (MQ=11)      ntttCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGGAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTAGCACCTTCCATAACGc                                                                                                                                                                                                                                             >  1:222638/2‑251 (MQ=17)                     ggTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTATGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                                                                                                                                                                                        >  1:79657/1‑241 (MQ=25)                              ccATTACCAGCATTGGGGTAAAAAGACACAATGACTTACACCATGTTTTACCCTGGCAGAATCTGCTACATGCGGGACAGGAAGCTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGCATTCACTGTGGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATATCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacca                                                                                                                                                                                                                                  >  1:286998/1‑238 (MQ=11)                                                        cTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTCCTCCATGCGGGCCAAGAAAATCTCTTTTCTAGTCTGCCGCCGCTTACTGCTGAATTAACTGTCGGCGCAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACACACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAAcc                                                                                                                                                                                             >  1:312102/1‑249 (MQ=17)                                                           gATGACTTCCACCATGTTTTGCCATGGCCGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGAtt                                                                                                                                                                                        >  1:272044/1‑251 (MQ=255)                                                                      ccnTGTNTNNCNNTGGNAGAATCTGCTCCATGCGGGACNNGAAAATCTCNTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGNTGACTCATGNNGANCCCTGNTNNNTNNNTCNAGATGACAAACATGATCTCNTANNAGGGACNTGTTCGCACCTTCCATAACGCTGTAg                                                                                                                                                                                                                                        <  2:79657/192‑1 (MQ=11)                                                                                                tCCATGCGGGACAAGAAAATCTCTTTTCTGATCTGACGGCGCTTACTGCCGAATTCACTGTCGACAAAGGTAAGTTGATCAATCATGATGCACACTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGCTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGCTATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTg                                                                                                                                                   >  1:148436/1‑251 (MQ=21)                                                                                                     gCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCCGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACATTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTaa                                                                                                                                               >  1:57979/1‑250 (MQ=255)                                                                                                                                                   aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTg                                                                                                                                       <  1:101427/212‑1 (MQ=255)                                                                                                                                                                 gAAGGTAAGTTGATGACTAATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                                                                                                                                                                                                    >  1:100458/1‑137 (MQ=255)                                                                                                                                                                     gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGATCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTaa                                                                               >  1:299541/1‑250 (MQ=255)                                                                                                                                                                         gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTcc                                                                                                               >  1:273357/1‑214 (MQ=255)                                                                                                                                                                               gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGCCCAACCTGATCTCCTATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCACAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAACGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGcac                                                                     >  1:198880/1‑250 (MQ=255)                                                                                                                                                                                            nnnnnnnnnntGNNNNNAGANGACAAACATGATCTCATAGCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAaattaat                                                                                                                                                        <  2:222638/144‑1 (MQ=38)                                                                                                                                                                                                        ncTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgttg                                           >  1:285164/2‑251 (MQ=255)                                                                                                                                                                                                            cAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTGTCATCCTATGAAAGTAATTTCTGTTAAAAGc                                                                                                                                           <  1:271129/151‑1 (MQ=255)                                                                                                                                                                                                                             ctcATATCAGGGACTTGTTCGCACCTTACATAACGCTGTAGCCACCAGAACAGATATTGCGGACCGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTaa                                                                               >  1:269802/1‑194 (MQ=255)                                                                                                                                                                                                                                  atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCa                  <  1:144976/250‑1 (MQ=255)                                                                                                                                                                                                                                        gACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTcc           >  1:107452/1‑251 (MQ=255)                                                                                                                                                                                                                                            tGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACATAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTGTCTTGATGAGAAGCTGATGCAAAATTCCGtn        <  1:36185/250‑2 (MQ=255)                                                                                                                                                                                                                                                 gCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTggn                                                                                                                                     <  1:87836/120‑2 (MQ=255)                                                                                                                                                                                                                                                 gCACCTTCCATAACGCTGTAGCCACCAGAACAGAGATTGCGGAGCGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTg                                                                                                                                                   <  1:303974/106‑1 (MQ=255)                                                                                                                                                                                                                                                  cACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATACTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTtata  >  1:10457/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   GCATGGTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATA  >  NC_000913/2289865‑2290355

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