Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,135,180	IS2 (–) +5 bp	coding (1346‑1350/2181 nt)	katG →	catalase/hydroperoxidase HPI

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4135180 =	0 (0.000)	40 (1.030)	36/490	0.2	100%	coding (1346/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 4135184	0 (0.000)	41 (1.060)	39/490	0.2	100%	coding (1350/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

AGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTG  >  NC_000913/4134934‑4135427                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACAtt                                                                                                                                                                                                                                                     <  1:251782/251‑1 (MQ=255)       gTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCTAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAAATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatt                                                                                                                                                                                                                                               <  1:295190‑M1/251‑7 (MQ=255)        tAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggattt                                                                                                                                                                                                                                              >  2:395598‑M1/1‑244 (MQ=255)         ccagagaaataactggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCt                                                                                                                                                                                                                                             >  2:271101‑M1/239‑251 (MQ=255)          aaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgc                                                                                                                                                                                                                                            <  1:334680‑M1/251‑10 (MQ=255)            gagaaataactggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGa                                                                                                                                                                                                                                           <  1:228189‑M1/15‑1 (MQ=255)              gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccct                                                                                                                                                                                                                                        >  1:327637‑M1/1‑238 (MQ=255)              gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccct                                                                                                                                                                                                                                        >  2:69300‑M1/1‑238 (MQ=255)                  aTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTAGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGCCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatat                                                                                                                                                                                                                                    >  2:162466‑M1/1‑234 (MQ=255)                      tGGTGACCGACCTGACGATGCGTTTTGATCCTGAGCTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCATGGTTCAAACTGACGCACAGGGATATGGGGCAGAAATATCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATATACAACCCGACCGAGCAGGACATTtggatttgcccctatatttcc                                                                                                                                                                                                                                >  1:296467‑M1/1‑230 (MQ=255)                       gtctggctgctacatcatgtttccgggcaacgatggagaccgtcatccccggttcaaagatctgctgaacacttgcgatcttttcctgtgtggtacgcagtctgcgtttctccggccctaagacatcaatcatctgttctccaatggctagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcaggtgtccagtctaACATTATCGATCTGAAATTCGCGATTg                                                                                                                                                                                                                               >  2:258288‑M1/225‑251 (MQ=255)                       ggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatctcttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTg                                                                                                                                                                                                                               >  1:236196‑M1/225‑251 (MQ=255)                       ggTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttcc                                                                                                                                                                                                                                <  2:332580‑M1/250‑22 (MQ=255)                          tggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGAt                                                                                                                                                                                                                                 >  1:51054‑M1/222‑246 (MQ=255)                          tggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGAt                                                                                                                                                                                                                                 <  2:51054‑M1/25‑1 (MQ=255)                            ccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCGGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacat                                                                                                                                                                                                                          <  2:216094‑M1/251‑28 (MQ=255)                              gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatct                                                                                                                                                                                                                        >  2:128321‑M1/1‑222 (MQ=255)                               gcgacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactaggattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTc                                                                                                                                                                                                                        >  2:64263‑M1/217‑250 (MQ=255)                                     cGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatc                                                                                                                                                                                                                  >  1:354732‑M1/1‑215 (MQ=255)                                      gCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcac                                                                                                                                                                                                                <  1:307009‑M1/251‑38 (MQ=255)                                               gggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTg                                                                                                                                                                                                       <  2:376454‑M1/51‑1 (MQ=255)                                                  tCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAAAGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTgggatttgcccctttatttccagacatctgttatcacttaacccantac                                                                                                                                                                                                    >  2:354404‑M1/1‑202 (MQ=255)                                                             agaagaTCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgc                                                                                                                                                                                         <  2:89972‑M1/251‑61 (MQ=255)                                                                  catccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCt                                                                                                                                                                                     >  2:322977‑M1/182‑250 (MQ=255)                                                                  catccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTg                                                                                                                                                                                    <  1:17530‑M1/70‑1 (MQ=255)                                                                   tctctcGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgc                                                                                                                                                                                      >  1:169413‑M1/1‑185 (MQ=255)                                                                             tttCCTCAACGATCCGAAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtg                                                                                                                                                                         >  1:341495‑M1/1‑175 (MQ=255)                                                                             caaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCtt                                                                                                                                                                         >  2:25322‑M1/171‑251 (MQ=255)                                                                              ttCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtgg                                                                                                                                                                        >  2:200580‑M1/1‑174 (MQ=255)                                                                              caaagctctgctgaacaattgcgatcttttcctgtggggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcaggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTc                                                                                                                                                                        <  2:4897‑M1/82‑1 (MQ=255)                                                                                ccTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggc                                                                                                                                                                       <  2:183787‑M1/250‑79 (MQ=255)                                                                                agctctgctgaacaattgcgatctgtgcctgtgtggtacgccgtctgcgtttcgccggccctaagacatcaatcatctgttctccaatgacgagtctaaaaactagtattaagactatcacttatttaagtgatataggttgtctggagattcaggaggccagtctaACATTTTCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCCGGTAGTGAGCTGGTATCGGGGGGCTGGGCATCTGCTTCTa                                                                                                                                                                      >  2:171934‑M1/168‑251 (MQ=255)                                                                                           cGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataaccc                                                                                                                                                           >  2:389119‑M1/1‑161 (MQ=255)                                                                                                  gTTCAACGAAGCCTTTGCCCGAGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaaccaattacaagcccgctgcagcagatattcccgtggcgagcgataaccaagccaa                                                                                                                                                     >  2:376865‑M1/1‑154 (MQ=255)                                                                                                             tgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCggtggt                                                                                                                                         <  1:378388‑M1/113‑1 (MQ=255)                                                                                                                              gcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCt                                                                                                                        >  1:376454‑M1/122‑251 (MQ=255)                                                                                                                                  aaCTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaac                                                                                                                    <  2:169413‑M1/251‑130 (MQ=255)                                                                                                                                           cctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCg                                                                                                           <  1:351664‑M1/143‑1 (MQ=255)                                                                                                                                                  catcactcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGcgcg                                                                                                    <  2:377795‑M1/150‑1 (MQ=255)                                                                                                                                                   atGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttct                                                                                                    <  2:194479‑M1/250‑146 (MQ=255)                                                                                                                                                     ggggCCGAAATCTCTCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTTATCGGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGTCCGAGCAGGACATTtggatttgcccctatatttccagacatcttttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgctctatgcggatgccattcgttataatgctcgaacgcctctgcaaggttcttt                                                                                                  <  1:56454‑M1/250‑148 (MQ=255)                                                                                                                                                       atcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggtgccagtctaACATTATCGATCTGAAATTCGCGAtt                                                                                                                                                                                                                                >  2:245571‑M1/97‑122 (MQ=255)                                                                                                                                                       atcatctgttctccaatgactagtctaaaaactagaattaagactatcccatatttatgtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGAtt                                                                                                                                                                                                                                <  1:245571‑M1/26‑1 (MQ=255)                                                                                                                                                                  tctaatgcatagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccattctaACATTATCGATCTGAGATTCACGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCGGGGCATCTGATTCTACCTTCCGTGGTGGCGACAAACGCGTTGGTGCCAACGGTGCGCGTCTGGGATTAATGCCTCCGCGCgactgggataggaacgc                                                                                    <  2:45414‑M1/166‑9 (MQ=255)                                                                                                                                                                     aatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAAcgccgc                                                                                 <  2:132710‑M1/169‑1 (MQ=255)                                                                                                                                                                       cATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctg                                                                               <  2:350027‑M1/251‑167 (MQ=255)                                                                                                                                                                            agtctaaaaactagtattaatactatcgattaattaattgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCGGTCTGTTAGTGTGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCTACAAACGCGGTGGTGCCACCGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGTTGTGAACGCCGCAACCGtt                                                                          >  2:96226‑M1/76‑251 (MQ=255)                                                                                                                                                                               ctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGt                                                                       <  1:100022‑M1/179‑1 (MQ=255)                                                                                                                                                                                taaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGCACGCCGCAGCCGTTCGTg                                                                      <  2:3649‑M1/180‑1 (MQ=255)                                                                                                                                                                                taaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTg                                                                      <  2:1523‑M1/180‑1 (MQ=255)                                                                                                                                                                                   aaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGctc                                                                   <  1:52685‑M1/183‑1 (MQ=255)                                                                                                                                                                                    aactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGTGCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTGCCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCTGCCGTtcgtgttct                                                                  >  2:94754‑M1/68‑247 (MQ=255)                                                                                                                                                                                      ctagtattaagactataacttatttaagtgatattggttgtctggagattcagggggccagtctaACCTTATCGATCTGAAATTCGCGATTGCGGATTCTAGACTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGAGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCtgg                                                                <  1:87942‑M1/186‑2 (MQ=255)                                                                                                                                                                                        agtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCt                                                              <  1:124197‑M1/188‑1 (MQ=255)                                                                                                                                                                                             gaagaTCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgta                                                         <  1:162466‑M1/251‑189 (MQ=255)                                                                                                                                                                                              aagaTCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtag                                                        <  2:303749‑M1/251‑190 (MQ=255)                                                                                                                                                                                                   tatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAAc                                                                                                                                                 <  1:24259‑M1/105‑1 (MQ=255)                                                                                                                                                                                                   tatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAAc                                                                                                                                                 >  2:24259‑M1/53‑157 (MQ=255)                                                                                                                                                                                                     gatctgGCAAGAACCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcac                                                    <  1:257161‑M1/248‑194 (MQ=255)                                                                                                                                                                                                      cacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATATGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGAATATGCTTCTACCTTCCGGGGTGGCGACAAACGCGGGGGTGCCAACGGTGAGCGTATGGCATTAATGCAGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAAt                                                <  2:226409‑M1/202‑1 (MQ=255)                                                                                                                                                                                                          tatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATcc                                              >  2:117794‑M1/46‑249 (MQ=255)                                                                                                                                                                                                          tatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCATGTTCTGGAGAAAATcc                                              <  1:206815‑M1/204‑1 (MQ=255)                                                                                                                                                                                                          ggCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatc                                            >  1:20292‑M1/1‑50 (MQ=255)                                                                                                                                                                                                           gCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcg                                           <  1:200580‑M1/251‑203 (MQ=255)                                                                                                                                                                                                              aGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttt                                        >  1:267548‑M1/1‑46 (MQ=255)                                                                                                                                                                                                                  tgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTc                                    <  1:258823‑M1/214‑1 (MQ=255)                                                                                                                                                                                                                    gCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacg                                   >  1:255617‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                             gCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcaattaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctg                         >  1:232614‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                              ccGATCTACAACCCGACCGCGCAGGACATTtggatttgcccatatatttcaagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttctaatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttctcggtttcacgcagctctctgc                        >  2:100929‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                   cTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttctcacttaacccattacaagcccgccgccgcagatatccccgcggccagcgataacccagcgcactatgcggatgccattcggtataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgctacagatggagtaacgctttatcgctttcacgaagctctctgctattc                   >  2:65979‑M1/1‑25 (MQ=255)                                                                                                                                                                                                                                   agattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTgg                   <  2:260895‑M1/231‑1 (MQ=255)                                                                                                                                                                                                                                      ttcagggcgccagtctaaCATTCTCGCTCTGAAACTCCCGATTGCGGCTTCCGCTCTGTCTGTTAGTGAGCTGGTCTCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTCGCGGCCACAAACGCGGCGGTGCCAACGGTGCGCGTCTGGCATTACTGCCGCAGCGCGACTGGGATCTCAACGCCGCACCCGTTCGTGCCCTGCCTGTCCTGGCGAAAATCCAGCAAGAGTCTGGTAAAACCTCCCTGGCgg                >  2:156073/18‑251 (MQ=255)                                                                                                                                                                                                                                       aaCCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggatctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgtt               >  2:317095‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                        cagggggccagtctaACATTATCGATCTGAAATTCGCGCTTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGat              <  2:33856‑M1/236‑1 (MQ=255)                                                                                                                                                                                                                                        aCCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatgttcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgtt               >  2:246207‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                        aCCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgtta              <  1:297977‑M1/251‑232 (MQ=255)                                                                                                                                                                                                                                          ccgaccgaGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagcgctctgctattccgttact            <  2:141439‑M1/251‑234 (MQ=255)                                                                                                                                                                                                                                           cgaccgaGCAGGAAATTtgtagttgccactatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggtgctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttactc           <  2:325664‑M1/251‑235 (MQ=255)                                                                                                                                                                                                                                            gggcttgtctaACATTATCGAGCTGACATTCGCGATTGCGGACTCTGGTCTGTCTGGTAGTGCGCTGGTATCGGTGGTCTGGGAATCTGCTTCTACCTCACGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTTTCATTAATGCTGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCa          >  1:386999‑M1/12‑251 (MQ=255)                                                                                                                                                                                                                                            gaccgaGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttactct          >  2:292237‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                                            gaccgaGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttactct          >  2:103677‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                                                cagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGACGCAGCAGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGt      <  2:222218‑M1/244‑1 (MQ=255)                                                                                                                                                                                                                                                    ctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTg  <  1:359706‑M1/248‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        AGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTG  >  NC_000913/4134934‑4135427

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.