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 BS168_reference_genome 49,252 G→C 21.2% Q208H (CAG→CAC yabE → putative cell wall shaping enzyme
RA BS168_reference_genome 96,240 G→A 7.9% intergenic (+19/‑152) trnJ‑Ala → / → rrnW‑16S tRNA‑Ala/ribosomal RNA‑16S
RA BS168_reference_genome 143,659 C→T 100% S100F (TCT→TTT)  rpsE → ribosomal protein S5
RA BS168_reference_genome 154,264 C→T 5.0% Y141Y (TAC→TAT rplM → ribosomal protein L13
RA BS168_reference_genome 165,749 Δ1 bp 100% intergenic (+42/‑5) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn
RA BS168_reference_genome 165,751 Δ1 bp 100% intergenic (+44/‑3) rrnI‑5S → / → trnI‑Asn ribosomal RNA‑5S/tRNA‑Asn
RA BS168_reference_genome 165,825:1 +C 100% noncoding (72/72 nt) trnI‑Asn → tRNA‑Asn
RA BS168_reference_genome 166,037:1 +T 100% intergenic (+4/‑27) trnI‑Gly → / → trnI‑Arg tRNA‑Gly/tRNA‑Arg
RA BS168_reference_genome 166,344 Δ1 bp 100% intergenic (+16/‑156) trnI‑Ala → / → rrnH‑16S tRNA‑Ala/ribosomal RNA‑16S
RA BS168_reference_genome 353,447 T→G 10.3% V197G (GTG→GGG)  ffoR → Fur‑regulated NADPH:ferredoxin oxidoreductase
RA BS168_reference_genome 376,018 Δ1 bp 100% intergenic (‑220/‑14) hxlA ← / → hxlR 3‑hexulose‑6‑phosphate synthase (HPS)/positive regulator of hxlAB expression (formaldehyde sensing)
RA BS168_reference_genome 377,502 C→T 5.7% H179Y (CAT→TAT)  srfAA → surfactin synthetase
RA BS168_reference_genome 389,065 T→G 14.5% L441R (CTT→CGT)  srfAB → surfactin synthetase
RA BS168_reference_genome 409,328 A→G 7.7% D208D (GAT→GAC tcyC ← cystine ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 432,699 T→G 20.1% F110V (TTT→GTT)  pbtN → petrobactin iron‑siderophore ABC transporter (permease)
RA BS168_reference_genome 457,443 A→T 5.4% I141F (ATT→TTT)  pxpA → oxoprolinase subunit A
RA BS168_reference_genome 463,401 G→T 17.6% E53* (GAA→TAA)  yczI → conserved protein of unknown function
RA BS168_reference_genome 507,401 T→G 14.7% V179G (GTT→GGT)  ydbO → putative cation efflux system
RA BS168_reference_genome 515,020 T→A 28.9% Y199F (TAT→TTT)  ydcA ← putative rhomboid protease
RA BS168_reference_genome 515,043 G→C 28.4% C191W (TGC→TGG ydcA ← putative rhomboid protease
RA BS168_reference_genome 515,053 G→C 54.8% S188W (TCG→TGG)  ydcA ← putative rhomboid protease
RA BS168_reference_genome 538,088 G→A 6.5% G10S (GGC→AGC)  conC → ICEBs1 mobile element: subunit of the conjugation machinery
RA BS168_reference_genome 557,865:1 +T 100% intergenic (+416/‑8) yddT → / → ydzN putative exported phage protein/hypothetical protein
RA BS168_reference_genome 608,185 G→T 20.7% intergenic (‑110/+61) ydgF ← / ← dinB putative amino acid permease/nuclease inhibitor
RA BS168_reference_genome 648,025 T→G 14.5% F29C (TTT→TGT)  tatCY → component of the twin‑arginine pre‑protein translocation pathway
RA BS168_reference_genome 670,425 T→G 11.5% A113A (GCT→GCG ydjE → putative sugar kinase (ribokinase family); prophage region 3
RA BS168_reference_genome 672,996 T→A 5.2% L335M (TTG→ATG)  ydjG → putative phage replication protein; prophage region 3
RA BS168_reference_genome 754,709 C→T 5.2% T77I (ACA→ATA)  yesE → conserved hypothetical protein
RA BS168_reference_genome 892,022 A→T 100% N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
RA BS168_reference_genome 944,551 C→T 100% P22L (CCT→CTT)  perR → transcriptional regulator (Fur family)
RA BS168_reference_genome 1,023,150 C→A 30.0% intergenic (+50/‑200) yhdF → / → bcaP putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter
RA BS168_reference_genome 1,023,175 C→A 24.5% intergenic (+75/‑175) yhdF → / → bcaP putative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter
RA BS168_reference_genome 1,025,078 A→C 19.2% T72P (ACG→CCG)  yhdH → putative sodium‑dependent transporter
RA BS168_reference_genome 1,073,325 C→G 5.9% G131G (GGG→GGC scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,083,515 A→T 8.2% I72N (ATT→AAT)  hmoB ← heme‑degrading monooxygenase
RA BS168_reference_genome 1,099,980 G→A 5.9% E274E (GAG→GAA lplJ → lipoate‑protein ligase
RA BS168_reference_genome 1,132,245 A→C 16.8% G380G (GGT→GGG yhjO ← putative permease (recent HGT island)
RA BS168_reference_genome 1,224,521 A→G 100% E356G (GAA→GGA)  oppD → oligopeptide ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 1,254,032 C→T 8.5% R654K (AGG→AAG)  yjcD ← putative ATP‑dependent DNA helicase (sporulation associated)
RA BS168_reference_genome 1,317,152:1 +GT 100% intergenic (+23/+5) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
RA BS168_reference_genome 1,317,153:1 +T 100% intergenic (+24/+4) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
RA BS168_reference_genome 1,349,495 G→T 100% A325A (GCC→GCA pitB ← putative low‑affinity inorganic phosphate transporter
RA BS168_reference_genome 1,474,890 C→T 100% L45L (TTG→TTA kre ← regulator of transcription factor ComK function via modulation of mRNA stability
RA BS168_reference_genome 1,495,070 T→C 9.5% E100G (GAA→GGA)  sppO ← spore protein cse15
RA BS168_reference_genome 1,511,930 G→T 29.9% S88* (TCA→TAA)  ykoA ← putative membrane protein of unknown function
RA BS168_reference_genome 1,517,908 G→A 5.0% R15H (CGC→CAC)  abhA → transcriptional regulator
RA BS168_reference_genome 1,565,798 C→A 19.8% intergenic (‑89/‑51) ylbA ← / → ylbB conserved hypothetical protein/putative enzyme
RA BS168_reference_genome 1,605,538 T→A 100% intergenic (+48/‑92) sigE → / → sigG RNA polymerase sporulation‑specific sigma‑29 factor (sigma‑E)/RNA polymerase sporulation‑specific sigma factor (sigma‑G)
RA BS168_reference_genome 1,616,570 C→A 15.9% Q102K (CAA→AAA)  ylyA → sporulation‑related RNA polymerase‑binding factor
RA BS168_reference_genome 1,690,758 C→A 100% R214S (CGT→AGT)  codY → transcriptional regulator, GTP and BCAA‑dependent
RA BS168_reference_genome 1,908,853 T→G 20.3% F454V (TTT→GTT)  yndD → spore germination receptor subunit
RA BS168_reference_genome 1,914,659 C→A 5.5% I10I (ATC→ATA yndK → putative phage/plasmid replication protein
RA BS168_reference_genome 1,914,857 T→A 9.1% I76I (ATT→ATA yndK → putative phage/plasmid replication protein
RA BS168_reference_genome 1,918,261 C→A 20.0% intergenic (‑5/‑145) lexA ← / → dinK transcriptional repressor of the SOS regulon/cell division inhibitor
RA BS168_reference_genome 1,940,723 A→C 22.4% K33N (AAA→AAC eglS → endo‑1,4‑beta‑glucanase
RA BS168_reference_genome 1,957,976 C→A 26.9% intergenic (‑57/+51) BSU_18275 ← / ← yngK conserved protein of unknown function/putative exported glycoside hydrolase
RA BS168_reference_genome 1,976,395 C→A 5.6% G2052* (GGA→TGA)  ppsC ← non‑ribosomal plipastatin synthetase C involved in synthesis of plipastatin
RA BS168_reference_genome 1,994,711 C→T 100% E1083K (GAA→AAA)  ppsA ← non‑ribosomal plipastatin synthetase A involved in synthesis of plipastatin
RA BS168_reference_genome 2,015,087 C→T 6.5% S103S (AGC→AGT gltC → transcriptional regulator (LysR family) (GltC‑glutamate)
RA BS168_reference_genome 2,021,779 T→A 8.4% N230I (AAC→ATC)  yoaB ← negatively charged metabolite transporter
RA BS168_reference_genome 2,043,841 T→C 5.4% L55S (TTA→TCA)  yozG → putative antitoxin of a toxin / antitoxin system; putative transcriptional regulator
RA BS168_reference_genome 2,086,527 G→A 5.0% A27V (GCA→GTA)  yozB ← putative integral inner membrane protein
RA BS168_reference_genome 2,097,080:1 +A 100% intergenic (‑104/+46) azoJ ← / ← yocK FMN‑dependent NADH‑azoreductase/putative general stress protein
RA BS168_reference_genome 2,119,350 G→T 22.1% intergenic (‑223/‑43) cdaS ← / → mdtK spore diadenylate cyclase/Na+‑driven efflux MatE family protein
RA BS168_reference_genome 2,271,424 T→C 100% R78R (AGA→AGG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,271,505 C→T 100% V51V (GTG→GTA uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
RA BS168_reference_genome 2,271,523 A→C 100% D45E (GAT→GAG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta
MC JC BS168_reference_genome 2,344,013 Δ39 bp 100% coding (2570‑2608/2745 nt) ponA → peptidoglycan glycosyltransferase (penicillin‑binding proteins 1A and 1B)
RA BS168_reference_genome 2,358,909 C→T 8.3% A4T (GCT→ACT)  bshBA ← N‑acetyl‑alpha‑D‑glucosaminyl L‑malate deacetylase 1; malate N‑acetylglucosamine N‑acetyl hydrolase
RA BS168_reference_genome 2,358,944 A→G 7.4% R127R (CGT→CGC mgsA ← methylglyoxal synthase
RA BS168_reference_genome 2,358,977 C→G 7.5% V116V (GTG→GTC mgsA ← methylglyoxal synthase
RA BS168_reference_genome 2,358,986 T→C 7.0% E113E (GAA→GAG mgsA ← methylglyoxal synthase
RA BS168_reference_genome 2,364,921 A→G 7.5% L58L (TTG→CTG)  qcrA ← menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)
RA BS168_reference_genome 2,364,928 A→T 7.2% V55V (GTT→GTA qcrA ← menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)
RA BS168_reference_genome 2,364,949 A→G 8.5% D48D (GAT→GAC qcrA ← menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)
RA BS168_reference_genome 2,379,904 A→G 8.2% H123H (CAT→CAC aroF ← chorismate synthase
RA BS168_reference_genome 2,379,919 T→C 8.9% A118A (GCA→GCG aroF ← chorismate synthase
RA BS168_reference_genome 2,379,922 A→T 11.4% G117G (GGT→GGA aroF ← chorismate synthase
RA BS168_reference_genome 2,387,878:1 +A 100% coding (711/735 nt) yphF ← putative lipoprotein of unknown function
RA BS168_reference_genome 2,436,674 A→T 6.7% Q210H (CAA→CAT ypuA → putative exported protein
RA BS168_reference_genome 2,441,380 T→A 6.3% S138C (AGT→TGT)  spoVAD ← stage V sporulation protein AD (uptake of pyridine‑2,6‑dicarboxylic acid)
RA BS168_reference_genome 2,476,636 A→G 5.2% Q265R (CAG→CGG)  yqjM → NADPH‑dependent flavin oxidoreductase (acting on cinnamaldehyde‑related compounds)
RA BS168_reference_genome 2,480,646 2 bp→AT 100% 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 100% 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 100% intergenic (+42/+83) zwf → / ← gndA glucose‑6‑phosphate 1‑dehydrogenase (NADP‑dependent)/NADP+‑dependent 6‑P‑gluconate dehydrogenase
RA BS168_reference_genome 2,528,711 T→G 27.8% N182T (AAC→ACC)  folD ← methylenetetrahydrofolate dehydrogenase; methenyltetrahydrofolate cyclohydrolase
RA BS168_reference_genome 2,568,536 A→T 7.0% intergenic (‑46/+37) yqgT ← / ← ltaSC putative d,l‑endopeptidase/enzyme responsible for polyglycerolphosphate LTA synthesis
RA BS168_reference_genome 2,581,726:1 +T 100% intergenic (‑109/+45) pstS ← / ← pbpA phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)
RA BS168_reference_genome 2,732,512 T→C 7.3% pseudogene (36/159 nt) yrdDc ← hypothetical protein; C‑terminal part of YrdD
RA BS168_reference_genome 2,772,038 G→T 5.1% A241S (GCC→TCC)  oatA → peptidoglycan O‑acetyltransferase
RA BS168_reference_genome 2,772,560 A→G 8.9% T415A (ACC→GCC)  oatA → peptidoglycan O‑acetyltransferase
RA BS168_reference_genome 2,818,348 G→T 20.6% P5Q (CCA→CAA)  yrzK ← conserved protein of unknown function
RA BS168_reference_genome 2,881,829 T→C 9.1% E321G (GAG→GGG)  lonA ← class III heat‑shock ATP‑dependent LonA protease
RA BS168_reference_genome 3,012,504 A→G 5.8% Y211Y (TAT→TAC argH ← argininosuccinate lyase
RA BS168_reference_genome 3,013,822 A→G 6.0% S175P (TCT→CCT)  argG ← argininosuccinate synthase
RA BS168_reference_genome 3,035,116 A→C 14.9% V107G (GTT→GGT)  dgcP ← diguanylate cyclase
RA BS168_reference_genome 3,097,107 G→A 5.3% G109D (GGC→GAC)  msmR → transcriptional regulator (LacI family)
RA BS168_reference_genome 3,114,257 A→G 6.4% L345L (TTG→CTG)  ytrF ← metabolite permease involved in resistance to cell wall inhibitors
RA BS168_reference_genome 3,162,506 G→T 6.2% M141I (ATG→ATT cotI → spore coat kinase
RA BS168_reference_genome 3,194,626 T→C 5.6% intergenic (‑99/+9) trnSL‑Ala1 ← / ← uppP tRNA‑Ala/undecaprenyl‑pyrophosphate phosphatase
RA BS168_reference_genome 3,363,614 G→A 8.3% L226L (CTC→CTT metN ← methionine ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 3,412,058 G→T 12.6% L125F (TTG→TTT rsiO → acid stress‑sensitive anti sigma factor
RA BS168_reference_genome 3,421,976 A→T 6.2% T69S (ACA→TCA)  yvsG → putative membrane hydrolase
RA BS168_reference_genome 3,478,630 T→G 19.5% R261S (AGA→AGC pgm ← phosphoglycerate mutase
RA BS168_reference_genome 3,531,060 C→G 24.1% G142G (GGC→GGG pnbA → para‑nitrobenzyl esterase (intracellular esterase B)
RA BS168_reference_genome 3,531,092 A→T 25.3% Y153F (TAT→TTT)  pnbA → para‑nitrobenzyl esterase (intracellular esterase B)
RA BS168_reference_genome 3,592,481 T→G 18.0% S305S (TCA→TCC yvoD ← putative integral inner membrane protein with nucleotide binding domain
RA BS168_reference_genome 3,755,776 T→G 15.9% S29R (AGC→CGC)  ywoC ← putative enzyme of isochorismatase family
RA BS168_reference_genome 3,766,421 T→G 20.3% S53S (TCA→TCC ywnA ← putative transcriptional regulator
RA BS168_reference_genome 3,770,058:1 +A 100% intergenic (‑250/+46) ywzF ← / ← csbD conserved protein of unknown function/stress response protein
RA BS168_reference_genome 3,774,612 C→A 27.7% intergenic (‑51/+43) ywmE ← / ← ywmD conserved protein of unknown function/putative exported protein
RA BS168_reference_genome 3,831,264 C→T 5.3% N168N (AAC→AAT ywiC → putative integral inner membrane protein
RA BS168_reference_genome 3,907,379 A→C 30.9% intergenic (‑65/‑465) ywcI ← / → vpr conserved protein of unknown function/extracellular serine protease
RA BS168_reference_genome 3,935,823 Δ1 bp 100% intergenic (+35/+1) ywbE → / ← rlmI conserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
RA BS168_reference_genome 4,019,488 A→C 17.1% L41R (CTT→CGT)  yxiK ← conserved protein of unknown function
RA BS168_reference_genome 4,020,611 A→C 16.5% S99S (TCT→TCG yxzG ← putative nucleic acid binding protein
RA BS168_reference_genome 4,028,640 C→A 7.1% G637* (GGA→TGA)  wapA ← cell wall‑associated tRNA nuclease precursor; intercellular growth inhibitor
RA BS168_reference_genome 4,153,115 G→T 100% L192I (CTT→ATT)  walK ← two‑component sensor histidine kinase [YycG]
RA BS168_reference_genome 4,155,390:1 +A 100% intergenic (‑181/+43) trnY‑Lys ← / ← purA tRNA‑Lys/adenylosuccinate synthetase
RA Exported 121 G→C 5.6% intergenic (–/‑2857)  / → eryR –//vntifkey=4
RA Exported 131 G→C 6.7% intergenic (–/‑2847)  / → eryR –//vntifkey=4
RA Exported 2,179 Δ1 bp 100% intergenic (–/‑799)  / → eryR –//vntifkey=4
JC Exported 2,655 Δ3 bp 100% intergenic (–/‑321)  / → eryR –//vntifkey=4
RA Exported 3,729 T→G 7.8% intergenic (+14/‑142) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,735 A→T 11.1% intergenic (+20/‑136) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,741 G→T 11.2% intergenic (+26/‑130) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,743 A→T 6.4% intergenic (+28/‑128) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,749 A→T 6.3% intergenic (+34/‑122) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,750 C→A 5.6% intergenic (+35/‑121) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,783 G→T 6.0% intergenic (+68/‑88) eryR → / → ampR /vntifkey=4//vntifkey=4
JC Exported 3,843 Δ5 bp 100% intergenic (+128/‑24) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 4,625 A→C 6.5% K252T (AAG→ACG)  ampR → /vntifkey=4
RA Exported 4,648 T→A 10.6% Y260N (TAC→AAC) ‡ ampR → /vntifkey=4
RA Exported 4,649 A→C 8.8% Y260S (TAC→TCC) ‡ ampR → /vntifkey=4
RA Exported 5,049 C→T 100% intergenic (+318/‑839) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,173 T→G 6.4% intergenic (+442/‑715) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,464 T→G 5.4% intergenic (+733/‑424) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 6,006 T→G 6.0% V40G (GTT→GGT)  lacZalpha → /vntifkey=4
RA Exported 6,031 A→C 6.7% A48A (GCA→GCC lacZalpha → /vntifkey=4
RA Exported 6,033 A→C 5.3% H49P (CAT→CCT)  lacZalpha → /vntifkey=4
RA Exported 6,047 A→C 5.4% S54R (AGC→CGC)  lacZalpha → /vntifkey=4
RA Exported 6,256 A→C 5.9% intergenic (+45/–) lacZalpha → /  /vntifkey=4/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160964–164898 165619–164903 6–4656 21 [15] [19] 20 [rrnI‑16S]–[rrnI‑5S] [rrnI‑16S],rrnI‑23S,[rrnI‑5S]
* * ÷ BS168_reference_genome 166444–166813 171227–170391 3579–4784 21 [18] [17] 24 rrnH‑16S–[rrnH‑5S] rrnH‑16S,rrnH‑23S,[rrnH‑5S]
* * ÷ BS168_reference_genome 171467–171679 176227–172978 1300–4761 21 [18] [17] 20 rrnG‑16S–[rrnG‑5S] rrnG‑16S,rrnG‑23S,[rrnG‑5S]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome 2152031 =73 (0.990)5 (0.080) 4/200 NT 7.5% pseudogene (406/420 nt) spsMn spore capsular polysaccharide biosynthesis enzyme fragment; N‑terminal part of SpsM
?BS168_reference_genome = 2286432 60 (0.940)intergenic (+24/‑16) sprA/spsMc serine‑type phage integrase; phage SPbeta/sporulation capsular polysaccharide biosynthesis enzyme C‑terminal end; restored upon phage SPbeta excision