Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	568,449	IS2 (–) +5 bp	coding (135‑139/333 nt)	emrE →	multidrug/betaine/choline efflux transporter EmrE

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	568449 =	0 (0.000)	26 (0.670)	26/474	0.6	100%	coding (135/333 nt)	emrE	multidrug/betaine/choline efflux transporter EmrE
	?	NC_000913	= 1469240	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 568453	0 (0.000)	28 (0.730)	27/472	0.5	100%	coding (139/333 nt)	emrE	multidrug/betaine/choline efflux transporter EmrE
	?	NC_000913	1467910 =	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCAC  >  NC_000913/1469002‑1469473                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          gCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAact                                                                                                                                                                                                                               >  2:116004/1‑251 (MQ=33)  ccAGAGAAATAACTGGCTTGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTACTCTAATAACTAGGATTAAGACTATCACTTATTTAAGTGATATTGGTTGCCTGGAGATTCAGGGGGCCAGTCGAGACCCCAGAactt                                                                                                                                                                                                                              <  2:183851/251‑1 (MQ=11)                 gctggctgCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTcc                                                                                                                                                                                                               <  1:341576/251‑1 (MQ=40)                           acacCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGaaa                                                                                                                                                                                                      >  2:292342/1‑250 (MQ=255)                             acCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAAtt                                                                                                                                                                                                   >  1:252113/1‑251 (MQ=255)                                           cAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTAc                                                                                                                                                                                     >  2:96238/1‑251 (MQ=255)                                                  ggAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGa                                                                                                                                                                              <  1:293270/251‑1 (MQ=255)                                                                                                                         gTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTGGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACAttt                                                                                                       <  1:97904/251‑1 (MQ=255)                                                                                                                                   cccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAgg                                                                                             >  1:316687/1‑251 (MQ=255)                                                                                                                                       aaGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGaaa                                                                                          <  2:252113/250‑1 (MQ=255)                                                                                                                                                        ttCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAAc                                                                        >  1:314774/1‑251 (MQ=255)                                                                                                                                                           tcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTATGTTTTTGAACTAt                                                                     >  1:359789/1‑251 (MQ=255)                                                                                                                                                                  aCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGgaaga                                                               >  1:45523/1‑250 (MQ=255)                                                                                                                                                                     aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAAc                                                             <  2:413096/249‑1 (MQ=255)                                                                                                                                                                     aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAAc                                                             <  2:413088/249‑1 (MQ=255)                                                                                                                                                                                     ttAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGa                                            >  2:141386/1‑250 (MQ=255)                                                                                                                                                                                              tCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAat                                  >  2:62881/1‑251 (MQ=255)                                                                                                                                                                                                        aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGtt                        <  1:75323/251‑1 (MQ=255)                                                                                                                                                                                                        aaGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCATCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGt                         <  2:256976/250‑1 (MQ=255)                                                                                                                                                                                                            gATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCaa                     <  1:141386/250‑1 (MQ=255)                                                                                                                                                                                                                 tGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACATGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAAccc                <  1:41833/250‑1 (MQ=255)                                                                                                                                                                                                                      gTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCTGGAAc           <  2:323714/250‑1 (MQ=255)                                                                                                                                                                                                                              aTTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCAc  >  1:9926/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          GCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCAC  >  NC_000913/1469002‑1469473

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