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)	38 (0.820)	36/488	0.3	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)	27 (0.580)	24/488	0.8	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 577690	0 (0.000)					coding (78/498 nt)	rrrD	lysozyme

AACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTC  >  NC_000913/566527‑567020                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aaCTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGT‑AAAACAACTTCCTGGTATGGGCGGTAGACCTCTAAAGCGGGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAGGGAGCAAGTGTATGACGGGCAAAGAGGCAATGATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGGCGTTGCCGCGGTAACAGGCTCAACCGTAATCt                                                                                                                                                                                                                                                     <  2:465401/251‑1 (MQ=255)                      gAATTAACCGTGGTGAGACGATACCTGAAACAGT‑AAAACAACTTCCTGTTATGGGCGGTAGAACTCTAAATCGTGCCCAGGCTCTGGCGACGATCGCCGAAATTACAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTCCCTGCGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTgg                                                                                                                                                                                                                                >  1:432393/1‑251 (MQ=255)                        aTTAACCGTGGTGAGACGATACCTGAACCAGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGac                                                                                                                                                                                                                              <  2:95034/251‑1 (MQ=255)                            aCCGTGGTGAGACGATACCTGAACCAGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACAcgc                                                                                                                                                                                                                          <  2:313447/251‑1 (MQ=255)                             ccGTGGTGAGACGATACCTGAACCAGT‑AAAACAACTTCCTGTTAAAGCCCGTCGACCTCTAAATCATGCACCGGCTCCGGCCACGATACAAGACATAAAAGCTAAGTTCGGACTGACAGGAGCAACTGTATGCCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAAAGCTTCTGTGCACAGGACGCTGACGCGGTAACAGGAGCAACCGTAATATGATCTTACCCAGCAATAGAGGacacccg                                                                                                                                                                                                                         >  1:449745/1‑248 (MQ=255)                                         gATACATGAACCACTCAAAACA‑CTACCTGTTATGGGCGGTAGACCTCTAAATCCTGCACCGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGCAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATGACTTGGGGCCGCATAAGAGCTTCTGTGCACAGGACGTTGACGCGGTAACAGGCGCAACCGTAATCCGTGCTTACCCAGCAATAGTGGAAACGCGGCTAAGCGAga                                                                                                                                                                                                              >  2:39250/1‑249 (MQ=255)                                                 aaCCAGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtca                                                                                                                                                                                                     <  2:120314/251‑1 (MQ=255)                                                   ccAGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                                                                                                                                                   <  1:292857/251‑1 (MQ=255)                                                     aGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAACTCGTGCACAGGCTCTGGCGAAGATCGCAGACATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtc                                                                                                                                                                                                  >  2:409784/1‑250 (MQ=255)                                                              aacTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGAc                                                                                                                                                                                         >  1:55371/1‑251 (MQ=255)                                                                           gggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGGGACTCACATGACaaa                                                                                                                                                                             <  1:158213/250‑1 (MQ=255)                                                                               ggTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACACAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGCGGTGACTCACATGACAAAAAc                                                                                                                                                                          >  2:112015/1‑249 (MQ=255)                                                                                        ctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCGCAGGCAGAGGTGACTCACATGACAAAAACAGTATCAAc                                                                                                                                                                 <  2:203702/249‑1 (MQ=255)                                                                                         ttaaTCGTGCACATGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGCGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATAACAAAAAAAGTATCAAcc                                                                                                                                                                >  2:103525/3‑249 (MQ=255)                                                                                           aaTCGTGCACAGGCTCTGGCGAAGATTGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTa                                                                                                                                                            <  1:109811/251‑1 (MQ=255)                                                                                                           tGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAAc                                                                                                                                             >  1:13261/1‑250 (MQ=255)                                                                                                                 aGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGCGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTc                                                                                                                                      >  1:228040/1‑251 (MQ=255)                                                                                                                 aGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTc                                                                                                                                      <  1:434631/251‑1 (MQ=255)                                                                                                                      gCAGAAATTAAAGCTACGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACCCGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTg                                                                                                                                 >  2:254610/1‑251 (MQ=255)                                                                                                                          aaaTTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAAtt                                                                                                                             <  2:83659/251‑1 (MQ=255)                                                                                                                               aaaGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGc                                                                                                                         <  1:100046/250‑1 (MQ=255)                                                                                                                                    tAAGTTCGGACTGAAAGGAGCAAGTGTATGTCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAg                                                                                                                                                              >  1:342862/1‑208 (MQ=255)                                                                                                                                    tAAGTTCGGACTGAAAGGAGCAAGTGTATGTCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAg                                                                                                                                                              <  2:342862/208‑1 (MQ=255)                                                                                                                                     aaGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAg                                                                                                                  <  2:55371/251‑1 (MQ=255)                                                                                                                                        ttCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGccc                                                                                                               >  1:415161/1‑251 (MQ=255)                                                                                                                                                gAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCtt                                                                                                       <  1:139952/251‑1 (MQ=255)                                                                                                                                                     gAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTg                                                                                                   >  1:382799/1‑250 (MQ=255)                                                                                                                                                     gAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAACTAGGAAGGCATACGAGCTTCTGCGCACAGCACGTTGACGCGGTACCAGGAGCAACCGCAAACTGCCCTTACCAAACACTAGTGGACACGCGGCTAAGTGAGCAAACACACAGTCAGCGGTGACTCACATGACAAAAACAGTATCAACCAGTACCAAACCCCGTAAACAGCATACGCCTGAATTTCGCAGTGAAGCCCTGAAGCTCGCTg                                                                                                   >  1:417983/1‑250 (MQ=12)                                                                                                                                                         aagtgtataacacgcaacggcgCAATTATTCACTACCCGAGGACGCATAAGAACTTCTGTGCACCGGACGTTGCCGCGGCCACAGGCGCAAACGTAATCTGATCTTACCCAGCAATAGTGGACCCGCGGCTAAGTGAGTAAAATCTCAGTAAGCGGTGACTCACAAGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCATGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGc                                                                                              >  2:281655/22‑251 (MQ=17)                                                                                                                                                                                 attattCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAACAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGCATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCgg                                                                      >  2:387757/1‑250 (MQ=255)                                                                                                                                                                                     ttCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGATGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCCGCCCGTGAAc                                                                  >  1:82686/1‑251 (MQ=255)                                                                                                                                                                                       cATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTc                                                                >  1:190371/1‑251 (MQ=255)                                                                                                                                                                                        aTTACCTGAGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCGGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCGGAACGCATCGGTGTGACTGCCGCAGCCAGTGAACTCa                                                               <  2:82686/251‑1 (MQ=255)                                                                                                                                                                                                 ggACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATg                                                      >  1:387735/1‑251 (MQ=255)                                                                                                                                                                                                                                 gCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAAc                        >  2:342270/1‑249 (MQ=39)                                                                                                                                                                                                                                            aGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTTGCCCAGTAAACAGCAAAATcagca           <  1:388507/251‑1 (MQ=18)                                                                                                                                                                                                                                                   aCCGTACTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTACACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGCAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAACTCAGCAGACGtct    >  2:3372/1‑251 (MQ=9)                                                                                                                                                                                                                                                     cGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGtcttc  <  1:409784/251‑1 (MQ=11)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    AACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGT‑AAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTC  >  NC_000913/566527‑567020

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