Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,199,869	+GA :: IS3 (+) +4 bp	coding (935‑938/1128 nt)	intE ←	putative integrase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	315229 =	NA (NA)	32 (0.550) +TC	29/488	0.9	100%	noncoding (1/1255 nt)	IS3	repeat region
	?	NC_000913	= 1199872	0 (0.000)					coding (935/1128 nt)	intE	putative integrase
*	?	NC_000913	1094245 =	NA (NA)	53 (0.900)	51/492	0.2	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	1199869 =	0 (0.000)					coding (938/1128 nt)	intE	putative integrase

CGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCAT  >  NC_000913/1093995‑1094490                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      cGTGCGAAAGCTTTAGAGGCAAGGATTCGAGAAAATGTTGAGCGTTTAACTGGCGATAATCATGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTCTCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCAtttt                                                                                                                                                                                                                                                       >  2:343059/1‑251 (MQ=255)          gCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTAcc                                                                                                                                                                                                                                              <  1:528420/251‑1 (MQ=255)            tttAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTAccc                                                                                                                                                                                                                                             >  1:218215/1‑250 (MQ=255)              tAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCAc                                                                                                                                                                                                                                           <  1:463993/250‑1 (MQ=255)                         aTTCGGGAAACTGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGCAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGaac                                                                                                                                                                                                                               >  2:546543/1‑249 (MQ=255)                         aTTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCGCTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGGGAATGTTTGTTAGTCA‑CCATAAAGCGGCATTGGGAACCTACCCACGTAATATGGaca                                                                                                                                                                                                                               <  2:113983/251‑1 (MQ=255)                               gAAAATGTTGCGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTCTTGGAGCTGTTGTGACTCAGAAAAATACGTTAAAAACAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGGGGTTGTGAGGGATGAGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGCTACTACCCACGAAATATGGACACAGGcc                                                                                                                                                                                                                         >  2:328853/1‑251 (MQ=255)                                        gAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAg                                                                                                                                                                                                                 <  1:312512/250‑1 (MQ=255)                                                aaCTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGt                                                                                                                                                                                                         <  2:263640/250‑1 (MQ=255)                                                   tGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTc                                                                                                                                                                                                     >  1:187291/1‑251 (MQ=255)                                                             ccTGAATTCGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACGTTACGAGGCGAAAGAGATCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTcc                                                                                                                                                                                           <  2:482864/251‑1 (MQ=255)                                                               tGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAACCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCgg                                                                                                                                                                                         >  2:205634/1‑251 (MQ=255)                                                                    aCGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTCCGAGGCGAAAGAGACCGGGCGTCATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTg                                                                                                                                                                                     >  1:88943/1‑250 (MQ=255)                                                                     cGCTATACCACAAAAAGTGACGCTTAGTATTGGCGCTGTTGTGACTCAGGCAAATGCGTTAAACCCAAATGCGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGCCACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTg                                                                                                                                                                                     >  2:395138/1‑249 (MQ=255)                                                                        tataCCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                                                                                                                                >  1:629/1‑251 (MQ=255)                                                                         ataCCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                                                                                                                                <  2:298449/250‑1 (MQ=255)                                                                           aCCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgcc                                                                                                                                                                             >  2:503132/1‑251 (MQ=255)                                                                                  aaaGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCa                                                                                                                                                                       >  1:543548/1‑250 (MQ=255)                                                                                      tGACGATTAGTATGGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CAATAAAGCGGCCTTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCGGCCACACCAACtat                                                                                                                                                                  >  1:210644/1‑249 (MQ=255)                                                                                         cGATTAGTATTGGCGCGGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGGCTGGCTGATAACGCACTTTACGAGGTGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CGAGAAAGCGGTATTTTGACCCGACCCTCGTAATAGGGACACAGGCCTAAGCGAGTTTCGTGTTTTCAAAGTGTTCCGGACTGAGGCCGCCACACCAACTGTgcc                                                                                                                                                               <  1:522912/251‑1 (MQ=38)                                                                                           aTTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgc                                                                                                                                                             >  1:533738/1‑251 (MQ=255)                                                                                                 aTTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTCTCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGCGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGAACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                                                                                                                                       >  1:335610/1‑251 (MQ=255)                                                                                                             gtgACTCAGGAAAATGCGTTAACCCAAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGGGAATTTTTGTGAGTCC‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATca                                                                                                                                           <  2:195716/251‑1 (MQ=255)                                                                                                                    aGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCg                                                                                                                                     >  1:104997/1‑250 (MQ=255)                                                                                                                       aaaaTGCGTTAAACCCAAATGAGATCTATCGAATGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCAACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTATGCCGCCGCCACCGATTGTAATCACATTCGatct                                                                                                                                 >  1:340270/1‑249 (MQ=255)                                                                                                                                               tCTATCGACTGGCTGATAACGCACTTTACGAGGCGCAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGCGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATCATTAAACACCGTTGCCCGCatt                                                                                                          >  2:290717/1‑250 (MQ=255)                                                                                                                                                  aTCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTAttt                                                                                                      >  2:424160/1‑251 (MQ=255)                                                                                                                                                         ggCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGc                                                                                                <  1:15326/250‑1 (MQ=255)                                                                                                                                                                   gCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAgtgt                                                                                     >  1:485078/1‑251 (MQ=255)                                                                                                                                                                                      agaCCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTcc                                                                   >  1:586082/1‑250 (MQ=255)                                                                                                                                                                                       gaCCGGGCGTAATAATGTGGTTGTGAGGGATGTGGTTAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATTTGGACACAGGCCTAAGCGATGTTCTTGTTTGCAAATTGTTCCTTACTGAGGCTGCCACACCAACTGTTCCGCCGCCGCCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCAc                                                                 <  1:42264/251‑1 (MQ=18)                                                                                                                                                                                       gaCCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCAc                                                                 >  1:384223/1‑251 (MQ=255)                                                                                                                                                                                            ggCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCC‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCa                                                            <  2:629/251‑1 (MQ=255)                                                                                                                                                                                                      ggtggtTGTGAGGGATGTGGTGAATTTTTGTGAGTCC‑CCATAAAGCGGCATTTTGATCCTACCCACGAAATAGGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGaaa                                                  <  2:218215/251‑1 (MQ=255)                                                                                                                                                                                                          gtgttgAGTGATGTGGTGAATTTTTGTGA‑TCAACCATAAAGCGGCATTTTTATCCTCCCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGAATGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTACACTTTCAGCGAATGAaagaag                                              <  2:398815/247‑1 (MQ=14)                                                                                                                                                                                                                ggggATGTGGTGAATTTTTGTGAGGCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTcc                                        <  1:431404/250‑1 (MQ=255)                                                                                                                                                                                                                aGGGAGGTGGTGAATGTTTGTGAGTCA‑CCATAAAGCGCCATTTTGATCCTACCCACGTAATATGGACACAGGCTGAAGCGAGGTTCTTGTTTACAAATTTTTCCGGTCTGTGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGAAATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATATATTCCACTGTCAGCGAATGAAAGAAGCTTTcc                                        <  2:335610/251‑1 (MQ=14)                                                                                                                                                                                                                  ggATGTGGTGAATTTTGTTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATAGGAACACAGGCCTAAGCGAGGTTCTTTTTTTCAAATTGTTCCTGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTACATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCAc                                      <  1:70846/251‑1 (MQ=14)                                                                                                                                                                                                                           gAATTTTTGTGGGTCA‑CCAAAAAGCGTCATTTTTATCCTACCCACGTAATATGGACTCAGGCCTAAGCGAGGTTCTTTTTTTCTAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCAt                              <  1:341833/250‑1 (MQ=11)                                                                                                                                                                                                                           gAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCAtt                             >  2:226202/1‑251 (MQ=255)                                                                                                                                                                                                                                  tgtgAGTCA‑CCATAAAGCGGAATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCAGCCACCCCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCAGGCTGATGAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGCAGCTTTCCACGCAGGCATTATCGt                        >  2:172657/1‑249 (MQ=255)                                                                                                                                                                                                                                        gtccccATAAAGCGGCCTTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCCCCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAgcagca                 <  1:393944/247‑1 (MQ=21)                                                                                                                                                                                                                                                  aGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTgcg       <  1:111569/251‑1 (MQ=33)                                                                                                                                                                                                                                                      gCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAACGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCa   <  1:395138/251‑1 (MQ=14)                                                                                                                                                                                                                                                       cATTTTGATCCTACCCCCGTAATATGGACACAGGCCTAAGGGAGGTTTTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACGTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCATCAACCTTTTGCGCTCAt  <  2:557768/251‑1 (MQ=9)                                                                                                                                                                                                                                                       cATTTTGATCCTACCCACGTAATAGGGACACAGGCCTAAGCTAGGTTCTTTTTTTCAAATTTGTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCAt  <  2:437831/251‑1 (MQ=9)                                                                                                                                                                                                                                                        aTTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGATCTTGTTTTTAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCAt  >  1:240868/1‑250 (MQ=14)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      CGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAAAGTGACGATTAGTATTGGCGCTGTTGTGACTCAGGAAAATGCGTTAAACCCAAATGAGATCTATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCA‑CCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCAT  >  NC_000913/1093995‑1094490

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