Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
MC JC | NC_000913 | 1,978,503 | Δ776 bp | insB1–insA | insB1, insA |
Missing coverage evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | start | end | size | ←reads | reads→ | gene | description | |||
* | * | ÷ | NC_000913 | 1978503 | 1979278 | 776 | 11 [0] | [0] 11 | insB1–insA | insB1,insA |
New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 1978502 | 0 (0.000) | 10 (0.320) | 7/114 | 1.5 | 100% | intergenic (‑305/+16) | flhD/insB1 | flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B |
? | NC_000913 | 1979279 = | 0 (0.000) | intergenic (‑64/‑474) | insA/uspC | IS1 repressor TnpA/universal stress protein |
GGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAG > NC_000913/1978436‑1978568 | gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt < 1:198132‑M1/68‑2 (MQ=255) tttCAATAATGCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact < 2:197638‑M1/68‑14 (MQ=255) tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact < 2:35692‑M1/68‑14 (MQ=255) gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggc < 2:197618‑M1/68‑26 (MQ=255) tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat > 2:521640‑M1/1‑33 (MQ=255) tCACACAAAACACTCAATTACTTAACAAAAATGtataagtcatacttttgttttgggtgtatttcaat > 2:947852‑M1/1‑33 (MQ=255) acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct > 2:1126902‑M1/1‑31 (MQ=255) cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt < 1:1152730‑M1/68‑48 (MQ=255) cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt > 1:1220591‑M1/1‑21 (MQ=255) ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc < 1:199056‑M1/68‑53 (MQ=255) tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag > 1:892153‑M1/1‑2 (MQ=255) | GGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAG > NC_000913/1978436‑1978568 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ 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. |