Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,649,932	IS2 (+) +5 bp	coding (139‑143/192 nt)	ydfD →	lysis protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1467910 =	NA (NA)	18 (0.590)	18/490	0.6	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	1649932 =	0 (0.000)					coding (139/192 nt)	ydfD	lysis protein
*	?	NC_000913	= 1649936	0 (0.000)	16 (0.530)	14/490	0.8	100%	coding (143/192 nt)	ydfD	lysis protein
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

AGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTTAG  >  NC_000913/2070038‑2070517                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       agaAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtcctcc                                                                                                                                                                                                                                       >  1:83188/1‑251 (MQ=35)           ggctggctGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTa                                                                                                                                                                                                                              <  1:71562/250‑1 (MQ=40)                          cATGTTGCCGGGCAACGAGGGAGACCGTCCTCCCCGGTGCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGATATTCAGGGGGCCAGTCTAGTGGCGAAACATGCTCCCGGGTTTTGAGTCTCATTGATGGTc                                                                                                                                                                                                              <  2:141494/251‑1 (MQ=14)                                       aaCGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACtt                                                                                                                                                                                                  >  2:265696/1‑250 (MQ=255)                                                                              aaCAATTGCGATCTTTTCTTGTGTGGTACGCCGTCTGCGTTTCTCCGGTCCTAAGACATCAATCAGCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCTTCCTCCCGTGTTTTAATTCTCATTGATGGTCAGTAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGAt                                                                                                                                                          <  1:176568/251‑1 (MQ=255)                                                                                                      ggTACGCCGTCTTCGTTTCTCCGGCCCTAAGACATCAATCCTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCTCTTAGTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGATTCTGCGCTGGAGCgtgt                                                                                                                                  <  1:258832/251‑1 (MQ=255)                                                                                                                            ggCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAACCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGAAGTCTGGAGATTAAGGGGGCCAGTCTAGTGGCGCAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGCTAGATGAGACTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCatat                                                                                                            >  1:154428/1‑251 (MQ=37)                                                                                                                                cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGTTATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTcc                                                                                                         >  2:205889/1‑250 (MQ=255)                                                                                                                                  aaGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGTTATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTcc                                                                                                         <  1:205889/248‑1 (MQ=255)                                                                                                                                       atcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTa                                                                                                 <  1:262749/251‑1 (MQ=255)                                                                                                                                         caatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTAc                                                                                                <  1:256848/250‑1 (MQ=255)                                                                                                                                                            gacgaGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATGATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACGGTACGGTTCAAAGGCAATTGGCga                                                                            <  2:209772/247‑2 (MQ=255)                                                                                                                                                                aGTCTAAAAGTTAGTATTAAGAGTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGGGGTGAAGCATCCTCCCGTGTTTTAATTCTCATTGGTGGTAAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTAt                                                                         <  2:106530/250‑1 (MQ=39)                                                                                                                                                                        aaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCAGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCAGTGTTTTAATTCTCATTGATGGTAAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAGCtt                                                                >  1:192960/1‑251 (MQ=255)                                                                                                                                                                                 tAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGTGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCTGTATTGTTAACTTCATCACCa                                                        <  1:55490/250‑1 (MQ=255)                                                                                                                                                                                     aCTATCCCTTATTTAAGTGATATTGGTTGGCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAg                                                   <  2:236218/251‑1 (MQ=255)                                                                                                                                                                                     aCTATCACTTATTTAAGTGATATTGGTTGTCTGTAGATTCAGGGGGCCAGTCTAGTGGCGATGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTAGGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAg                                                   <  2:222260/251‑1 (MQ=255)                                                                                                                                                                                                   aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTTATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCgg                                      <  2:83188/250‑1 (MQ=255)                                                                                                                                                                                                     gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTACCTTATCAGCGCGcc                                                                                                                                                                                       >  2:55656/1‑103 (MQ=255)                                                                                                                                                                                                     gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGcc                                                                                                                                                                                       <  1:55656/103‑1 (MQ=255)                                                                                                                                                                                                     gTGATAGTGGTTGGCTGGAGATGCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGGACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACt                                   <  1:156696/251‑1 (MQ=255)                                                                                                                                                                                                                                        cTAGTGGCGAAGCTTCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATCATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCCCCGGACTCCCTATATCACTTAAATTGATACACCttttgag  >  2:92168/1‑246 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       AGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTTAG  >  NC_000913/2070038‑2070517

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