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 annotation gene description
RA AE015451_phaM1 144,577 C→T intergenic (+556/‑190) PP_0136 CDS (conserved exported protein of unknown function) → / → gltP CDS (glutamate_aspartate‑proton DAACS transporter) –/–
RA AE015451_phaM1 278,744 +A intergenic (‑239/+75) PP_0224 CDS (Monooxygenase_ DszC family) ← / ← sctC CDS (sulfur compound ABC transporter ‑ ATP‑binding subunit) –/–
RA AE015451_phaM1 307,866 +C intergenic (+881/‑988) hslO CDS (33 kDa chaperonin) → / → PP_0254 CDS (conserved protein of unknown function) –/–
RA AE015451_phaM1 336,125 +T coding (210/237 nt) PP_0278 CDS (conserved protein of unknown function) →
RA AE015451_phaM1 353,067 A→G intergenic (+34/+84) hisF CDS (Imidazole glycerol phosphate synthase subunit HisF) → / ← cbcV CDS (choline _ betaine _ carnitine ABC transporter ‑ ATP binding subunit) –/–
MC JC AE015451_phaM1 386,865 Δ1,248 bp coding (5‑1252/1254 nt) glyA‑I CDS (serine hydroxymethyltransferase) →
RA AE015451_phaM1 499,203 A→G intergenic (+47/‑249) uhpA CDS (DNA‑binding response regulator UhpA) → / → PP_0411 CDS (Polyamine ABC transporter_ ATP‑binding protein) –/–
RA AE015451_phaM1 539,294 A→C T558P (ACC→CCC)  rpoB CDS (DNA‑directed RNA polymerase subunit beta) →
RA AE015451_phaM1 542,514 C→G L252V (CTG→GTG)  rpoC CDS (DNA‑directed RNA polymerase subunit beta') →
MC JC AE015451_phaM1 780,595 Δ1,248 bp coding (4‑1251/1254 nt) glyA‑II CDS (serine hydroxymethyltransferase) ←
RA AE015451_phaM1 1,070,248 +GA intergenic (+48/+94) aroP‑I CDS (aromatic amino acid transport protein) → / ← PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein) –/–
MC JC AE015451_phaM1 1,124,024 Δ5,919 bp [gcvT‑I CDS (aminomethyltransferase)][gcvH‑I CDS (glycine cleavage system H protein 1)] [gcvT‑I CDS (aminomethyltransferase)], tdcG‑II CDS (L‑serine dehydratase), gcvP‑I CDS (glycine dehydrogenase), [gcvH‑I CDS (glycine cleavage system H protein 1)]
RA AE015451_phaM1 1,419,549 +C intergenic (‑38/+98) PP_1242 CDS (conserved protein of unknown function) ← / ← PP_1244 CDS (conserved protein of unknown function) –/–
RA AE015451_phaM1 1,499,480 2 bp→AC intergenic (+57/‑106) trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain) –/–
RA AE015451_phaM1 1,499,498 1 bp→CG intergenic (+75/‑89) trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain) –/–
RA AE015451_phaM1 1,499,508 +C intergenic (+85/‑79) trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain) –/–
RA AE015451_phaM1 1,713,060 A→G T2A (ACC→GCC)  tsaB CDS (endopeptidase) →
MC JC AE015451_phaM1 1,738,086 Δ40,059 bp [PP_1532 CDS (putative Phage integrase)][PP_1586 CDS (putative Killer protein)] 60 genes
RA AE015451_phaM1 1,932,640 +C intergenic (‑107/+43) rluA CDS (23S rRNA (746) and tRNA (32) pseudouridine synthase) ← / ← minE CDS (Cell division topological specificity factor) –/–
RA AE015451_phaM1 2,586,629 Δ1 bp intergenic (+535/‑7) PP_t39 tRNA (tRNA‑Leu) → / → PP_2266 CDS (DNA‑directed RNA polymerase_ bacteriophage and mitochondrial type) –/–
MC JC AE015451_phaM1 2,586,636 Δ39,187 bp [PP_2266 CDS (DNA‑directed RNA polymerase_ bacteriophage and mitochondrial type)][PP_2297 CDS (putative integrative genetic element Ppu40_ integrase)] 42 genes
MC JC AE015451_phaM1 2,798,045 Δ693 bp coding (1‑693/693 nt) endX CDS (extracellular DNA endonuclease) ←
MC JC AE015451_phaM1 3,366,695 Δ15,777 bp [PP_2964 CDS (putative transposase)][PP_2984 CDS (conserved protein of unknown function)] 20 genes
MC JC AE015451_phaM1 3,412,992 Δ35,612 bp [PP_3026 CDS (putative Phage recombinase)]PP_5557 CDS (protein of unknown function) 47 genes
MC JC AE015451_phaM1 3,760,510 Δ22,673 bp PP_3325 CDS (putative Outer membrane ferric siderophore receptor)[nikD CDS (nickel ABC transporter ATP‑binding subunit)] 22 genes
MC JC AE015451_phaM1 3,819,910 Δ966 bp coding (1‑966/966 nt) endA CDS (endonuclease I) →
RA AE015451_phaM1 3,845,705 +G intergenic (‑72/+26) PP_3394 CDS (putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase) ← / ← PP_3395 CDS (Transcriptional regulator_ LysR family) –/–
RA AE015451_phaM1 3,919,448 G→C A14P (GCC→CCC)  mobA CDS (Molybdenum cofactor guanylyltransferase) →
RA AE015451_phaM1 3,951,159 T→C intergenic (‑123/+72) PP_3486 CDS (putative cytochrome c) ← / ← PP_3487 CDS (conserved exported protein of unknown function) –/–
RA AE015451_phaM1 3,951,161 A→T intergenic (‑125/+70) PP_3486 CDS (putative cytochrome c) ← / ← PP_3487 CDS (conserved exported protein of unknown function) –/–
MC JC AE015451_phaM1 4,372,651 Δ54,766 bp PP_3849 CDS (putative calcium‑binding protein_ hemolysin‑type)[PP_3920 CDS (putative Phage integrase)] 78 genes
RA AE015451_phaM1 4,586,030 2 bp→TC intergenic (+113/+101) PP_4061 CDS (conserved protein of unknown function) → / ← PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase) –/–
RA AE015451_phaM1 4,586,035 +G intergenic (+118/+97) PP_4061 CDS (conserved protein of unknown function) → / ← PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase) –/–
RA AE015451_phaM1 4,586,057 +C intergenic (+140/+75) PP_4061 CDS (conserved protein of unknown function) → / ← PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase) –/–
RA AE015451_phaM1 4,740,811 (T)7→6 intergenic (+48/+705) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,740,816 T→G intergenic (+53/+700) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,740,819 2 bp→CT intergenic (+56/+696) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,228 +C intergenic (+465/+288) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,230 A→G intergenic (+467/+286) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,232 C→T intergenic (+469/+284) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,235 C→T intergenic (+472/+281) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,238 C→G intergenic (+475/+278) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,240 2 bp→CT intergenic (+477/+275) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,246 Δ1 bp intergenic (+483/+270) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,254 +T intergenic (+491/+262) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,258 A→G intergenic (+495/+258) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,262 Δ1 bp intergenic (+499/+254) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,263 A→C intergenic (+500/+253) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,272 A→G intergenic (+509/+244) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,274 A→C intergenic (+511/+242) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,276 G→C intergenic (+513/+240) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,278 C→T intergenic (+515/+238) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,283 +T intergenic (+520/+233) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,289 C→G intergenic (+526/+227) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,292 C→T intergenic (+529/+224) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,294 C→G intergenic (+531/+222) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,296 A→C intergenic (+533/+220) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
RA AE015451_phaM1 4,741,298 C→T intergenic (+535/+218) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–
MC JC AE015451_phaM1 4,919,155 Δ69,176 bp [PP_4328 CDS (conserved protein of unknown function)][ycgR CDS (Flagellar brake protein YcgR)] 68 genes
MC JC AE015451_phaM1 5,393,059 Δ5,661 bp [hsdR CDS (Type I restriction enzyme EcoAI R protein)][hsdS CDS (Type I restriction‑modification system_ S subunit)] [hsdR CDS (Type I restriction enzyme EcoAI R protein)], hsdM CDS (Type I restriction enzyme EcoEI M protein), [hsdS CDS (Type I restriction‑modification system_ S subunit)]
RA AE015451_phaM1 5,555,608 +GCC intergenic (‑84/+87) PP_4887 CDS (conserved protein of unknown function) ← / ← PP_4888 CDS (putative methyl‑accepting chemotaxis transducer) –/–
RA AE015451_phaM1 5,659,798 G→A R281H (CGC→CAC)  metK CDS (methionine adenosyltransferase) →
RA AE015451_phaM1 5,674,751 +G intergenic (+113/+27) rhlE‑II CDS (ATP‑dependent RNA helicase) → / ← yceI CDS (polyisoprene‑binding protein_ acidic stress response factor) –/–
RA AE015451_phaM1 5,674,753 G→C intergenic (+115/+25) rhlE‑II CDS (ATP‑dependent RNA helicase) → / ← yceI CDS (polyisoprene‑binding protein_ acidic stress response factor) –/–
RA AE015451_phaM1 5,681,415 +CGGG intergenic (‑133/+71) PP_4986 CDS (Channel protein_ hemolysin III family) ← / ← PP_4987 CDS (putative Chemotaxis protein) –/–
RA AE015451_phaM1 5,701,209 C→A intergenic (+1955/‑2050) PP_5002 CDS (conserved protein of unknown function) → / → PP_5009 CDS (conserved protein of unknown function) –/–
MC JC AE015451_phaM1 5,918,836 Δ4,582 bp [gcvP‑II CDS (glycine dehydrogenase)][gcvT‑II CDS (aminomethyltransferase)] [gcvP‑II CDS (glycine dehydrogenase)], gcvH‑II CDS (glycine cleavage system H protein 2), [gcvT‑II CDS (aminomethyltransferase)]
RA AE015451_phaM1 6,011,161 T→C intergenic (+90/+20) xpt CDS (Xanthine phosphoribosyltransferase) → / ← PP_5266 CDS (putative acetyl‑CoA hydrolase family protein) –/–
MC JC AE015451_phaM1 6,159,202 Δ6,596 bp PP_5404 CDS (Transposition protein_ TnsC‑related protein)PP_5407 CDS (Transposition protein_ TnsD‑related protein) PP_5404 CDS (Transposition protein_ TnsC‑related protein), PP_5405 CDS (transposase), PP_5406 CDS (transposase), PP_5407 CDS (Transposition protein_ TnsD‑related protein)
RA pS4413_qbsFH 87 (CCGCG)1→2 intergenic (–/‑1866)  / → aadA CDS –/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ AE015451_phaM1 2529185–2530862 2530883 22–1699 380 [338] [0] 978 PP_2218 CDS (transposase)–[PP_2219 CDS (conserved protein of unknown function)] PP_2218 CDS (transposase),[PP_2219 CDS (conserved protein of unknown function)]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451_phaM1 = 1780066257 (0.320)169 (0.230) 69/256 2.7 44.4% intergenic (+25/+81) PP_1587 CDS (NhaP‑type Na+_H+ and K+_H+ antiporters)/dapC CDS (N‑succinyl‑L_L‑diaminopimelate aminotransferase alternative) –/–
?AE015451_phaM1 1780129 = 179 (0.240)intergenic (+88/+18) PP_1587 CDS (NhaP‑type Na+_H+ and K+_H+ antiporters)/dapC CDS (N‑succinyl‑L_L‑diaminopimelate aminotransferase alternative) –/–
* ? AE015451_phaM1 2301453 =0 (0.000)956 (1.240) 112/264 0.4 100% intergenic (‑97/‑234) sbcD CDS (Exonuclease SbcD)/PP_5497 CDS (conserved protein of unknown function with rhs repeat‑associated core domain) –/–
?AE015451_phaM1 2925429 = NA (NA)intergenic (‑177/+388) PP_2563 CDS (putative Antibiotic biosynthesis protein)/PP_2564 CDS (putative Quercetin 2_3‑dioxygenase) –/–
* ? AE015451_phaM1 = 23014540 (0.000)816 (1.050) 131/266 0.1 100% intergenic (‑98/‑233) sbcD CDS (Exonuclease SbcD)/PP_5497 CDS (conserved protein of unknown function with rhs repeat‑associated core domain) –/–
?AE015451_phaM1 5545321 = NA (NA)intergenic (‑461/+416) rpsF CDS (30S ribosomal protein S6)/rlmB CDS (23S rRNA (guanosine‑2'‑O‑)‑methyltransferase RlmB) –/–
* ? AE015451_phaM1 = 2530859NA (NA)967 (1.250) 111/264 0.4 100% coding (98/369 nt) PP_2219 CDS (conserved protein of unknown function)
?AE015451_phaM1 2530884 = 0 (0.000)coding (123/369 nt) PP_2219 CDS (conserved protein of unknown function)
* ? AE015451_phaM1 4496471 =28 (0.040)902 (1.140) 125/270 0.2 97.0% coding (781/960 nt) PP_3987 CDS (Site‑specific recombinase_ phage integrase family)
?AE015451_phaM1 4536203 = NA (NA)intergenic (+730/‑85) PP_4022 CDS (conserved protein of unknown function)/PP_4024 CDS (transposase) –/–
* ? AE015451_phaM1 = 449647883 (0.100)872 (1.100) 137/270 0.1 91.3% coding (788/960 nt) PP_3987 CDS (Site‑specific recombinase_ phage integrase family)
?AE015451_phaM1 = 4538243 NA (NA)intergenic (+45/‑275) PP_4025 CDS (transposase)/PP_4027 CDS (conserved protein of unknown function) –/–
* ? AE015451_phaM1 = 5699642678 (0.860)9721 (1.310) 144/248 0.0 96.9% intergenic (+388/‑3617) PP_5002 CDS (conserved protein of unknown function)/PP_5009 CDS (conserved protein of unknown function) –/–
?pS4413_qbsFH 12 = 0 (0.000)intergenic (–/‑1941) –/aadA CDS –/–