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)	9 (0.890)	7/298	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)	16 (1.580)	12/298	0.0	100%	coding (1350/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

GCGTTCACATCCCAGTCGCGCTGCGGCATTAATGCCAGACGCGCACCGTTGGCACCACCGCGTTTGTCGCCACCACGGAAGGTAGAAGCAGATGCCCAGGCCACCGATACCAGCTCACTAACAGACAGACCAGAATCCGCAATCGCGAATTTCAGATCGATAATGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4135345‑4135180 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTA  >  NC_000913/4498181‑4498422                                                                                                                                                                                                                                                                                                                                                                                                                          GCGTTCACATCCCAGTCGCGCTGCGGCATTAATGCCAGACGCGCACCGTTGGCACCACCGCGTTTGTCTCCACCACGTAAGGTAGAAGCAGATGCCCAGGCCACCGATACCAGCTCACTAACAGACAGACCAGAATCCGCAATCGCGAATTTCAGATCTATAATGTTAGACTGGCCCCCTTAATTTCCAGACATCCAAAATCACTTAAATAAGTGATAGTCTTAAAACTAGTCTTTAGACTTGTCATTGGA                                                                                                                                                               <  1:201632/251‑1                               AATGCCAGACGCGCACCGGTGGCACCACCGCGTTTGTCGCCACCACTGAAGGTAGAAGCAGATGCCCAGGCCACCGATACCAGCTCACTAACAAACAGACCAGAATCCGCAATCGCGAATGTCAGATCGATAATGTTAGACTGGCCCCCTGAATCTCCAGACAACCAAGATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTATTCATTGGAGAACAGATGATTGATGTCTTA                                                                                                                                          <  1:251782/242‑1                                                                    CGCCACCACGGAAGGTAGAAGCAGATGCCCAGGCCACCGATACCAGCTCACTAACAGACAGACCAGAATCCGCAATCGCGAATTTCAGATCGATAATGTTAGACTTGCCCCCTGAATCTCCAGACAACCAATATCACTGAAATAAGTGATATTCTTAGTACGAGTTTTg                                                                                                                                                                              <  1:72979/169‑2                                                                     GCCACCACGGAAGGTAGAAGCAGATGCCCAGGCCACAGATACCAGCTCACTAACAGACAGACCAGAATACGCAATCGCGAATTTCAGATCGATAATGTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAAT                                                                                           >  1:162991/1‑251                                                                                       AGCAGATGCCCAGGCCACCGATACCAGCTCACTAACAGACAGACCAGAATCCGCAATCGCGAATTTCCGATCGATAATGTTAGACTGGCCC                                                                                                                                                                                                                                         >  1:250338/1‑91                                                                                                                                    AGAATCCGCAATCGCGAATTTCAGATCGATAATGTTAGACTGGCCCCCTGAATCTCAAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGC                            >  1:146549/1‑251                                                                                                                                         nCGCAATCGCGAATTTCAGATCGATAATGTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGCGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGG                                  >  1:47638/2‑240                                                                                                                                            CAATCGCGAATTTCAGATCGATAATGTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCACAATTGTTCAGCAGAGCTTTGAACCGGGGATGAC                                                             >  1:91143/1‑210                                                                                                                                                               GATAATGTTAGACTGGCCCC‑TGACTCTCCAGACAACCACTATCACTTAAATAAGTGATCGTATTAATACTCGTTATTCGACTAGTCATTGGGGAACAGATGATTGTTGTCTTAAGGCCGGAGACACGAAGACGGCGTACCACACAGGAAACGATCGCACTTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTA  >  1:218668/1‑249                                                                                                                                                                                                                                                                                                                                                                                                                          GCGTTCACATCCCAGTCGCGCTGCGGCATTAATGCCAGACGCGCACCGTTGGCACCACCGCGTTTGTCGCCACCACGGAAGGTAGAAGCAGATGCCCAGGCCACCGATACCAGCTCACTAACAGACAGACCAGAATCCGCAATCGCGAATTTCAGATCGATAATGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4135345‑4135180 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTA  >  NC_000913/4498181‑4498422

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