Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	2,238,416	Δ19 bp	coding (856‑874/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	= 2238415	1 (0.020)	30 (0.690)	28/486	0.5	96.8%	coding (875/1521 nt)	mglA	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
	?	NC_000913	2238435 =	1 (0.020)					coding (855/1521 nt)	mglA	D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit

TTTTGTAGTTGCGAATATTGGAAATTAACGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACA  >  NC_000913/2238176‑2238644                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      ttttGTAGTTGCGAATATTGGAAATTAACGATTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgag                                                                                                                                                                                                                             <  2:381775‑M1/250‑11 (MQ=255)           gctaatatTGGAAATTAACGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaa                                                                                                                                                                                                                   <  1:269635‑M1/247‑21 (MQ=255)           gCGAATATTGGAAATTAACGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaa                                                                                                                                                                                                                   <  1:188729‑M1/250‑21 (MQ=255)            cGAATATTGGAAATTAACGAGTTAAAACAAATATACAGATAGGCATAAATTGCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATACATGGTTTATGGCTTCGTTGGAATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTAGCCACCACAAGACCGGCAATACAGAatcaaacgggacatcgcgaatg                                                                                                                                                                                                                 >  1:52534‑M1/1‑229 (MQ=255)                           aaCGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcag                                                                                                                                                                                                   <  1:281428‑M1/250‑37 (MQ=255)                                       ccAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCCCCTCAGTTACCATTGCAAATCCATGGTTTATGGCTTCGTGGGCATTATGGTTATTGATCTTTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCACAACAATACCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTg                                                                                                                                                                                      <  2:352317/251‑1 (MQ=255)                                                                                  cAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtca                                                                                                                                            <  1:6602‑M1/250‑92 (MQ=255)                                                                                             gCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATATGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatccaacgagacatcgcgaatcgacggctggcgcagtgacgtccggttacgtaactcgaggatgacttcgcccggtttgttttctttgtcagggcaacgctgg                                                                                                                                >  2:48502‑M1/1‑148 (MQ=255)                                                                                                          ataTGGGTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCGACGTTATGGGGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCCCCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtggcgtcaggttacgtacctcgaggaggactccgaccggtttgttttctttgtcagggaaacgatgggtaagagaacgcc                                                                                                                   >  2:303056‑M1/2‑135 (MQ=255)                                                                                                            aTGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacc                                                                                                                                                                                <  1:318603‑M1/188‑56 (MQ=255)                                                                                                                      tGGCATTATGGTTATTGATATGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgccttagcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaacaatcatgg                                                                                                       >  2:345347‑M1/1‑123 (MQ=255)                                                                                                                                gTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatctt                                                                                             >  1:74587‑M1/1‑113 (MQ=255)                                                                                                                                   aTTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtccggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacactggttaagagaacgcccaaccatcatggcgatgatcttgtc                                                                                          >  1:162362‑M1/1‑110 (MQ=255)                                                                                                                                   aTTGACCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTGTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtc                                                                                          <  2:74587‑M1/251‑142 (MQ=255)                                                                                                                                     tGATCTGTTTGCCGTGCAACGTTATGGGGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGTTACGTTTCTCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtcc                                                                                         <  2:324342‑M1/250‑143 (MQ=255)                                                                                                                                          tGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtc                                                                                   >  1:467764‑M1/1‑103 (MQ=255)                                                                                                                                                         gTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctgccagcggt                                                                    <  1:236374‑M1/251‑164 (MQ=255)                                                                                                                                                           tATGGTGCCAGCAGATTTCCTGCGAATCCCAAATAACGTCTCAGCAAAAGCGGTACGTTTCGCACCCATCAGACCGGCAATACCGAatccaacgcgacatcgcgaatcgacggctggcgcattgacgtcaggttacgtacctcgtggatgacttcgcccggtttgttttctttgtcagtgaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtccggccagcggtt                                                                   <  2:198553‑M1/250‑165 (MQ=255)                                                                                                                                                               gTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCGCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcattgcgatgatcttgtccatcgtcagtcctgccagcggttcggtg                                                              <  1:128509‑M1/251‑170 (MQ=255)                                                                                                                                                                  ccAGCCGATTTCTCGAGAATGCAAAAAAACGTCTCACCAATATCGGTACGTTTCGCCCCAACCAGACCGGCAATACCGAatcaaacgagacatcgcgtatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgaattcgcccggtttgatttctttgtcagggaaacgctggttaagagaacgcacaaccatcatggcgatgataatgtccatagtcagtcatgccggcggttcggtgacg                                                           >  2:284622‑M1/1‑79 (MQ=255)                                                                                                                                                                         aTTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCCCCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctgccagcggttcggtggcgatccact                                                    <  1:366027‑M1/251‑180 (MQ=255)                                                                                                                                                                              tcGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatgg                                                                                                       <  1:351943‑M1/195‑129 (MQ=255)                                                                                                                                                                              tcGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatgg                                                                                                       >  2:351943‑M1/1‑67 (MQ=255)                                                                                                                                                                                gcgAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcacacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctg                                                                            >  2:331432‑M1/1‑65 (MQ=255)                                                                                                                                                                                gcgAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctg                                                                            <  1:331432‑M1/220‑156 (MQ=255)                                                                                                                                                                                  gAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctgccagcggttcggtggcgatccactgaccgtcgc                                           >  1:79433‑M1/1‑63 (MQ=255)                                                                                                                                                                                         aaaTAACGTCTCAACAATAGCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAagcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcgtgtccatcgtcagtcctgccagcggttcggtggcgatccactgaccgtcgcgcaatac                                    <  2:337243‑M1/251‑196 (MQ=255)                                                                                                                                                                                                    caacaaTATCGGTACGCTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctgccagcggttcggtggcgatccactgaccgtcgcgcaatacggtaacttca                          <  1:344170‑M1/250‑206 (MQ=255)                                                                                                                                                                                                      acaaTATCGGTTCGTTTCGCCCCCACCAGACCGGCAATACCGAatcaaacgagacatcgcgaatcggcggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctgccagcgtttcggtggcgatccactgaccgtcgcgcaatacggtaacttcatca                       <  2:50389‑M1/251‑209 (MQ=255)                                                                                                                                                                                                           atcgtgccgtttcgaCCCCACCAGACCGGCAATACCGAatcaaacgcgacatcgcgaatcgacggctggcgcagttactgcaagttacgtacctcgaggaagacgtcgcccggtttgttttctttgtcagggaaacgctggttaagagatcgcccaaccatcagggcgatgatcttgtccatcgtcagtcctgccagcggttcggtggcgatccactgaccgtcgcgcaatacggtaacttcatcacataa                  <  2:431810‑M1/244‑214 (MQ=255)                                                                                                                                                                                                                           cccACCAGACCGGAAATCCCGAatcaaacgagacatcgcgaatcgacggctggcgcagtgacgtcaggttacgtacctcgaggatgacttcgcccggtttgttttctttgtcagggaaacgctggttaagagaacgcccaaccatcatggcgatgatcttgtccatcgtcagtcctgccagcggttcggtggcgatccactgaccgtcgcgcaatacggtaacttcatcacataactggaagatttcttcc  <  2:192108‑M1/251‑230 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                      TTTTGTAGTTGCGAATATTGGAAATTAACGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACA  >  NC_000913/2238176‑2238644

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