Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,737,137	IS2 (+) +5 bp	coding (37‑41/825 nt)	bax ←	putative glycoside hydrolase Bax

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3187426	NA (NA)	77 (1.220)	76/482	0.1	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	3737137 =	0 (0.000)					coding (41/825 nt)	bax	putative glycoside hydrolase Bax
*	?	NC_000913	= 3737141	2 (0.030)	61 (0.960)	54/484	0.3	96.8%	coding (37/825 nt)	bax	putative glycoside hydrolase Bax
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

AACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAT  >  NC_000913/4497933‑4498428                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        aaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGC‑TTTGATGGTCCCTTG‑t                                                                                                                                                                                                                                                         >  1:535081/1‑249 (MQ=255)                             ggAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                                                                                                                                                                           <  2:506062/251‑1 (MQ=255)                                     tCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCATGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCACGACAACCAATATCACt                                                                                                                                                                                                                   >  2:428149/1‑251 (MQ=255)                                        aGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTa                                                                                                                                                                                                                 <  1:521333/250‑1 (MQ=255)                                           attattTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATa                                                                                                                                                                                                             <  1:530517/251‑1 (MQ=255)                                                               aGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                                                                                                                                                                           >  1:3023/1‑249 (MQ=255)                                                                  cAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                                                                                                                                                                                      >  1:174568/1‑251 (MQ=255)                                                                                 gATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAt                                                                                                                                                                         <  1:30007/249‑1 (MQ=255)                                                                                  atgaatttaatatttacCGCTACCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAACAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGAGGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCACCTGAATCTCCAGACAACCAATATCACTTAAATACGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                                                                                                                      >  2:25698/5‑251 (MQ=38)                                                                                  aTGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAGGGGGGGC‑GGAGGGAATTTGA‑GGGCCCTTGGAAGACGGGCCCCTGGAATCTC‑AGGAAAACAAAATCACTTAAAAAAGGGAAAGTCTTAATAATAGGTTTTAggataagtaatggg                                                                                                                                                                     >  1:148214/1‑238 (MQ=255)                                                                                         tactatttacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACa                                                                                                                                                                <  2:161178/250‑1 (MQ=255)                                                                                                tacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattg                                                                                                                                                        >  2:82158/1‑251 (MQ=255)                                                                                                tacCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt                                                                                                                                                         <  2:95911/250‑1 (MQ=255)                                                                                                   cGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAGATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATg                                                                                                                                                     >  2:339164/1‑251 (MQ=255)                                                                                                             aGGGGGGATTATCAGAATCAGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGGAAGGAGTGCTGCAGGGA‑GTTGATGGTGCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTCAGACTAGTCATTGTAGAACAGATGATTGATGTCTTAgggac                                                                                                                                           >  2:427960/1‑249 (MQ=37)                                                                                                                 gggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGgag                                                                                                                                       <  1:558102/251‑1 (MQ=255)                                                                                                                  ggATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCgg                                                                                                                                         >  1:330707/1‑248 (MQ=255)                                                                                                                                cgcgGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGt                                                                                                                        >  1:543674/1‑251 (MQ=255)                                                                                                                                         tCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcacac                                                                                                                 >  1:86703/1‑249 (MQ=255)                                                                                                                                            ttatttatATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGa                                                                                                             >  1:554020/1‑250 (MQ=255)                                                                                                                                                gtatataCAATCAATAAAAAAGTAACGTCATAAAAAATCCTGAACTCACAGGGTGATGCGGGATAACTGACGAAAGGAGTGCTGCAAGGA‑TTGGATGGGTGCTTG‑GAGACTGGCACCATGAATCTAAAGATAACCAATATCACTTACCTAAGTGATAGTGTTAATACTAGTTTTTAGAGTAGTCATTAGAGAGCAGAAGATTAGTGTGTTAGGGCCGGATAAACGCAGACGGCGTACCACACAGGAAAAgc                                                                                                        >  2:186630/2‑250 (MQ=16)                                                                                                                                                 tatataCATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCa                                                                                                                                >  1:129882/1‑226 (MQ=255)                                                                                                                                                 tatataCATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCa                                                                                                                                <  2:129882/226‑1 (MQ=255)                                                                                                                                                       cATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAAt                                                                                                 >  1:461382/1‑251 (MQ=255)                                                                                                                                                           cAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATACCTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGtt                                                                                             >  1:66327/1‑251 (MQ=255)                                                                                                                                                                    agaaaCCTCATAAAAACTCCTGAACTCAACGGGCGATGCGGGATAACTGACATAAGGAGTGCTGCATGGA‑TTTGATGTGCCCTTG‑TAGACTGACCGCCTGAATCTTCAGACACACAATATCACTTAAATAAGTGATAGTCTTAATACTAGATTTTAGACTCGTAATTTGAGAACAGATGATTGATGTCTTATGTCTGGAGAAACGCAGACGGAGGATCAAACATGACATAATCGCAATTGTTcagaaaagg                                                                                    >  2:54387/4‑245 (MQ=12)                                                                                                                                                                     gTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCt                                                                                   >  2:37815/1‑251 (MQ=255)                                                                                                                                                                       aaCCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCttt                                                                                 >  1:32959/1‑251 (MQ=255)                                                                                                                                                                           tCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                                                                                                                                 <  1:83308/183‑1 (MQ=255)                                                                                                                                                                           tCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                                                                                                                                 >  2:83308/1‑183 (MQ=255)                                                                                                                                                                                       cTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACgg                                                                 >  1:195146/1‑251 (MQ=255)                                                                                                                                                                                           aCTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCCGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGtctc                                                             >  1:178060/1‑251 (MQ=255)                                                                                                                                                                                                  gggggATGCGGGATAACTGACGTAAGGAGTACTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAAACAACCAATATCACTTAAGTAAGGGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGGTGATTGATGTCt                                                                                                                                                  >  2:101502/5‑159 (MQ=18)                                                                                                                                                                                                  gggTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCt                                                                                                                                                  <  1:101502/159‑1 (MQ=255)                                                                                                                                                                                                                       gTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACTCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCCCGTTGCCCGGCAACATGGTGTagca                                 <  1:69594/251‑1 (MQ=255)                                                                                                                                                                                                                                gctgcAGGGA‑TTTGAAGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAtccagcca                         <  2:359502/250‑1 (MQ=255)                                                                                                                                                                                                                                        gA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCACCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTtct                <  2:3023/251‑1 (MQ=37)                                                                                                                                                                                                                                             tGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTgg            >  1:159321/1‑251 (MQ=35)                                                                                                                                                                                                                                             tGATGGTCCCTTG‑TAGACTGGCACCATGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGACGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCACCATGGTGTAGCAGCCAGCCGGATATTTCTCTgg            >  2:428015/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCACTATCACTTAAATAAGTGCTAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGCCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGt         >  2:574324/1‑250 (MQ=18)                                                                                                                                                                                                                                                 ggTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGt         >  1:486947/1‑250 (MQ=33)                                                                                                                                                                                                                                                  gTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTa        >  1:220644/1‑250 (MQ=33)                                                                                                                                                                                                                                                        tG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAt  >  1:569959/1‑250 (MQ=11)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        AACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGA‑TTTGATGGTCCCTTG‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAAT  >  NC_000913/4497933‑4498428

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