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)11 (1.100) 9/322 0.1 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)16 (1.600) 13/322 0.0 100% coding (1350/2181 nt) katG catalase/hydroperoxidase HPI
?NC_000913 = 4499511 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

AACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCT  >  NC_000913/4497933‑4498368
                                                                                                                                                                                                                                                        |                                                                                                                                                                                           
aaCAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAg                                                                                                                                                                                           <  1:122672/251‑1 (MQ=255)
                 gCAAATGCCCGGGAATTATTCTAGGCTTATTTTCTTTGTTACAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAAtct                                                                                                                                                                          >  1:93738/1‑251 (MQ=255)
                                                                                                      cccTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                                                                     >  1:189328/1‑251 (MQ=255)
                                                                                                        cTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                                                                                                                              >  1:12090/1‑192 (MQ=255)
                                                                                                                                                    taCATTCAATAAAAAAGTAACCTCATAAAAAATCCTTTACTCACTGGGTGATGCGGGATAACTTACGTTGGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGGCTGGCCCCCTGAAtc                                                                                                                                                                           <  1:160047/119‑1 (MQ=255)
                                                                                                                                                    taCATNNAANNNAANNGTAACCTCATAAAAAATCCTGTANTCACAGGGTGATGCGGGATAACTGACGTANGNANTGCTNNAGGNATTTGATGGTNCCTTGTAGACTGGCCCCCTGAAtc                                                                                                                                                                           >  2:160047/1‑119 (MQ=255)
                                                                                                                                                                      aaCCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATNGTCNNNNNAc                                                                                                                                >  2:128339/1‑144 (MQ=255)
                                                                                                                                                                                        gTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                                                                                                          <  1:100855/148‑1 (MQ=255)
                                                                                                                                                                                         tACTCACAGGGAGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCATTGTAGACTGGCCCCCTGACTCTCCAGACAACCAATCTCACTTAAATAAGTGATAGTCTTCATCCTAGTTTTTAGAATAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGtct  >  1:215508/1‑251 (MQ=17)
                                                                                                                                                                                                     gNNGCGGGANNACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                                                                                                                              <  2:12090/99‑1 (MQ=255)
                                                                                                                                                                                                                 cTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGNTTGNNNNNtt                                                                                     >  2:136676/1‑144 (MQ=255)
                                                                                                                                                                                                                                        aTTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAgg                                                  <  1:168208/156‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                           
AACAGTAATCATCCCGGGCAAATGCCCGGGAATTATTCTAGGATTATTTTCTTTGTTAAAAAAGACAAACGGTATTAACTGATGTATTTACTATTTACCGCTCCCTGCAGGGGGGATTATCAGAATCCGCGGTTATTCATTATTTATATACATTCAATAAAAAAGTAACCTCATAAAAAATCCTGTACTCACAGGGTGATGCGGGATAACTGACGTAAGGAGTGCTGCAGGGATTTGATGGTCCCTTGTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCT  >  NC_000913/4497933‑4498368

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