Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	577,688	IS3 (+) +3 bp	coding (76‑78/498 nt)	rrrD →	lysozyme

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 316483	NA (NA)	46 (0.890)	44/488	0.2	100%	noncoding (1255/1255 nt)	IS3	repeat region
	?	NC_000913	577688 =	0 (0.000)					coding (76/498 nt)	rrrD	lysozyme
*	?	NC_000913	566777 =	NA (NA)	48 (0.930)	47/488	0.2	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 577690	0 (0.000)					coding (78/498 nt)	rrrD	lysozyme

GTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCT  >  NC_000913/316236‑316733                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           gTATGAGCGAACAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAAAAAGAAACTCGCTTAGGCCTGTGTCCCTATTACGTGGGTAGGATCATa                                                                                                                                                                                                                                                          >  2:268230/1‑250 (MQ=2)           aaaaGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGACTGTGTCCATATTACGTGGGTAGGATCATAAAGGCCTCggg                                                                                                                                                                                                                                               >  2:161751/1‑251 (MQ=18)                  tgcggcgaCGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGTAACAATTGGAAAACAAGAACCTCGCGTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTg                                                                                                                                                                                                                                         <  2:250260/243‑1 (MQ=255)                    ctgctACGATACTGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGGAAAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                                                                                                                                                                                                                                      >  2:374365/1‑251 (MQ=255)                      gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                                                                                                                                                                     <  1:142330/250‑1 (MQ=40)                                tGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACCAGAAACTCGCCTAGGCCTGTGTCCATATTACGTGGGTATGAACATCAAGACCTAGGGTCCTTGGCCACCCTTTTGGcc                                                                                                                                                                                                                          >  2:445466/1‑250 (MQ=11)                                 gCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCAt                                                                                                                                                                                                                         >  1:104766/1‑251 (MQ=255)                                   cTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGGAAACAAGAACCTCGCCTTGGCCTGTGTCCATATTAAGTGGGTAGGATCATCACGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGt                                                                                                                                                                                                                       >  2:26436/1‑251 (MQ=255)                                     gCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGt                                                                                                                                                                                                                       <  1:268230/249‑1 (MQ=255)                                         ggAAAGCTTCTTTCATTCGATGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGACGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAAGCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGGGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat                                                                                                                                                                                                                 <  2:71085/251‑1 (MQ=255)                                          gAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATACATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGACTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGAACATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCat                                                                                                                                                                                                                 >  2:528445/1‑250 (MQ=255)                                                    ttCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATAc                                                                                                                                                                                                       <  1:412209/250‑1 (MQ=255)                                                       aTTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTTTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTa                                                                                                                                                                                                   <  2:110654/251‑1 (MQ=255)                                                       aTTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGt                                                                                                                                                                                                    >  1:183862/1‑250 (MQ=255)                                                           gCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAg                                                                                                                                                                                               <  1:116697/251‑1 (MQ=255)                                                                    ggAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGttt                                                                                                                                                                                      >  1:146995/1‑251 (MQ=255)                                                                         gTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGTCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGTGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAAc                                                                                                                                                                                  <  1:310555/250‑1 (MQ=255)                                                                           aTCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGGCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCgg                                                                                                                                                                               >  1:382052/1‑250 (MQ=255)                                                                               aTGGAGACCACTTTATCAGCCGGGAAATAATGCGGGCCACGGCGTTCAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACCATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCCTATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCCAACGTGAGc                                                                                                                                                                           >  1:376326/1‑251 (MQ=18)                                                                                            tATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGAt                                                                                                                                                               <  1:161751/250‑1 (MQ=255)                                                                                             aTCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAATTCCTCGGGTCGTTGACCACTCTTTTGTCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTCCAGCCGTGAGCAAATCGACCGAtt                                                                                                                                                              <  2:512952/250‑1 (MQ=255)                                                                                                    gggAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGGTCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATAcc                                                                                                                                                                                                      >  1:288095/1‑203 (MQ=255)                                                                                                    gggAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGGCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACGCTTTTGGCATGTCAGCATATTCACATAcc                                                                                                                                                                                                      <  2:288095/203‑1 (MQ=255)                                                                                                           aatGCGGGCAACGGTGTTTAATTATATCGAATGTGAATACCAACCCAGCCAGCGGCACAGTTGGACTGCCGGCACCCGACCGGAACCATTAGCCCACACGACCCTCGCTTAGGCCTGTGTCCCTAATAAGTGGGTAGGATCCTACAGTCCTCGACTCGATGTCCACTCTTTTGGCATGACCGCACATTCACATACCGTACCAGCACAAGTTTCAACCCTGAGCAAATCGACCACTTTTGTTGAAGGTTCa                                                                                                                                                >  1:510446/1‑250 (MQ=255)                                                                                                                       ggTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGAACAGTTGGTGTGGCGGACTCAGTCCGGACCAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATAATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTACAACGTGGAAACg                                                                                                                                    >  2:533274/1‑250 (MQ=255)                                                                                                                                    atCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCg                                                                                                                       <  1:374365/250‑1 (MQ=255)                                                                                                                                      cGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCt                                                                                                                    <  1:132848/251‑1 (MQ=255)                                                                                                                                        aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCtt                                                                                                                   >  1:169525/1‑250 (MQ=255)                                                                                                                                                tACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAAttct                                                                                                           <  2:296643/250‑1 (MQ=255)                                                                                                                                                     tCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGCACAATTTGAAAACAAGAACCTCGCCTAGGCCTGTGTCCCTATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGGATGTCAGCATATTCACCTACCGTTCCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTATTGACGGTTCAAACGTGGAAAAGAAATAGATATCCGTCTTTAAAttctttaatc                                                                                                     >  2:201027/1‑246 (MQ=37)                                                                                                                                                           gcggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATTCCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGCAACGAAATAGAAATCCGTCTTAAAATTCTTCATAcctatgg                                                                                               >  2:538694/1‑245 (MQ=255)                                                                                                                                                                ggcACAGTTGGTGTGGCGGCCTCCGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTaa                                                                                           >  2:165462/1‑250 (MQ=255)                                                                                                                                                                   acaGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAACACACGAAACTCGCTTAGGCCTGTGTCCATATTACGCGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCt                                                                                             >  1:479909/1‑245 (MQ=255)                                                                                                                                                                   acaGTTGGTGGGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATAGTCACATACCGTACCAGCAAAAGTTTCAACCGGGAGCAAATCGACCGATTTGGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAGATTCTTCATACATCTGGCt                                                                                             <  2:479909/245‑1 (MQ=255)                                                                                                                                                                            gtgGCGGCCTCAGTCCGGAACAATGTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATTAcgtcgt                                                                              <  2:56504/251‑1 (MQ=255)                                                                                                                                                                                 ggcCTCAGTCAGGAACAATTTGAAAACAAGAACCTCGCTCAGCCCTGTGTCCATATTACGAGCGTAGGATCATAAACTCCTCGGGTCGTTGTCCCCTCTTTTGGCATGTCAGCATATTCACATCCCGTACCAGCAAAAGTTTCAACCGTGCGCAAATCGACCGTATTTGTTGACGGTTCAAACGAGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGttt                                                                         >  1:337683/1‑251 (MQ=255)                                                                                                                                                                                      cAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCGGGCTAAATGATGACGTCGTAGTTTGAct                                                                     <  1:184787/250‑1 (MQ=255)                                                                                                                                                                                      cAGTCCGGAACAATTTGAAAACAAGAACCGCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGGACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAGAATGCTTCATACATCTGGCTAAATGATGACGTCTTAGTTTGAct                                                                     <  2:67067/250‑1 (MQ=255)                                                                                                                                                                                             gAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATc                                                             <  1:172404/251‑1 (MQ=255)                                                                                                                                                                                                     tGAAAACAAGAACCTCGCTTAGGCCTGTGGTCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGGCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGTTAATCGGCCGATTGTGTTGACGGTTCAAACGTGGAAAGGAAATAGGTATCCGTCTTAAAATTCTGCATACATTGGGGTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGAt                                                      <  1:323656/250‑1 (MQ=255)                                                                                                                                                                                                               aaCCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGacac                                           <  1:7696/251‑1 (MQ=255)                                                                                                                                                                                                                    cGCTTAGGCCTGTGTCCATATTACGTTGGTAGGATCATAAAGTCCTCGGGTCGTTCTCCACTCTTTTGGCATGTCAGCATAGTCACATACCGGACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGCTGATTTACTCGACACTTCaa                                      <  2:241642/251‑1 (MQ=255)                                                                                                                                                                                                                    cGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTc                                        <  2:146995/249‑1 (MQ=255)                                                                                                                                                                                                                        tAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCTTCGGGTAGTTGTACACTCTTTTGGCATGTCAGCATATTAACATACCGTACCAGCAAAAGTTTCCACCGTGAGCAACTCGACCGATTTTGTTGACGGTTCCAACGTGGAAACGAAGTAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATgac                                  >  1:539368/1‑249 (MQ=255)                                                                                                                                                                                                                          ggCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCat                                >  1:24308/1‑251 (MQ=255)                                                                                                                                                                                                                                 gtCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTaa                          >  1:313176/1‑250 (MQ=255)                                                                                                                                                                                                                                     atatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAg                     <  1:165462/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATa         <  2:313176/250‑1 (MQ=255)                                                                                                                                                                                                                                                       cATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGtc   <  1:528445/251‑1 (MQ=255)                                                                                                                                                                                                                                                       cATAAAGTCCTCGGGTCGGTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGgtgt    <  1:445466/250‑1 (MQ=255)                                                                                                                                                                                                                                                        aTAAAGCCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGtct  >  1:416112/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           GTATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTGTCT  >  NC_000913/316236‑316733

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