breseq version 0.35.4 revision f352f80f4bc9
mutation predictions | marginal predictions | summary statistics | genome diff | command line log |
Predicted mutations | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NZ_CP022686 | 457,894 | Δ1 bp | 100% | pseudogene (288/997 nt) | CIW80_RS02445 ← | DUF968 domain‑containing protein |
RA | NZ_CP022686 | 560,880 | C→T | 100% | intergenic (+62/+53) | CIW80_RS03085 → / ← pilV | hypothetical protein/shufflon system plasmid conjugative transfer pilus tip adhesin PilV |
RA | NZ_CP022686 | 601,295 | T→G | 23.4% | V1607G (GTC→GGC) | irp1 → | yersiniabactin polyketide synthase HMWP1 |
MC JC | NZ_CP022686 | 628,478 | Δ1,286 bp | 100% | [clbP] | [clbP] | |
RA | NZ_CP022686 | 756,115 | A→G | 100% | G181G (GGT→GGC) | galE ← | UDP‑glucose 4‑epimerase GalE |
RA | NZ_CP022686 | 794,460 | C→T | 100% | A568V (GCG→GTG) | dgcE → | diguanylate cyclase |
RA | NZ_CP022686 | 829,452 | G→T | 21.6% | intergenic (‑41/+186) | gatY ← / ← fbaB | tagatose‑bisphosphate aldolase subunit GatY/class I fructose‑bisphosphate aldolase |
RA | NZ_CP022686 | 917,427:1 | +TT | 100% | intergenic (+130/‑91) | yeiP → / → yeiQ | elongation factor P‑like protein YeiP/mannitol dehydrogenase family protein |
RA | NZ_CP022686 | 1,000,401 | A→C | 19.9% | V955G (GTG→GGG) | yfaL ← | AIDA‑I family autotransporter adhesin YfaL/EhaC |
RA | NZ_CP022686 | 1,005,085 | C→A | 13.2% | T380T (ACC→ACA) | nrdA → | ribonucleoside‑diphosphate reductase subunit alpha |
RA | NZ_CP022686 | 1,174,272 | A→C | 15.5% | intergenic (+1/‑43) | CIW80_RS05855 → / → CIW80_RS05860 | tRNA‑Val/tRNA‑Val |
RA | NZ_CP022686 | 1,174,279 | A→T | 13.0% | intergenic (+8/‑36) | CIW80_RS05855 → / → CIW80_RS05860 | tRNA‑Val/tRNA‑Val |
RA | NZ_CP022686 | 1,360,724 | G→T | 13.9% | intergenic (‑59/+265) | kgtP ← / ← rrf | alpha‑ketoglutarate permease/5S ribosomal RNA |
RA | NZ_CP022686 | 1,360,750 | G→T | 13.9% | intergenic (‑85/+239) | kgtP ← / ← rrf | alpha‑ketoglutarate permease/5S ribosomal RNA |
RA | NZ_CP022686 | 1,504,361 | G→T | 19.6% | V266V (GTG→GTT) | CIW80_RS07520 → | YgcG family protein |
RA | NZ_CP022686 | 1,674,244 | T→G | 17.4% | T182P (ACC→CCC) | ygfB ← | UPF0149 family protein YgfB |
RA | NZ_CP022686 | 1,954,923 | T→G | 15.3% | T92P (ACG→CCG) | ttdR ← | DNA‑binding transcriptional activator TtdR |
RA | NZ_CP022686 | 2,049,089 | T→G | 16.9% | N256T (AAC→ACC) | deaD ← | ATP‑dependent RNA helicase DeaD |
RA | NZ_CP022686 | 2,049,094 | G→T | 16.9% | T254T (ACC→ACA) | deaD ← | ATP‑dependent RNA helicase DeaD |
RA | NZ_CP022686 | 2,179,878 | T→G | 12.4% | K35T (AAG→ACG) | rplQ ← | 50S ribosomal protein L17 |
RA | NZ_CP022686 | 2,185,075 | A→C | 10.2% | V46G (GTT→GGT) | rpsE ← | 30S ribosomal protein S5 |
RA | NZ_CP022686 | 2,214,419 | T→C | 100% | K43R (AAA→AGA) | rpsL ← | 30S ribosomal protein S12 |
RA | NZ_CP022686 | 2,320,775 | C→T | 100% | G274S (GGC→AGC) | CIW80_RS11675 ← | DNA‑processing protein DprA |
RA | NZ_CP022686 | 2,394,287 | T→G | 16.9% | N115T (AAC→ACC) | chuS ← | hematinate‑forming heme oxygenase ChuS |
RA | NZ_CP022686 | 2,408,861 | Δ1 bp | 100% | pseudogene (1097/3115 nt) | mdtF → | multidrug efflux pump RND permease MdtF |
JC | NZ_CP022686 | 2,590,563 | +TTGG | 100% | coding (109/153 nt) | CIW80_RS26455 → | hypothetical protein |
RA | NZ_CP022686 | 2,736,578 | C→A | 8.9% | T180T (ACC→ACA) | rffA → | dTDP‑4‑amino‑4,6‑dideoxygalactose transaminase |
RA | NZ_CP022686 | 2,808,135 | C→A | 19.3% | P444T (CCT→ACT) | ubiB → | ubiquinone biosynthesis regulatory protein kinase UbiB |
RA | NZ_CP022686 | 2,868,740 | T→G | 25.9% | Q22P (CAG→CCG) | fdoH ← | formate dehydrogenase O subunit beta |
JC | NZ_CP022686 | 3,196,873 | +GCGTAATCGGTT | 100% | coding (944/1122 nt) | ampC ← | BlaEC family class C beta‑lactamase |
RA | NZ_CP022686 | 3,348,357 | T→G | 10.0% | intergenic (‑173/‑27) | CIW80_RS16685 ← / → CIW80_RS26495 | inovirus Gp2 family protein/hypothetical protein |
RA | NZ_CP022686 | 3,424,795 | A→C | 100% | intergenic (‑41/+27) | yjiJ ← / ← yjiK | MFS transporter/YjiK family protein |
RA | NZ_CP022686 | 3,451,834 | C→A | 12.9% | H275N (CAT→AAT) | lgoD → | L‑galactonate oxidoreductase |
RA | NZ_CP022686 | 3,516,093 | T→G | 8.9% | V219G (GTG→GGG) | nhaA → | Na+/H+ antiporter NhaA |
MC JC | NZ_CP022686 | 3,559,837 | Δ2,343 bp | 100% | coding (10‑2352/2352 nt) | polB ← | DNA polymerase II |
RA | NZ_CP022686 | 3,656,188 | A→G | 100% | S549P (TCC→CCC) | htrE ← | outer membrane usher protein |
RA | NZ_CP022686 | 3,770,509 | A→T | 10.8% | intergenic (‑245/+218) | CIW80_RS26265 ← / ← CIW80_RS18890 | ATP‑binding protein/hypothetical protein |
RA | NZ_CP022686 | 3,792,805 | A→G | 100% | Y421C (TAC→TGC) | fucP → | L‑fucose:H+ symporter permease |
RA | NZ_CP022686 | 3,797,481 | T→G | 13.9% | L130V (TTA→GTA) | CIW80_RS19065 → | LacI family DNA‑binding transcriptional regulator |
MC JC | NZ_CP022686 | 3,844,426 | Δ1,036 bp | 100% | coding (2‑1037/1056 nt) | dinB → | DNA polymerase IV |
RA | NZ_CP022686 | 3,868,846 | T→G | 11.1% | T383P (ACG→CCG) | ecpC ← | fimbrial usher EcpC |
RA | NZ_CP022686 | 3,935,324 | T→G | 15.3% | N705T (AAC→ACC) | lacZ ← | beta‑galactosidase |
RA | NZ_CP022686 | 4,270,807 | C→A | 9.7% | R484R (CGC→CGA) | sucA → | 2‑oxoglutarate dehydrogenase E1 component |
RA | NZ_CP022686 | 4,368,341 | T→G | 17.5% | D47A (GAC→GCC) | ybiY ← | glycyl‑radical enzyme activating protein |
RA | NZ_CP022686 | 4,368,372 | T→G | 16.7% | S37R (AGC→CGC) | ybiY ← | glycyl‑radical enzyme activating protein |
RA | NZ_CP022686 | 4,461,732 | T→G | 10.2% | intergenic (+86/‑84) | aroA → / → ycaL | 3‑phosphoshikimate 1‑carboxyvinyltransferase/metallopeptidase YcaL |
RA | NZ_CP022686 | 4,678,047 | A→C | 31.1% | intergenic (+45/‑76) | agn43 → / → CIW80_RS23450 | autotransporter adhesin Ag43/hypothetical protein |
MC JC | NZ_CP022686 | 4,844,868 | Δ1,662 bp | 100% | [umuD]–[umuC] | [umuD], [umuC] | |
RA | NZ_CP022686 | 4,870,912 | Δ1 bp | 100% | intergenic (‑108/‑120) | dhaK ← / → dhaR | dihydroxyacetone kinase subunit DhaK/dihydroxyacetone kinase operon transcriptional regulator DhaR |
RA | pMUT1‑hok‑empty | 1,314:1 | +A | 100% | intergenic (–/+403) | – / ← putative HTH domain protein | –/– |
RA | pMUT1‑hok‑empty | 1,376 | Δ1 bp | 100% | intergenic (–/+341) | – / ← putative HTH domain protein | –/– |
RA | pMUT1‑hok‑empty | 4,123:1 | +T | 100% | intergenic (‑341/–) | putative rom (RNA I inhibition modulator protein) ← / – | –/– |
RA | pMUT1‑hok‑empty | 4,127 | C→A | 100% | intergenic (‑345/–) | putative rom (RNA I inhibition modulator protein) ← / – | –/– |
RA | pMUT1‑hok‑empty | 4,132 | 4 bp→TA | 100% | intergenic (‑350/–) | putative rom (RNA I inhibition modulator protein) ← / – | –/– |
Unassigned missing coverage evidence | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | start | end | size | ←reads | reads→ | gene | description | |||
* | * | ÷ | NZ_MW240712 | 337–2618 | 3173 | 556–2837 | 4187 [4166] | [367] NA | [JNG15_RS00005]–mobA | [JNG15_RS00005],JNG15_RS00010,JNG15_RS00015,mobA |
* | * | ÷ | pMUT1‑hok‑empty | 1453–1574 | 3986–1722 | 149–2534 | 2180 [2044] | [2100] 2150 | putative HTH domain protein–putative rom (RNA I inhibition modulator protein) | putative HTH domain protein,uncharacterized protein,hypothetical protein,putative rom (RNA I inhibition modulator protein) |
Unassigned new junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NZ_CP022686 | 262466 = | 0 (0.000) | 108 (0.990) | 51/250 | NT | 100% | coding (320/1647 nt) | fadK | medium‑chain fatty‑acid‑‑CoA ligase |
? | NZ_CP022686 | = 1524130 | NA (NA) | intergenic (+8/+232) | CIW80_RS07590/syd | IS3‑like element ISEc52 family transposase/SecY‑interacting protein | |||||
* | ? | NZ_CP022686 | = 262468 | 0 (0.000) | 114 (1.050) | 51/250 | NT | 100% | coding (322/1647 nt) | fadK | medium‑chain fatty‑acid‑‑CoA ligase |
? | NZ_CP022686 | = 1761697 | NA (NA) | intergenic (‑72/‑70) | CIW80_RS08710/CIW80_RS08720 | IS3‑like element ISEc52 family transposase/IS3‑like element IS600 family transposase | |||||
* | ? | pMUT1‑hok‑empty | 456 = | 1675 (0.470) | 2341 (0.660) | 184/250 | NT | 66.4% | intergenic (–/+1261) | –/putative HTH domain protein | –/– |
? | pMUT1‑hok‑empty | 4136 = | 691 (0.200) | intergenic (‑354/–) | putative rom (RNA I inhibition modulator protein)/– | –/– | |||||
* | ? | pMUT1‑hok‑empty | = 491 | 2418 (0.680) | 1724 (0.610) | 144/200 | NT | 56.8% | intergenic (–/+1226) | –/putative HTH domain protein | –/– |
? | pMUT1‑hok‑empty | = 4098 | 691 (0.240) | intergenic (‑316/–) | putative rom (RNA I inhibition modulator protein)/– | –/– |