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 39 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,536,1630GA49.8% 337.7 / inf 448intergenic (‑12/+132)rrf/PA14_RS253605S ribosomal RNA/23S ribosomal RNA
*NC_0084636,206,0690TG48.7% 628.9 / inf 729intergenic (‑174/+23)PA14_RS28345/PA14_RS28350type VI secretion system tip protein VgrG/Hcp family type VI secretion system effector
*NC_008463114,7940AG44.6% 925.6 / inf 838K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
*NC_0084635,461,0390AG39.4% 695.5 / inf 512G1272G (GGT→GGCPA14_RS25000MBG domain‑containing protein
*NC_0084635,461,0090AG34.3% 761.9 / inf 461N1282N (AAT→AACPA14_RS25000MBG domain‑containing protein
*NC_0084635,461,1350AG33.0% 627.1 / inf 590D1240D (GAT→GACPA14_RS25000MBG domain‑containing protein
*NC_008463114,8650TC32.8% 1580.9 / inf 756Y529H (TAC→CAC) vgrG1btype VI secretion system spike protein VgrG1b
*NC_0084634,616,4540TG31.4% 316.2 / 25.3 161intergenic (‑174/+28)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084634,616,4500TG27.5% 325.6 / 14.2 167intergenic (‑170/+32)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084632,954,0030TG25.7% 845.7 / 54.9 376Q75P (CAG→CCG) PA14_RS13685DUF6162 family protein
*NC_0084635,910,2560AC24.8% 712.6 / 35.7 308intergenic (+55/‑90)aceE/aceFpyruvate dehydrogenase (acetyl‑transferring), homodimeric type/dihydrolipoyllysine‑residue acetyltransferase
*NC_0084635,460,2830AG23.6% 1039.8 / 282.4 432N1524N (AAT→AACPA14_RS25000MBG domain‑containing protein
*NC_0084634,590,2130TC23.3% 1419.2 / inf 516T254T (ACT→ACCPA14_RS21015site‑specific integrase
*NC_0084634,616,4480GC22.8% 362.1 / 33.5 164intergenic (‑168/+34)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084635,125,7990TG22.8% 363.8 / 15.6 206intergenic (‑54/+32)PA14_RS23455/PA14_RS23460glutathione S‑transferase N‑terminal domain‑containing protein/cytochrome c1
*NC_0084631,763,5340TG22.8% 395.9 / 16.8 259L13R (CTC→CGC) PA14_RS08240RHS repeat protein
*NC_0084633,515,8640TA22.6% 624.9 / 37.3 295intergenic (+618/+266)PA14_RS16005/PA14_RS16010AAA family ATPase/DUF4354 family protein
*NC_0084635,248,9500TG22.4% 882.1 / 10.8 351L666R (CTG→CGG) tpsA4two‑partner secretion system exoprotein TpsA4
*NC_0084632,621,4640AT22.4% 913.7 / 71.9 403intergenic (‑9/+29)aat/trxBleucyl/phenylalanyl‑tRNA‑‑protein transferase/thioredoxin‑disulfide reductase
*NC_0084636,352,9930AT22.2% 1089.6 / 74.3 434intergenic (+73/‑46)dgcB/PA14_RS29035dimethylglycine demethylation protein DgcB/electron transfer flavoprotein subunit alpha/FixB family protein

Marginal new junction evidence (lowest skew 10 of 504 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_008463 = 5935541482 (0.950)35 (0.080) 31/258 23.5 7.9% coding (251/981 nt) PA14_RS27150 hypothetical protein
?NC_008463 5935938 = 407 (0.940)coding (648/981 nt) PA14_RS27150 hypothetical protein
* ? NC_008463 = 729560471 (0.930)39 (0.080) 32/298 25.6 7.7% coding (157/1092 nt) PA14_RS03435 anhydro‑N‑acetylmuramic acid kinase
?NC_008463 729571 = 462 (0.920)coding (146/1092 nt) PA14_RS03435 anhydro‑N‑acetylmuramic acid kinase
* ? NC_008463 = 229722466 (0.920)22 (0.050) 22/270 28.0 4.6% coding (676/882 nt) PA14_RS01030 triphosphoribosyl‑dephospho‑CoA synthase
?NC_008463 230146 = 488 (1.070)coding (217/300 nt) PA14_RS01035 malonate decarboxylase subunit delta
* ? NC_008463 = 729041527 (1.040)30 (0.060) 22/300 29.8 5.3% coding (676/1092 nt) PA14_RS03435 anhydro‑N‑acetylmuramic acid kinase
?NC_008463 729055 = 536 (1.060)coding (662/1092 nt) PA14_RS03435 anhydro‑N‑acetylmuramic acid kinase
* ? NC_008463 = 473590623 (1.230)19 (0.040) 17/260 29.8 3.4% coding (4727/7431 nt) chpA chemotaxis signal transduction system protein ChpA
?NC_008463 474032 = 548 (1.250)coding (5169/7431 nt) chpA chemotaxis signal transduction system protein ChpA
* ? NC_008463 = 6439354561 (1.110)18 (0.040) 16/260 30.4 3.7% coding (973/1290 nt) PA14_RS29425 MFS transporter
?NC_008463 6439844 = 445 (1.020)coding (483/1290 nt) PA14_RS29425 MFS transporter
* ? NC_008463 = 3304506479 (0.950)15 (0.030) 15/268 31.5 3.4% coding (1023/1449 nt) PA14_RS15105 aldehyde dehydrogenase family protein
?NC_008463 3304825 = 417 (0.920)coding (704/1449 nt) PA14_RS15105 aldehyde dehydrogenase family protein
* ? NC_008463 = 3143295468 (0.930)15 (0.030) 14/268 32.0 3.4% coding (731/783 nt) PA14_RS14385 TIM barrel protein
?NC_008463 3143692 = 443 (0.980)coding (334/783 nt) PA14_RS14385 TIM barrel protein
* ? NC_008463 = 1530568687 (1.360)15 (0.030) 14/272 32.3 2.7% coding (1099/1728 nt) PA14_RS07160 acyl‑CoA dehydrogenase family protein
?NC_008463 1531091 = 455 (0.990)coding (1622/1728 nt) PA14_RS07160 acyl‑CoA dehydrogenase family protein
* ? NC_008463 = 419388492 (0.970)14 (0.030) 13/278 33.3 2.8% coding (1484/1734 nt) ggt gamma‑glutamyltransferase
?NC_008463 419804 = 503 (1.070)coding (1068/1734 nt) ggt gamma‑glutamyltransferase