Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,238,374	IS1 (–) +9 bp	coding (908‑916/1521 nt)	mglA ←	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2238374 =	0 (0.000)	38 (0.800)	37/486	0.3	100%	coding (916/1521 nt)	mglA	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
	?	NC_000913	3583428 =	NA (NA)					noncoding (1/768 nt)	IS1	repeat region
*	?	NC_000913	= 2238382	0 (0.000)	50 (1.050)	45/486	0.1	100%	coding (908/1521 nt)	mglA	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
	?	NC_000913	= 3584195	NA (NA)					noncoding (768/768 nt)	IS1	repeat region

ATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCG  >  NC_000913/3583200‑3583672                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     aTTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAt                                                                                                                                                                                                                                >  1:387895/1‑251 (MQ=255)                         tATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTg                                                                                                                                                                                                          >  1:169153/1‑249 (MQ=255)                           acTTACTTTGTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATGTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACGTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                                                                                                                                                                                                      <  2:292585/250‑1 (MQ=255)                                       aaGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGTTGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                                                                                                                                          >  2:179221/1‑251 (MQ=255)                                        agagCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                                                                                                                                         >  1:241015/1‑251 (MQ=255)                                                                    gATTACGATTATTTAGCCTCTCGATTCCAGAAACTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTc                                                                                                                                                             <  2:50354/251‑1 (MQ=255)                                                                       tACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAATCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtt                                                                                                                                                                                                              <  2:282328/199‑1 (MQ=255)                                                                       tACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGCTCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtt                                                                                                                                                                                                              >  1:282328/1‑199 (MQ=255)                                                                               tttAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGCAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCACCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTggcg                                                                                                                                                  >  1:70263/1‑248 (MQ=255)                                                                                  aGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                                                                                                                                               >  1:199112/1‑251 (MQ=255)                                                                                          gATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACACTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAACGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTCACGGGGGGGTGCgg                                                                                                                                       >  2:500324/1‑250 (MQ=255)                                                                                              ccAGAACCTTCAGTTACAAACCTTTTTGCCGTAAATGAATTTTATAGACATACTGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACg                                                                                                                                   >  2:172557/1‑251 (MQ=255)                                                                                                       tCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc                                                                                                                           <  2:169153/250‑1 (MQ=255)                                                                                                        cAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc                                                                                                                         >  1:158030/1‑251 (MQ=255)                                                                                                        cAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAc                                                                                                                         <  2:145492/251‑1 (MQ=255)                                                                                                                          ccGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACACTGATCTACAAGAACGTTTAGTGGTCGGGGATTATAGCATTTCAATATTTACTGATGACATTAAAGGTGATGCTGCCAACTTACTGATTTATTGTATGAGGGTGTGTTGGAGGTGCTCCAGTGGCTTCTTTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCTCAACGGCAAAAGCACCGCCGGTCATCAGCGCTa                                                                                                       <  2:285682/251‑1 (MQ=21)                                                                                                                              aaaTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTAtctc                                                                                                   >  1:355946/1‑251 (MQ=255)                                                                                                                                aTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGGGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTg                                                                                                 >  1:488978/1‑251 (MQ=255)                                                                                                                                                   tgatgaTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAAc                                                                               >  1:190898/1‑250 (MQ=255)                                                                                                                                                    gatgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGCCATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACAt                                                                             >  2:320920/1‑251 (MQ=255)                                                                                                                                                               aaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCGGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCa                                                                   <  1:270805/250‑1 (MQ=255)                                                                                                                                                                     tGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCACCTGCAGTTCACt                                                            >  1:61275/1‑251 (MQ=255)                                                                                                                                                                        gcTACAAGAAAGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTGAAGATGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGAGTAACGGCAAAAGCACCGCAGGACATCAGCGCGATCTCTGCTCTCACTGCCGAAAAACATGGCAACTGCAGTTCACTTa                                                          >  2:161137/2‑250 (MQ=17)                                                                                                                                                                               gAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCtt                                                  >  2:249428/1‑251 (MQ=255)                                                                                                                                                                                aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAATGATTTAGTGTATGTTGGTGTTTTTGAGGGGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTtc                                                 >  1:499651/1‑251 (MQ=255)                                                                                                                                                                                   ggTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCa                                              <  1:324690/251‑1 (MQ=255)                                                                                                                                                                                      tagggATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAGCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAAc                                            <  1:135415/246‑1 (MQ=255)                                                                                                                                                                                                     ataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGa                             >  1:113556/1‑250 (MQ=255)                                                                                                                                                                                                     ataGCATGTCAATATTTACTTATGACATTAAAGTTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGGTGTGCTACAGTGGCTTTTGTTTCTATCAGCTGTACATCATGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAACAGCACAGCCGTACCTCAGCGCTATCTCTGCTCTCATTTCCTTAAACAATGGCAACTGCAGTTCACTTCCACCGCTTATCACTCCGGTACGCACCAGa                             >  2:15186/1‑250 (MQ=14)                                                                                                                                                                                                      taGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACCCCGCTTCTCAACCCGGTACGCACCAGaa                            >  1:503845/1‑250 (MQ=255)                                                                                                                                                                                                                  atTTACTTATGACATTAAAGGTGATTCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATg               <  1:490185/251‑1 (MQ=38)                                                                                                                                                                                                                  atTTACTTATGACATTAAAGGTGATGCTACCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAAGGCAAAAGCACCGCAGGACATCAGCGATATCTCTGCTCTCACTGCCGTAAAACATGACAACTGCAGTTAACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATg               >  1:313009/1‑251 (MQ=255)                                                                                                                                                                                                                    ttACTTAGGGCATTAAAGTTGATTCTGCCAACTTACTGATTTAGGGTATGTTGGTTTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTACCATCTGTCCCTCCTGTGCAGCTACTTACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCCCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATgg              <  2:499651/250‑1 (MQ=11)                                                                                                                                                                                                                         tATGACATTAAAGGGGATGCTGCCAACTTACGGATTTAGTGTATGATGGTTTTTTTGAGGTGTGCCAGTGGCTTCGGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACGGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCTGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATg         <  2:113556/250‑1 (MQ=11)                                                                                                                                                                                                                               aTTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCg  >  1:217348/1‑251 (MQ=31)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCG  >  NC_000913/3583200‑3583672

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.