Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,651,851	IS2 (–) +5 bp	intergenic (+1866/+873)	ydfD → / ← ynfP	lysis protein/protein YnfP

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1651548	NA (NA)	23 (0.530)	22/486	0.9	100%	noncoding (706/706 nt)	IS2	repeat region
	?	NC_000913	= 1651855	0 (0.000)					intergenic (+1870/+873)	ydfD/ynfP	lysis protein/protein YnfP
*	?	NC_000913	1651851 =	0 (0.000)	41 (0.950)	37/486	0.2	100%	intergenic (+1866/+877)	ydfD/ynfP	lysis protein/protein YnfP
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATAC  >  NC_000913/2070027‑2070507                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             gCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGTGGGCCAGTCTAGTGGCg                                                                                                                                                                                                                                        <  1:366597/251‑1 (MQ=14)                    atagctggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGGACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAACACTAGTATTAAGACTATCACTTATTTATGTGATATTGGTTGTCTGGAGATTCAGGGGGCAAGTCTAGGGGCGAAGCATCCTCCCGTGtttt                                                                                                                                                                                                                     >  2:159538/4‑251 (MQ=255)                         tggctgCTACACCATGTTGCCGGGAAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtc                                                                                                                                                                                                                >  1:178925/1‑251 (MQ=255)                                tACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCTTTgt                                                                                                                                                                                                         >  2:346616/1‑250 (MQ=255)                                 acacCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTg                                                                                                                                                                                                          <  1:75734/249‑1 (MQ=255)                                         ttGCCGGGCAACGAGGTAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATATTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAAGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTTATTCTCCTTGTTGGTCaggg                                                                                                                                                                                                >  2:438192/1‑250 (MQ=255)                                           gtcGGGCAACGAGGGTGACCGTCATTCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTTTGGTACGCCGTCTTCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCGCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCaggagg                                                                                                                                                                                              <  1:63836/249‑1 (MQ=255)                                                         gagaCCGGCAGCCCCGGTTCAAAGCTCTGCTGAACAATGGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATTACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATAGTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATGCACATTGATGGTCAGGAGTTAACTGATCAgcgc                                                                                                                                                                                <  2:39708/251‑1 (MQ=17)                                                               gTCATCCCCGGTTCAAATCTCTGCTGAACAATTGCTATCTGTTCCTGTGTGGTACGCCGTCGGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGATATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGga                                                                                                                                                                          <  2:445503/251‑1 (MQ=255)                                                                      cctgTTCAAATCCCTGCTGAACAATTTCGATTTTTACCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACGAGTCTAAATACTAGTATTAAGACTATCACTTCTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCTTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAAtt                                                                                                                                                                   <  1:281693/248‑1 (MQ=18)                                                                                gcgctgctgAACAATTGCGATCTTTTCCTTCGTGGTACGCCCGCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCCCTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGgtgttgt                                                                                                                                                         <  2:264861/248‑4 (MQ=21)                                                                                                          ccTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGGGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGGGTCTGCGCTgg                                                                                                                               >  2:297729/1‑251 (MQ=255)                                                                                                                tgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGAATACTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGGAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGCGTCTGAGCTGGAGCg                                                                                                                           >  1:85733/1‑249 (MQ=255)                                                                                                                  gTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGtt                                                                                                                       <  1:297729/251‑1 (MQ=255)                                                                                                                         gTCTGCGTTTCGCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGTTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTTGGACTCGTGATAGATGAGTCTGCGCTGGAGCGTGTTAAGGTAg                                                                                                                <  1:153606/251‑1 (MQ=255)                                                                                                                                 ttCTCCGGCCCTAAGACATCAATCAGCTGTTCTCCAATGACTAGTCTAAAAACTTGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAgg                                                                                                         <  2:134782/250‑1 (MQ=255)                                                                                                                                     ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTc                                                                                                       <  1:435219/248‑1 (MQ=255)                                                                                                                                     ccGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCat                                                                                                    <  1:266930/251‑1 (MQ=255)                                                                                                                                          ccTAAGACATCAATCATCTGGTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTcc                                                                                               <  2:165878/251‑1 (MQ=255)                                                                                                                                            tAAGACATCAATCATCTATTAACCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTAATATTGGATATCTGGAGATTAAGGGTGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGTGCAGGAGGTAACTTATCAGAGCGCAGGAAATAATTATGGTGTGGGACTGTTATTAGCTGAGTCTGCGATGGAGCGTTTTGGGGTAGTGAAAGGTCCATATTACGt                                                                                             >  2:102490/1‑251 (MQ=21)                                                                                                                                                      atcatcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCAATTATTTAAGTGATATTGGTTGTCTGGAGATTCAAGCGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGCTGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGGTAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGGGAAAGGTCCATATTCCGTACTGTACGGt                                                                                   >  1:373873/1‑251 (MQ=255)                                                                                                                                                             gTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAg                                                                             <  1:448268/250‑1 (MQ=255)                                                                                                                                                                      tGACTAGTCTAAAAAATAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCg                                                                   >  1:104260/1‑251 (MQ=255)                                                                                                                                                                       gACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCgg                                                                  <  1:346616/251‑1 (MQ=255)                                                                                                                                                                           aggcgaaaaaCTAGTATTAATACTATCACTGATTTAAGTGTTATTGGTTGTCGGGATATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCTCGTGTTTTAATTCTCATTGATGTTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATATATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGTCAATTTGCGGTAtt                                                              <  2:447136/246‑1 (MQ=255)                                                                                                                                                                              atAAAAACTAGTATTAATACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGGGGTAGTGAAAGGTCAATGTTCCGTACTGTACGGTTCACAGGAAATTGGCgggcttgtt                                                           >  2:29745/2‑244 (MQ=38)                                                                                                                                                                                    aCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTACGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTaaattt                                                     >  2:352129/1‑247 (MQ=255)                                                                                                                                                                                              aGACTATCACTTATTTAAGTGATATTGGTTGTCTGTAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCaaa                                            >  1:52366/1‑250 (MQ=255)                                                                                                                                                                                                                    tatTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATACCCAAAAAAGGGGAATCTCCGGACTCCCt                     >  2:347213/1‑251 (MQ=255)                                                                                                                                                                                                                                ggAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTAGCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATGCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTaaa         <  1:438192/251‑1 (MQ=255)                                                                                                                                                                                                                                        aGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGGATTGTTAACTTCATCACCCAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTtataa  >  2:290670/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             GCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATAC  >  NC_000913/2070027‑2070507

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