Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,353,837	IS2 (+) +5 bp	coding (877‑881/1518 nt)	lysU ←	lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3186096 =	NA (NA)	23 (0.730)	23/490	0.4	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 4353841	0 (0.000)					coding (877/1518 nt)	lysU	lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase
*	?	NC_000913	4353837 =	0 (0.000)	39 (1.230)	35/492	0.1	100%	coding (881/1518 nt)	lysU	lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

ATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGG  >  NC_000913/3185865‑3186333                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              aTCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTAGTAGACTGGCCCCCTGAAtc                                                                                                                                                                                                                             <  1:19410/250‑1 (MQ=255)            aaTTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                                                                                                                                                                 <  2:2488/251‑1 (MQ=255)               ttCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                                                                                                                                                               <  1:32773/250‑1 (MQ=255)                 cAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                                                                                                                                                                            >  2:77209/1‑251 (MQ=255)                        cGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCGAAATCTCCAGACAACCAATATCACTTaaa                                                                                                                                                                                                     <  2:22827/251‑1 (MQ=255)                             ggAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAg                                                                                                                                                                                                 <  2:12551/250‑1 (MQ=255)                                aTCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGAt                                                                                                                                                                                             >  1:225499/1‑251 (MQ=255)                                                                                               cGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGGTATTAGACTGGCCCCCGGAATCTCCAGACAACCAATATCAGTTAAATAAGTGATAGTCTTAATTCTAGGTTGTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGga                                                                                                                              <  1:162139/251‑1 (MQ=255)                                                                                                  gTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATAAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGaa                                                                                                                           <  1:169103/251‑1 (MQ=255)                                                                                                          aTTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATTCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGAc                                                                                                                   <  2:103624/251‑1 (MQ=255)                                                                                                                                         aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATTTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGt                                                                                    <  2:107204/251‑1 (MQ=255)                                                                                                                                               tAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTcagc                                                                               >  1:108276/1‑250 (MQ=255)                                                                                                                                               tAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTcagc                                                                               <  1:15028/250‑1 (MQ=255)                                                                                                                                               tAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATAATCCGTTATTAGACTGGCCCCCTGAATATCCAGACAACCAATATCACTTAAATAAGAGATAGCCTTAATACTAGTTTTAAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGCCCGGAGAAACGCAGACGGCGTACGACACAGGAACAGATCGCAATTGTTcagc                                                                               >  2:80285/1‑250 (MQ=21)                                                                                                                                                aaTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACAACAATAAGCGAATTGGTTAAAAATATTTATTGATCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTAATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTCCCACACAGGAAAAGATCGCAATTGATcagc                                                                               >  1:204989/1‑249 (MQ=18)                                                                                                                                                           gATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAAcc                                                                  >  1:213723/1‑251 (MQ=255)                                                                                                                                                                  tGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATg                                                           >  1:1217/1‑251 (MQ=255)                                                                                                                                                                                           aaTAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTGCCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAAc                                  <  2:15655/251‑1 (MQ=255)                                                                                                                                                                                                                           cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTgg  <  1:165738/251‑1 (MQ=34)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              ATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGG  >  NC_000913/3185865‑3186333

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