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 77 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,820,7280GA52.6% ‑0.8 / 115.8 78D4270N (GAT→AAT) pksLpolyketide synthase of type I (bacillaene synthesis)
*BS168_reference_genome2,003,2180GT51.1% 16.7 / 54.5 45intergenic (+36/+58)iseA/trnSL‑Arg1inhibitor of DL‑endopeptidases involved in cell‑separation/tRNA‑Arg
*BS168_reference_genome608,2141.A50.0% 47.2 / 84.0 46intergenic (‑139/+32)ydgF/dinBputative amino acid permease/nuclease inhibitor
*Exported4,1370CT49.2% 1330.6 / inf 8122L89L (CTC→CTTampR/vntifkey=4
*BS168_reference_genome2,003,2090TG48.9% 8.3 / 36.0 47intergenic (+27/+67)iseA/trnSL‑Arg1inhibitor of DL‑endopeptidases involved in cell‑separation/tRNA‑Arg
*BS168_reference_genome911,9590TG47.8% ‑2.8 / 13.9 23intergenic (+31/+5)icaC/yfiRacetyl‑glucosamine metabolite exporter component involved in biofilm formation/transcriptional regulator (TetR/AcrR family)
*Exported4,1140AG40.5% 4115.8 / inf 6732I82V (ATT→GTT) ampR/vntifkey=4
*BS168_reference_genome515,0390GT37.0% 20.5 / 10.4 30Q193K (CAA→AAA) ydcAputative rhomboid protease
*BS168_reference_genome223,3630CG36.6% 71.0 / 13.3 44K209N (AAG→AACybdMputative protein kinase
*BS168_reference_genome2,748,6900TG35.8% 64.2 / 16.3 67K43T (AAG→ACG) csnchitosanase
*BS168_reference_genome3,335,7560AC35.5% 32.2 / 11.2 31F173V (TTT→GTT) pucExanthine dehydrogenase, iron‑sulfur subunit
*BS168_reference_genome1,023,1020CA35.1% 81.8 / 13.6 57intergenic (+2/‑248)yhdF/bcaPputative NAD(P)‑dependent dehydrogenase/branched‑chain amino acid transporter or methylthioribose‑degradation metabolite exporter
*BS168_reference_genome3,333,2770TG34.9% 50.0 / 20.9 43F39C (TTT→TGT) pucLurate oxidase with peroxide reductase N‑terminal domain
*BS168_reference_genome2,003,2140AT34.8% 12.9 / 24.9 46intergenic (+32/+62)iseA/trnSL‑Arg1inhibitor of DL‑endopeptidases involved in cell‑separation/tRNA‑Arg
*BS168_reference_genome463,3940CA34.4% 79.7 / 23.1 64R50R (CGC→CGAyczIconserved protein of unknown function
*BS168_reference_genome1,511,9300GT34.3% 25.9 / 10.1 35S88* (TCA→TAA) ykoAputative membrane protein of unknown function
*BS168_reference_genome3,918,0510AT34.1% 55.1 / 12.5 44Y71N (TAC→AAC) qoxAcytochrome aa3‑600 quinol oxidase (subunit II)
*BS168_reference_genome608,1850GT33.9% 106.4 / 10.4 61intergenic (‑110/+61)ydgF/dinBputative amino acid permease/nuclease inhibitor
*BS168_reference_genome238,2320TG33.3% 46.8 / 23.1 54pseudogene (217/285 nt)ybfEcconserved protein of unknown function (C‑terminal end)
*BS168_reference_genome824,6920TG32.6% 61.3 / 17.5 46K9T (AAA→ACA) fecEiron‑dicitrate ABC transporter (permease)

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 = 3485387106 (1.140)39 (0.020) 28/242 1.2 2.4% coding (80/1395 nt) araE arabinose‑related compounds permease
?Exported = 5546 3032 (0.700)intergenic (+815/‑342) ampR/lacZalpha /vntifkey=4//vntifkey=4
* ? BS168_reference_genome 3484069 =92 (0.990)25 (0.010)
+34 bp
13/184 2.3 2.2% intergenic (‑295/+3) cggR/araE transcriptional regulator of gapA/arabinose‑related compounds permease
?Exported 6207 = 2936 (0.650)coding (320/324 nt) lacZalpha /vntifkey=4
* ? BS168_reference_genome = 261143491 (0.980)7 (0.000) 7/252 4.9 0.53% intergenic (‑97/+22) cdd/unpK cytidine/deoxycytidine deaminase/undecaprenol kinase
?Exported 5613 = 2524 (0.560)intergenic (+882/‑275) ampR/lacZalpha /vntifkey=4//vntifkey=4
* ? BS168_reference_genome = 9501103 (1.110)6 (0.000) 5/252 5.6 0.42% intergenic (+42/‑309) gyrA/rrnO‑16S DNA gyrase (subunit A)/ribosomal RNA‑16S
?Exported 6212 = 2766 (0.620)intergenic (+1/–) lacZalpha/– /vntifkey=4/–
* ? BS168_reference_genome 2653329 =85 (0.910)5 (0.060) 3/242 6.4 8.2% pseudogene (337/471 nt) sigKn RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (N‑terminal half)
?BS168_reference_genome = 2701361 31 (0.350)pseudogene (24/417 nt) sigKc RNA polymerase sporulation‑specific sigma‑K factor precursor (Sigma‑27) (C‑terminal fragment)
* ? BS168_reference_genome 4171249 =105 (1.130)3 (0.000) 3/250 6.5 0.11% coding (1205/1308 nt) yybO putative permease
?Exported 1541 = 5293 (1.190)intergenic (–/‑1437) –/eryR –//vntifkey=4
* ? BS168_reference_genome 1829200 =122 (1.310)4 (0.000) 3/250 6.5 0.21% coding (7648/12789 nt) pksM trans AT polyketide synthase of type I involved in bacillaene synthesis
?Exported 5197 = 3639 (0.820)intergenic (+466/‑691) ampR/lacZalpha /vntifkey=4//vntifkey=4