Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,580,649	IS5 (–) +4 bp	coding (811‑814/1395 nt)	hsdS ←	type I restriction enzyme EcoKI specificity subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	25 (0.560)	22/482	0.7	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 4580652	0 (0.000)					coding (811/1395 nt)	hsdS	type I restriction enzyme EcoKI specificity subunit
*	?	NC_000913	= 3653230	NA (NA)	35 (0.780)	33/482	0.3	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	4580649 =	0 (0.000)					coding (814/1395 nt)	hsdS	type I restriction enzyme EcoKI specificity subunit

GTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTCCAAAGACAATGTTATAAGCAAGT  >  NC_000913/3652987‑3653478                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         gTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTg                                                                                                                                                                                                                                                   <  1:425110/251‑1 (MQ=11)   ttCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAATAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATTAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGa                                                                                                                                                                                                                                                  <  1:51267/250‑1 (MQ=11)     ccAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGCCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGacac                                                                                                                                                                                                                                               >  2:184227/1‑251 (MQ=14)       aGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTa                                                                                                                                                                                                                                             <  2:147857/251‑1 (MQ=32)        gCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAt                                                                                                                                                                                                                                            <  2:457119/251‑1 (MQ=33)            aTAAGGTCGCCGGCCATTTCCAGTCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttttt                                                                                                                                                                                                                                        <  1:278852/251‑1 (MQ=21)                   cGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTTCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCGGTTCTATGGCTCCAGATGGCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGGGACACTATTTTTGTACAt                                                                                                                                                                                                                                  <  2:313429/250‑1 (MQ=255)                   cGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATa                                                                                                                                                                                                                                 <  2:78526/251‑1 (MQ=38)                    gCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATa                                                                                                                                                                                                                                 >  1:362790/1‑250 (MQ=38)                           aTTACCAGCCTTGGGGTAAAACGGCTCGATGACTTACACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGCATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCCGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCTTAGGTGACACTATTTTTTTACATAAAGAacac                                                                                                                                                                                                                         >  2:244862/1‑251 (MQ=255)                              aCCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAccc                                                                                                                                                                                                                       <  1:244862/250‑1 (MQ=255)                                  gCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccc                                                                                                                                                                                                                  <  2:178496/251‑1 (MQ=255)                                    cTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGTGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                                                                                                                                                                                                                >  1:2984/1‑251 (MQ=255)                                                 ggCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGCCTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATaacg                                                                                                                                                                                                   >  2:431444/1‑249 (MQ=255)                                                 ggCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAg                                                                                                                                                                                                   <  2:143564/251‑1 (MQ=255)                                                              ccaccaTGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCt                                                                                                                                                                                      <  1:431444/251‑1 (MQ=255)                                                               caccaTGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAGGGACAATATCtt                                                                                                                                                                                     <  2:37011/251‑1 (MQ=255)                                                                         tGCCATGGCAGAATCTGCTACATGCGGGACAAGAACATCTCTTTTCTGGTATGACGGCGCTTACTGCTGAATTCACTGCCGGCGAAGGTAAGCTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAAAATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAACGATGGACAATATCTTGacaaacaa                                                                                                                                                                            >  1:223804/1‑250 (MQ=18)                                                                                  aGAATCTGCTCCATGCGGGACAAGAAAATCTCTTGTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTTATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACTCTATTTTTTTACATAAAGAACACCCTACCACCAAATTATAGAAAGATGGACAATATCTTGACAAACAAAGCGAAGGtg                                                                                                                                                                  <  2:44329/251‑2 (MQ=255)                                                                                   gAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGAGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGttt                                                                                                                                                                 <  2:341920/251‑1 (MQ=255)                                                                                    aaTCTGCTCCATGCGGGACAAGAAAAGCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAAGCGAAGGTTTa                                                                                                                                                                <  1:352995/251‑1 (MQ=255)                                                                                        tGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAtttt                                                                                                                                                            >  1:368752/1‑251 (MQ=255)                                                                                         gCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTc                                                                                                                                                           <  1:5264/251‑1 (MQ=255)                                                                                            ccATGCGGGACAAGAAAATCTCTTTTCTTGTCTGACGGCGCTTACTGCGGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGTCACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAAc                                                                                                                                                        <  1:203188/251‑1 (MQ=255)                                                                                              aTGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGCACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACat                                                                                                                                                      >  2:96073/1‑251 (MQ=255)                                                                                                       aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGa                                                                                                                                                                                                  >  2:123687/1‑198 (MQ=255)                                                                                                       aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGa                                                                                                                                                                                                  <  1:123687/198‑1 (MQ=255)                                                                                                            aaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAAt                                                                                                                                        <  2:171903/251‑1 (MQ=255)                                                                                                                   tttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                                                                                                                                 >  1:385946/1‑251 (MQ=255)                                                                                                                         gTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGCGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTCATGACAACCTTTTAc                                                                                                                           >  2:161351/1‑251 (MQ=255)                                                                                                                             gACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTGTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATGTTCAACATTATTAACCATTAATGACAGCCTTTTACGAg                                                                                                                        <  2:43695/250‑1 (MQ=255)                                                                                                                                  cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCGTTTACGAGCAAAAg                                                                                                                  <  2:66272/251‑1 (MQ=255)                                                                                                                                        ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGaa                                                                                                             >  2:195218/1‑250 (MQ=255)                                                                                                                                            tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGc                                                                                                        >  2:79812/1‑251 (MQ=255)                                                                                                                                                 tCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCg                                                                                                    <  1:16354/250‑1 (MQ=255)                                                                                                                                                   ccTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACa                                                                                                 <  2:400048/250‑1 (MQ=255)                                                                                                                                                           cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCGTCCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGTCAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTaa                                                                                         <  2:160819/251‑1 (MQ=255)                                                                                                                                                           cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTaa                                                                                         <  1:10013/251‑1 (MQ=255)                                                                                                                                                                    gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAAttttt                                                                                   >  1:322454/1‑248 (MQ=255)                                                                                                                                                                              cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGAt                                                                       >  2:363837/1‑250 (MQ=255)                                                                                                                                                                              cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATGGCTTCGACAATAATTAAATTTTTAGAGGCCAGGAt                                                                       <  1:157860/250‑1 (MQ=255)                                                                                                                                                                               atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATaa                                                                     <  2:34134/251‑1 (MQ=255)                                                                                                                                                                                                aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAg                                                                                                                                                                     >  1:304268/1‑138 (MQ=255)                                                                                                                                                                                                aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAg                                                                                                                                                                     <  2:304268/138‑1 (MQ=255)                                                                                                                                                                                                  ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAAGATCTTGACAAACAAATCGAAGGTTTATTTTCAACAGTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGa                                                  <  1:56569/251‑1 (MQ=255)                                                                                                                                                                                                                               cAGGGACTTGTTCGCACCTTCCTTAGGTGACCCTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTag                                                                                 >  2:43414/1‑191 (MQ=255)                                                                                                                                                                                                                               cAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTag                                                                                 <  1:43414/191‑1 (MQ=255)                                                                                                                                                                                                                                    aCTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTCCAAAGACa                 >  2:289787/1‑250 (MQ=255)                                                                                                                                                                                                                                                  tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTCCAAAGACAATGTTATAAGCaa    <  1:5843/249‑1 (MQ=255)                                                                                                                                                                                                                                                  tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTCCAAAGACAATGTTATAAGCAAg   <  2:349797/250‑1 (MQ=255)                                                                                                                                                                                                                                                   ccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTCCAAAGACAATGTTATAAGCAAGt  <  1:425759/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         GTTTCCAGCGGATAAGGTCGCCGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCGGAAATTGCTTCGACAATAATTAAATTTTTAGAGGCCAGGATAACCAATCTTTCATACAATGACATATTAAAATATCAGCAAGAATTCCAAAGACAATGTTATAAGCAAGT  >  NC_000913/3652987‑3653478

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