Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	606,417	IS4 (–) +12 bp	coding (956‑967/1119 nt)	ybdK ←	putative glutamate‑‑cysteine ligase 2

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	606417 =	0 (0.000)	25 (0.870)	23/484	0.3	100%	coding (967/1119 nt)	ybdK	putative glutamate‑‑cysteine ligase 2
	?	NC_000913	= 4503515	NA (NA)					noncoding (1/1426 nt)	IS4	repeat region
*	?	NC_000913	= 606428	0 (0.000)	17 (0.600)	17/482	0.5	100%	coding (956/1119 nt)	ybdK	putative glutamate‑‑cysteine ligase 2
	?	NC_000913	4502090 =	NA (NA)					noncoding (1426/1426 nt)	IS4	repeat region

CTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGCCAGTAACCTT  >  NC_000913/4501844‑4502338                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         tctaaagggtgtttcataacattgttaaatgtaagttaaaacaccattgtcagggatattcttctgtaaggcaattcccggcttagtcaccggcccagatctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcatgcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAAt                                                                                                                                                                                                                                                       <  2:107232‑M2/4‑1 (MQ=255)    aaagggtgtttaataacattgttaaatgtaagttaaaacaccattgtcagggatattcttcagtaaggcaattcccggcttagtcaccggcccagatctcacaatgctttttcaccagcacaatcagcgagccgccatcggcgacgaaagcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCca                                                                                                                                                                                                                                                   <  2:134282‑M2/8‑1 (MQ=255)    aaaaaTGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGcc                                                                                                                                                                                                                                                    >  1:190879/1‑250 (MQ=255)        ggtgtttcataacattgttaaatgtaagttaaaacaccattgtcagggatattcttctgtaaggcaattcccggcttagtcaccggcccagatctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCcagtca                                                                                                                                                                                                                                               <  2:268076‑M2/12‑1 (MQ=255)                         tCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactgactg                                                                                                                                                                                                                               <  1:10521/250‑1 (MQ=255)                                 aagttaaaacaccattgtcagggatattcttcggtaaggcaatccccggcttagtcaccggcccagatctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcgccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCttttt                                                                                                                                                                                                                      <  2:235804‑M2/37‑1 (MQ=255)                                       ccAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGc                                                                                                                                                                                                                >  1:31545/1‑251 (MQ=255)                                                      aGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCATTATTCTCGTTGGCAATCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCTTCGTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTcc                                                                                                                                                                                                  <  1:300119/250‑1 (MQ=255)                                                      aGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTcc                                                                                                                                                                                                  <  1:265440/250‑1 (MQ=255)                                                                      ggCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAATCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTa                                                                                                                                                                                   <  1:3921/249‑1 (MQ=255)                                                                             gCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCg                                                                                                                                                                          <  1:281177/251‑1 (MQ=255)                                                                                ggcttagtcaccggcccagatctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGa                                                                                                                                                                       <  1:167680‑M2/84‑1 (MQ=255)                                                                                             cTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCct                                                                                                                                                          >  1:162639/1‑251 (MQ=255)                                                                                                tAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAACATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCCGATAATTTAACAGTTGATGCTGAAACGACATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAACAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCAAGGGCCTCTCCTTTACCACTCTCTGGAAGTCCCTTTCCcttct                                                                                                                                                       >  1:307532/1‑251 (MQ=255)                                                                                                cagatctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCcttct                                                                                                                                                       >  1:22665‑M2/152‑251 (MQ=255)                                                                                                   atctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTa                                                                                                                                                                                                                                           >  1:142156‑M2/149‑164 (MQ=255)                                                                                                   atctcacaatgctttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgaccgggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGGCAGTTa                                                                                                                                                                                                                                           <  2:142156‑M2/16‑1 (MQ=255)                                                                                                         tGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGCTAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTa                                                                                                                                               >  2:211992/1‑250 (MQ=255)                                                                                                               tttttcaccagcccaatcagcgagccgccatcggcgacgaaatcgcgccttagctgcgcttcattcagaccgctgacgccctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGTGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCTGTAATTTcac                                                                                                                                        >  2:101878‑M2/137‑251 (MQ=255)                                                                                                                         gcccaatcagcgagccgccatcggcgacgaaatcgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGGCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTccc                                                                                                                              <  2:179519‑M2/125‑1 (MQ=255)                                                                                                                               cAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACtt                                                                                                                        >  1:51500/1‑251 (MQ=255)                                                                                                                                   cgagcagccatcggcgacgaaaacgcgcattagctgcgcttcattcagaccgctgacgaccaggcgatgcagggactcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGACAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCATCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGGCGGTACTTTCACAAGTTGTCCCATACTTGCa                                                                                                                     >  2:206216‑M2/117‑250 (MQ=255)                                                                                                                                                gggCGTTGGCGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAg                                                                                                       <  1:182658/250‑1 (MQ=255)                                                                                                                                                    gTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTACTGCCAGGCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTGCCTTCTTGTCGGGAATTTCACAAGTTGTCCCCTACTTGCAAGATCGCGCATCAGCTcc                                                                                                   >  1:162295/1‑251 (MQ=255)                                                                                                                                                      tGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCg                                                                                                  <  2:66171/250‑1 (MQ=255)                                                                                                                                                          cgcgcattagctgcgcttcattcagaccgctgacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGtat                                                                                              <  1:144021‑M2/157‑1 (MQ=255)                                                                                                                                                                                gCTCAGCATTCTCGTTGGCAAACAGATCATTTAACAGTTGATGCTGAAACGCAATAGGCATCTTGGTTCTGTAATGCAAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTACCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGc                                                                        >  2:132326/1‑250 (MQ=255)                                                                                                                                                                                   cAGTATTATCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAATCGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATg                                                                    <  1:18635/251‑1 (MQ=255)                                                                                                                                                                                     gTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGtc                                                                  <  2:51500/251‑1 (MQ=255)                                                                                                                                                                                         tctcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGtcatc                                                               >  1:266374/1‑250 (MQ=255)                                                                                                                                                                                         gacgacctggcgatgcagggcctcaatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGtcat                                                                <  1:98977‑M2/187‑1 (MQ=255)                                                                                                                                                                                            cGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCCTCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCgcgc                                                                                                            >  2:174974/1‑202 (MQ=255)                                                                                                                                                                                            cGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCgcgc                                                                                                            <  1:174974/202‑1 (MQ=255)                                                                                                                                                                                            cGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCa                                                           <  2:85614/251‑1 (MQ=255)                                                                                                                                                                                                                 aatcgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATg                                      <  1:215342‑M2/213‑1 (MQ=255)                                                                                                                                                                                                                   aGTTGATGATGAAACGAAATAGGCATCGTGGTGCTGTAATGCCAGTCAGTTAAGCAACTGCCTGGCTCTTTTTCGGGGCTGTGGGGTAGTTCCAGTGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGAt                                    <  2:52297/251‑1 (MQ=255)                                                                                                                                                                                                                    cgcgctcgatgcaccaattttatgtgcggaaggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGAt                                    <  1:153430‑M2/215‑1 (MQ=255)                                                                                                                                                                                                                                    aaTAGGCGTCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTg                   >  2:273431/1‑251 (MQ=255)                                                                                                                                                                                                                                       aGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGTTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCTCCATTCCGCATGATTCTGAGAAACTCAGTTGAtt                <  1:132326/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aggggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCTCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCAGCTCCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGCCAGTAAc     <  1:206216‑M2/246‑1 (MQ=255)                                                                                                                                                                                                                                                     ggTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGCCAGTAACCtt  >  1:52692‑M2/3‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         CTAAAAAATGCGTTAATGTTGAGTTCCTCGAACCCTTACCAATTAAGGGACATAGCCAGATACGCAACTGGCTGATGCAAGCAGTCTAATGTCTTTAGATCCGCTGGCCTGCGGGCATGACAATGTCAAACAGTCAAAGATTGTGGCGTTGACGGTCAGTCTTGGGGACAAAAAAGCTCAGTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCCTTTACCACTCTCGGGAAGGCCCTTTCCCTTCTTGTCGGTAATTTCACAAGTTGTCCCATACTTGCAAGATCGCGCATCAGCTCCGGTATACGTCCCGGTGAAGCGCCCTGCAATGTCATCAGCATTCTCATCACCATTCCGCATGATTCTGAGAAACTCAGTTGATTCGGCCAGTAACCTT  >  NC_000913/4501844‑4502338

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