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)	40 (0.610)	39/488	0.7	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)	55 (0.830)	51/488	0.3	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 577690	0 (0.000)					coding (78/498 nt)	rrrD	lysozyme

TATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAA  >  NC_000913/316237‑316723                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  tATGAGCGCAAAAGGGTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCAATGGAGAACACTTTATCAGACGGGAAATAATGCGGGCAACGGGGTTTAATTATATCGAATGTGATTACAATCGGAGGCGGCGGCACAGTTGGTGTGGCGGACTCAGTCCGGAAACATTTGAAAACAAGAAACTCGCTTAGGCCTGTGTCCATATTACGTGGGAAGGATCATaa                                                                                                                                                                                                                                                >  1:429849/1‑250 (MQ=11)     aGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTcc                                                                                                                                                                                                                                           <  1:540332/251‑1 (MQ=33)           aaaGGATGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCGGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTTGGATCATAAAGTCCTCgggg                                                                                                                                                                                                                                     >  1:222249/1‑250 (MQ=255)                  gctgctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTc                                                                                                                                                                                                                              <  1:375861/251‑1 (MQ=39)                        aCGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCActc                                                                                                                                                                                                                         <  1:657673/250‑1 (MQ=255)                        aCGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGGCCTCGGGTCGTTGTCCActc                                                                                                                                                                                                                         <  2:690752/250‑1 (MQ=255)                             aaTGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTgg                                                                                                                                                                                                                   >  1:591360/1‑251 (MQ=255)                              attccTGCGTGGAAAGCTTCTGTCATTCGCTTAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCTGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCTCAGTTGGTGTGGCTGCCTCAGTCCGGAACAATTTGAAAACAAGAACCGCGCTTAGGCCTGTGTCCATATTACGTGGGTAGTATCATAAAGTCCTCGGGTGGTTGTCCACTCTTTTGGc                                                                                                                                                                                                                  <  1:481609/248‑1 (MQ=14)                              aTGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGc                                                                                                                                                                                                                  >  1:115913/1‑251 (MQ=255)                                gCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCCGTCCGGAACAATTTGAAAAAAAGAACCTCGGTTAGGCCTGTGTCCCTATTACGCGGGGAGGATCATAAAGGCCTCGGGGCGTTGTACACTCTTTTGGCAt                                                                                                                                                                                                                >  2:534755/1‑251 (MQ=14)                                    ctgTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTTTCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGTGGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATGTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGTTCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGt                                                                                                                                                                                                              <  1:675139/247‑1 (MQ=255)                                      gTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTACTTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGc                                                                                                                                                                                                          >  1:501055/1‑251 (MQ=255)                                        ggAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGttcgcat                                                                                                                                                                                                        >  2:80156/1‑245 (MQ=255)                                          aaaGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCatat                                                                                                                                                                                                      >  1:16030/1‑251 (MQ=255)                                           aaGCTTCCTTCATTCCCTGAAAGCGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGCGGCGGCGGCACCGTTGGTGTGGCGGCCTCAGTCCAGAAAAAATTGAAAAAAAGAACCTCTCTTAGGCCTGTGTCTATATTACGTGGGTAGGATCTTAAAGTCCTCGGGTCGTTTTCCACTCTTTTGGCATGTCAGCACAtt                                                                                                                                                                                                     >  2:67757/1‑251 (MQ=14)                                                   ttCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACAt                                                                                                                                                                                                <  1:411791/248‑1 (MQ=255)                                                       ttCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTAc                                                                                                                                                                                         >  1:487469/1‑251 (MQ=255)                                                                aGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCCTATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCCGCaaa                                                                                                                                                                                  >  1:535201/1‑249 (MQ=255)                                                                      aTGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCaa                                                                                                                                                                          <  2:665085/251‑1 (MQ=255)                                                                      aTGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCa                                                                                                                                                                           >  1:660440/1‑250 (MQ=255)                                                                         tATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCg                                                                                                                                                                       >  1:680068/1‑251 (MQ=255)                                                                                   gaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATc                                                                                                                                                             >  1:420356/1‑251 (MQ=255)                                                                                               aGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGAGGCGGCGGCACAGTTGGTGTGGCGGCCTCATTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGCACATAAAGTCCTCGGGTGGTTGGCCACTCTTTTGGAATGTCAGCATATACACATAACGGACCAGCAAAAGTTTCAACCGTGAGCAAATCAACCGATTttgt                                                                                                                                                 >  2:121760/1‑251 (MQ=18)                                                                                                   gggAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGAc                                                                                                                                             <  1:645201/251‑1 (MQ=255)                                                                                                   gggAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTGTCAACCGTGAGCAAATCGACCGATTTTGTTGAc                                                                                                                                             <  2:115913/251‑1 (MQ=255)                                                                                                                ggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTgg                                                                                                                                >  1:216871/1‑251 (MQ=255)                                                                                                                ggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTgg                                                                                                                                <  1:80156/251‑1 (MQ=255)                                                                                                                 gCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGa                                                                                                                               <  1:323726/251‑1 (MQ=255)                                                                                                                    aCGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAAc                                                                                                                            <  1:534755/251‑1 (MQ=255)                                                                                                                      ggTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGa                                                                                                                          <  2:16030/251‑1 (MQ=255)                                                                                                                         gtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAAAAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAAt                                                                                                                       <  2:28372/251‑1 (MQ=255)                                                                                                                                   atCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCg                                                                                                              <  2:353180/250‑1 (MQ=255)                                                                                                                                      gAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCGCGCTTATGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑Ctt                                                                                                         <  1:372204/251‑1 (MQ=255)                                                                                                                                                  aaTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCAt                                                                                             >  1:539267/1‑251 (MQ=255)                                                                                                                                                   aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGCTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATa                                                                                            <  1:77843/251‑1 (MQ=255)                                                                                                                                                                         gtgtGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGAcgtc                                                                      >  1:285665/1‑251 (MQ=255)                                                                                                                                                                             ggcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTa                                                                   <  1:147946/250‑1 (MQ=255)                                                                                                                                                                              ggggcCTCAGTCCGGAACAAGTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGt                                                                 <  2:535201/249‑1 (MQ=255)                                                                                                                                                                                 gcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTg                                                              >  1:626082/1‑251 (MQ=255)                                                                                                                                                                                              aCAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCg                                                 >  1:454326/1‑251 (MQ=255)                                                                                                                                                                                                                  ccGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGGTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAAGTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTGATTCGACACTTCa                             <  1:469698/250‑1 (MQ=255)                                                                                                                                                                                                                     cTTAGGCCTGTGTCCCTATTACGTGGGTAGGATCATCAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCCGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATg                          >  2:86795/1‑251 (MQ=255)                                                                                                                                                                                                                              gtgtCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTa                 <  2:660440/251‑1 (MQ=255)                                                                                                                                                                                                                                  ccATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATg             >  1:543373/1‑251 (MQ=255)                                                                                                                                                                                                                                    atatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCt                                                                                       <  2:524195/175‑1 (MQ=255)                                                                                                                                                                                                                                    atatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCt                                                                                       >  1:524195/1‑175 (MQ=255)                                                                                                                                                                                                                                    atatTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATg             >  1:491740/1‑249 (MQ=255)                                                                                                                                                                                                                                         aCGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGcc      <  2:216871/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGcc      >  2:538018/1‑251 (MQ=255)                                                                                                                                                                                                                                         aCGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGcc      >  1:434557/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGcct     <  2:501055/251‑2 (MQ=255)                                                                                                                                                                                                                                           gTGGTTAGGATCATAAAGTCCGCGGGTCGGTGTCCACGCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGGGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATATAGA‑GCGTGCTTAAAATTCTTCATACATCTGGGTAAATTATGACGTCGTAGTTTGACTTCTTCATCCGGTAATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCaa    <  1:647303/251‑1 (MQ=255)                                                                                                                                                                                                                                             gggTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCaaaa  <  2:543373/251‑1 (MQ=255)                                                                                                                                                                                                                                             gggTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCaaa   <  1:503455/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  TATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGT‑CTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAA  >  NC_000913/316237‑316723

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