Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,315,674	IS2 (+) +5 bp	coding (2187‑2191/2673 nt)	rcsD →	RcsD phosphotransferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2068941 =	NA (NA)	11 (0.940)	11/304	0.1	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	2315674 =	0 (0.000)					coding (2187/2673 nt)	rcsD	RcsD phosphotransferase
*	?	NC_000913	= 2070271	NA (NA)	16 (1.370)	14/304	0.0	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 2315678	0 (0.000)					coding (2191/2673 nt)	rcsD	RcsD phosphotransferase

CGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAG  >  NC_000913/2070078‑2070432                                                                                                                                                                                                  |                                                                                                                                                                  cGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGTCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATAGTGGTTGTCTGGAGTTTCAGGGGGCCTGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATc                                                                                                          <  1:135757/251‑1 (MQ=255)   aGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGGTTGAATTCTCATTGATGGTCAGGAGGTAACTTATCAg                                                                                                        <  1:246098/251‑1 (MQ=255)    gngAGACCGTCATCAACGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAgc                                                                                                       >  1:216414/1‑251 (MQ=255)                           aaGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGGGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTAGCAGCGCGCCGGAGATAATTATGGtgtg                                                                                <  1:212337/251‑1 (MQ=255)                               tctGCTGAACAATGGCGATCTTTTCCTGTGTGGAACGCCATCTGCGTTTCTCCGGCACTCAGACCTCAATCATCTCTTCTACAATGAATAGTCTACAAAATAGAATTAAGACTATCCATTATTTAAGTGATAATGCATGTCTGGAGATTCAGGGGGCACGTCTAGTGGCGAAGCATACTCCCGTGTTTTAATTATCATTGATGGTCAGGAGGTAACTTATCAGCGCGCATGCGATAATTATGGTGTggg                                                                              >  1:181726/1‑249 (MQ=255)                                       aCAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCg                                                                                                                                                             <  1:117034/162‑1 (MQ=32)                                          aTTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCATTCATCTGTTCTCCAATGAATAGTCTAAAAACTACTATTAAGACTATCACTTATTTAAAAGAAATTGTTTGTCTGGCGATTCAGGGGGCCAGTCTAGTGGAGAAGCATCCTCCCGTGTTTTATTTCTCATTGATGGTCAGGAGGCAACTTCTCAg                                                                                                        >  1:153841/1‑212 (MQ=17)                                             gCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCTGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAg                                                              <  1:20040/251‑1 (MQ=255)                                                                 aCGCCGTCTTCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCagg                                                                                                                      <  1:150071/175‑1 (MQ=255)                                                                  cGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAg                                         <  1:89196/251‑1 (MQ=255)                                                                                ctcCGGCACTAAGACATCAATCATCTGTTCTCCAATGACTTGTATAAAAACTAGTATTAAGCCTATCACTTATTTAAGTGATATCGGTTGTCAGGAGATTCAGGGGGCCAGTCTAATGGCGAAGCATCATCCAGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATAAGCGCGCCGGAGATACTTATGGTGTgg                                                                               >  1:179818/1‑199 (MQ=21)                                                                                      nccTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATACTCCCGTGTTTTACTTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGgtg                                                                                  >  1:51911/2‑190 (MQ=255)                                                                                      cccTAAGTCATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCTCTTATTTAAGTGATAGTGGTTGTCTTGATATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCTTTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTAGGGTGTGGGGCTGTTGATAGATGAGGCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTc                     <  1:278786/251‑1 (MQ=37)                                                                                                          nntctcCANNGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAg                                                                                                                                                          <  2:4688/96‑1 (MQ=255)                                                                                                                                  cTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGACAGTCTAGTGGCGACGCATCCTCCCGTGTTTTAATTCTCATTGCTGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTACATATTCCGTACTGTACGGTTCACAg  >  1:180882/1‑226 (MQ=255)                                                                                                                                                                      gTTGTCTGGAGATTCAGGGGGCCAGTCTAGTNGCGNAGCANCCNCCCGTGTTTTAATTCNCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTNNGNNNTNGGACTGtt                                                                        <  2:187450/120‑1 (MQ=255)                                                                                                                                                                      gTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAATGATGGTCAGGAGGTAACTTATCAACGCGCCGGGGATAATTCTGGTGTGGGGCTGtt                                                                        >  1:187450/1‑120 (MQ=38)                                                                                                                                                                                           ccAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTAc                >  1:225674/1‑155 (MQ=255)                                                                                                                                                                                                  |                                                                                                                                                                  CGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAG  >  NC_000913/2070078‑2070432

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