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 23,520 Δ1 bp 6.2% coding (250/624 nt) dgk ← deoxyguanosine kinase
RA BS168_reference_genome 56,607 G→A 9.5% A86T (GCA→ACA)  glmU → bifunctional glucosamine‑1‑phosphate N‑acetyltransferase/UDP‑N‑acetylglucosamine pyrophosphorylase
RA BS168_reference_genome 77,835 A→G 6.7% L284L (CTA→CTG ftsH → ATP‑dependent cytoplasmic membrane protease
RA BS168_reference_genome 89,128 T→A 5.7% V134V (GTT→GTA lysS → lysyl‑tRNA synthetase
RA BS168_reference_genome 96,240 G→A 7.6% 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 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 182,667 A→T 8.0% D219E (GAT→GAA feuA ← ferri‑bacillibactin‑binding lipoprotein
RA BS168_reference_genome 233,974 T→C 7.5% intergenic (‑79/+20) glpQA ← / ← glpT secreted glycerophosphoryl diester phosphodiesterase/sn‑glycerol‑3‑phosphate permease
RA BS168_reference_genome 335,553 T→A 8.3% L75L (CTT→CTA ycgG → conserved protein of unknown function (similar to phage protein)
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 378,514 T→G 16.0% V516G (GTT→GGT)  srfAA → surfactin synthetase
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,164 C→A 29.5% intergenic (‑89/+82) ydgF ← / ← dinB putative amino acid permease/nuclease inhibitor
RA BS168_reference_genome 608,185 G→T 43.1% intergenic (‑110/+61) ydgF ← / ← dinB putative amino acid permease/nuclease inhibitor
RA BS168_reference_genome 682,074 C→T 8.3% A97T (GCT→ACT)  ydjP ← putative aminoacrylate hydrolase
RA BS168_reference_genome 770,241 G→T 7.7% R3M (AGG→ATG)  rhgW → rhamnogalacturonan endolyase
RA BS168_reference_genome 771,300 C→T 100% A356V (GCG→GTG)  rhgW → rhamnogalacturonan endolyase
RA BS168_reference_genome 892,022 A→T 100% N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon
RA BS168_reference_genome 961,032:1 +A 100% intergenic (‑46/‑362) katA ← / → ssuB vegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein)
RA BS168_reference_genome 1,012,003 G→T 7.3% A419E (GCA→GAA)  lytF ← gamma‑D‑glutamate‑meso‑diaminopimelate muropeptidase (major autolysin, vegetative)
RA BS168_reference_genome 1,023,150 C→A 34.8% 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,029,038 T→G 25.1% T184P (ACC→CCC)  asiMA ← negative regulator of the activity of sigma‑M
RA BS168_reference_genome 1,042,924 C→T 10.3% S441S (TCG→TCA nhaC ← Na+/H+ antiporter
RA BS168_reference_genome 1,126,886 G→A 9.1% F189F (TTC→TTT ntdC ← biosynthesis of neotrehalosadiamine (3,3'‑diamino‑3,3'‑dideoxy‑alpha,beta‑trehalose), glucose‑6‑P 3‑dehydrogenase; kanosamine biosynthesis (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,306,435 A→G 8.7% H317R (CAT→CGT)  uxuA → D‑mannonate dehydratase
RA BS168_reference_genome 1,317,152:1 +G 88.7% intergenic (+23/+5) phrA → / ← yjpA secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme
RA BS168_reference_genome 1,317,152:2 +T 84.4% 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,320,221 G→C 8.4% A178P (GCA→CCA)  xpdC → phage PBSX; manganese‑containing peroxidase
RA BS168_reference_genome 1,323,881 C→T 100% P27L (CCG→CTG)  xkdD → phage PBSX; conserved hypothetical protein
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,498,659 T→C 7.3% V72A (GTA→GCA)  mobB → molybdopterin‑guanine dinucleotide biosynthesis protein B
RA BS168_reference_genome 1,511,930 G→T 32.9% S88* (TCA→TAA)  ykoA ← putative membrane protein of unknown function
RA BS168_reference_genome 1,527,317 G→A 8.0% Q29Q (CAG→CAA ykyA → putative lipoprotein
RA BS168_reference_genome 1,542,428 G→T 11.5% L181L (CTG→CTT ylaA → conserved protein of unknown function
RA BS168_reference_genome 1,561,301 C→T 8.3% P279L (CCG→CTG)  ctaC → cytochrome caa3 oxidase (subunit II)
RA BS168_reference_genome 1,565,793 G→T 16.2% intergenic (‑84/‑56) ylbA ← / → ylbB conserved hypothetical protein/putative enzyme
RA BS168_reference_genome 1,614,239 A→T 6.2% I295F (ATT→TTT)  ileS → isoleucyl‑tRNA synthetase
RA BS168_reference_genome 1,616,570 C→A 15.5% Q102K (CAA→AAA)  ylyA → sporulation‑related RNA polymerase‑binding factor
RA BS168_reference_genome 1,690,759 G→A 100% R214H (CGT→CAT)  codY → transcriptional regulator, GTP and BCAA‑dependent
RA BS168_reference_genome 1,999,145 A→T 30.1% I224K (ATA→AAA)  dacC ← D‑alanyl‑D‑alanine carboxypeptidase
RA BS168_reference_genome 2,048,992 A→G 9.9% I154V (ATC→GTC)  penP → beta‑lactamase precursor (putative hydrolase)
RA BS168_reference_genome 2,067,071 A→T 7.0% F653L (TTT→TTA yobI ← putative phage repair NTPase with transmembrane helices; putative defective prophage 6
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,165,244 A→T 8.3% V124E (GTG→GAG)  nrdEB ← SPbeta phage ribonucleoside reductase alpha (large) subunit
RA BS168_reference_genome 2,190,712 A→T 6.6% V25D (GTC→GAC)  pghZ ← gamma‑polyglutamate hydrolase; phage SPbeta
RA BS168_reference_genome 2,198,642 C→T 10.5% A42T (GCT→ACT)  yoqH ← conserved hypothetical protein; phage SPbeta
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,298,184 C→A 9.0% G112C (GGT→TGT)  ypjQ ← putative phosphatidylglycerophosphatase
RA BS168_reference_genome 2,308,345 C→T 9.4% intergenic (‑6/‑150) degR ← / → ypzA activator of degradative enzymes (aprE, nprE, sacB) production or activity/putative spore coat protein
RA BS168_reference_genome 2,323,088 G→A 9.5% A171V (GCC→GTC)  kdgA ← 2‑keto‑3‑deoxygluconate‑6‑phosphate aldolase
RA BS168_reference_genome 2,332,900 T→A 7.9% T38S (ACA→TCA)  cotD ← spore coat protein (inner)
RA BS168_reference_genome 2,345,064 C→T 6.1% G117G (GGG→GGA nth ← endonuclease III, apurinic apyrimidic DNA lyase
RA BS168_reference_genome 2,366,434 C→T 5.5% P395P (CCG→CCA ypiA ← conserved protein of unknown function
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,408,513 T→G 25.4% K317N (AAA→AAC gmmB ← subunit of the platform for genome maintenance machinery
RA BS168_reference_genome 2,418,340 T→C 7.8% I274V (ATC→GTC)  resC ← factor required for cytochrome c synthesis
RA BS168_reference_genome 2,439,190 C→T 6.6% E222K (GAG→AAG)  spoVAF ← stage V sporulation protein AF
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,497,144 T→A 7.0% T309T (ACA→ACT bkdB ← branched‑chain alpha‑keto acid dehydrogenase E2 subunit (lipoamide acyltransferase)
RA BS168_reference_genome 2,575,387 A→G 7.2% C112R (TGC→CGC)  yqgM ← putative glycosyltransferase
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,652,544 T→G 13.1% D85A (GAC→GCC)  nucB ← sporulation‑specific matrix degrading exported DNase
RA BS168_reference_genome 2,805,209 G→A 21.4% G185S (GGC→AGC)  glnP → glutamine ABC transporter (permease)
RA BS168_reference_genome 2,820,393 G→A 6.6% A28V (GCA→GTA)  dtd ← gly‑tRNA(Ala) deacylase / D‑Tyr‑tRNATyr deacylase
RA BS168_reference_genome 2,865,943 A→T 8.2% F221Y (TTC→TAC)  folC ← folyl‑polyglutamate synthase
RA BS168_reference_genome 2,910,929 A→G 6.4% intergenic (‑183/+187) lysCA ← / ← uvrC aspartokinase II alpha subunit (aa 1‑>408)/excinuclease ABC (subunit C)
RA BS168_reference_genome 3,019,150 C→T 8.2% Q103Q (CAG→CAA ytfI ← conserved sporulation‑related protein
RA BS168_reference_genome 3,019,221 T→C 8.8% S80G (AGC→GGC)  ytfI ← conserved sporulation‑related protein
RA BS168_reference_genome 3,099,285 T→C 7.1% S55P (TCA→CCA)  msmF → carbohydrate ABC transporter (permease)
RA BS168_reference_genome 3,139,266 T→C 79.0% intergenic (‑40/‑165) rpmEB ← / → ythA ribosomal protein L31/putative cytochrome bd menaquinol oxidase subunit I
RA BS168_reference_genome 3,290,436 G→A 7.4% T350T (ACC→ACT dhbC ← isochorismate synthase (siderophore‑specific)
JC BS168_reference_genome 3,383,490 (ATTATAATT)2→1 15.0% intergenic (+15/‑67) yusZ → / → mrgA putative short‑chain acyl dehydrogenase/metalloregulation DNA‑binding stress protein
RA BS168_reference_genome 3,451,374 G→A 6.9% G115G (GGC→GGT smpB ← tmRNA‑binding protein
RA BS168_reference_genome 3,516,325 A→C 26.0% S187A (TCC→GCC)  epsM ← putative O‑acetyltransferase involved in biofilm matrix formation
RA BS168_reference_genome 3,550,926 C→T 9.7% A192A (GCG→GCA mdxK ← maltose phosphorylase
RA BS168_reference_genome 3,644,663 C→T 7.0% D212N (GAC→AAC)  degU ← two‑component response regulator
RA BS168_reference_genome 3,667,987 Δ1 bp 7.3% coding (1925/2703 nt) ggaB ← poly(glucosyl N‑acetylgalactosamine 1‑phosphate) glucosyltransferase
RA BS168_reference_genome 3,728,415 T→C 6.3% N16S (AAC→AGC)  BSU_36215 ← hypothetical protein
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,782,538 G→A 7.9% F125F (TTC→TTT atpD ← ATP synthase (subunit beta, component F1)
RA BS168_reference_genome 3,820,834 T→C 6.7% K216K (AAA→AAG ywjA ← putative ABC lipid transporter (ATP‑binding protein)
RA BS168_reference_genome 3,863,729 G→T 28.7% Q178K (CAA→AAA)  lipL ← amidotransferase of the octanoyl moiety from octanoyl‑GcvH to the E2 subunit of pyruvate dehydrogenase
RA BS168_reference_genome 3,913,926 C→A 16.7% intergenic (‑413/‑83) rodA ← / → ywcE glycosyltransferase involved in extension of the lateral walls of the cell/holin‑like protein required for proper spore morphogenesis and germination
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,958,235 C→A 8.5% A282E (GCG→GAG)  ilvK → branched‑chain amino acid aminotransferase
RA BS168_reference_genome 4,011,659 T→G 13.5% L419R (CTT→CGT)  citH → secondary transporter of divalent metal ions/citrate complexes
RA BS168_reference_genome 4,031,276 T→C 7.8% R124G (AGA→GGA)  yxxF ← putative transporter
RA BS168_reference_genome 4,068,272 T→A 6.9% K89* (AAG→TAG)  yxeA ← conserved protein of unknown function
RA BS168_reference_genome 4,084,116 T→C 8.0% I90V (ATT→GTT)  iolA ← methylmalonate‑semialdehyde dehydrogenase
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,165,287 T→A 5.2% L112F (TTA→TTT gdpP ← phosphodiesterase acting on cyclic dinucleotides
RA BS168_reference_genome 4,196,097 T→A 5.5% D79V (GAT→GTT)  maa ← maltose O‑acetyltransferase
RA pSIJ663_pHT315_araE_xy 121 G→C 5.7% intergenic (–/‑2857)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 131 G→C 6.3% intergenic (–/‑2847)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 137 A→C 5.4% intergenic (–/‑2841)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 2,159 C→T 5.1% intergenic (–/‑819)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 2,206 C→A 9.6% intergenic (–/‑772)  / → eryR –/–
JC pSIJ663_pHT315_araE_xy 2,210 +TCTCG 12.9% intergenic (–/‑768)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 2,239 A→G 7.4% intergenic (–/‑739)  / → eryR –/–
JC pSIJ663_pHT315_araE_xy 2,655 Δ3 bp 100% intergenic (–/‑321)  / → eryR –/–
RA pSIJ663_pHT315_araE_xy 3,729 T→G 7.2% intergenic (+14/‑142) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,735 A→T 9.9% intergenic (+20/‑136) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,741 G→T 10.1% intergenic (+26/‑130) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,743 A→T 5.6% intergenic (+28/‑128) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,749 A→T 5.8% intergenic (+34/‑122) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 3,783 G→T 5.7% intergenic (+68/‑88) eryR → / → ampR –/–
JC pSIJ663_pHT315_araE_xy 3,843 Δ5 bp 100% intergenic (+128/‑24) eryR → / → ampR –/–
RA pSIJ663_pHT315_araE_xy 4,625 A→C 7.0% K252T (AAG→ACG)  ampR →
RA pSIJ663_pHT315_araE_xy 4,648 T→A 10.6% Y260N (TAC→AAC) ‡ ampR →
RA pSIJ663_pHT315_araE_xy 4,649 A→C 8.3% Y260S (TAC→TCC) ‡ ampR →
RA pSIJ663_pHT315_araE_xy 5,049 C→T 100% intergenic (+318/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,173 T→G 6.4% intergenic (+442/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,302 A→C 5.1% intergenic (+571/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,464 T→G 6.1% intergenic (+733/–) ampR → /  –/–
RA pSIJ663_pHT315_araE_xy 5,467 T→G 6.2% intergenic (+736/–) ampR → /  –/–

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ BS168_reference_genome 160977–161238 165598–164903 3666–4622 13 [9] [10] 17 [rrnI‑16S]–[rrnI‑5S] [rrnI‑16S],rrnI‑23S,[rrnI‑5S]
* * ÷ BS168_reference_genome 166444–166703 171220–170392 3690–4777 11 [10] [10] 11 rrnH‑16S–[rrnH‑5S] rrnH‑16S,rrnH‑23S,[rrnH‑5S]
* * ÷ BS168_reference_genome 171471–171679 176232–172533 855–4762 12 [10] [10] 11 rrnG‑16S–[rrnG‑5S] rrnG‑16S,rrnG‑23S,[rrnG‑5S]
* * ÷ BS168_reference_genome 3485377–3485468 3486758 1291–1382 13 [7] [3] 46 [araE]–[araR] [araE],[araR]

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome = 265333346 (0.960)3 (0.070) 3/222 NT 9.7% pseudogene (341/471 nt) sigKn RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
?BS168_reference_genome 2701367 = 12 (0.260)pseudogene (30/417 nt) sigKc RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)
* ? BS168_reference_genome 3486759 =3 (0.060)20 (0.010)
+50 bp
15/132 NT 94.7% intergenic (+1/+48) araR/yvbT transcriptional repressor of the ara regulon (LacI family)/putative alkanal monooxygenase
?pSIJ663_pHT315_araE_xy 9781 = 1 (0.000)intergenic (+5050/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 555114 (0.000)6619 (1.020) 193/222 NT 99.8% intergenic (+820/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 9889 = NA (NA)intergenic (+5158/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 56145 (0.000)11 (0.000) 10/226 NT 69.0% intergenic (+883/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 9532 = 5 (0.000)intergenic (+4801/–) ampR/– –/–