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

AATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATT‑CATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCC  >  NC_000913/4497939‑4498404                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                aaTCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTGTTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATT‑CATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGc                                                                                                                                                                                                                         <  1:256989/251‑1 (MQ=255)                                 tAGGAGTATTTTCTTTGTTAAAAAAGACAAACGGTAGTAACTGATGTATTTCATATTAACCCCTCACTGCAGGGGGGCTTATCAGAATACGCGGTAATTACATCATT‑ATATACATTCAATAAAAAAGCACCATCATCAAAAATCATGTTCTCACAGGGTGCCGCGGGATAACTGACGCAAGAAGCGCTGCAGGGACTTGATGGACCCTTGTAGACTGGACCACTGCATCTCCAGACAACCAATCTCaatt                                                                                                                                                                                          >  1:328380/1‑247 (MQ=255)                                                                cGGTATTAACTGATTTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGTTTATT‑CATTATTTATAGACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCATGGATTTGATGGTCCCTTTTAGACTGGCCCTCTGAATCTCCAGACAACCAGTATCACTTAAATAAGTGATAGTCTTAATACTATTTTTTAg                                                                                                                                                          <  1:110330/251‑1 (MQ=39)                                                                            aTGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATT‑CATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                                                                                              <  1:71873/251‑1 (MQ=255)                                                                              gTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATT‑CATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag                                                                                                                                            <  1:100412/251‑1 (MQ=255)                                                                                                                              gTTATT‑CATTATTTATATACATTCAATCAAAAAGTACCCTCATAAAAAATCCTGTACTCCCAGGGTGATGCGGGATAACTGACGTAAGGAGAGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGTCAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATg                                                                                                                             >  1:266043/1‑218 (MQ=255)                                                                                                                                                            aaaaGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTcagcag                                                               <  1:305005/251‑1 (MQ=255)                                                                                                                                                                    ccTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTg                                                        >  1:186065/1‑250 (MQ=255)                                                                                                                                                                                       cTNNCAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                                                                                                                                                   <  2:4451/92‑1 (MQ=255)                                                                                                                                                                                       cTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                                                                                                                                                   >  1:4451/1‑92 (MQ=255)                                                                                                                                                                                       cTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTcc                                    >  1:21919/1‑251 (MQ=255)                                                                                                                                                                                                 gATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGccc                          >  1:7249/1‑251 (MQ=255)                                                                                                                                                                                                        gATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTCAATCAGCGATAGTCTTAAt                                                                                                                                                                      >  1:172645/1‑104 (MQ=25)                                                                                                                                                                                                               gACGTAAGGAGTGCTGCAGGGATTTGATGTTCCCTTGTAGACTTGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCGTAATACTAGTTTTGAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCTGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACGGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGgtgn             <  1:40976/250‑2 (MQ=14)                                                                                                                                                                                                                      ggAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTCGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGc                                                            >  1:293341/1‑196 (MQ=255)                                                                                                                                                                                                                         gTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCagccagcc  >  1:276340/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTTTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTcn                                <  1:79157/199‑2 (MQ=17)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                AATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATT‑CATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCC  >  NC_000913/4497939‑4498404

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.