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)	43 (1.510)	40/482	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)	33 (1.160)	31/480	0.1	97.0%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	3429988 =	1 (0.040)					coding (36/258 nt)	yrdB	DUF1488 domain‑containing protein YrdB

CAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACG  >  NC_000913/2169941‑2170409                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               cAAACATTAAATCGGCATTAATATCCTGTAATGAACGATCGGCAATGGAACCATGCATCGAGAGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCCTTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCCATAGTg                                                                                                                                                                                                                            >  2:184178/1‑251 (MQ=255)                                                gAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATg                                                                                                                                                                             <  1:216518/251‑1 (MQ=255)                                                                              aCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGCAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGACAGAGGTGACTCACATGACAAAAACAGTATCCACCAGTAAAAAAccc                                                                                                                                               >  1:23745/1‑251 (MQ=255)                                                                                        cacaGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGc                                                                                                                                     >  2:77366/1‑251 (MQ=255)                                                                                               gacCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTTAAAAACCCCGTAAACAGCATTCGcc                                                                                                                              <  1:77366/251‑1 (MQ=255)                                                                                                                       ggAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTAGAGAATACCGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGa                                                                                                                                                                                     >  1:35286/1‑172 (MQ=255)                                                                                                                       ggAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATGAATCCTTCGGCGATAAGCATGTTCGTACGTCCGCTGGCGAGAAGAACGGTCTGAGCTTACCCAGCAAGAGTGGACACGCGGCTAAGTGAGTAAACTCTCATTCAGAGGTGa                                                                                                                                                                                     <  2:35286/172‑1 (MQ=37)                                                                                                                                                     ggTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGc                                                                        >  2:18888/1‑251 (MQ=255)                                                                                                                                                     ggTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGACAGAGGTGACTCACATGACAAAAACAGGATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGGGAAGCCCTGAAGCTTGCTGATCGCATCGGTTGTACTGCCg                                                                         <  2:23745/250‑1 (MQ=255)                                                                                                                                                        aaTCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGc                                                                     <  1:18888/251‑1 (MQ=255)                                                                                                                                                              cGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGAGAGAAAACTCTCAGTCAGAGGAGACTCACATGAAAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCACCTGAATTTCGCAGTGAAGCCATGAAGCTTGCTGAACGCATCGGTGTTACGGCCGCAGCACGTg                                                                >  2:255575/1‑250 (MQ=18)                                                                                                                                                                    attgttGGCGGTCATTAATCCTGCGGCGATAAGCATTGTCGGACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACGCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCACCCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTTAAGCCCTGAAGCTTGCTTAAGGCATCGGTGTTATTGCAGCAGCGCGTGGACTCAg                                                         <  2:207633/247‑1 (MQ=17)                                                                                                                                                                     tttttGGCGGTCATTAATCCATTGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTAAGTCAGAGGTGACTCAAATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTAGCAGGAATTTCGCCGTGAAGCCCTGACGCTTGATGAACACATCGGTGTTACTGCAGCAGCCCGTGAACTCaac                                                        >  2:138495/1‑249 (MQ=17)                                                                                                                                                                                   tAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATcac                                           <  1:184178/250‑1 (MQ=255)                                                                                                                                                                                              cGATAAGCATTGTCGGACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATGTCGCATTGAAGCCCTGAAGCGGGCTGAACGCATCGGTGTTAGTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACt                               <  1:23955/251‑1 (MQ=255)                                                                                                                                                                                                aTAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGAGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCATGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTgg                             >  2:298997/1‑251 (MQ=255)                                                                                                                                                                                                                           agaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCACCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACACCTGGCGCAGTAAACAGCAAAATCAGCAGACg  >  2:278557/1‑251 (MQ=18)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACG  >  NC_000913/2169941‑2170409

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