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

GAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCG  >  NC_000913/4499276‑4499699
                                                                                                                                                                                                                                           |                                                                                                                                                                                            
nagagTAACGGAATAGCAGAGAGCTTCGTGAAAACGATACAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCACAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATg                                                                                                                                                                               >  1:58150/2‑251 (MQ=16)
                                     tAAAGCGTGACTACATCAGTATCATGCCCAGACCAGACGGGGTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACa                                                                                                                                           >  1:281424/1‑250 (MQ=255)
                                                    tcaggatcatgacacaaacAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGaa                                                                                                                                                                             >  1:258248/19‑201 (MQ=11)
                                                                                            gAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCAGAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGAGAACAGATGTCTGGAAATGTAGGGGCAAAGCCCCTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTGTCAGTACATGCGACTGTATTACAGGTAGCCATATATGGTTACCTGTGATCAGGAAGAATAAGtn                                                                                   <  1:64088/251‑2 (MQ=18)
                                                                                                       agGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGtttn                                                                        <  1:79107/251‑2 (MQ=255)
                                                                                                                 gCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTg                                                              <  1:135600/251‑1 (MQ=255)
                                                                                                                         aCGAATGGCATCCGCATAGTGCGCTGGGGTATCGCTCGCCACGGGAATATCTGCGGCAGCTGGCTTGTAATGGGTTAAGTGATATCAGATGTCTGGAAATATAGGGGCAAATTCACTTTTTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                                                                                        <  1:299455/169‑1 (MQ=18)
                                                                                                                                                       aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAgg                         <  1:286451/250‑1 (MQ=255)
                                                                                                                                                                           cTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                                                                                         >  1:323557/1‑118 (MQ=255)
                                                                                                                                                                             gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCg  <  1:162374/251‑1 (MQ=255)
                                                                                                                                                                                                             ananNTNNCTGNAAATATAGGGGCAAATCCACTNNGTGCTGATANGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGNCTGCCATATNNGANTACCTNTNNNCNNNAANAATAAGTGGCTGAGGTTTCNAANNGTCTGCNGGTATATTCCTCTCTCCCTTAAAAAc                                  <  2:71578/187‑1 (MQ=255)
                                                                                                                                                                                                             aCAGATGTCTGGAAATATAGGGGCAAATCCACTTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGACGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAAc                                  >  1:71578/1‑187 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                            
GAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACTGTCTGCCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCG  >  NC_000913/4499276‑4499699

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