| 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 | 15 [0] | [0] 15 | 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) | 14 (0.560) | 12/116 | 0.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:881594‑M1/68‑2 (MQ=255) cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTTCTTAACATAAATGtataag < 2:776194‑M1/68‑7 (MQ=255) aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg < 2:271234‑M1/68‑18 (MQ=255) ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt < 1:289577‑M1/68‑19 (MQ=255) gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt < 1:633230‑M1/68‑26 (MQ=255) gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt < 2:464411‑M1/68‑26 (MQ=255) tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat < 1:732062‑M1/68‑30 (MQ=255) tCACACAAAACACTCAATTACTTAACATAAATGtataagtcat > 1:710566‑M1/1‑33 (MQ=255) tCACACAAAACACTCAATTACTTAACATAAATGtataagtcat < 2:710571‑M1/43‑11 (MQ=255) acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt > 1:26350‑M1/1‑29 (MQ=255) aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc > 2:326648‑M1/1‑14 (MQ=255) aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc > 2:637327‑M1/1‑14 (MQ=255) ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta > 1:694682‑M1/1‑12 (MQ=255) tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac < 2:495799‑M1/68‑58 (MQ=255) tGtattagtcatacttttgatttgggtgaatttcaaactcataaaaagtttttcgctacgctagcaag > 1:959847‑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 < 8 ≤ 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. |