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 =	6 (0.130)	25 (0.520)	24/484	0.9	80.4%	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	6 (0.130)	28 (0.590)	25/480	0.8	82.1%	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

ACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACC  >  NC_000913/2288698‑2289162                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    aCTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                                                                                                                         <  2:219737/250‑1 (MQ=255)   tGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGCATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGGCCCATAAACCAGCTTGATATCTTTTTTCTGGTCAACACCCTCAACTGTCTCACCATAATCAAAACTTGATGCCCATTCTATGGGTGCATATAAGGAAAGTGCTAATAAGCAGGGAAATTCTtctc                                                                                                                                                                                                                      >  2:125678/1‑251 (MQ=255)      gTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCACCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGTTGCCAATGCTATGGGTGAATATAAGGAAGGGGCGAATAAGCGGGGAAATTCTTCTCGGc                                                                                                                                                                                                                   >  2:477033/1‑251 (MQ=255)                             cAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                                                                                                                             <  1:253927/250‑1 (MQ=255)                                   aCTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                                                                                                                                      <  2:26521/251‑1 (MQ=255)                                         ccccTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGcc                                                                                                                                                                                 <  1:233288/249‑1 (MQ=255)                                           ccTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCCACTGGATATTTTTTTTCAGGACAACAACATCAACTGTGTCACCCTAATCAACAGTTGATGTCAATGCTATGGGAGAATATAAGGGAGGTGCGAATAAGCGGGGGAAATCTTTTTGGCTGACTTAGTCATTTCATttcttactttttgaggcgat                                                                                                                                                                              >  2:206175/1‑236 (MQ=255)                                               cGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAttttt                                                                                                                                                                          <  1:103754/251‑1 (MQ=255)                                                       tCAGTATTATTACTTACACCAAAGTCTTGTATATGCTTTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCTATTTTTTCTCCCg                                                                                                                                                                   <  2:73380/250‑1 (MQ=255)                                                         aGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGGTATCTTTTTCCAGGACAAAACCATCAACTGTCTCACCATACTCAACAGTTGATGCAAATGCTATTGGTGAATATAAGGAAGGTGCGAATAAGGGGGGGAATTCTTCTCGGCTGACTCTGTCATTTCATTTCTTCATGTTTGAGCCGCTTTTTTCTTCCGTa                                                                                                                                                                 >  1:268491/1‑250 (MQ=255)                                                            attattACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt                                                                                                                                                              <  1:54785/250‑1 (MQ=255)                                                              tattaCTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTTGTATTATTGGCGGTCCCATAAACCATCTGGATTTCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTTATACCAATGCTATGGGTGAATATAAGGAAGTTGCGAATAAGCGGTGAAATTCTTCGCGGCTGACTCAGTCATTTCATTTCTTCCTGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                                                                                                                            <  1:243736/250‑1 (MQ=38)                                                                 tacttacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATTAACCCGCTGGATATCTTTTTTCCGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCCATGCTATGGGGGAATATAAGGAAGGGGCGAATAAGCGGGGGAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTTAGCCGGTTTTTTCTCCCGTAAATGCCttt                                                                                                                                                        >  2:102564/1‑250 (MQ=37)                                                                  acttacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGa                                                                                                                                                       <  1:282364/251‑1 (MQ=255)                                                                     tacACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCTGTAAATGCCTTGaa                                                                                                                                                      <  1:119541/249‑1 (MQ=255)                                                                                      tatGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCGGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGtt                                                                                                                                   >  1:346795/1‑251 (MQ=255)                                                                                          cTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCCGCTGGGTATCTTTTTTCAGGACAAAACCATCCACTGTCTCACCATACTAAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGGATAAGCGGGGAAATTCTTCTCGGCTGACTCAGGCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGGATCCGCCTATTTAGACCGTttctt                                                                                                                               >  1:78792/1‑251 (MQ=25)                                                                                                             attTTCGTATTATTGGCGGTCCCATAAACCCGCTGGCTATCTTTTTCCAGGACAAAACCATCAACTGTCTCACCATACTCAACAGTTGTTGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAACGGGGGAATTCTTCTCGGCTGACTCAGTTATTTCATTTTTTCATGTTTTAGCCGATTTTTTCTCCCGTCAATGCCTTGAATCATCCTATTTTGACCGTTTCTTTTCCCTTTAAGGCGttttc                                                                                                            >  2:22308/1‑248 (MQ=255)                                                                                                                           ggCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTg                                                                                              <  1:362814/251‑1 (MQ=255)                                                                                                                                           aGCTGGCTATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTg                                                                               >  2:366151/1‑250 (MQ=255)                                                                                                                                                              ggACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGccc                                                           <  1:102564/251‑1 (MQ=255)                                                                                                                                                                      ccATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAgg                                                   >  1:426647/1‑251 (MQ=255)                                                                                                                                                                            aCTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGcc                                             >  2:273175/1‑251 (MQ=255)                                                                                                                                                                               gTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGc                                          <  1:206175/251‑1 (MQ=255)                                                                                                                                                                                                   tGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGtt                       <  1:10938/250‑1 (MQ=40)                                                                                                                                                                                                       gCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTa                  <  1:125678/251‑1 (MQ=38)                                                                                                                                                                                                                  gggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCa        <  1:350714/250‑1 (MQ=33)                                                                                                                                                                                                                       aaTATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAAcc  >  1:22436/1‑251 (MQ=31)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    ACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACC  >  NC_000913/2288698‑2289162

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