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 BS168_reference_genome 143,659 C→T S100F (TCT→TTT)  rpsE → ribosomal protein S5
RA BS168_reference_genome 165,749 Δ1 bp intergenic (+42/‑5) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn
RA BS168_reference_genome 165,751 Δ1 bp intergenic (+44/‑3) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn
RA BS168_reference_genome 165,826 +C intergenic (+1/‑4) trnI‑Asn → / → trnI‑Thr tRNA‑Asn/tRNA‑Thr
RA BS168_reference_genome 166,037 +T intergenic (+4/‑27) trnI‑Gly → / → trnI‑Arg tRNA‑Gly/tRNA‑Arg
RA BS168_reference_genome 166,345 Δ1 bp intergenic (+17/‑155) trnI‑Ala → / → rrnH‑16S tRNA‑Ala/ribosomal RNA‑16S
RA BS168_reference_genome 376,025 (G)8→7 intergenic (‑227/‑7) hxlA ← / → hxlR 3‑hexulose‑6‑phosphate synthase (HPS)/positive regulator of hxlAB expression (formaldehyde sensing)
RA BS168_reference_genome 557,870 (T)5→6 intergenic (+421/‑3) yddT → / → ydzN putative exported phage protein/hypothetical protein
RA BS168_reference_genome 892,022 A→T N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
RA BS168_reference_genome 944,623 C→T A46V (GCT→GTT)  perR → transcriptional regulator (Fur family)
RA BS168_reference_genome 1,317,152 +GT intergenic (+23/+5) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
RA BS168_reference_genome 1,317,158 (T)5→6 coding (257/258 nt) yjpA ← putative enzyme
RA BS168_reference_genome 1,474,890 C→T L45L (TTG→TTA kre ← regulator of transcription factor ComK function via modulation of mRNA stability
RA BS168_reference_genome 1,685,128 C→G R490G (CGC→GGC)  topA → DNA topoisomerase I
RA BS168_reference_genome 1,693,526 (A)7→8 coding (1031/1611 nt) fliF → flagellar basal‑body M‑ring protein
RA BS168_reference_genome 2,097,087 (A)7→8 intergenic (‑111/+39) azoJ ← / ← yocK FMN‑dependent NADH‑azoreductase/putative general stress protein
RA BS168_reference_genome 2,271,424 T→C R78R (AGA→AGG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,271,505 C→T V51V (GTG→GTA uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,271,523 A→C D45E (GAT→GAG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,387,885 (A)7→8 coding (704/735 nt) yphF ← putative lipoprotein of unknown function
RA BS168_reference_genome 2,480,646 2 bp→AT intergenic (+21/+103) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,480,654 Δ1 bp intergenic (+29/+96) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,480,667 Δ1 bp intergenic (+42/+83) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,581,732 (T)6→7 intergenic (‑115/+39) pstS ← / ← pbpA phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)
RA BS168_reference_genome 2,891,153 G→T A322E (GCG→GAG)  leuB ← 3‑isopropylmalate dehydrogenase
RA BS168_reference_genome 3,075,913 G→A E183K (GAA→AAA)  ytfP → putative NAD(FAD)‑utilizing dehydrogenase
RA BS168_reference_genome 3,770,066 (A)8→9 intergenic (‑258/+38) ywzF ← / ← csbD conserved protein of unknown function/stress response protein
RA BS168_reference_genome 3,935,825 Δ1 bp coding (1190/1191 nt) rlmI ← large subunit ribosomal RNA m5C1962 methyltransferase
RA BS168_reference_genome 4,153,115 G→T L192I (CTT→ATT)  walK ← two‑component sensor histidine kinase [YycG]
RA BS168_reference_genome 4,155,395 (A)5→6 intergenic (‑186/+38) trnY‑Lys ← / ← purA tRNA‑Lys/adenylosuccinate synthetase
JC pSIJ663_pHT315_araE_xy 2,655 Δ3 bp intergenic (–/‑321)  / → eryR –/–
JC pSIJ663_pHT315_araE_xy 3,843 Δ5 bp intergenic (+128/‑24) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 5,049 C→T intergenic (+318/–) ampR → /  –/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160970–161238 165598–164903 3666–4629 13 [9] [11] 23 [rrnI‑16S]–[rrnI‑5S] [rrnI‑16S],rrnI‑23S,[rrnI‑5S]
* * ÷ BS168_reference_genome 166431–166702 171222–170391 3690–4792 15 [11] [12] 13 rrnH‑16S–[rrnH‑5S] rrnH‑16S,rrnH‑23S,[rrnH‑5S]
* * ÷ BS168_reference_genome 171470–171679 176234–171682 4–4765 13 [10] [12] 15 rrnG‑16S–[rrnG‑5S] rrnG‑16S,rrnG‑23S,[rrnG‑5S]
* * ÷ BS168_reference_genome 3485351–3485467 3486758 1292–1408 14 [11] [0] 47 [araE]–[araR] [araE],[araR]
* * ÷ pSIJ663_pHT315_araE_xy 5552 11147–9804 4253–5596 7503 [44] [4225] 4273 ampR/– –/–

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome 3486759 =0 (0.000)27 (0.010)
+50 bp
19/154 0.4 96.8% intergenic (+1/+48) araR/yvbT transcriptional repressor of the ara regulon (LacI family)/putative alkanal monooxygenase
?pSIJ663_pHT315_araE_xy 9781 = 3 (0.000)intergenic (+5050/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 555144 (0.010)7226 (1.060) 205/244 0.0 99.4% intergenic (+820/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 9889 = NA (NA)intergenic (+5158/–) ampR/– –/–