breseq  version 0.35.4  revision f352f80f4bc9
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 70 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,003,2140AT58.1% 2.5 / 15.6 31intergenic (+32/+62)iseA/trnSL‑Arg1inhibitor of DL‑endopeptidases involved in cell‑separation/tRNA‑Arg
*BS168_reference_genome608,2141.A56.7% 44.2 / 45.9 30intergenic (‑139/+32)ydgF/dinBputative amino acid permease/nuclease inhibitor
*BS168_reference_genome3,918,0570CA51.4% 14.1 / 32.5 39V69L (GTA→TTA) ‡qoxAcytochrome aa3‑600 quinol oxidase (subunit II)
*BS168_reference_genome1,291,9230TG47.1% 31.8 / 37.5 68F194V (TTT→GTT) cypEcytochrome P450 CYP109B1, monooxygenase
*BS168_reference_genome2,003,2090TG46.9% 2.0 / 20.9 32intergenic (+27/+67)iseA/trnSL‑Arg1inhibitor of DL‑endopeptidases involved in cell‑separation/tRNA‑Arg
*BS168_reference_genome3,774,6220GC44.7% 39.8 / 34.6 48intergenic (‑61/+33)ywmE/ywmDconserved protein of unknown function/putative exported protein
*BS168_reference_genome238,2350CG41.7% 68.8 / 22.1 70pseudogene (214/285 nt)ybfEcconserved protein of unknown function (C‑terminal end)
*BS168_reference_genome1,185,5610TG40.6% 29.6 / 10.2 32A187A (GCT→GCGyitQconserved protein of unknown function
*BS168_reference_genome3,766,7300CA40.0% 51.7 / 11.7 50R565I (AGA→ATA) ureCurease (alpha subunit)
*BS168_reference_genome238,2320TG39.7% 48.1 / 35.9 68pseudogene (217/285 nt)ybfEcconserved protein of unknown function (C‑terminal end)
*BS168_reference_genome3,485,0040GT38.6% ‑3.9 / 10.2 63Q155K (CAG→AAG) araEarabinose‑related compounds permease
*BS168_reference_genome3,774,6210AC37.8% 55.7 / 13.8 47intergenic (‑60/+34)ywmE/ywmDconserved protein of unknown function/putative exported protein
*BS168_reference_genome4,107,4090TG37.3% 46.8 / 14.7 59K173N (AAA→AACqdoRtranscriptional regulator of qdoI
*BS168_reference_genome3,918,0560AT36.8% 37.4 / 10.8 39V69E (GTA→GAA) ‡qoxAcytochrome aa3‑600 quinol oxidase (subunit II)
*BS168_reference_genome911,9590TG36.4% 30.1 / 12.0 33intergenic (+31/+5)icaC/yfiRacetyl‑glucosamine metabolite exporter component involved in biofilm formation/transcriptional regulator (TetR/AcrR family)
*BS168_reference_genome2,620,9450CG36.1% 43.9 / 11.3 36L77V (CTT→GTT) yqeWputative Na+/anion cotransporter
*BS168_reference_genome3,333,2770TG34.8% 63.5 / 16.9 46F39C (TTT→TGT) pucLurate oxidase with peroxide reductase N‑terminal domain
*BS168_reference_genome3,774,6150TA34.8% 57.1 / 16.9 46intergenic (‑54/+40)ywmE/ywmDconserved protein of unknown function/putative exported protein
*BS168_reference_genome463,3940CA34.3% 97.5 / 17.9 70R50R (CGC→CGAyczIconserved protein of unknown function
*BS168_reference_genome1,641,6480TG34.2% 49.6 / 10.2 39intergenic (+53/‑24)yloC/remAconserved protein of unknown function/essential sporulation DNA binding protein; regulator of biofilm formation

Marginal new junction evidence (sorted from low to high skew)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome 9457 =107 (1.300)1134 (14.480) 4/238 5.5 91.4% coding (2464/2466 nt) gyrA DNA gyrase (subunit A)
?BS168_reference_genome 3484075 = 112 (1.430)coding (1392/1395 nt) araE arabinose‑related compounds permease
* ? BS168_reference_genome = 70063169 (0.840)5 (0.000) 3/244 6.0 0.14% coding (400/1296 nt) purB adenylosuccinate lyase
?pSIJ663_pHT315_araE_xy 673 = 7321 (2.270)intergenic (–/‑2305) –/eryR –/–
* ? BS168_reference_genome = 4128586130 (1.580)4 (0.000)
+T
3/248 6.1 0.23% coding (468/1926 nt) fbp fructose‑1,6‑bisphosphatase
?pSIJ663_pHT315_araE_xy = 5098 3442 (1.040)intergenic (+367/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 50483456 (1.050)157 (0.050) 16/248 9.1 4.1% intergenic (+317/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 5052 = 3855 (1.180)intergenic (+321/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 56147 (0.000)15 (0.000) 13/244 9.9 68.5% intergenic (+883/–) ampR/– –/–
?pSIJ663_pHT315_araE_xy 9532 = 7 (0.000)intergenic (+4801/–) ampR/– –/–
* ? pSIJ663_pHT315_araE_xy = 18386854 (2.080)9 (0.000) 8/246 11.8 0.17% intergenic (–/‑1140) –/eryR –/–
?pSIJ663_pHT315_araE_xy 2572 = 3960 (1.220)intergenic (–/‑406) –/eryR –/–