New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4135180 =	3 (0.370)	11 (1.340)	6/334	0.2	78.6%	coding (1346/2181 nt)	katG	catalase/hydroperoxidase HPI
	?	NC_000913	4498181 =	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCT  >  NC_000913/4134940‑4135404                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  gTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGCTTTGATCCTGCGTTCGAGAAGATCTCTAGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGAAGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGACGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACAtttggatt                                                                                                                                                                                                                         >  1:107494/1‑245 (MQ=255)     aCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCACTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAACCTTTGCACGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCAGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACAtttggatttgcc                                                                                                                                                                                                                     >  1:155949/1‑241 (MQ=255)          cAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGCGAAGATCTATCGTAGTTTCCTCAACGTTACGCGGGCGTACAACGAAGCCTTTGCCCGTGCCTGGGTCACACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAACGAAGATCTGATCTGGCAAGATACGCTGACGCAGCCGCTCTACAACCCGACCGAGCCGGCCAtttgtctttgaccctat                                                                                                                                                                                                                >  1:114064/1‑236 (MQ=255)             aactggctggctgctacaccatgttgccgggcaacgagggacaccgtcatcaccggttccagactctgctgcacacttgcgatcttttcctgtgtggtccgccgtctgcgtttctacggccctaagacatccataatctgttctccactgactagtctcaaaactagtattaagactatcacttatttaagagatcttggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTcgcg                                                                                                                                                                                                             <  1:67552‑M1/23‑1 (MQ=255)                         gctactccatgttgccgggcaacgagggagaacgtcataccaggttcaaagctctgctgaacacttgctatcggttcctgtgtggtacgccgcctgcgtatctccgggacaaaggcatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttcagtgatattggttgtctggaaattcagggggccagtctaACATAATCGATCTGAAATTCGCGATTGCGGATTc                                                                                                                                                                                                  >  1:3922‑M1/217‑250 (MQ=255)                                                                        aaagctctgctgaacaagtgcgatcgtgtccggtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtagtaagactatctcttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCTATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGcttn                                                                                                                                                 >  1:14843‑M1/170‑250 (MQ=255)                                                                                      caattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCgtggtg                                                                                                                                   >  1:35853‑M1/156‑251 (MQ=255)                                                                                                                    ntctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTgcgc                                                                                                     <  1:102405‑M1/126‑1 (MQ=255)                                                                                                                              ctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtcgaACATTATCGATCTGAAATTCGCGATTGCGGATTCTggtctn                                                                                                                                                                                           >  1:30585‑M1/116‑155 (MQ=255)                                                                                                                                  ggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctgtagagtcagggggccagtctaACATTAGCGATCGGAGATTCGCGATTTCGGATGCTGGTCTGTC‑GTTGAGTGAGCTGGTTGCGGTGGCCTGGGCATCTGGTTCTACCTTCCGTGGTGGCGACAAACGCGGCGGTGCCAACGGTGCGCGTCTGGCATTAATg                                                                                       >  1:93900‑M1/112‑251 (MQ=255)                                                                                                                                                                              aactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTTTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGctct                                           <  1:166886‑M1/184‑1 (MQ=255)                                                                                                                                                                                                                    nttgtctggagattcagggggccagtataACAGTATCGCTCTGATATTAGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGCGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAgaaagagtatggcaa       <  1:49108‑M1/220‑8 (MQ=255)                                                                                                                                                                                                                      tgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGAGTCTAGTCTGTCTGTT‑AGTGAGCTGGTATCGGGGGCCTGGGCATCTGCTTCTACCTTCCGGGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGTGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGc    >  1:71159‑M1/28‑250 (MQ=255)                                                                                                                                                                                                                       ntctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGACTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCt  <  1:31022‑M1/225‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  GTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTT‑AGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCT  >  NC_000913/4134940‑4135404

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