Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	2,239,575	Δ5 bp	coding (770‑774/999 nt)	mglB ←	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	= 2239574	3 (0.080)	29 (0.750)	28/480	0.4	90.6%	coding (775/999 nt)	mglB	D‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein
	?	NC_000913	2239580 =	3 (0.080)					coding (769/999 nt)	mglB	D‑galactose/methyl‑galactoside ABC transporter periplasmic binding protein

AGACTGCGCCCAATCAGTCTTATTTCTTGCTGAATTCAGCCAGGTTGTCTTTATCTACGCCAACATAAGGTACGCGGACCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTT‑GTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGC  >  NC_000913/2239331‑2239817                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     aGACTGCGCCCAATCAGTCTTATTTCTTGCTGAATTCAGCCAGGTTGTCTTTATCTACGCCAACATAAGGTACGCGGACCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacac                                                                                                                                                                                                                                               >  1:155302‑M1/1‑244 (MQ=255)            aaTCAGTCTTATTTCTTGCTGAATTCAGCCAGGTTGTCTTTATCTACGCCAACATAAGGTACGCGGACACCTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGGACGGTGCCC‑GCCAGTGCACCGGCTTTCACCAGCGCAAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctgg                                                                                                                                                                                                                                    >  2:410805‑M1/1‑233 (MQ=255)                aGTCTTATTTCTTGCTGAATTCAGCCAGGTTGTCTTTATCTACGCCAACATAAGGTACGCGGACCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggactt                                                                                                                                                                                                                                >  1:15698‑M1/1‑229 (MQ=255)                             tctGAATTCAGCCAGGTGGTTTTTTACTCCCCCAACATAAGGTACGCGACCCTCTTTTTTGCCGTTTTTCCAGTTGGTGCCATCA‑GCCGCACTTTTACCACCGGCCAGGTTTTCC‑GCCAAATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTCCAGTCCGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgcttt                                                                                                                                                                                                                    <  1:324683‑M1/249‑35 (MQ=255)                             gCTGAATTCAGCCAGGTTGTCTTTATCTACGCCAACATAAGGTACGCGGACCAATTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCAAGGTTTTTC‑GCCAGGTCCATGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCAACcaaacaccggaatgctggacctgttgtgtgctttc                                                                                                                                                                                                                   >  2:300938‑M1/1‑216 (MQ=255)                                                           gCCAACATAAGGTACGCGGACCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCCTCcaaccaccggaatgctggacttgttgtgtgctttaggcgcttcaaccgcgcccattgccatcccc                                                                                                                                                                                     >  2:82999‑M1/1‑186 (MQ=255)                                                                 aTAAGGTACGCGGCCCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATAAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCccaacaccggaatgctggacttgttgtgtgctttcagcgcttccaccgcccccattgccatcccttcgttt                                                                                                                                                                               >  2:387900‑M1/1‑180 (MQ=255)                                                                          gcgGACCCCTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggactgctggact                                                                                                                                                                                                                                 >  2:377268‑M1/1‑171 (MQ=255)                                                                             gACCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCCCCTTTACCATCGGCCAGGGTTTTC‑GCCAGATCAACGGACGCTTTCGCCTGGTTGTTACCATCGTTCAGTACGCTGCCC‑GCAAGTGCTCCGGATTTCACCAGAGCCAGCTCTTCTGGCAACGCCTC‑A‑G‑‑‑ACCACCGGACTGCTGGAATT‑GTTGTGTGCTTACAGCGATTCAACCACGCCCATTGCCATCGCATCgttgatagcggatccc                                                                                                                                                             >  2:230158/1‑239 (MQ=255)                                                                                  cTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGGGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttc                                                                                                                                                              >  1:338140‑M1/1‑163 (MQ=255)                                                                                   tttGTTGTCGGTTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGCCGCATTAGCCTGGTTGTTAGAATCGTTCAGTACGCTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCAAGCGCTTCTGGCAGCGCATCcaaacaccgaaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgtggttggcgctaaccacttcg                                                                                                                                                             >  1:272558‑M1/1‑162 (MQ=255)                                                                                      gttgaCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACAATCTGCCAGGATGTTC‑GACAGATCAAAGGTAGCTTTAACCTGGTTGTTAGAGTCGTTCAGTACGGTGACC‑GCCAATGAACCGGATTTCACCAGCACCAACGCATCTGGCAGCGCAT‑‑‑‑‑CCAAACACCGGAAAGCTGGACTTAGACAT‑TGCTTTCAGAGAGTCAACCGCACACATTTCCATCTCCTCGTTGTTGGCGATCACCACtgcgttt                                                                                                                                                    >  2:221644/1‑245 (MQ=255)                                                                                       ttgtCGTTTTCCCAGTTGGTCCCA‑CACGACGCATCTTTACCATCG‑CCAAGTTTTTCCGCCAGATAAAAGGTACCTCTCGCCTTGTTGTTATCTTCGTTCATTTCG‑TGCCCCGCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCAT‑‑‑‑‑CCAAACACCGGAATGCTGGACTT‑GTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGAtt                                                                                                                                                    <  1:151465/250‑1 (MQ=255)                                                                                          tCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgt                                                                                                                                                      >  2:82868‑M1/1‑155 (MQ=255)                                                                                                                                     gCCAGGTTTTTC‑GCCAGATCAAAGGTAGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGCTTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcg                                                                                                            >  1:345178‑M1/1‑113 (MQ=255)                                                                                                                                      ccAGGTTTTTC‑GCAAGATCAAAGGACGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGACC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgattcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagaaagccaggagtccatcttatctttcgc                                                                                                           >  2:132785‑M1/1‑112 (MQ=255)                                                                                                                                                          aaaGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgctcattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacatt                                                                                        >  1:159011‑M1/1‑93 (MQ=255)                                                                                                                                                                    ggacgtggtctcttccaatctttcagtccggtgccc‑gccagtgcACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatcta                                                                              <  1:1487‑M1/230‑169 (MQ=255)                                                                                                                                                                                ttagcataatacggcacggtgccc‑gccagtgCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcggccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgt                                                                   <  1:410805‑M1/235‑180 (MQ=255)                                                                                                                                                                                        gtccagtACGGGGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttg                                                          <  1:295469‑M1/248‑189 (MQ=255)                                                                                                                                                                                         ttCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatcta                                                                              >  2:1487‑M1/1‑62 (MQ=255)                                                                                                                                                                                                cgttcccc‑gccagTGCACCTGATTTCACCAGCGCCAGCGCTTCTGTCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctt                                                   <  1:405843‑M1/245‑196 (MQ=255)                                                                                                                                                                                                  ggttcc‑gccaGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactttaactgttcagttttgatgcctttat                                                <  1:157553‑M1/247‑199 (MQ=255)                                                                                                                                                                                                   gcccc‑gcCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatc                                               <  2:407689‑M1/249‑200 (MQ=255)                                                                                                                                                                                                               gCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattct                                     <  2:155302‑M1/250‑210 (MQ=255)                                                                                                                                                                                                                  ccGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattcttta                                  >  2:136395‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                     gATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgccttagcggtgtcccccattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattctttaatcc                              >  2:260567‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                         tCACCAGCGCCAGCGCTTCTGGCAGCGCATCcaaccaccggaatgctggacttgttgtgtgctatcagcgcttcaacagcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcggcccagacagccaggcgtacatcttatctttcgccttcgcggtgtaccacattgcggtatctaactgaaactgttcagttttgatgccttattcgttcaattctttaatcacgt                           >  2:229509‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                           cccagcgccagcgcTGCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccagcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattctttaatcacgtaag                        <  2:162346‑M1/250‑223 (MQ=255)                                                                                                                                                                                                                                 gccagcgcTTCTGGCAGCGCATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattctttaatcacgtaagtggtac                  <  1:198440‑M1/251‑229 (MQ=255)                                                                                                                                                                                                                                 gccagcgcTTCTGGCAGCGCATCcaaacaccggaatgatggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattctttaatcacgtaagtggtac                  >  1:109426‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                                 gcgcATCcaaacaccggaatgctggacttgttgtgtgctttcagcgcttcaaccgcgcccattgccatcgcatcgttgttggcgataaccacttcgattttgttggcgttcgggccagacagccaggcgtccatcttatctttcgcctgagcggtgtcccacattgcggtatctaactgtaactgttcagttttgatgcctttatcgttcaattctttaatcacgtaagtggtacgtgcttctgcatccgg  <  1:82868‑M1/251‑245 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     AGACTGCGCCCAATCAGTCTTATTTCTTGCTGAATTCAGCCAGGTTGTCTTTATCTACGCCAACATAAGGTACGCGGACCACTTTGTTGTCGATTTTCCAGTTGGTGCCATCA‑GCCGCACCTTTACCATCGGCCAGGTTTTTC‑GCCAGATCAAAGGTCGCTTTCGCCTGGTTGTTAGCATCGTTCAGTACGGTGCCC‑GCCAGTGCACCGGATTTCACCAGCGCCAGCGCTTCTGGCAGCGCATCGACGCCAAACACCGGAATGCTGGACTT‑GTTGTGTGCTTTCAGCGCTTCAACCGCGCCCATTGCCATCGCATCGTTGTTGGCGATAACCACTTCGATTTTGTTGGCGTTCGGGCCAGACAGCCAGGCGTCCATCTTATCTTTCGCCTGAGCGGTGTCCCACATTGCGGTATCTAACTGTAACTGTTCAGTTTTGATGCCTTTATCGTTCAATTCTTTAATCACGTAAGTGGTACGTGCTTCTGC  >  NC_000913/2239331‑2239817

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