Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,470,820	IS2 (+) +5 bp	intergenic (‑177/+253)	yodE ← / ← yfdK	toxic protein YodE/putative tail fiber assembly protein YfdK

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2068941 =	NA (NA)	37 (0.910)	36/496	0.2	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	2470820 =	0 (0.000)					intergenic (‑177/+257)	yodE/yfdK	toxic protein YodE/putative tail fiber assembly protein YfdK
*	?	NC_000913	= 2470824	0 (0.000)	26 (0.640)	25/496	0.5	100%	intergenic (‑181/+253)	yodE/yfdK	toxic protein YodE/putative tail fiber assembly protein YfdK
	?	NC_000913	= 2997691	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

ACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGAC  >  NC_000913/2997451‑2997931                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 acgtcAGAGAAATAACGGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTGCAAAGCTCTGCTGAACAATTGCGATCTTTTCTTGGGTTGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAAACATCTGTTCTCCAATGACTAGTCTAAAAAGTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGctc                                                                                                                                                                                                                                         <  2:223872/246‑1 (MQ=11)             tAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGgtg                                                                                                                                                                                                                            >  1:94442/1‑251 (MQ=38)                   gctggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATc                                                                                                                                                                                                                      <  1:221324/251‑1 (MQ=40)                                    gTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAATCTCTGCTGAACAATTGCGATCTTTTCCTGTGTTGTACGCCGTCTGCGTTTCTCCGGCCTTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATAGTTGTTGTCTGGAGATTCAGGGTGACAGTCTAAACTTGCTCTTTTCTTCTGGTGTTGATCTTTGTAAGATTTGTAgg                                                                                                                                                                                                     <  2:327383/251‑1 (MQ=14)                                    gTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAgg                                                                                                                                                                                                     >  1:1339/1‑251 (MQ=255)                                              aaCGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAACAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGCGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAAc                                                                                                                                                                                            >  1:151037/1‑250 (MQ=255)                                                                                    gAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATGAAGTAAAGCCGCATCTGCTTCTCTGTAATGTTTATc                                                                                                                                                     <  2:451282/251‑1 (MQ=255)                                                                                                   tttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAAcc                                                                                                                                      >  1:371550/1‑251 (MQ=255)                                                                                                     tCCTGTGTGGTTCGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCGTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCt                                                                                                                                     <  2:366510/250‑1 (MQ=255)                                                                                                       ctgagtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAAATAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGAAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTggg                                                                                                                                  >  2:116527/5‑251 (MQ=255)                                                                                                             ggTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATtct                                                                                                                            <  1:152529/251‑1 (MQ=255)                                                                                                                    cGTCTGCGCTTATCCAGCCCTAAAACATTAAACATATGATCTCCAATGACTAGTCTAAAAACTAGTATTAAAACTATCACTTCTTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGACAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGACAGATTTGTAGGATAAATAACCTTAAGTAAAGCCGCATCGGCTTCTCTGTCATGTTTATCTTTTTCCAGATAACCTGGGGCTTCTCCCCCt                                                                                                                      >  2:391493/1‑250 (MQ=25)                                                                                                                           gTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGTGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAg                                                                                                              <  2:92341/251‑1 (MQ=255)                                                                                                                                        tAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAA‑CTAGTATTAAGACCATCACTTATTTAAATAATATTGGTTTTCTGGAGCTTCAAGGAGCCCGGCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACCTTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATATTTTTCCAGCTAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTAttt                                                                                                >  1:447696/1‑250 (MQ=37)                                                                                                                                         aaGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTc                                                                                                <  2:122408/251‑1 (MQ=255)                                                                                                                                                     atcTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCtt                                                                                    >  2:320130/1‑251 (MQ=255)                                                                                                                                                        tGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTACAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTATCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTAt                                                                                 >  1:246431/1‑251 (MQ=255)                                                                                                                                                          ttCTCCAATGACTAGTCTAAAAACTAGTAATAAGACTATCACTTATTTAAGTGATATTGGTTGTTGGGAGATTCAGGGTGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTGTGAAAGATTTGTAGTGTATATAACATTAAGGAAAGCCGCATCTGCTTCTCTGTAATGTGTATCTTTTTCCAGATAACCGGGGGATTCTCCCCCTAAATCAAGCTTTTACGTTATTCTTTCATAACCTTTATc                                                                                <  2:225494/250‑1 (MQ=37)                                                                                                                                                             accAATGACTAGTCTAAAAACTAGGATTAAGACTATCACTTATTTAAGTGATATTGGTTGGCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAATCGGTAACGTTATTCTTTCATAACCTTTATCTATg                                                                            <  2:432233/250‑1 (MQ=255)                                                                                                                                                                   gactagacggacaaCTAATATTCAGACAATCCCTTATTTAAGAGATATTGACTGACTGGAGATCCAGGAAACCAGTCTAAAATAGCTCTTGTCTTATGGTGTTTATCTTTGAAAGGTTTGTCGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCTAGATAACCTGGTGATTATCACCCAAAATCAAGCTGTAACGTTATTCTTTCATTGCCTTTATCTATGAggtgg                                                                      >  2:300601/13‑251 (MQ=255)                                                                                                                                                                    aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAggtggt                                                                     >  1:351571/1‑251 (MQ=255)                                                                                                                                                                               aaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCgtt                                                          <  1:392444/251‑1 (MQ=255)                                                                                                                                                                                   aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCgttgtt                                                       <  1:8918/250‑1 (MQ=255)                                                                                                                                                                                               aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTCAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGCAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCAt                                          >  1:405945/1‑251 (MQ=255)                                                                                                                                                                                                        ttAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCAt                                 >  1:30893/1‑251 (MQ=255)                                                                                                                                                                                                          aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTg                               <  1:411780/251‑1 (MQ=255)                                                                                                                                                                                                           aGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAATCTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGGGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGa                              <  1:171370/251‑1 (MQ=255)                                                                                                                                                                                                             tGATATTGGTTGTCTGGACCGTCAGGAGGCCAGTCTCACCTTCATCTTTTCTTCTGGTGTTTATCTTTGACAGATTTGTAGGATCAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTCATGTTTATCTTTTTCAAGATAACCTGGGGATTCTCCCCCTAACTCAAACTGTAACGTTATTCTTTCATAACCTTTATCCATGAGGTGGTTAAACATCGTTGTTGAAAGAGTTTCATCCAGTTCAtagcct                            >  2:277383/1‑246 (MQ=255)                                                                                                                                                                                                                             gagaTTCAGGGGGCCATGCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACCTTAAGTAAAGCCGCATCGTCTTCGCTGTAATGTTTATCTTTTTCCAGATATCCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCgtagta            <  1:214601/251‑1 (MQ=255)                                                                                                                                                                                                                             gagaTTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCgtagt             >  1:349281/1‑250 (MQ=255)                                                                                                                                                                                                                                   cAGGGGGCCATTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCATTTCATTGACTTTAGCTTCGCGTATTAAATCa       <  2:321045/250‑1 (MQ=255)                                                                                                                                                                                                                                        ggCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGAc  >  1:340880/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 ACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGAC  >  NC_000913/2997451‑2997931

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.