Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	789,812	IS2 (–) +5 bp	coding (164‑168/1149 nt)	galK ←	galactokinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	72 (0.890)	60/490	0.3	97.2%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	789812 =	2 (0.030)					coding (168/1149 nt)	galK	galactokinase
*	?	NC_000913	= 382590	NA (NA)	102 (1.270)	89/486	0.0	98.0%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	= 789816	2 (0.030)					coding (164/1149 nt)	galK	galactokinase

GTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTT  >  NC_000913/382343‑382828                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               gTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAAtt                                                                                                                                                                                                                                             <  1:59082/251‑1 (MQ=11)   gCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTa                                                                                                                                                                                                                                            <  1:321612/250‑1 (MQ=11)          ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAg                                                                                                                                                                                                                                    <  2:665953/251‑1 (MQ=33)                 gTAACGGAATAGCAGAGAGCTTCGTGAACACGATAAATCGTGACTACATCAGTATCATTCCCAAACCAGACGGGTTAACGGCAGCAAAGATCCTTGCAGAGGCGTTCGAGCAGTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGttt                                                                                                                                                                                                                             <  1:289618/251‑1 (MQ=255)                         aTAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                                                                                                                                     <  2:463648/251‑1 (MQ=40)                           aGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCa                                                                                                                                                                                                                   <  1:665254/251‑1 (MQ=255)                              agagagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATAAGGCAGTTTGCGCCAATCACt                                                                                                                                                                                                                 >  1:509559/1‑250 (MQ=255)                                     ttCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                                                                                                                                                         >  1:471542/1‑251 (MQ=255)                                                     aGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGCATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATACCGGCAcc                                                                                                                                                                                          >  1:413701/1‑250 (MQ=255)                                                       cGTGACTACATCAGTATCCTGCCCAAACCAGACAGGTTAACGGCAGCACAGAACCTTACCGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTATGTGATACCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGACAccctg                                                                                                                                                                                       >  1:280339/1‑248 (MQ=17)                                                          ggttacatcaggaccaggtccAAACCAGACGGGTTAACGGCAGCAAAGACACTTGCAGAGGCGTTAGAGCAGTATAACGAATGGCATCCGCCGAGTGCGCTGGGGTATCGTTAGCCAAGGGAATATTTGCGGCATCGGGCTTGTAATGGGTTAAGTGATAACAGATGCCTGGAAATACAGGGGCAAATCCAATAACTCATCAGGAAGTTGGAGCAACTAACTCCCGGACCCTGAATAACGGAACCATGACa                                                                                                                                                                                    <  2:755008/232‑1 (MQ=2)                                                          gACTACATCAGTATCATGCACAAACCAGACAGGTTCACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGACATCCGCATAGTGCGCTGGGTTATAGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAGGAGATAACAGATGACTGGAAATATAGGGGCAAATCCAATTAACTATCAGGCAGTTTGCGCAAATCACTACCGGACCCTCAATAACGGCACCATGAc                                                                                                                                                                                     >  1:750860/1‑250 (MQ=17)                                                                aTCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACg                                                                                                                                                                              <  2:284364/251‑1 (MQ=255)                                                                   aGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGt                                                                                                                                                                           <  1:528734/251‑1 (MQ=255)                                                                       tCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGTCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTc                                                                                                                                                                       >  1:127420/1‑251 (MQ=255)                                                                         aTGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCCTCCGCATAGTGCGCTGGGTTATAGCTCGCCACGGGAATATCTGCGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCCAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCCTGCCAAATACGGGTATTcgc                                                                                                                                                                     >  2:371886/1‑251 (MQ=255)                                                                            cccAAACCAGACGGGTTAACTGCAGCAAAGATCCTTGCAGAGGCGTTCGAGCATGATAACGAATGGCATCCGCATAGTGCGCTGGTTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGcc                                                                                                                                                                  <  1:270955/251‑1 (MQ=255)                                                                                 aCCAGACGTGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa                                                                                                                                                             <  1:314971/251‑1 (MQ=255)                                                                                   cAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaaca                                                                                                                                                           <  2:817094/251‑1 (MQ=255)                                                                                    aGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCg                                                                                                                                                                 >  1:415937/1‑244 (MQ=255)                                                                                    aGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCTCGCCg                                                                                                                                                                 <  2:415937/244‑1 (MQ=255)                                                                                     gacggtgtAACGTCAGCAAATAACCTTGCAGAGGCGTGCGAGTATTATAACGAAGTGCATCCGCATAGTGCGCGGGGTTAACGCTCGCCACTGGAATATCTTCGTCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATTTCTGGAAATATAGGGGCAAATCCAATTTCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCTCCATGACAAATACGGGTATTCGCGCCGATaacaac                                                                                                                                                         <  2:513389/244‑1 (MQ=18)                                                                                       cGGGTTAACGGCCGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAAAGAATGGCATCAGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGAGGTAACAGATGTCTGGAAATCTAGGGGCACATCCCATTACCTATCAGGCAGTTTGAGCCAATCACTACCGGACCCTGAATAACGGCCCCCTGACAAATACGGGTATTCGCGCCGATAACAACtt                                                                                                                                                       >  2:189524/1‑251 (MQ=18)                                                                                                 gcagcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAAGATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACACCTTCACCTGCAGt                                                                                                                                             <  2:507378/251‑1 (MQ=255)                                                                                                   agcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGttt                                                                                                                                           <  1:468712/251‑1 (MQ=255)                                                                                                      aaagaaccattatgagtcgtgcGAGCATTAGAACGAATGGCAGCCGTATGGTGCGCTGGGTTATCGCGCGCCGCGGGAATATCTGCGGCAGCGGGCTTGTAAGGGGTTTTGTGAGAACAGATGTCTGGAAATATAGGGGCAAAGCCAATTTCCTATCAGGCAGTTTGCGCCAATCACTACCGGGCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                                                                                                                                        <  2:280339/230‑1 (MQ=255)                                                                                                       aaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTTCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcatc                                                                                                                                       <  1:356675/251‑1 (MQ=255)                                                                                                         gAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcatca                                                                                                                                      >  1:797357/1‑250 (MQ=255)                                                                                                          aaCCTTGCAGAGGCGTTCGATCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaat                                                                                                                                    <  1:298143/251‑1 (MQ=255)                                                                                                                 caggggtgttctggCATTATAACGCATGGCATCCGCATAGTGCGCTGGTTTATCGCGCGCCTCTGGAATATCTGCGGCAGCGGGCTGGTAATGGGTTAAGTGATTACAGATGTTTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCTCTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCGCCTGTAGTTTCATCAATAATTAcc                                                                                                                             <  2:555957/238‑1 (MQ=255)                                                                                                                 cAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAcc                                                                                                                             <  1:751776/251‑1 (MQ=255)                                                                                                                 cAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAcc                                                                                                                             <  1:174480/251‑1 (MQ=255)                                                                                                                   gagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTg                                                                                                                           <  1:631769/251‑1 (MQ=255)                                                                                                                      gCGTTCGAGCATTATAACGAATGTCATCCGCATAGTGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTGTGCGCCAAACACTACCGAACCCTGTATAAc                                                                                                                                                                                                <  2:582090/179‑1 (MQ=17)                                                                                                                      gCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGCATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGAGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                                                                                                                                                >  1:582090/1‑179 (MQ=255)                                                                                                                              gCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCt                                                                                                                <  1:325218/251‑1 (MQ=255)                                                                                                                                   ctaACGAATAGCATCCGCATAGTGCGCTGGGTTCTCGCTCGCCACGGGAATCTCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGtt                                                                                                                                                                                                                              >  2:39226/2‑136 (MQ=255)                                                                                                                                   ataACGAATGGCATCCGCATAGTGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGtt                                                                                                                                                                                                                              <  1:39226/136‑1 (MQ=255)                                                                                                                                            ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTctct                                                                                                  >  1:766491/1‑251 (MQ=255)                                                                                                                                               aTCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCttt                                                                                                >  1:82666/1‑250 (MQ=255)                                                                                                                                                   gCATAGTGCGCTGGGTTATCGCTCGCCAGGGGAATATCTGCGGCAGCGGTCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACAt                                                                                           <  1:14253/251‑1 (MQ=255)                                                                                                                                                           cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                                   <  1:138589/251‑1 (MQ=255)                                                                                                                                                           cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                                   <  1:531406/251‑1 (MQ=255)                                                                                                                                                                       gttcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGTTTAAGGGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTAGCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCTCCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATGAAGCTTCAGATAGTTTTTc                                                                       <  2:425205/249‑1 (MQ=255)                                                                                                                                                                          tgtcACGGGAATATCGGCGGCATCTGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATATGGGCAAATCCAATTACCTATCAGGCAGTTTGCTCCAATCGCTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGAGAACAACTACACCTGCAGTTTTATCAATAATTTCCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTCAGCTGCAGATAGTTTTTCAg                                                                     <  2:129668/247‑1 (MQ=38)                                                                                                                                                                           gcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGACAAATTCAATTAACTATCAGGCAGTTTGAGCCAATCACTACCGGCCCCTGAATAACGGAAGCATGACAAATACGGGCATTCGCGCCGATAACAACTACACCTGCAGTTTCATCAATAATTACCTGTACTGCACTCTGAATATTTTTCTCTTTCACATTAAGCTTCAGATAGTTTTTCAGtt                                                                   >  1:834050/1‑251 (MQ=37)                                                                                                                                                                             cACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTcc                                                                 >  1:652148/1‑251 (MQ=255)                                                                                                                                                                              aCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCCGGCAGTTTGCGCCAATCACTACCGGAACCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAAATTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATCTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCt                                                                >  2:822178/1‑251 (MQ=255)                                                                                                                                                                                    atatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACa                                                          <  1:170500/251‑1 (MQ=255)                                                                                                                                                                                        cTGCGGCAGCGTGCTTGCAATGGGGTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCGCTTTAACATTAAGCTTCAGATAGTTTTTCAGt                                                                    <  2:127344/237‑1 (MQ=255)                                                                                                                                                                                        cTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGt                                                                    >  1:127344/1‑237 (MQ=255)                                                                                                                                                                                        cTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGa                                                        >  1:625325/1‑249 (MQ=255)                                                                                                                                                                                          gCGGCAGCGGGGTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAgg                                                     <  2:797357/250‑1 (MQ=255)                                                                                                                                                                                           gggCAGCGGGCTTGTAATGGGTTAAGTGTTAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGtt                                                   <  2:471542/250‑1 (MQ=255)                                                                                                                                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTg                                                  >  1:772139/1‑251 (MQ=255)                                                                                                                                                                                                        gtactGGGGTAAGAGATACCAGATGTCTGGCAATCTAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACAGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAAAAACTTCACCTGCAGTTTCATCAATAATTCCCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCTGATGGTTTTTCCGTTCCTGTGCCATATGGTTGGTAATAAGCGGt                                      >  2:16445/5‑251 (MQ=255)                                                                                                                                                                                                        gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGt                                      <  1:815831/251‑1 (MQ=255)                                                                                                                                                                                                         aaaagggTTAAGTGATAACAGATGTCTGGAAATAGAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGtt                                     <  2:558492/247‑1 (MQ=255)                                                                                                                                                                                                          aaTGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTg                                    >  1:69285/1‑251 (MQ=255)                                                                                                                                                                                                                      gATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccacca                        <  1:706131/251‑1 (MQ=255)                                                                                                                                                                                                                          aCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGc                    >  1:357233/1‑251 (MQ=255)                                                                                                                                                                                                                            aGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCgg                  >  1:325599/1‑251 (MQ=255)                                                                                                                                                                                                                                   tGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACTAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCg           >  2:144796/1‑251 (MQ=255)                                                                                                                                                                                                                                          ataGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCCGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgggt    >  1:147297/1‑251 (MQ=255)                                                                                                                                                                                                                                           gaGGGTCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGTGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGGCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGtt   <  1:804467/250‑1 (MQ=255)                                                                                                                                                                                                                                            aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGttt  <  1:189524/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTT  >  NC_000913/382343‑382828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.