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)12 (1.250) 11/314 0.0 100% coding (1346/2181 nt) katG catalase/hydroperoxidase HPI
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 41351840 (0.000)10 (1.050) 8/314 0.1 100% coding (1350/2181 nt) katG catalase/hydroperoxidase HPI
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

TGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGA  >  NC_000913/4134952‑4135375
                                                                                                                                                                                                                                        |                                                                                                                                                                                               
aactggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttgtcctgggtggtacgccgtcggcggttctccggccctaagacatcaatcatctgttcgccaatgacgagtcgaaaaactagtattaagacgatcacttatgtaagtgatattggttgtctggagattcagggggccaggctaACATGAGCGATCTGAAATtcgg                                                                                                                                                                                >  1:166336‑M1/229‑249 (MQ=255)
                                            aagaTCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGTTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGATCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccg                                                                                                                                   <  1:35904‑M1/251‑63 (MQ=255)
                                                       gttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttatgtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCt                                                                                                                         >  1:160911‑M1/174‑250 (MQ=255)
                                                               ctctgctgaaccattgcgatcttttcatgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGc                                                                                                                             <  1:172498‑M1/73‑1 (MQ=255)
                                                                        acaattgcgatcttttcctgtgtggtacgccgtctgcgtttcgccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactagcacttatttaagtgatattggttgtctggagagtcaggtggccagtctaACATTATCGATCTAAAATTCGTGATTTCGCATTCTCGTCTGTCTGTTAGt                                                                                                                                                    >  1:182681‑M1/157‑206 (MQ=255)
                                                                                 atcttttcctgtgtgggacgccgtctgcgtttctccggccctaagacatcagtcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagggatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAAGTCGCGATTGCGGATTCTGGTCTGGCTGATAGTGAGCTGGTATCGGTGGCCTGGGc                                                                                                                             >  1:182054‑M1/148‑220 (MQ=255)
                                                                                           cTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTAGCTACATCGGGCCGGAAGAGCCGACAGCAGATCTGATCTGGCACGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgcc                                                                                    >  1:126387‑M1/1‑142 (MQ=255)
                                                                                                 ncgTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgctaacccagcgcactatgcggattccattcgt                                                                              >  1:56749‑M1/2‑136 (MQ=255)
                                                                                                                        acaGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgc                                                        >  1:11841‑M1/1‑113 (MQ=255)
                                                                                                                                      catctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTgcg                                                                         <  1:119851‑M1/125‑1 (MQ=255)
                                                                                                                                               cGCTACAGCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcactgaacccattacaagcccgctgccgcagatattcccgtggcgagcgagaacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacc                                 <  1:92552‑M1/250‑161 (MQ=255)
                                                                                                                                                ccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATACGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAAcgcc                               >  1:219022‑M1/85‑251 (MQ=255)
                                                                                                                                                      actagtctaaaaactagtattaagactatcacgtatgtaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGtctgtctg                                                                                                                                                         >  1:208305‑M1/79‑123 (MQ=255)
                                                                                                                                                           ccGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttggg                    <  1:80873‑M1/251‑174 (MQ=255)
                                                                                                                                                               aaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTGAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCGACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGct                >  1:177088‑M1/70‑251 (MQ=255)
                                                                                                                                                                  actagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACAGTATCGATCGGAAATTCGCGAGTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGGGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAgc                                                   >  1:34142‑M1/67‑213 (MQ=255)
                                                                                                                                                                             gatctgatctGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatcggttatcacttaacccattacaagcccgctgccgcagctattcccgtagcgagcgctaacccagcgcactatgcggatgccattcgttatactgctcgaacgcctctgcaagcttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtc  >  1:12243‑M1/1‑60 (MQ=255)
                                                                                                                                                                                  gatctgGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtc  <  1:6519‑M1/246‑192 (MQ=255)
                                                                                                                                                                                  gatctgGCAAGATCCGCTGCCGCAGCCGATCTACAACCAGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgtta                                                                                                                                                                >  1:101393‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                             gATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcn                                                             <  1:24582‑M1/160‑133 (MQ=255)
                                                                                                                                                                                                              ggagattcagggtgccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGGCAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTc      >  1:117565‑M1/23‑214 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                               
TGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGA  >  NC_000913/4134952‑4135375

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