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.070)17 (1.200) 14/310 0.1 94.3% coding (1346/2181 nt) katG catalase/hydroperoxidase HPI
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 41351841 (0.070)14 (0.990) 12/310 0.1 93.1% coding (1350/2181 nt) katG catalase/hydroperoxidase HPI
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

CCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTG  >  NC_000913/4134945‑4135427
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
gagaaataactggctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGa                                                                                                                                                                                                                                           <  1:244886‑M1/15‑1 (MQ=255)
           ggTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTAGTTTCATCAACGATCTGCAGGCGTTAAACGAAGCCTTTGCCCGTGACTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATATGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttcc                                                                                                                                                                                                                                >  1:229275‑M1/1‑229 (MQ=255)
            nctggctgctacaccatgttgccgggcaacgagggagaccgtcatccccggttcaaagctctgctgaaccattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactcgtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGc                                                                                                                                                                                                                              <  1:95121‑M1/28‑1 (MQ=255)
                              cGTTTTGCTCCTGAGTTAGAGAAGATCTCTCGTCGTTTCCTCACCGATCCGCAGGCGGTAAACGAAGCCTTTGCCCGTGCCTGGGTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGACGTGCCGAAAGAAGATCTGATCTGGCAAGATACGCTGCCTCAGCCGATCTACAACCCGCCCGAGCAGGAAATTtggatttgcccctatatttccatacatctgttatcactta                                                                                                                                                                                                             >  1:250121‑M1/1‑210 (MQ=255)
                                   gggcaacgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGAGCTGAAATGCGCGATTGCGGATTCTGGTCTGTCTTTTAGTg                                                                                                                                                                                                       >  1:111668‑M1/201‑251 (MQ=255)
                                         cgagggagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttctcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGAt                                                                                                                                                                                                                          <  1:186659‑M1/32‑1 (MQ=255)
                                          gAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaag                                                                                                                                                                                                 <  1:8908‑M1/250‑53 (MQ=255)
                                             ngagaccgtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaAc                                                                                                                                                                                                                                                        <  1:196853‑M1/2‑1 (MQ=255)
                                                    gtcatccccggttcaaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGcc                                                                                                                                                                                      <  1:16479‑M1/68‑1 (MQ=255)
                                                                 caaagctctgctgaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactcgtctaaaaactcgtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCtt                                                                                                                                                                         <  1:227007‑M1/81‑1 (MQ=255)
                                                                             gaacaattgcgatcttttcctgtgtggtacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACATTCCgtg                                                                                                                                                             <  1:102584‑M1/93‑1 (MQ=255)
                                                                                     cGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTTCATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGTCCGAGCAGGACATTATCGATCTGTAAGTCGCGATTGCGGAGTCTGGTCTGTCTGTTAGTGAGCTGGTAGCGGTGGCCTGGGCATCTGCTTCTAcc                                                                                                                                                                    <  1:275707/236‑1 (MQ=255)
                                                                                       nnnnnnttnnnnngtgntacgccgtctgcgtttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctgcagattcagggggccagtctaAc                                                                                                                                                                                                                                                        >  2:196853‑M1/149‑150 (MQ=255)
                                                                                                nnnnnngcnnnnnCCTNGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatt                                                                                                                                                                                                                                               <  2:196023‑M1/144‑7 (MQ=255)
                                                                                                            gCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataan                                                                                                                                                             <  1:78177‑M1/220‑89 (MQ=255)
                                                                                                                      aaCTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCTGAAGTGTCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtgtatttttccctatatttccatgcagctgttatcacttaacccattacaatcccgctgccgcagatgttcccgtgtcgagcgatagcccatcgcactatgcggatgccattcgttataatgctcgaac                                                                                                                    <  1:158849‑M1/251‑130 (MQ=255)
                                                                                                                                             ccGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACACCCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggctgccattcgttataatgctcgaacgcctctgcaaggttctttgctgc                                                                                             >  1:243124‑M1/1‑99 (MQ=255)
                                                                                                                                                        caatgactagtctaaaaactagtattaagaatatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATAGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGCGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGc                                                                                                                                       <  1:120641‑M1/115‑1 (MQ=255)
                                                                                                                                                                  tctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAACATCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGACATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGCGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCg                                                                        <  1:312402‑M1/178‑1 (MQ=255)
                                                                                                                                                                   ctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGAGTCTGGTCTGTCTGGTAGTGAGCTGGTATCGGtgn                                                                                                                                                                                        >  1:95783‑M1/73‑137 (MQ=255)
                                                                                                                                                                          gCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctc                                                                                                               >  1:19064‑M1/1‑70 (MQ=255)
                                                                                                                                                                           ccGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTgggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcanat                                                                                                                                                                                  >  2:192368‑M1/1‑69 (MQ=255)
                                                                                                                                                                                     nnnnnnnnntnttnnagggntattggttgtctggagattcagggngccagtctaACATTATCGAGCTGAAATTCGCNANTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCNTCTGCTTCTANCTGNNGTGGNGNCNACAAACGCGGTGGGGCCAACGGTGCGCGTNTGGCATTAATGCCGCAGCGCGACNNGGNNNTGNNcgccg                                                                                  >  2:244886‑M1/55‑222 (MQ=255)
                                                                                                                                                                                         gatctgGCAAGATCCGCTGCCGCAGCCGATCTACAACCAGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacacattacaagcccgctgccgcagatattcacgtggcgagc                                                                                                                                                                   >  1:139790‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                     tCCGCTGCCGCATCCGATCTACAACCCGACCGAGCAGGACATTtggatttgaacctatatttcccgacatatgttatcacttaccccattacaagcccgctgccgcagatattcccgtgtcgagcgataacccagagcactatgc                                                                                                                                               >  1:222973‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                     tCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTtggatttgcccctntatntccanacntctgttatcacttaanccattacaagcccgctgccgcagatattcccgtggcgagcgatanncnnncncactatgc                                                                                                                                               <  2:222973‑M1/145‑103 (MQ=255)
                                                                                                                                                                                                      ccGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACCTTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgacattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcac                                    >  1:127227‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                                      gagattcagggggccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCCTTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTGCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTg                    <  1:23869‑M1/230‑1 (MQ=255)
                                                                                                                                                                                                                             cccGACCGAGCAGGACATTtggatttgcccctatatttccagacatctgttatcacttaacccattacaagcccgctgccgcagatattcccgtggcgagcgataacccagcgcactatgcggatgccattcgttataatgctcgaacgcctctgcaaggttctttgctgccgttaacccgtctggtttgggcatgatactgatgtagtcacgctttatcgttttcacgaagctctctgctattccgttan             <  1:195022‑M1/251‑233 (MQ=255)
                                                                                                                                                                                                                                  gccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCAt         <  1:8365‑M1/241‑1 (MQ=255)
                                                                                                                                                                                                                                   ccagtctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGCCAAACGCGGTGGTGCCAAc                                                                                                                                   <  1:321894‑M1/119‑1 (MQ=255)
                                                                                                                                                                                                                                        ctaACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAAAGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATTTCATAGTGCTg  <  1:46752‑M1/248‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
CCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGGCATCTGCTTCTACCTTCCGTGGTGGCGACAAACGCGGTGGTGCCAACGGTGCGCGTCTGGCATTAATGCCGCAGCGCGACTGGGATGTGAACGCCGCAGCCGTTCGTGCTCTGCCTGTTCTGGAGAAAATCCAGAAAGAGTCTGGTAAAGCCTCGCTGGCGGATATCATAGTGCTG  >  NC_000913/4134945‑4135427

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.