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 | 121,760 | T→C | intergenic (+87/‑159) | rlmG → / → rpoB | 23S rRNA m2G1835 methyltransferase/RNA polymerase (beta subunit) |
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 | 497,435 | C→T | L88L (CTG→CTA) | dctB ← | membrane bound C4‑dicarboxylate binding sensor protein |
RA | BS168_reference_genome | 529,425 | G→A | intergenic (+3/+80) | trnS‑Leu2 → / ← ydcL | tRNA‑Leu/ICEBs1 mobile element: integrase |
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 |
JC | BS168_reference_genome | 961,043 | (46 bp)1→2 | intergenic (‑57/‑351) | katA ← / → ssuB | vegetative catalase 1/aliphatic sulfonate ABC transporter (ATP‑binding protein) |
RA | BS168_reference_genome | 1,221,661 | C→A | T23K (ACA→AAA) | oppB → | oligopeptide ABC transporter (permease) |
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,445,938 | G→A | pseudogene (301/699 nt) | ykvQ → | putative sporulation‑specific glycosylase (HGT island) |
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,580,038 | (T)9→8 | intergenic (+43/‑83) | bshC → / → mraZ | malate glucosamine cysteine ligase/inhibitor of RsmH and transcriptional regulator |
RA | BS168_reference_genome | 1,704,223 | C→T | Q5* (CAA→TAA) | fliZ → | flagellar regulatory 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,373,376 | C→A | R44L (CGT→CTT) | trpB ← | tryptophan synthase (beta subunit) |
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 | 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 | 160977–161667 | 165598–164903 | 3237–4622 | 16 [11] | [10] 14 | [rrnI‑16S]–[rrnI‑5S] | [rrnI‑16S],rrnI‑23S,[rrnI‑5S] |
* | * | ÷ | BS168_reference_genome | 166461–166703 | 171222–170066 | 3364–4762 | 12 [10] | [11] 12 | rrnH‑16S–[rrnH‑5S] | rrnH‑16S,rrnH‑23S,[rrnH‑5S] |
* | * | ÷ | BS168_reference_genome | 171470–171679 | 176227–172979 | 1301–4758 | 18 [9] | [10] 15 | rrnG‑16S–[rrnG‑5S] | rrnG‑16S,rrnG‑23S,[rrnG‑5S] |
* | * | ÷ | BS168_reference_genome | 529439 | 549887 | 20449 | 14 [4] | [10] 14 | ydcL–yddM | ydcL,immA,immR,xis,ydzL,ydcO,BSU_04849,helP,conQ,nicK,ydcS,ydcT,yddA,conB,conC,conD,conE,yddF,conG,cwlT,yddI,yddJ,yddK,rapI,phrI,yddM |
* | * | ÷ | BS168_reference_genome | 3485349–3485465 | 3486758 | 1294–1410 | 12 [10] | [1] 66 | [araE]–[araR] | [araE],[araR] |
* | * | ÷ | pSIJ663_pHT315_araE_xy | 5552 | 11149–9806 | 4255–5598 | 4274 [44] | [2491] 2576 | ampR/– | –/– |
Unassigned new junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | BS168_reference_genome | 3486759 = | 1 (0.020) | 35 (0.030) +50 bp |
22/150 | 0.3 | 96.7% | 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 | = 5551 | 44 (0.010) | 4089 (1.030) | 195/240 | 0.0 | 99.0% | intergenic (+820/–) | ampR/– | –/– |
? | pSIJ663_pHT315_araE_xy | 9889 = | NA (NA) | intergenic (+5158/–) | ampR/– | –/– |