Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
MC JC | NC_000913 | 3,090,347 | Δ708 bp | coding (1‑708/708 nt) | endA → | DNA‑specific endonuclease I |
Missing coverage evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | start | end | size | ←reads | reads→ | gene | description | |||
* | * | ÷ | NC_000913 | 3090347 | 3091054 | 708 | 26 [0] | [0] 25 | endA | DNA‑specific endonuclease I |
New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 3090346 | 0 (0.000) | 24 (0.790) | 22/98 | 0.1 | 100% | intergenic (+94/‑1) | yggI/endA | conserved protein YggI/DNA‑specific endonuclease I |
? | NC_000913 | 3091055 = | 0 (0.000) | intergenic (+1/‑79) | endA/rsmE | DNA‑specific endonuclease I/16S rRNA m(3)U1498 methyltransferase |
AATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACC > NC_000913/3091007‑3091104 | caaaacagctttcgctacgttgctggctcgttttaacacggagaagtgcc > 1:1385854‑M2/49‑50 (MQ=255) aaacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa < 1:2650877‑M2/4‑1 (MQ=255) aacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa > 1:1519790‑M2/46‑49 (MQ=255) cagctttcgctacgttgctggctcgttttaacacggagaagtgCCTacac > 1:1937881‑M2/44‑50 (MQ=255) agctttcgctacgttgctggctcgttttaacacggagaagtgCCTACACt < 1:220887‑M2/8‑1 (MQ=255) ctttcgctacgttgctggctcgttttaacacggagaagtgCCTACACTAg < 1:58697‑M2/10‑1 (MQ=255) cgctacgttgctggctcgttttaacacggagaagtgCCTACACTAGCgg > 1:877734‑M2/37‑49 (MQ=255) ctacgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGAt < 1:198590‑M2/16‑1 (MQ=255) cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCt > 1:1416272‑M2/32‑50 (MQ=255) gttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCtt > 1:1426068‑M2/31‑50 (MQ=255) ttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCttt < 1:2061011‑M2/21‑1 (MQ=255) ttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCttt > 1:2521843‑M2/30‑50 (MQ=255) gctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCtttt < 1:907123‑M2/22‑1 (MQ=255) gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGtt > 1:2099428‑M2/24‑49 (MQ=255) gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGtt < 1:1691872‑M2/26‑1 (MQ=255) gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGt < 1:1859657‑M2/25‑1 (MQ=255) gttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAccc > 1:517117‑M2/20‑50 (MQ=255) ttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAaccct > 1:760155‑M2/19‑49 (MQ=255) ttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc < 1:239187‑M2/32‑1 (MQ=255) tttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc > 1:2149456‑M2/18‑49 (MQ=255) acacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTAccc < 1:2187623‑M2/37‑1 (MQ=255) agaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAcgcg > 1:871117‑M2/8‑50 (MQ=255) aagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTa < 1:321280‑M2/45‑1 (MQ=255) gCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAAc > 1:2468744‑M2/2‑50 (MQ=255) ccTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAAcc > 1:1375104/1‑50 (MQ=255) | AATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACC > NC_000913/3091007‑3091104 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ 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. |