Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,238,948	Δ1 bp :: IS186 (–) +6 bp :: Δ2 bp	coding (337‑342/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	2238948 =	0 (0.000)	15 (0.590)	15/482	0.6	100%	coding (342/1521 nt)	mglA	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
	?	NC_000913	2514275 =	NA (NA)					noncoding (3/1345 nt)	IS186	repeat region
*	?	NC_000913	= 2238953	0 (0.000)	18 (0.710)	18/480	0.4	100%	coding (337/1521 nt)	mglA	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
	?	NC_000913	= 2515616	NA (NA)					noncoding (1344/1345 nt)	IS186	repeat region

GAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACT  >  NC_000913/2514028‑2514513                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               gAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCAt                                                                                                                                                                                                                                             >  1:252673/1‑251 (MQ=255)                       ccGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtgtg                                                                                                                                                                                                                       >  1:280473/1‑251 (MQ=255)                         gtgtgctgacggcttattctcTGTGTTTCTGGTTTCCGTCGCTAACTTCTCATCAATCGGTATTATCGCAGGTGCAGTTAATGGCCTGAAAGAGGAGCAAGGTTACGTGGTTTCTCGTTTCGGTCTTAAGCTGGTTTTCTGCTCTACCCT‑TGTGAGTGTGTGGTCTGCGTCAATCACAGCACTGGTGCTGTAAGACCATACATAAAAAATCCGAGGATAATTCCCATAATCGCTACCTTAAGggtggtggt                                                                                                                                                                                                                     <  2:22415/234‑1 (MQ=255)                         gcgTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                                                                                                                                                                                     <  1:8911/251‑1 (MQ=255)                                        catcTCTGTGTTACTGGTTTCCTTCGCTACCATATCTTCAATCAGGATTATCGCAGGTGCAGTTACAGGCCTGAATGCAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCTTCATTCGCAGCACTGGTGCTGTCAGACCAAACATGATAAAGCCGGGGCTACTTCCCATACGCGCTAACGTAAGGGTTGTGGTATTACGCATGCTATg                                                                                                                                                                                                      >  2:14327/1‑251 (MQ=255)                                                tgtTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                                                                                                                              >  1:78297/1‑251 (MQ=255)                                                        gTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                                                                                                                      >  2:14719/1‑251 (MQ=255)                                                            ccTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                                                                                                                   <  1:69074/250‑1 (MQ=255)                                                               tCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                                                                                                                                <  1:14719/250‑1 (MQ=255)                                                                cGCTAACTTCTCTTCAAGCGGGATTATCGCAGGGGCAGTTAAAGGCCTGAATGAAGAGCAAGGGAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACGTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                                                                                                               <  1:34360/250‑1 (MQ=255)                                                                   tAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGAGGTATTAGGCCTGATATGATTTAACGTGCAGTTGAATTACtctcc                                                                                                                                                                           >  2:72017/1‑250 (MQ=255)                                                                                                   ccgTTAAAGGCCTGAATGAAGCGAAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGCTATGATATAACGTGCCGATGCATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                                                                                                                           >  2:46756/3‑251 (MQ=255)                                                                                                         aaGGCCTTAATGAAGAGCAAGGTAACGTGTTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGTATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTTCGCCTGATATGATTTAACGTGCCGATGAATTACTCTCTCGATAACTGGTCAGCAATTCTGTCCCATATTGGTAAGc                                                                                                                                     <  1:157266/251‑1 (MQ=255)                                                                                                               tgaatgaaGGGCAGGGTAACGTGGTTTCTCGCTTCGGTCTGAATCTGGTTTACGGCTCTTCCCT‑GGTGAGTGTGCTGACTGCGTCAAGCGCAGCTCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                                                                                                                                <  2:16207/250‑1 (MQ=25)                                                                                                                  atgaaTAGCAAGGTAACGTGGTTTCTCGTTTCTGTCTTAAGCTGGTTTACGGCTCTACCCT‑GGTTAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGTCCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTTCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAAc                                                                                                                            <  1:279437/251‑1 (MQ=25)                                                                                                                                  cGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGCAGAACTGGATACTTCGGCACg                                                                                                            >  2:91558/1‑251 (MQ=255)                                                                                                                                                   gTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCTGGGCTCTa                                                                                            >  1:67679/1‑250 (MQ=255)                                                                                                                                                         aGCTGGGTGACGTCTCTACCCTTGGT‑AGTGTGCTGTCCGCGGCAATCGCAGCACGGGTGCTGTAATACCAGACATAAAAAAGCCGGGGATAATTCCCATAAGCTCTAACTTAAGTGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACGCTCACGATAGCTGGTCAGCCATTCTGGCCCATATTGTTAAGCCCCAATAACTGGATTCGTCGGCACGTAATGCCGGGGTTCTACCccgcc                                                                                     <  2:280473/251‑1 (MQ=17)                                                                                                                                                                ttACGGCTCTTCCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATAtt                                                                                                                                            <  2:178286/189‑1 (MQ=255)                                                                                                                                                                ttACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATAtt                                                                                                                                            >  1:178286/1‑189 (MQ=255)                                                                                                                                                                  aCGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGaa                                                                             >  1:112674/1‑250 (MQ=255)                                                                                                                                                                   cGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAAt                                                                           <  1:242719/251‑1 (MQ=255)                                                                                                                                                                   cGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAAt                                                                           >  1:178438/1‑251 (MQ=255)                                                                                                                                                                       tctACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTCAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGt                                                                       >  1:244997/1‑251 (MQ=255)                                                                                                                                                                                           ggccTGCGGCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATGCCCATTAGCGCTAACTTAAGGGTTGTGGTATGACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTAc                                                    <  2:244997/248‑1 (MQ=255)                                                                                                                                                                                                   tCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTgggg                                            >  1:155658/1‑251 (MQ=255)                                                                                                                                                                                                      aTCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGggctg                                         <  1:222117/251‑1 (MQ=255)                                                                                                                                                                                                                       ggaataCCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTGAAGGGTTGTGGTATTACTCCTGATATGATTTAACTTGCCGATGAATTACTCTCACGATAACGGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCTGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCggggg                        <  2:67679/246‑1 (MQ=255)                                                                                                                                                                                                                         aaGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCgggggg                       >  1:231675/1‑250 (MQ=255)                                                                                                                                                                                                                          aGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGAt                     >  1:148985/1‑251 (MQ=255)                                                                                                                                                                                                                                      aaaaaGCCGGGGATAATTCCCCTAAGCGCTAACTTAAGGGTTGGGGTATTACGCTTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGGCAGCAAGTCGGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATTTCATTACGTGa         <  2:148985/251‑1 (MQ=11)                                                                                                                                                                                                                                             cGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACt  >  1:266583/1‑251 (MQ=33)                                                                                                                                                                                                                                                  aTAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCt                                                <  1:62597/200‑1 (MQ=32)                                                                                                                                                                                                                                                  aTAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCt                                                >  2:62597/1‑200 (MQ=32)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               GAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCT‑GGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACT  >  NC_000913/2514028‑2514513

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