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

TTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGAT  >  NC_000913/4130503‑4130779                                                                                                                                    |                                                                                                                                                  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)                            cGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgc                                                                                                        <  4:61015‑M1/148‑44 (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)                                       aCAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctac                                                                                             >  4:336089‑M1/1‑94 (MQ=255)                                        cAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctaca                                                                                            >  4:294708‑M1/1‑93 (MQ=255)                                              cGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                                    <  1:302775/118‑1 (MQ=255)                                              cGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                                    >  2:302775/1‑118 (MQ=255)                                                           gTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacg                                                                        >  2:33669‑M1/1‑74 (MQ=255)                                                             gAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGAt                                                                       >  4:82473/1‑148 (MQ=255)                                                                aGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                                    <  4:380586/148‑1 (MQ=255)                                                                       tgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAc                                                            <  1:35450/149‑1 (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)                                                                                           tCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactgg                                        >  1:129143‑M1/1‑42 (MQ=255)                                                                                           tCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTg                                         >  3:47961/1‑148 (MQ=255)                                                                                              ggATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGcc                                      <  4:351685/148‑1 (MQ=255)                                                                                                aTGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgc                                   <  1:362413‑M1/149‑113 (MQ=255)                                                                                                        aGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTg                            >  3:357276/1‑148 (MQ=255)                                                                                                                   gcgcCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTATTGGTGCGCGTAtt                 >  3:130772/1‑148 (MQ=255)                                                                                                                       cTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgcgtattgtcaccagccacgatgat            <  2:193554‑M1/149‑136 (MQ=255)                                                                                                                       cTGGCGGCTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgcgtattgtcaccagccacgatgat            >  1:193554‑M1/1‑14 (MQ=255)                                                                                                                         ggcggcTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagcta                                                                                              >  4:164094‑M1/1‑12 (MQ=255)                                                                                                                         ggcggcTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagcta                                                                                              <  3:164094‑M1/65‑54 (MQ=255)                                                                                                                              cTCCCGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgcgtattgtcaccagccacgatgatgtctgt      <  4:249778‑M1/148‑142 (MQ=255)                                                                                                                               tCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTAtt                 >  1:133170/1‑136 (MQ=255)                                                                                                                               tCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTAtt                 <  2:133170/136‑1 (MQ=255)                                                                                                                                cccGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgcgtattgtcaccagccacgatgatgtctgttt    >  3:70859‑M1/1‑5 (MQ=255)                                                                                                                                cccGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACga   <  1:254720/149‑1 (MQ=255)                                                                                                                                 ccGTactttacagccggtttctggcagcgaattcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgcgtattgtcaccagccacgatgatgtctgtttga  >  1:176329‑M1/1‑4 (MQ=255)                                                                                                                                 ccGTactttacagccggtttctggcagcgaaatcgacctgcaactgcgctgtagctacacgccggatcaaggttccacgcgcattgaacgcgtgctggcctccggtactggtgcgcgtattgtcaccagccacgatgatgtctgtttg   <  4:130772‑M1/148‑145 (MQ=255)                                                                                                                                 ccGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACga   <  4:47961/148‑1 (MQ=255)                                                                                                                                    |                                                                                                                                                  TTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGAT  >  NC_000913/4130503‑4130779

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