Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
MC JC | NC_000913 | 3,453,928 | Δ15,926 bp | gspA–[chiA] | 16 genes gspA, gspC, gspD, gspE, gspF, gspG, gspH, gspI, gspJ, gspK, gspL, gspM, gspO, bfr, bfd, [chiA] |
Missing coverage evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | start | end | size | ←reads | reads→ | gene | description | |||
* | * | ÷ | NC_000913 | 3453928 | 3469853 | 15926 | 42 [0] | [0] 41 | gspA–[chiA] | gspA,gspC,gspD,gspE,gspF,gspG,gspH,gspI,gspJ,gspK,gspL,gspM,gspO,bfr,bfd,[chiA] |
New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 3453927 | 0 (0.000) | 38 (1.330) | 32/92 | 0.0 | 100% | coding (1/420 nt) | gspB | putative general secretion pathway protein B |
? | NC_000913 | 3469854 = | 0 (0.000) | intergenic (‑1/+291) | chiA/tufA | endochitinase/translation elongation factor Tu 1 |
CACCAGCCCCATACCAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGAA > NC_000913/3469806‑3469903 | gtgacctgtttctttctgctcacgcgcggcgatatagaattcaaacatta > 1:348419‑M2/49‑50 (MQ=255) tgacctgtttctttctgctcacgcgcggcgatatagaattcaaacattat < 1:2037696‑M2/3‑1 (MQ=255) gacctgtttctttctgctcacgcgcggcgatatagaattcaaacattat > 1:359023‑M2/47‑49 (MQ=255) acctgtttctttctgctcacgcgcggcgatatagaattcaaacattata < 1:1035780‑M2/4‑1 (MQ=255) acctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATaa < 1:2470732‑M2/5‑1 (MQ=255) cctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATAAt < 1:224651‑M2/6‑1 (MQ=255) tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATcc < 1:1008142‑M2/8‑1 (MQ=255) tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCt < 1:1424369‑M2/10‑1 (MQ=255) tctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCtt < 1:192830‑M2/11‑1 (MQ=255) ctttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgt > 1:1136261‑M2/38‑50 (MQ=255) tttctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgt > 1:1810072‑M2/37‑49 (MQ=255) tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtgtg < 1:1455139‑M2/14‑1 (MQ=255) tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGAc > 1:1522825‑M2/35‑50 (MQ=255) tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGAc > 1:1169340‑M2/35‑50 (MQ=255) ctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACg > 1:2283186‑M2/34‑50 (MQ=255) tcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaa > 1:1621373‑M2/30‑49 (MQ=255) cacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaa > 1:425721‑M2/29‑49 (MQ=255) cacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaa > 1:91798‑M2/29‑49 (MQ=255) cgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTg > 1:1454945‑M2/25‑50 (MQ=255) gcggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGc < 1:2682069‑M2/27‑1 (MQ=255) cggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGc > 1:1513457‑M2/23‑49 (MQ=255) ggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaa < 1:274385‑M2/29‑1 (MQ=255) gatatcgaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaac < 1:1143559‑M2/32‑1 (MQ=255) gatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaa > 1:921993‑M2/19‑49 (MQ=255) atatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaaca < 1:462520‑M2/33‑1 (MQ=255) atagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa > 1:215447‑M2/16‑49 (MQ=255) atagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa > 1:2613327‑M2/16‑49 (MQ=255) tagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaaaa < 1:2166755‑M2/36‑1 (MQ=255) tagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaaa < 1:93488‑M2/35‑1 (MQ=255) agaattcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAg > 1:1049277‑M2/14‑50 (MQ=255) agaattcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAg > 1:391930‑M2/14‑50 (MQ=255) agaattcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAg > 1:1570638‑M2/14‑50 (MQ=255) aattcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTc > 1:1686225‑M2/12‑50 (MQ=255) ttcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCt < 1:2080550‑M2/41‑1 (MQ=255) ttcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCt < 1:111006‑M2/41‑1 (MQ=255) tcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTa > 1:243012‑M2/9‑50 (MQ=255) caaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc > 1:1320774‑M2/8‑50 (MQ=255) acatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCa < 1:2442914‑M2/46‑1 (MQ=255) atTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAg > 1:1353347‑M2/3‑49 (MQ=255) tTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAgg > 1:860493‑M2/2‑49 (MQ=255) tataATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGaa > 1:1143331/1‑50 (MQ=255) | CACCAGCCCCATACCAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAGGAA > NC_000913/3469806‑3469903 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ 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. |