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)	43 (0.970)	41/486	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)	70 (1.580)	62/486	0.0	100%	coding (1350/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GGTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCG  >  NC_000913/4499263‑4499758                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        ggTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAcc                                                                                                                                                                                                                                                       >  1:286000/1‑251 (MQ=11)      ggAGTCCAGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACTGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTGGTAAGGAGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTtgtg                                                                                                                                                                                                                                                  <  1:406034/251‑1 (MQ=255)           ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGat                                                                                                                                                                                                                                             <  2:181533/251‑1 (MQ=34)              gagaTTAACGGAATAGCAGAGAGCGTCGTGAAAACGATAAATCGTGACTACATCAGTTTCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCTTTCGAGCGTTATAACGAATGGCATCCGCATAGTGCGCTGGGTTAGCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATTGGTTAAGTGAGAACAGATGTCTGGAAATATAGGGGCAAATTCACCTTGTGCTGATATg                                                                                                                                                                                                                                          <  2:282395/251‑1 (MQ=11)                 agTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAAt                                                                                                                                                                                                                                       <  2:276310/251‑1 (MQ=37)                           tAGCAGAGAGCTTCGTGAAAACGATAAATCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTACTGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGGGCTGATATGAATACGGATAttt                                                                                                                                                                                                                             >  2:344527/1‑251 (MQ=255)                            aGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATAttt                                                                                                                                                                                                                             <  1:292760/250‑1 (MQ=255)                               agagagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAGCAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                                                                                                                                                         <  1:340396/251‑1 (MQ=255)                                           gAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACttt                                                                                                                                                                                                              <  1:390330/250‑1 (MQ=255)                                             aaaCGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCGTGTGCTGATATGAATACGGATATTTTCTACTGAGTACTTTTCa                                                                                                                                                                                                           <  2:212981/251‑1 (MQ=255)                                              aaCGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGAAACCGCATAGAGCGCTGGGTTATAGCTCGCCACGGGAATATCTGCAGCAGCGGGCTTGTAATGGGCTAAGTGATAACGGATGTCTGGAAATATAGGGGAAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                                                                                                                                                            >  1:449435/1‑249 (MQ=17)                                                cGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGt                                                                                                                                                                                                         >  1:86745/1‑250 (MQ=255)                                                 gATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTAc                                                                                                                                                                                                       >  1:345587/1‑251 (MQ=255)                                                     aaGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGGACAAGc                                                                                                                                                                                                   <  2:56159/251‑1 (MQ=255)                                                       gcatgactaactCAGTATCATGCCAGAAACAGACAGGTTAACGGCAGCAAAAAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCAGCATAGTGCGCTGGGCTCTCGCTCGACACACTATTATCTGCGGCAGCGGGCTTGTAATTGTTTAAGTGTTAATAGATGTCTGTAAGTATAGGGGACAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGATTACTTTTCAGTACAcgcgt                                                                                                                                                                                                 >  2:64142/12‑250 (MQ=12)                                                             cTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCCGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTAt                                                                                                                                                                                           >  1:116071/1‑251 (MQ=255)                                                                 aTCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCAGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                                                                                                                                                       <  1:363349/251‑1 (MQ=255)                                                                            gCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGGGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGTGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCTCCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGAGTGTATTGCTGTCCGCCatat                                                                                                                                                                            <  2:215347/251‑1 (MQ=17)                                                                              ccAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCatatat                                                                                                                                                                          >  2:352251/1‑251 (MQ=255)                                                                                     aGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTAcc                                                                                                                                                                   >  2:37807/1‑251 (MQ=255)                                                                                         gggTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGAAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                                                                                                                                                               >  2:196001/1‑251 (MQ=255)                                                                                          ggTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                                                                                                                                                               >  1:88184/1‑250 (MQ=255)                                                                                                            aCCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGAGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGc                                                                                                                                            >  1:51975/1‑251 (MQ=255)                                                                                                                  cAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGt                                                                                                                                      >  1:302070/1‑251 (MQ=255)                                                                                                                       gCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCaa                                                                                                                                 <  1:92745/251‑1 (MQ=255)                                                                                                                                    ataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGt                                                                                                                    <  1:257667/251‑1 (MQ=255)                                                                                                                                              gCATCCGCTTAGTGCGCGGGGTTACCGCTAGCCACGCGAATAGCTACGGCAGCGGGCTTGTAATGGGTTAAGAGATAACAGACGGCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAACACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAATATGTCTGCCGGTATATTCctct                                                                                                          >  2:379643/1‑251 (MQ=255)                                                                                                                                                  ccGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTcc                                                                                                      <  2:302070/251‑1 (MQ=255)                                                                                                                                                         gTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTaaaa                                                                                               <  2:345587/251‑1 (MQ=255)                                                                                                                                                          tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTaaaa                                                                                               >  1:142453/1‑250 (MQ=255)                                                                                                                                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCt                                                                                                                          >  2:190331/1‑219 (MQ=255)                                                                                                                                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCt                                                                                                                          <  1:190331/219‑1 (MQ=255)                                                                                                                                                                 ggTTATCGCGCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTACGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAAc                                                                                       >  2:291910/1‑251 (MQ=255)                                                                                                                                                                     aTCGCGCGCCACGGAAATATCTGTGGCAGCGGGCTTGTAATGGTTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATAGTCCTCTCTCCCTTAAAAACCATAACAgg                                                                                    <  1:203836/250‑1 (MQ=255)                                                                                                                                                                               aCGGGAATATCTGCGGCAGCGTGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCCCCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGCGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCtt                                                                         <  1:373085/251‑1 (MQ=255)                                                                                                                                                                                                aGCNGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTTATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCGGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTGCGTCTGGCAGCGGATTg                                                        <  2:410472/251‑1 (MQ=255)                                                                                                                                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGaa                                                      <  1:123900/251‑1 (MQ=255)                                                                                                                                                                                                     gCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATgg                                                   <  1:162521/251‑1 (MQ=255)                                                                                                                                                                                                         gTAATGGGTTTCGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGt                                               >  2:285722/1‑251 (MQ=255)                                                                                                                                                                                                                                                       ccACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCAACTGAAAAATACCTGCg  >  2:370951/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        GGTGCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCG  >  NC_000913/4499263‑4499758

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