Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	4,130,635	Δ13 bp	27.0%	coding (801‑813/2433 nt)	metL →	Bifunctional aspartokinase/homoserine dehydrogenase 2

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4130634	27 (0.620)	10 (0.240)	10/246	NT	27.0%	coding (800/2433 nt)	metL	Bifunctional aspartokinase/homoserine dehydrogenase 2
	?	NC_000913	4130648 =	28 (0.660)					coding (814/2433 nt)	metL	Bifunctional aspartokinase/homoserine dehydrogenase 2

TTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGA  >  NC_000913/4130503‑4130792                                                                                                                                                  |                                                                                                                                                 ttCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACtt                                                                                                                                               <  1:233982/149‑1 (MQ=255)      gcgTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAg                                                                                                                                          >  1:148542/1‑149 (MQ=255)        gTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGcc                                                                                                                                        >  2:160677/1‑149 (MQ=255)              aTCTTGAGCGACGTCGCCGGGGTATACAGTACCGCCACGTGCAAAGTGAAAGATGCCTGCCTGCTGCCGTTTCTGCGTCTCGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGttt                                                                                                                                   >  4:208930/1‑148 (MQ=255)              aTCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTc                                                                                                                                  >  1:68608/1‑149 (MQ=255)                cTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTg                                                                                                                                <  1:374344/149‑1 (MQ=255)                              ccGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAAGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGAc                                                                                                                  <  1:160677/149‑1 (MQ=255)                                   gTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCa                                                                                                             >  2:35450/1‑149 (MQ=255)                                      taCAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                                                             >  1:147616/1‑114 (MQ=255)                                      taCAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                                                             <  2:147616/114‑1 (MQ=255)                                         cgccggggtatacagtgccgacccgcgtaaagtgaaagatgcctgcctgctgccgttgctgcgtctggatgaggccagcgaactggcgcgcctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTgc                                                                                                        <  4:61015‑M2/43‑1 (MQ=255)                                              cGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                                                 <  1:302775/118‑1 (MQ=255)                                              cGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                                                 >  2:302775/1‑118 (MQ=255)                                                    acagtgccgacccgcgtaaagtgaaagatgcctgcctgctgccgttgctgcgtctggatgaggccagcgaactggcgcgcctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTac                                                                                             >  4:336089‑M2/95‑148 (MQ=255)                                                     cagtgccgacccgcgtaaagtgaaagatgcctgcctgctgccgttgctgcgtctggatgaggccagcgaactggcgcgcctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTaca                                                                                            >  4:294708‑M2/94‑148 (MQ=255)                                                             gAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGAt                                                                                    >  4:82473/1‑148 (MQ=255)                                                                aGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                                                 <  4:380586/148‑1 (MQ=255)                                                                       tgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAc                                                                         <  1:35450/149‑1 (MQ=255)                                                                        gtgaaagatgcctgcctgctgccgttgctgcgtctggatgaggccagcgaactggcgcgcctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAcg                                                                        >  2:33669‑M2/75‑149 (MQ=255)                                                                        gcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAc                                                                         >  4:302079/1‑148 (MQ=255)                                                                            gctgcCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCa                                                                    >  1:365885/1‑149 (MQ=255)                                                                                  gTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAc                                                              <  2:365885/149‑1 (MQ=255)                                                                                   ttGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                                                   <  4:82010/111‑1 (MQ=255)                                                                                   ttGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                                                   >  3:82010/1‑111 (MQ=255)                                                                                    tgctgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAcgc                                                            <  2:209727/149‑1 (MQ=255)                                                                                           tCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTg                                                      >  3:47961/1‑148 (MQ=255)                                                                                              ggATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGcc                                                   <  4:351685/148‑1 (MQ=255)                                                                                                        tctggatgaggccagcgaactggcgcgcctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTgg                                        >  1:129143‑M2/43‑149 (MQ=255)                                                                                                        aGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTg                                         >  3:357276/1‑148 (MQ=255)                                                                                                             atgaggccagcgaactggcgcgcctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTgcgc                                   <  1:362413‑M2/112‑1 (MQ=255)                                                                                                                   gcgcCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTATTGGTGCGCGTAtt                              >  3:130772/1‑148 (MQ=255)                                                                                                                               tCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTAtt                              >  1:133170/1‑136 (MQ=255)                                                                                                                               tCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTAtt                              <  2:133170/136‑1 (MQ=255)                                                                                                                                cccGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACga                <  1:254720/149‑1 (MQ=255)                                                                                                                                 ccGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACga                <  4:47961/148‑1 (MQ=255)                                                                                                                                    ctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACgatgat            >  1:193554‑M2/15‑149 (MQ=255)                                                                                                                                    ctggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACgatgat            <  2:193554‑M2/135‑1 (MQ=255)                                                                                                                                      ggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                                              >  4:164094‑M2/13‑65 (MQ=255)                                                                                                                                      ggcggctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                                              <  3:164094‑M2/53‑1 (MQ=255)                                                                                                                                           ctcccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGt      <  4:249778‑M2/141‑1 (MQ=255)                                                                                                                                           cGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGtt     >  1:344780/1‑149 (MQ=255)                                                                                                                                           cGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGtt     <  2:148542/149‑1 (MQ=255)                                                                                                                                           cGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGt      <  3:302079/148‑1 (MQ=255)                                                                                                                                             cccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGttt    >  3:70859‑M2/6‑148 (MQ=255)                                                                                                                                              ccgtACTTTACAGCCGGTTTCTGGCAGCGAATTCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTttga  >  1:176329‑M2/5‑149 (MQ=255)                                                                                                                                              ccgtACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTttg   <  4:130772‑M2/144‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGA  >  NC_000913/4130503‑4130792

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