Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,089,842	IS1 (+) +9 bp	intergenic (‑129/‑146)	yefM ← / → hisL	YefM antitoxin of the YoeB‑YefM toxin‑antitoxin pair and DNA binding transcriptional repressor/his operon leader peptide

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1978503 =	NA (NA)	15 (0.590)	15/484	0.6	93.5%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	2089842 =	1 (0.040)					intergenic (‑129/‑154)	yefM/hisL	YefM antitoxin of the YoeB‑YefM toxin‑antitoxin pair and DNA binding transcriptional repressor/his operon leader peptide
*	?	NC_000913	= 2089850	0 (0.000)	23 (0.900)	22/484	0.2	100%	intergenic (‑137/‑146)	yefM/hisL	YefM antitoxin of the YoeB‑YefM toxin‑antitoxin pair and DNA binding transcriptional repressor/his operon leader peptide
	?	NC_000913	3583428 =	NA (NA)					noncoding (1/768 nt)	IS1	repeat region

AATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATG  >  NC_000913/3583212‑3583669                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  aaTCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATACGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCATGCTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGAGATGCTGCCAACTTACTGATTTAGTGTatgat                                                                                                                                                                                                                 >  1:46463/1‑251 (MQ=255)         tGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGTTGCTGCCAACt                                                                                                                                                                                                                                    >  2:156767/1‑224 (MQ=255)         tGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACt                                                                                                                                                                                                                                    <  1:156767/224‑1 (MQ=255)                                cAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGt                                                                                                                                                                                  >  1:196422/1‑251 (MQ=255)                                    ggatggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCt                                                                                                                                                                              >  1:285303/1‑251 (MQ=255)                                       tggatTCGAAGTTCATGCATTACGATTCTTTAGCCTCTTGATTCCAGACCCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATc                                                                                                                                                                           >  1:231198/1‑251 (MQ=255)                                                     aTGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                                                                                                                                                             >  1:18525/1‑251 (MQ=255)                                                                  attTTGCCTCTCGATTCCATAACCTTCAGTTAAAAAACGTTTTGCCGTAAATGAATATTATAGAAATAATGATGATTTCAGAAACCCTGATCTAGAAGAACGGTTAGTGATCGTGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTAGGATGGTGTTTTTGAGGTGCTCCTGTGGCTTTTGTTTCTATCAGCTTTCCAACCTCTTCAGCTGCTGAc                                                                                                                                                <  2:90615/251‑1 (MQ=37)                                                                                        ccTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                                                                                                                                                                        >  1:286721/1‑141 (MQ=255)                                                                                        ccTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                                                                                                                                                                        <  2:286721/141‑1 (MQ=255)                                                                                                              acGTAAAAGAATTTTATAGAAATAATGATGACTACATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGAATTTAAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGCCGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCtc                                                                                                    <  2:97729/250‑2 (MQ=21)                                                                                                                              taGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGc                                                                                    >  1:89272/1‑251 (MQ=255)                                                                                                                                       tgatgaTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTCTCTCTGCTCTCACTGCCGTAAAAc                                                                            >  2:156054/1‑250 (MQ=255)                                                                                                                                                         tGATCTAAAAGCACGGGTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAACGGGGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTAACTCCTGTTCAACTACTGCCGGGGAGGTGTGTAACGGCAAAAGCACCGACGGACGACAGCGCTATCTCTGATCTCACTCCCGCACAACACGGCAAATGccattcaat                                                         >  2:274775/1‑243 (MQ=17)                                                                                                                                                          gaccgtcAAGAACGGTTAGTGAACGGGGATTAATGCCTTTCAATATTCACTTATGTCATTAAAGGTGATGCTACCAGCTTACTGATTTAGTTTATGATGGTGTTTTTGAGGTGCTCCAGTGGTCTCTGTTTTTATCAGATGTCCCTCCTGTTCAGCTATCGAGTGTGTGGTGCGCAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTTCTCTCACGGCCGTAAAACATGGCAACTGCAGTTCTCtt                                                        <  2:31742/245‑1 (MQ=14)                                                                                                                                                                     aCGGTTAGTGATAGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCAAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTtct                                             >  2:66745/1‑251 (MQ=255)                                                                                                                                                                                         ataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGa                          >  1:117569/1‑250 (MQ=255)                                                                                                                                                                                                                gACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATg  <  1:156054/251‑1 (MQ=33)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  AATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATG  >  NC_000913/3583212‑3583669

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