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 19,965 G→C 7.7% intergenic (+19/‑3) pdxS → / → pdxT glutamine amidotransferase for pyridoxal phosphate synthesis; pyridoxal 5'‑phosphate synthase complex, synthase subunit/glutamine amidotransferase for pyridoxal phosphate synthesis; pyridoxal 5'‑phosphate synthase complex, glutamine amidotransferase subunit PdxT
RA BS168_reference_genome 22,700 T→G 17.7% Q150H (CAA→CAC dck ← deoxyadenosine/deoxycytidine kinase
RA BS168_reference_genome 27,319 G→T 5.5% R169I (AGA→ATA)  dnaX → DNA polymerase III subunit tau subunit
RA BS168_reference_genome 66,187 G→T 21.0% L457F (TTG→TTT yabM → putative exporter
RA BS168_reference_genome 67,130 G→T 5.2% E242D (GAG→GAT yabN → putative fusion methylase and nucleotide pyrophosphohydrolase
RA BS168_reference_genome 95,851 T→A 27.4% intergenic (+5/‑10) trnJ‑Gly → / → trnJ‑Leu2 tRNA‑Gly/tRNA‑Leu
RA BS168_reference_genome 96,240 G→A 7.4% intergenic (+19/‑152) trnJ‑Ala → / → rrnW‑16S tRNA‑Ala/ribosomal RNA‑16S
RA BS168_reference_genome 126,612 A→T 5.3% K351* (AAA→TAA)  rpoC → RNA polymerase (beta' subunit)
RA BS168_reference_genome 143,659 C→T 100% S100F (TCT→TTT)  rpsE → ribosomal protein S5
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 196,501 T→G 20.5% F97V (TTT→GTT)  cdaA → diadenylate cyclase
RA BS168_reference_genome 199,191 A→C 13.2% N232T (AAC→ACC)  glmM → phosphoglucosamine mutase
RA BS168_reference_genome 236,174 G→T 10.4% R70S (AGG→AGT ybfA → putative transcriptional regulator with acetyltransferase domain
RA BS168_reference_genome 272,231 T→G 25.9% A144A (GCT→GCG gudD → glucarate dehydratase
RA BS168_reference_genome 360,188 A→C 15.9% L83R (CTG→CGG)  nasC ← assimilatory nitrate reductase (catalytic subunit)
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 557,865:1 +T 100% intergenic (+416/‑8) yddT → / → ydzN putative exported phage protein/hypothetical protein
RA BS168_reference_genome 606,089 A→G 6.5% P97P (CCT→CCC vmlR ← ATP‑binding cassette efflux transporter
RA BS168_reference_genome 608,164 C→A 17.4% intergenic (‑89/+82) ydgF ← / ← dinB putative amino acid permease/nuclease inhibitor
RA BS168_reference_genome 779,310 A→G 10.6% N826S (AAC→AGC)  yetA → putative enzyme
RA BS168_reference_genome 892,022 A→T 100% N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
RA BS168_reference_genome 944,586 T→C 100% S34P (TCT→CCT)  perR → transcriptional regulator (Fur family)
RA BS168_reference_genome 962,064 A→C 18.2% N224T (AAC→ACC)  ssuB → aliphatic sulfonate ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 977,627 G→A 5.1% G187R (GGA→AGA)  yhbJ → putative membrane fusion protein; putative exporter subunit (benzoate transcriptome)
RA BS168_reference_genome 1,023,150 C→A 28.1% 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,051,605 T→A 5.6% N190Y (AAT→TAT)  spaD ← spore coat associated protein
RA BS168_reference_genome 1,073,527 G→C 100% S64C (TCT→TGT)  scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
RA BS168_reference_genome 1,139,565 A→G 8.3% E1082E (GAA→GAG addB → ATP‑dependent deoxyribonuclease (subunit B)
RA BS168_reference_genome 1,174,075 A→T 5.6% D248V (GAT→GTT)  slpS → 2‑phospho‑3‑sulfolactate synthase
RA BS168_reference_genome 1,176,993 G→A 5.6% A87A (GCC→GCT yitG ← putative efflux transporter
RA BS168_reference_genome 1,222,981 C→G 5.5% P150R (CCT→CGT)  oppC → oligopeptide ABC transporter (permease)
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,385,573 C→A 5.0% P12P (CCG→CCT metE ← cobalamin‑independent methionine synthase
RA BS168_reference_genome 1,424,639 T→G 8.1% intergenic (‑205/+128) mtnK ← / ← mtnU methylthioribose kinase (methionine salvage pathway, promiscuous)/ketoglutaramate omega‑amidase
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,602,723 C→T 5.5% R1147R (CGC→CGT bpr → bacillopeptidase F
RA BS168_reference_genome 1,690,758 C→A 100% R214S (CGT→AGT)  codY → transcriptional regulator, GTP and BCAA‑dependent
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,145,824 C→A 17.1% M1F (TTG→TTT) †
C77F (TGC→TTC) 
BSU_19745 ←
yoyG ←
conserved hypothetical protein
putative toxin of a type I toxin family (sporulation operon)
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
RA BS168_reference_genome 2,274,225 T→C 6.8% K100E (AAA→GAA)  yokL ← putative N‑acetyltransferase; phage SPbeta
RA BS168_reference_genome 2,301,634 C→T 6.8% G269R (GGA→AGA)  ilvD ← dihydroxy‑acid dehydratase
RA BS168_reference_genome 2,357,654 T→A 7.6% T186T (ACA→ACT bshA ← N‑acetyl‑alpha‑D‑glucosaminyl L‑malate synthase; malate glycosyltransferase for bacillithiol synthesis
RA BS168_reference_genome 2,360,680 G→A 6.0% A17A (GCG→GCA ypjC → putative integral inner membrane protein
RA BS168_reference_genome 2,362,105 A→G 6.7% I78I (ATT→ATC ypjB ← sporulation septum‑associated protein
RA BS168_reference_genome 2,364,760 A→G 6.2% S111S (AGT→AGC qcrA ← menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)
RA BS168_reference_genome 2,367,169 G→A 10.2% A150A (GCC→GCT ypiA ← conserved protein of unknown function
RA BS168_reference_genome 2,367,180 C→T 10.0% A147T (GCT→ACT)  ypiA ← conserved protein of unknown function
RA BS168_reference_genome 2,367,213 C→T 8.8% V136I (GTT→ATT)  ypiA ← conserved protein of unknown function
RA BS168_reference_genome 2,367,214 A→G 9.0% P135P (CCT→CCC ypiA ← conserved protein of unknown function
RA BS168_reference_genome 2,367,241 C→T 6.9% Q126Q (CAG→CAA ypiA ← conserved protein of unknown function
RA BS168_reference_genome 2,368,719 G→A 6.3% A174A (GCC→GCT aroA ← 3‑phosphoshikimate 1‑carboxyvinyltransferase (5‑enolpyruvoylshikimate‑3‑phosphate synthase)
RA BS168_reference_genome 2,382,141 A→T 7.6% P275P (CCT→CCA hepT ← heptaprenyl diphosphate synthase component II
RA BS168_reference_genome 2,382,153 G→A 7.2% A271A (GCC→GCT hepT ← heptaprenyl diphosphate synthase component II
RA BS168_reference_genome 2,382,177 G→A 6.6% N263N (AAC→AAT hepT ← heptaprenyl diphosphate synthase component II
RA BS168_reference_genome 2,382,180 T→C 5.5% G262G (GGA→GGG hepT ← heptaprenyl diphosphate synthase component II
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,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,508,042 A→G 6.9% G54G (GGT→GGC mmgF ← 2‑methylisocitrate lyase
RA BS168_reference_genome 2,512,448 G→A 6.4% T130I (ACA→ATA)  mmgB ← 3‑hydroxybutyryl‑CoA dehydrogenase
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,610,685 C→T 6.9% V88I (GTC→ATC)  era ← maturation of 16S RNA and assembly of 30S ribosomal subunit GTPase
RA BS168_reference_genome 2,635,354 C→A 7.0% M86I (ATG→ATT gpr ← spore germination protease
RA BS168_reference_genome 2,648,641 A→G 5.6% F4S (TTT→TCT)  yqeF ← putative lipoprotein; putative esterase
RA BS168_reference_genome 2,717,858 G→A 6.6% A296T (GCT→ACT)  blt → efflux transporter (promiscuous, acetylated polyamines efflux)
RA BS168_reference_genome 2,728,535 A→T 9.1% intergenic (‑53/+112) brnQ ← / ← azlD low‑affinity branched‑chain amino acid transporter/branched‑chain amino acid / methionine exporter
RA BS168_reference_genome 2,756,537 G→A 5.9% G76R (GGA→AGA)  adhA → putative aldehyde dehydrogenase; carbonyl stress response
RA BS168_reference_genome 2,818,348 G→T 24.4% P5Q (CCA→CAA)  yrzK ← conserved protein of unknown function
RA BS168_reference_genome 2,856,358 T→C 5.5% K161R (AAA→AGA)  spoIVFB ← membrane metalloprotease [proteolysis of pro‑sigma(K)]
RA BS168_reference_genome 2,871,202 C→T 5.2% A516A (GCG→GCA spoVID ← morphogenetic spore protein (stage VI sporulation)
RA BS168_reference_genome 3,016,944 G→T 10.0% T231N (ACC→AAC)  ytxK ← putative nucleic acid methyltransferase
RA BS168_reference_genome 3,059,481 Δ1 bp 9.2% intergenic (‑80/+66) ytnP ← / ← trmB putative quorum‑quenching lactonase/tRNA (guanine‑N(7)‑)‑methyltransferase
RA BS168_reference_genome 3,059,482 Δ1 bp 100% intergenic (‑81/+65) ytnP ← / ← trmB putative quorum‑quenching lactonase/tRNA (guanine‑N(7)‑)‑methyltransferase
RA BS168_reference_genome 3,116,898 C→A 10.1% L23F (TTG→TTT ytrD ← ABC transporter, permease component involved in resistance to cell wall inhibitors
RA BS168_reference_genome 3,236,265 T→A 16.1% intergenic (+46/+157) yuxK → / ← yufK putative thiol‑disulfide oxidoreductase/putative integral membrane protein
RA BS168_reference_genome 3,349,360 A→T 7.3% M115K (ATG→AAG)  frlN ← fructose‑amino acid permease
RA BS168_reference_genome 3,368,850 C→A 25.0% E389* (GAA→TAA)  fadA ← acetyl‑CoA C‑acyltransferase
RA BS168_reference_genome 3,441,791 A→G 6.2% V580A (GTC→GCC)  copA ← copper transporter ATPase
RA BS168_reference_genome 3,511,543 C→G 10.4% A255G (GCT→GGT)  lutP → L‑lactate permease
RA BS168_reference_genome 3,658,834 G→T 36.9% pseudogene (59/486 nt) tuaAn ← putative undecaprenyl‑phosphate N‑acetylgalactosaminyl‑1‑phosphate transferase; N‑terminal part of TuaA
RA BS168_reference_genome 3,701,855 C→A 6.8% P149T (CCT→ACT)  rbsR → ribose operon repressor (LacI family, D‑ribose)
RA BS168_reference_genome 3,755,184 Δ1 bp 6.7% coding (108/1359 nt) ywoD ← putative efflux transporter
RA BS168_reference_genome 3,755,185 Δ1 bp 7.5% coding (107/1359 nt) ywoD ← putative efflux transporter
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,827,054 T→C 7.6% K964K (AAA→AAG narG ← nitrate reductase (alpha subunit)
RA BS168_reference_genome 3,882,167 A→T 6.0% intergenic (‑153/+24) rocG ← / ← bslB glutamate dehydrogenase/biofilm protein, member of the processed secretome contributing to biofilm hydrophobicity
RA BS168_reference_genome 3,907,393 G→C 30.3% intergenic (‑79/‑451) 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 3,954,693 A→G 6.1% K304R (AAG→AGG)  dltB → putative D‑alanine transferase from DCP to undecaprenylphosphate for lipoteichoic acid and wall teichoic acid
RA BS168_reference_genome 3,961,614 Δ1 bp 5.9% coding (261/309 nt) licB ← phosphotransferase system (PTS) lichenan‑specific enzyme IIB component
RA BS168_reference_genome 4,066,077 A→G 5.8% intergenic (+115/‑133) cotNE → / → yxeD inner spore coat protein/conserved protein of unknown function
RA BS168_reference_genome 4,141,515 A→G 6.3% intergenic (+41/+196) phrG → / ← argI secreted regulator of the activity of phosphatase RapG/arginase
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 BS168_reference_genome 4,201,544 A→G 6.3% D510D (GAT→GAC yyaE ← putative oxidoreductase (molybdoenzyme)
RA Exported 121 G→C 5.6% intergenic (–/‑2857)  / → eryR –//vntifkey=4
RA Exported 131 G→C 6.8% intergenic (–/‑2847)  / → eryR –//vntifkey=4
RA Exported 137 A→C 5.3% intergenic (–/‑2841)  / → eryR –//vntifkey=4
RA Exported 1,922 C→T 100% intergenic (–/‑1056)  / → eryR –//vntifkey=4
JC Exported 2,655 Δ3 bp 100% intergenic (–/‑321)  / → eryR –//vntifkey=4
RA Exported 3,729 T→G 6.9% intergenic (+14/‑142) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,735 A→T 9.7% intergenic (+20/‑136) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,741 G→T 9.8% intergenic (+26/‑130) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,743 A→T 5.6% intergenic (+28/‑128) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,749 A→T 5.7% intergenic (+34/‑122) eryR → / → ampR /vntifkey=4//vntifkey=4
RA Exported 3,783 G→T 5.2% 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.2% K252T (AAG→ACG)  ampR → /vntifkey=4
RA Exported 4,648 T→A 11.3% Y260N (TAC→AAC) ‡ ampR → /vntifkey=4
RA Exported 4,649 A→C 9.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.7% intergenic (+442/‑715) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,302 A→C 5.8% intergenic (+571/‑586) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 5,464 T→G 5.3% intergenic (+733/‑424) ampR → / → lacZalpha /vntifkey=4//vntifkey=4
RA Exported 6,006 T→G 6.2% V40G (GTT→GGT)  lacZalpha → /vntifkey=4
RA Exported 6,031 A→C 6.3% A48A (GCA→GCC lacZalpha → /vntifkey=4
RA Exported 6,256 A→C 6.2% intergenic (+45/–) lacZalpha → /  /vntifkey=4/–
RA Exported 6,277 A→C 5.7% intergenic (+66/–) lacZalpha → /  /vntifkey=4/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160970–161238 165598–164903 3666–4629 16 [15] [13] 18 [rrnI‑16S]–[rrnI‑5S] [rrnI‑16S],rrnI‑23S,[rrnI‑5S]
* * ÷ BS168_reference_genome 166441–166703 176232–172979 6277–9792 17 [15] [15] 16 rrnH‑16S–[rrnG‑5S] rrnH‑16S,rrnH‑23S,rrnH‑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 =58 (0.910)9 (0.160) 9/200 NT 15.5% pseudogene (406/420 nt) spsMn spore capsular polysaccharide biosynthesis enzyme fragment; N‑terminal part of SpsM
?BS168_reference_genome = 2286432 48 (0.870)intergenic (+24/‑16) sprA/spsMc serine‑type phage integrase; phage SPbeta/sporulation capsular polysaccharide biosynthesis enzyme C‑terminal end; restored upon phage SPbeta excision
* ? BS168_reference_genome = 215204661 (0.960)8 (0.150) 8/200 NT 13.7% intergenic (+1/+40) spsMn/sprB spore capsular polysaccharide biosynthesis enzyme fragment; N‑terminal part of SpsM/recombination directionality factor; phage SPbeta
?BS168_reference_genome 2286449 = 48 (0.870)pseudogene (2/606 nt) spsMc sporulation capsular polysaccharide biosynthesis enzyme C‑terminal end; restored upon phage SPbeta excision
* ? BS168_reference_genome 2653329 =56 (0.880)30 (0.490) 22/222 NT 45.3% pseudogene (337/471 nt) sigKn RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
?BS168_reference_genome = 2701361 19 (0.310)pseudogene (24/417 nt) sigKc RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)
* ? BS168_reference_genome = 265333356 (0.880)6 (0.100) 6/222 NT 14.2% pseudogene (341/471 nt) sigKn RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
?BS168_reference_genome 2701367 = 19 (0.310)pseudogene (30/417 nt) sigKc RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)