Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,429,988	IS3 (–) +3 bp	coding (34‑36/258 nt)	yrdB ←	DUF1488 domain‑containing protein YrdB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1094245 =	NA (NA)	30 (1.080)	29/486	0.1	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3429990	0 (0.000)					coding (34/258 nt)	yrdB	DUF1488 domain‑containing protein YrdB
*	?	NC_000913	2170171 =	NA (NA)	21 (0.760)	20/484	0.4	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	3429988 =	0 (0.000)					coding (36/258 nt)	yrdB	DUF1488 domain‑containing protein YrdB

ATACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGC  >  NC_000913/2169925‑2170404                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          aTACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATc                                                                                                                                                                                                                                       >  1:110731/1‑251 (MQ=255)    ccATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCtt                                                                                                                                                                                                                                     >  1:5241/1‑250 (MQ=255)        cAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAccc                                                                                                                                                                                                                                 >  2:167863/1‑250 (MQ=255)         aGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGATCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACCGACGACCAGCGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGTAGACTTTTGGTCATCACCGTGATGTTTTTGTCGGTCATTAATCCTTCGGCGATACTCCTTGTCGTACTTCCTCTGTCGAGAATAACGGGCTGATCTTAccc                                                                                                                                                                                                                                 >  1:133850/1‑249 (MQ=255)           cACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTTCTTCCGCTGTCGAGAGTAACGGTCTGAGCGTACCCAGCa                                                                                                                                                                                                                             <  2:38951/251‑1 (MQ=255)             ccGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATGGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAAt                                                                                                                                                                                                                           <  1:117843/251‑1 (MQ=255)                        tAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACAcg                                                                                                                                                                                                                 <  2:90694/250‑1 (MQ=255)                                 aTTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGGGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCCTATCTTACCCAGCAATAGTGGACCCGCGGCTAAGTg                                                                                                                                                                                                       >  2:272407/1‑251 (MQ=255)                                            gTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTTCCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACtctc                                                                                                                                                                                            >  2:10370/1‑251 (MQ=255)                                                  aaCGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGAGAGTAAACTCtcagtc                                                                                                                                                                                        >  1:155330/1‑249 (MQ=255)                                                                gAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATg                                                                                                                                                                        >  1:123797/1‑251 (MQ=255)                                                                  aCCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTTGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                                                                                                                      <  1:166937/251‑1 (MQ=255)                                                                           cGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGt                                                                                                                                                             <  1:167863/251‑1 (MQ=255)                                                                                       aTGGCGGAACGTGCACCCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTa                                                                                                                                                  >  1:70583/1‑250 (MQ=255)                                                                                              aCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAccc                                                                                                                                          >  1:146437/1‑251 (MQ=255)                                                                                                    ccaccaCAGACGACCAGAGTAATGTCTTTATTGTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa                                                                                                                                    <  2:110731/251‑1 (MQ=255)                                                                                                             acgacCAGAGTAATGGCGTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCg                                                                                                                           <  2:225099/251‑1 (MQ=255)                                                                                                                                         aaaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGc                                                                                               >  2:129011/1‑251 (MQ=255)                                                                                                                                                             agaCTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGtt                                                                           >  1:257788/1‑251 (MQ=255)                                                                                                                                                                              cGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGa                                                          <  2:146437/251‑1 (MQ=255)                                                                                                                                                                               gTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGa                                                          >  2:46092/1‑250 (MQ=255)                                                                                                                                                                                        ttGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTATACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACTCATTGTTGTTACTGCCGCAGCCCGTGAACTCAGCCt                                                 <  1:172698/250‑1 (MQ=255)                                                                                                                                                                                                  ttaatgcttaggcgaaaacaactGTCGCACTTCCGCTCTCGAGAATAACCGTCTGGTCTTACACAGCAACAGTGGACACGCGGCTAAGCGAGTAAACTCTCAGTCAGAGGTGAATCACATGACAAAAACAGTATCAACCAGTCCACACCCCCGTAAACAGCATACGCCTGAATTTCGCAGTGAAGCCATGCAGCTTGCTGAACGCATCGGTGTTACTGCCGCCGCCCGTACACTCAGCCTGCATGAATcac                                      >  2:253774/23‑251 (MQ=2)                                                                                                                                                                                                                             tACTTCCGCTGTCGAGAATAACGGTCTGATCTTCCCCAGCAATAGTGGACACGCGGCAAAGTGAGTAAACTCTAAGCCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTACACAGc           >  1:18895/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTCGAGAATAACGGTCTGATCTTACCCAGAAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTAAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTGGCTGAACGAATCGGTGTTACTGCCGCAGACCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTACACAGCAAAATcagc  >  1:187300/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          ATACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGC  >  NC_000913/2169925‑2170404

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