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)	16 (0.390)	16/486	1.0	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)	12 (0.290)	12/486	1.3	100%	coding (143/192 nt)	ydfD	lysis protein
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

TACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTT  >  NC_000913/2070030‑2070515                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     tACGCCAGAGAAATAAATGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGATTCCCCGGCCCTAAGACATCCATCATCTGTTCACCAATGACTAGTCTAAAAACTAGAATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGCGCTTCAGGGGGCCAGTCTAGTGGCGa                                                                                                                                                                                                                                               >  2:89969/1‑249 (MQ=11)      cAGAGAAATAACTTGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGGTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCCCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGTCGAAGCAtcc                                                                                                                                                                                                                                        <  2:282353/251‑1 (MQ=16)                   gggtggctgCTACACCATGTTGCCGGTCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCATGTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCAGCGGTTCTCCAACGACTAGTCTAAAAACTAGTATTGAGACTATCTCTTATTGAAGTGATATTGGTTGCCTGGGGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTa                                                                                                                                                                                                                            <  1:373306/247‑1 (MQ=14)                                acCATGGTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAGATCGCTGCTGAACAATTGCGATCTTTTCCGGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTATGTCATCAATCAGCTGTTCTCCAATGATTAGTCTAATAACGAGTATTAAGACTATCACTTTTTTAAGTGATATTGGTTGTCTGGAGATTCATGGGTCCTGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATgg                                                                                                                                                                                                              <  1:209359/251‑1 (MQ=14)                                                          ccGTCATCCCCGGTTCAAATCTCTGCTGAACTATTGCGATCGTGTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCTCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCg                                                                                                                                                                                    <  1:420650/251‑1 (MQ=255)                                                                                      aaCAATTGCGATCTTTTCCTTTGTGGTACGCCGTCGGCGTTTCTCCTGCCCTAAGACATCAATCATCTGTTCTCCAATTACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGGTGTCTGGAGATTCAGTGGGCCAGTCTAGTGGCGAAGCATTCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATAATGGTTTGTTACTGTTGAt                                                                                                                                                        <  2:272864/251‑1 (MQ=18)                                                                                                    tttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTg                                                                                                                                             <  1:41665/248‑1 (MQ=255)                                                                                                         tgtgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGCTTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTgg                                                                                                                                       >  2:215631/1‑249 (MQ=255)                                                                                                              ggTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCgtgt                                                                                                                                <  1:296902/251‑1 (MQ=255)                                                                                                                        cTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTg                                                                                                                      >  1:55725/1‑251 (MQ=255)                                                                                                                                   cGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGCGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCat                                                                                                            >  1:177338/1‑250 (MQ=255)                                                                                                                                   cGGCCCGAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGTCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTGATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATTAGTTTGCGCTGGAGCGTGTTGAGGTATTGAAAGGGCCata                                                                                                           <  2:204814/251‑1 (MQ=25)                                                                                                                                               atcaatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGt                                                                                                >  1:126108/1‑250 (MQ=255)                                                                                                                                                          gTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAgg                                                                                    >  1:71495/1‑251 (MQ=255)                                                                                                                                                                cAATGACTAGTCTAAAAACTAGTATTCCGACTATCACTTATTTAAGCGATATTGGTTGTCTGGAGCTTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTg                                                                              >  2:384648/1‑251 (MQ=255)                                                                                                                                                                   tGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCg                                                                           >  2:127570/1‑251 (MQ=255)                                                                                                                                                                             aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGtt                                                                   <  1:215631/249‑1 (MQ=255)                                                                                                                                                                                       aTTAAGACTATCACTTATTTAAGTGAGATGGGTTGTCTGGAGATTCAGGGGGCCAGTCTATTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATGGATGGTCAGGATGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGGGTTGAGTTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCGCAGGGAATTGGCGGTATTGTTAACTACATCAc                                                        <  2:362094/250‑1 (MQ=255)                                                                                                                                                                                             aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATGCCGTACTGTACGGTTCACAGGCAATTGGCTGTATTGTTAACTTCATCACCaaaaa                                                  <  2:55725/250‑1 (MQ=255)                                                                                                                                                                                                         ttAAGTGATATTGGTTGTCTGGTGATTCAGGGGGCCAGTCTAGTGGCGAAGTATCCTCCCGTGTTGTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCGGCTCTGGATCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGCTGAAtcttc                                      <  2:220628/250‑3 (MQ=255)                                                                                                                                                                                                          tAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCCGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGGGAATCTCCgg                                    >  1:91035/1‑251 (MQ=255)                                                                                                                                                                                                            aGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGGGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAAtt                                                                               >  1:16697/1‑206 (MQ=255)                                                                                                                                                                                                            aGTGAGATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTTCGGTTCACAGGCAAtt                                                                               <  2:16697/206‑1 (MQ=255)                                                                                                                                                                                                                     ggTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTatat                         <  1:130552/251‑1 (MQ=255)                                                                                                                                                                                                                            gggAGATTCAGGGGTCCAGTCTAGTGGCGAAGCAGCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACGTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGGTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTACGTATGGTACGGTTCACAGGCAAGTGGCGGTATTGTTAACTTCATCCCCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTaa                  <  2:398776/250‑1 (MQ=255)                                                                                                                                                                                                                             ggAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTTATGGTCAGGAGGTAACTTAGCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGGCTGCGCTGGAGCGTTTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTaa                  <  2:310371/250‑1 (MQ=255)                                                                                                                                                                                                                                            cAGTATAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGCGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGAGAATCTCCGGACTCCCTATATCACTTCAATTGCTACAACttttt  >  2:92870/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     TACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTT  >  NC_000913/2070030‑2070515

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