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 46 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,349,4950GT42.4% 135.3 / 290.8 196A325A (GCC→GCApitBputative low‑affinity inorganic phosphate transporter
*BS168_reference_genome3,174,6120TC37.9% 0.2 / 20.9 93noncoding (2307/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,174,6670TC37.9% 0.2 / 44.7 142noncoding (2252/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,174,5990AC36.9% 0.3 / 24.2 85noncoding (2320/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,174,6180AC35.4% 2.9 / 21.0 102noncoding (2301/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,174,5810TC32.4% 0.0 / 11.5 75noncoding (2338/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,174,5640AC30.6% 0.3 / 11.7 64noncoding (2355/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome1,346,7050AG30.1% 500.3 / 64.4 299intergenic (+27/‑26)xepA/xhlAphage PBSX; lytic exoenzyme/phage PBSX; putative enzyme
*BS168_reference_genome2,989,8260TC28.6% 201.6 / 23.7 130intergenic (‑261/+74)accD/maeBacetyl‑CoA carboxylase (carboxyltransferase beta subunit)/NADP‑dependent malic enzyme (conversion of malate into pyruvate, anabolic)
*BS168_reference_genome3,174,6470AC28.2% 10.3 / 21.7 133noncoding (2272/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome2,989,8290GC28.2% 149.7 / 15.1 128intergenic (‑264/+71)accD/maeBacetyl‑CoA carboxylase (carboxyltransferase beta subunit)/NADP‑dependent malic enzyme (conversion of malate into pyruvate, anabolic)
*BS168_reference_genome1,253,1560TG27.6% 577.7 / 69.8 262G18G (GGT→GGGyjcZputative type I toxin
*BS168_reference_genome1,253,1800CG27.4% 466.7 / 69.4 235G26G (GGC→GGGyjcZputative type I toxin
*BS168_reference_genome100,3950TG27.1% 102.2 / 28.7 71noncoding (2286/2928 nt)rrnW‑23Sribosomal RNA‑23S
*BS168_reference_genome2,625,6460AC26.8% 682.1 / 56.3 311G89G (GGT→GGGdnaJco‑factor of molecular chaperone
*BS168_reference_genome3,174,6250AG26.7% 7.0 / 18.3 102noncoding (2294/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,588,4210TC26.7% 599.0 / 58.1 323I314V (ATC→GTC) hisZhistidyl‑tRNA synthetase‑like subunit of ATP phophoribosyltransferase
*BS168_reference_genome3,174,6540AC25.9% 11.9 / 17.4 140noncoding (2265/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome3,174,6090GC25.5% 6.6 / 16.8 94noncoding (2310/2928 nt)rrnB‑23Sribosomal RNA‑23S
*BS168_reference_genome100,3610AG24.8% 226.1 / 21.5 129noncoding (2252/2928 nt)rrnW‑23Sribosomal RNA‑23S

Marginal new junction evidence (lowest skew 10 of 13 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? BS168_reference_genome = 2344012116 (0.350)76 (0.240) 66/282 5.7 40.8% coding (2569/2745 nt) ponA peptidoglycan glycosyltransferase (penicillin‑binding proteins 1A and 1B)
?BS168_reference_genome 2344052 = 111 (0.360)coding (2609/2745 nt) ponA peptidoglycan glycosyltransferase (penicillin‑binding proteins 1A and 1B)
* ? BS168_reference_genome 946697 =NA (NA)276 (0.840)
+G
7/296 18.0 52.2% noncoding (2/1555 nt) rrnD‑16S ribosomal RNA‑16S
?BS168_reference_genome 3178705 = 254 (0.770)intergenic (‑65/‑601) rrnB‑16S/thiT ribosomal RNA‑16S/thiamin permease
* ? BS168_reference_genome = 944791272 (0.820)6 (0.020) 6/298 18.6 2.1% coding (305/438 nt) perR transcriptional regulator (Fur family)
?BS168_reference_genome 944924 = 285 (0.860)coding (438/438 nt) perR transcriptional regulator (Fur family)
* ? BS168_reference_genome 35232 =305 (0.920)457 (1.460)
+TGAGAGTG
4/282 19.4 63.3% intergenic (+129/‑5) rrnA‑23S/rrnA‑5S ribosomal RNA‑23S/ribosomal RNA‑5S
?BS168_reference_genome = 951387 254 (0.770)intergenic (+42/‑70) rrnD‑23S/rrnD‑5S ribosomal RNA‑23S/ribosomal RNA‑5S
* ? BS168_reference_genome 3933706 =442 (1.340)3 (0.010) 3/280 20.0 0.71% coding (130/678 nt) cidB metabolic regulator controlling activity of murein hydrolases
?BS168_reference_genome = 3933702 421 (1.360)coding (126/678 nt) cidB metabolic regulator controlling activity of murein hydrolases
* ? BS168_reference_genome 3683219 =456 (1.380)3 (0.010) 3/282 20.0 0.73% coding (1047/1146 nt) tagB teichoic acid primase, CDP‑glycerol:N‑acetyl‑beta‑d‑mannosaminyl‑1, 4‑N‑acetyl‑d‑glucosaminyldiphosphoundecaprenyl glycerophosphotransferase
?BS168_reference_genome 3683258 = 390 (1.250)coding (1086/1146 nt) tagB teichoic acid primase, CDP‑glycerol:N‑acetyl‑beta‑d‑mannosaminyl‑1, 4‑N‑acetyl‑d‑glucosaminyldiphosphoundecaprenyl glycerophosphotransferase
* ? BS168_reference_genome = 2458690420 (1.270)3 (0.010) 3/294 20.3 0.72% coding (377/558 nt) nudF isopentenyl pyrophosphate and dimethylallyl pyrophosphate diphosphatase (moonlighting ADP‑ribose pyrophosphatase)
?BS168_reference_genome 2458695 = 413 (1.270)coding (372/558 nt) nudF isopentenyl pyrophosphate and dimethylallyl pyrophosphate diphosphatase (moonlighting ADP‑ribose pyrophosphatase)
* ? BS168_reference_genome = 428434365 (1.100)3 (0.010) 3/296 20.3 0.82% coding (1188/1422 nt) yclK two‑component sensor histidine kinase [YclJ]
?BS168_reference_genome 428440 = 363 (1.110)coding (1194/1422 nt) yclK two‑component sensor histidine kinase [YclJ]
* ? BS168_reference_genome = 1650535234 (0.710)3 (0.010) 3/298 20.4 1.3% coding (152/765 nt) prpC multitarget phosphorylated protein phosphatase
?BS168_reference_genome 1650653 = 226 (0.680)coding (270/765 nt) prpC multitarget phosphorylated protein phosphatase
* ? Exported = 39999100 (1.070)13 (0.000) 10/294 24.5 0.15% coding (129/861 nt) ampR /vntifkey=4
?Exported 4255 = 8763 (1.040)coding (385/861 nt) ampR /vntifkey=4