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

GGAAATTAACGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATT  >  NC_000913/2238195‑2238666                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        ttttgtagttgcgaatattggaaattaacgatttaaaaccaatatccagataggcataaattcccgttgagcggcgctcctcagttaccagtgcaaatccatggtttatggcttcgttggcattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACgag                                                                                                                                                                                                                                <  2:381775‑M2/10‑1 (MQ=255)           gctaatattggaaattaacgagttaaaaccaatatccagataggcataaattcccgttgagcggcgctcctcagttaccagtgcaaatccatggtttatggcttcgttggcattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGaa                                                                                                                                                                                                                      <  1:269635‑M2/20‑1 (MQ=255)           gcgaatattggaaattaacgagttaaaaccaatatccagataggcataaattcccgttgagcggcgctcctcagttaccagtgcaaatccatggtttatggcttcgttggcattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGaa                                                                                                                                                                                                                      <  1:188729‑M2/20‑1 (MQ=255)            cgaatattggaaattaacgagttaaaacaaatatacagataggcataaattgccgttgagcggcgctcctcagttaccagtgcaaatacatggtttatggcttcgttggaattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttagccaccacaagaccggcaatacagaATCAAACGGGACATCGCGAatg                                                                                                                                                                                                                    >  1:52534‑M2/230‑250 (MQ=255)                    ccAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCCCCTCAGTTACCATTGCAAATCCATGGTTTATGGCTTCGTGGGCATTATGGTTATTGATCTTTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCACAACAATACCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTg                                                                                                                                                                                                            <  2:352317/251‑1 (MQ=255)                           aacgagttaaaaccaatatccagataggcataaattcccgttgagcggcgctcctcagttaccagtgcaaatccatggtttatggcttcgttggcattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAg                                                                                                                                                                                                      <  1:281428‑M2/36‑1 (MQ=255)                                                                                  cagttaccagtgcaaatccatggtttatggcttcgttggcattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCa                                                                                                                                               <  1:6602‑M2/91‑1 (MQ=255)                                                                                             gcaaatccatggtttatggcttcgttggcattatggttattgatatgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCCAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCCGGTTACGTAACTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGCAACGCTgg                                                                                                                                   >  2:48502‑M2/149‑251 (MQ=255)                                                                                                          atatgggttcgttggcattatggttattgatctgtttgccgtgcgacgttatggggccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgccccccccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGGCGTCAGGTTACGTACCTCGAGGAGGACTCCGACCGGTTTGTTTTCTTTGTCAGGGAAACGATGGGTAAGAGAACGcc                                                                                                                      >  2:303056‑M2/136‑251 (MQ=255)                                                                                                            atggcttcgttggcattatggttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTAcc                                                                                                                                                                                   <  1:318603‑M2/55‑1 (MQ=255)                                                                                                                      tggcattatggttattgatatgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGCCTTAGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACAATCATgg                                                                                                          >  2:345347‑M2/124‑251 (MQ=255)                                                                                                                                gttattgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCtt                                                                                                >  1:74587‑M2/114‑251 (MQ=255)                                                                                                                                   attgatctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCCGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACACTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTc                                                                                             >  1:162362‑M2/111‑251 (MQ=255)                                                                                                                                   attgacctgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtgtcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTc                                                                                             <  2:74587‑M2/141‑1 (MQ=255)                                                                                                                                     tgatctgtttgccgtgcaacgttatggggccagccgatttctcgcgaataccaaataacgtctcaacaatatcgttacgtttctcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTcc                                                                                            <  2:324342‑M2/142‑1 (MQ=255)                                                                                                                                          tgtttgccgtgcaacgttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTc                                                                                      >  1:467764‑M2/104‑251 (MQ=255)                                                                                                                                                         gttatggtgccagccgatttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGt                                                                       <  1:236374‑M2/163‑1 (MQ=255)                                                                                                                                                           tatggtgccagcagatttcctgcgaatcccaaataacgtctcagcaaaagcggtacgtttcgcacccatcagaccggcaataccgaATCCAACGCGACATCGCGAATCGACGGCTGGCGCATTGACGTCAGGTTACGTACCTCGTGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGTGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCGGCCAGCGGtt                                                                      <  2:198553‑M2/164‑1 (MQ=255)                                                                                                                                                               gtgccagccgatttctcgcgaataccaaataacgtcgcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATTGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTg                                                                 <  1:128509‑M2/169‑1 (MQ=255)                                                                                                                                                                  ccagccgatttctcgagaatgcaaaaaaacgtctcaccaatatcggtacgtttcgccccaaccagaccggcaataccgaATCAAACGAGACATCGCGTATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGAATTCGCCCGGTTTGATTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCACAACCATCATGGCGATGATAATGTCCATAGTCAGTCATGCCGGCGGTTCggtgacg                                                              >  2:284622‑M2/80‑248 (MQ=255)                                                                                                                                                                         atttctcgcgaataccaaataacgtctcaacaatatcggtacgtttcgccccccccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACt                                                       <  1:366027‑M2/179‑1 (MQ=255)                                                                                                                                                                              tcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATgg                                                                                                          <  1:351943‑M2/128‑1 (MQ=255)                                                                                                                                                                              tcgcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATgg                                                                                                          >  2:351943‑M2/68‑195 (MQ=255)                                                                                                                                                                                gcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCACACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTg                                                                               >  2:331432‑M2/66‑220 (MQ=255)                                                                                                                                                                                gcgaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTg                                                                               <  1:331432‑M2/155‑1 (MQ=255)                                                                                                                                                                                  gaataccaaataacgtctcaacaatatcggtacgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTcgc                                              >  1:79433‑M2/64‑251 (MQ=255)                                                                                                                                                                                         aaataacgtctcaacaatagcggtacgtttcgcccccaccagaccggcaataccgaAGCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCGTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATAc                                       <  2:337243‑M2/195‑1 (MQ=255)                                                                                                                                                                                                    caacaatatcggtacgcttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTtca                             <  1:344170‑M2/205‑1 (MQ=255)                                                                                                                                                                                                      acaatatcggttcgtttcgcccccaccagaccggcaataccgaATCAAACGAGACATCGCGAATCGGCGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGTTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTtcatca                          <  2:50389‑M2/208‑1 (MQ=255)                                                                                                                                                                                                           atcgtgccgtttcgaccccaccagaccggcaataccgaATCAAACGCGACATCGCGAATCGACGGCTGGCGCAGTTACTGCAAGTTACGTACCTCGAGGAAGACGTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGATCGCCCAACCATCAGGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATaa                     <  2:431810‑M2/213‑1 (MQ=255)                                                                                                                                                                                                                           cccaccagaccggaaatcccgaATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTcc     <  2:192108‑M2/229‑1 (MQ=255)                                                                                                                                                                                                                              accagaccagcaagaccgaaTCACACGCGACATCGCGAATCGCCGGCTGGCGCAGAGACGTCAGGTTACGTACCAAGCGGATGACTTCGCCAGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATACTGGCCATGATCTTGTCCATCGTCAGTACTGCCAGCGGTTCGATGGCGATCCACTGACCGTCGCGCTATACGGTACCATCATCTCATAATAGGATGATTTCTTCCAtt  >  2:48601/20‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        GGAAATTAACGAGTTAAAACCAATATCCAGATAGGCATAAATTCCCGTTGAGCGGCGCTCCTCAGTTACCAGTGCAAATCCATGGTTTATGGCTTCGTTGGCATTATGGTTATTGATCTGTTTGCCGTGCAACGTTATGGTGCCAGCCGATTTCTCGCGAATACCAAATAACGTCTCAACAATATCGGTACGTTTCGCCCCCACCAGACCGGCAATACCGAGGATCTCCCCTTTATGCAGATCAAACGAGACATCGCGAATCGACGGCTGGCGCAGTGACGTCAGGTTACGTACCTCGAGGATGACTTCGCCCGGTTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATT  >  NC_000913/2238195‑2238666

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