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 =	1 (0.080)	12 (0.920)	11/290	0.1	92.0%	coding (1346/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
*	?	NC_000913	= 4135184	1 (0.080)	27 (2.070)	21/290	0.0	96.3%	coding (1350/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGAT  >  NC_000913/4134957‑4135415                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             gctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGc                                                                                                                                                                                                                   <  1:258718‑M1/28‑1 (MQ=255)      cgaccTGACGCTGCGTTTTGATCCTGAGTTCGCGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGACGCCTTTGCCCGTGACTGGTTCAAACTGACGCACAGGGATCTGGGGACGACATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCCAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGAtt                                                                                                                                                                                                              >  1:45689/1‑251 (MQ=255)                               gAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagc                                                                                                                                                                                     <  1:209547‑M1/251‑54 (MQ=255)                                   tCGAGAATATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCGGATCTGGCAAGAGCCGCTGCCGCAGCCGAGCTACAACCCGACCGAGCAGTACATTtgaactggcccctatattgccagacatctgttatcacttaacccattacaagcccg                                                                                                                                                                                  <  1:291210‑M1/250‑57 (MQ=255)                                           catccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTg                                                                                                                                                                        >  1:69737‑M1/182‑251 (MQ=255)                                                                      cAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccag                                                                                                                                              <  1:104755‑M1/251‑93 (MQ=255)                                                                      aattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGtggtgg                                                                                                                                             >  1:107604‑M1/155‑251 (MQ=255)                                                                              cNACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgcc                                                                                                                                                                             >  1:253579‑M1/1‑151 (MQ=255)                                                                              cAACGAAGCCTTTGCCCGTGCCTGGCTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCACGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagccatctgttatcacttaacccattacaagcccgctgccgcagat                                                                                                                                                                       >  1:123831‑M1/1‑151 (MQ=255)                                                                                              acgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccaggctaACATTAGCTATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTAGCGGTGGCCTGGGCATCTGCTGCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACgg                                                                                                                     >  1:305556‑M1/131‑251 (MQ=255)                                                                                                     ggTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatg                                                                                                                <  1:129242‑M1/250‑123 (MQ=255)                                                                                                        tCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgc                                                                                                                                                                                                                                 <  1:142940‑M1/134‑10 (MQ=255)                                                                                                                   cacaGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAACGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgtta                                                                                                                                                                                                         >  1:52257‑M1/1‑114 (MQ=255)                                                                                                                       gggATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttcacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgaggatgccattcgttataat                                                                                                                 >  1:237949‑M1/1‑110 (MQ=255)                                                                                                                                            ccaatgactagtctaaaaactagtattaaggctatcacttatttaagtgatattgggtgtctggagattcagggggccagtctaACATTATCGATCTTAAAAGCGCGATTGCTGATTC‑GGTGCTGTCTGTTAGTGAGCTTGTATCGGTGGCCTGTGCATCTGCTTCGACCGTCCGTGGTGGCGACAAACGCAGTGGTTCCAACGGTGCTCGTCTGGCATTAATGCCGCAGCGCGACTGATATGTGAAcgc                                                                        >  1:317116‑M1/85‑250 (MQ=255)                                                                                                                                                      gCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgc                                                                                                                                    <  1:158847‑M1/181‑103 (MQ=255)                                                                                                                                                                       aagaTCTGATCTGGCAAGATCCGCTGNCGCAGCNNATNTACAACCCGACCGAGCAGGNCANNngnnntngccccnatnnttccagacanctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgancgataacccnn                                                                                                                                              >  2:76560‑M1/1‑62 (MQ=255)                                                                                                                                                                       aagaTCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccag                                                                                                                                              <  1:76560‑M1/154‑93 (MQ=255)                                                                                                                                                                              gatctgGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggagttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcatatattcccgtggcgagcgagaatccagcgcactaggcgggtgccattcgttataatgctcgaacgcctctgcaaggttctttgctgcc                                                                                 <  1:151710‑M1/208‑154 (MQ=255)                                                                                                                                                                               cacttatttaagtgatattggttgtctggagattccgggggccagtctaACATTATCGCTCTGAAATTCGCGATTGCAGATTCTgg                                                                                                                                                                                                          <  1:277503‑M1/37‑1 (MQ=255)                                                                                                                                                                                 cttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCt                                                                                                                                                              <  1:99360‑M1/80‑1 (MQ=255)                                                                                                                                                                                  tgGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGACCAGGACCTTtagatttgcctctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattaacgaggagagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggttagggcatgatactgatggagtcacgctttat                                  >  1:301581‑M1/1‑51 (MQ=255)                                                                                                                                                                                   ggCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataattctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttggg                                                             <  1:211511‑M1/223‑174 (MQ=255)                                                                                                                                                                                   ggCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataaccc                                                                                                                                                >  1:133300‑M1/1‑50 (MQ=255)                                                                                                                                                                                    gCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGCGCAGGACATTtgactttgcccctatatttccagacatctgttatcacttaacacataacaagcccgcagccgaagatctcaccgtggagcgcgataaaccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaccccgtttggtttgggcatgatactgatgtagtcacgctttatcg                                >  1:221554‑M1/1‑49 (MQ=255)                                                                                                                                                                                    gCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacacgcccgctgccgcagatattcccgtggcgagcgataacccagcgcactacgcggatgccattcgttataatgctcgaacgcctctg                                                                                                  >  1:317218‑M1/1‑49 (MQ=255)                                                                                                                                                                                     cAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggt                                                                                            >  1:253039‑M1/1‑48 (MQ=255)                                                                                                                                                                                             gcggccgCAGCCGATCTACAACCCGACCGAGCAGGACATTttgatgtgcccctatatttccagacatctgttatcactgaacccattacaagcccgctgccgcagatattcccgt                                                                                                                                                               <  1:210496‑M1/112‑76 (MQ=255)                                                                                                                                                                                                gCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagccctctgccgcagatattcccgtggc                                                                                                                                                            <  1:92045‑M1/115‑79 (MQ=255)                                                                                                                                                                                                   ggtgtctggagattcagggggccagtctaTCATTATCGATCTGAAATGCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATc                                                                                                                                                                  >  1:213791‑M1/30‑105 (MQ=255)                                                                                                                                                                                                      gtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCGGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGGGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCg                                                            >  1:100888‑M1/27‑204 (MQ=255)                                                                                                                                                                                                        cGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatn                                                                                                                                                                                                       <  1:83241‑M1/64‑36 (MQ=255)                                                                                                                                                                                                               cAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagtcatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtcgggtttgggcatgatactgatgtagtcan                                         <  1:41453‑M1/215‑194 (MQ=255)                                                                                                                                                                                                                tcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTAAGTGAGATGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGGGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTgg                                           >  1:325560‑M1/17‑211 (MQ=255)                                                                                                                                                                                                                  cccGACCGAGCAGGACATTtggatttgcccctatatttcaagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggagagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttac  >  1:320587‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                  accgaccgagCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttac  >  1:92760‑M1/2‑19 (MQ=255)                                                                                                                                                                                                                        ccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGGGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGctct                                                      >  1:131202‑M1/9‑192 (MQ=255)                                                                                                                                                                                                                           gtctaACATTATCGATCTGAAATTCGCGATTGCggattatg                                                                                                                                                                                                           <  1:159818‑M1/36‑4 (MQ=255)                                                                                                                                                                                                                           gCAGGACATTtggatttgcccctatatttacagacatctgttatcacttatccccttaaaagaccgctgccgccgatattcccgtggcgagcgataaccaagcgcactatgcggatgccattcgttataatgctcgaacgcctctgccaggttctttgctgccgt                                                                               >  1:282299‑M1/1‑10 (MQ=255)                                                                                                                                                                                                                              ggACATTtggatttgcccctatctttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatac                                                     >  1:92334‑M1/1‑7 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             GTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGT‑CTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGAT  >  NC_000913/4134957‑4135415

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.