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 48 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084636,206,0690TG51.0% 591.5 / inf 770intergenic (‑174/+23)PA14_RS28345/PA14_RS28350type VI secretion system tip protein VgrG/Hcp family type VI secretion system effector
*NC_0084635,536,1630GA49.5% 309.3 / inf 414intergenic (‑12/+132)rrf/PA14_RS253605S ribosomal RNA/23S ribosomal RNA
*NC_008463114,7940AG40.5% 1200.1 / inf 851K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
*NC_0084634,616,4540TG39.6% 248.3 / 30.8 158intergenic (‑174/+28)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084635,461,1350AG37.0% 476.1 / inf 514D1240D (GAT→GACPA14_RS25000MBG domain‑containing protein
*NC_0084634,616,4500TG37.0% 256.9 / 23.3 171intergenic (‑170/+32)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084635,461,0390AG34.1% 908.8 / inf 571G1272G (GGT→GGCPA14_RS25000MBG domain‑containing protein
*NC_0084635,461,0090AG33.5% 878.9 / inf 541N1282N (AAT→AACPA14_RS25000MBG domain‑containing protein
*NC_0084635,910,2560AC32.2% 544.9 / 58.6 295intergenic (+55/‑90)aceE/aceFpyruvate dehydrogenase (acetyl‑transferring), homodimeric type/dihydrolipoyllysine‑residue acetyltransferase
*NC_008463114,8650TC30.7% 1715.1 / inf 754Y529H (TAC→CAC) vgrG1btype VI secretion system spike protein VgrG1b
*NC_0084634,616,4390CG29.9% 235.9 / 56.3 157intergenic (‑159/+43)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084631,763,5340TG26.6% 421.3 / 29.9 283L13R (CTC→CGC) PA14_RS08240RHS repeat protein
*NC_0084636,206,2270GA25.7% 137.8 / 62.4 71H128H (CAC→CATPA14_RS28350Hcp family type VI secretion system effector
*NC_0084632,621,4640AT25.4% 1007.9 / 86.7 476intergenic (‑9/+29)aat/trxBleucyl/phenylalanyl‑tRNA‑‑protein transferase/thioredoxin‑disulfide reductase
*NC_0084635,248,9500TG25.2% 925.8 / 22.1 354L666R (CTG→CGG) tpsA4two‑partner secretion system exoprotein TpsA4
*NC_0084632,954,0030TG25.2% 1015.6 / 36.7 428Q75P (CAG→CCG) PA14_RS13685DUF6162 family protein
*NC_0084633,515,8760TA24.7% 663.9 / 49.0 286intergenic (+630/+254)PA14_RS16005/PA14_RS16010AAA family ATPase/DUF4354 family protein
*NC_0084635,460,2830AG24.6% 1114.9 / inf 481N1524N (AAT→AACPA14_RS25000MBG domain‑containing protein
*NC_0084634,616,4570CG24.5% 297.1 / 32.6 159intergenic (‑177/+25)hda/PA14_RS21180DnaA regulatory inactivator Hda/DUF2066 domain‑containing protein
*NC_0084633,630,8250GT24.2% 370.6 / 30.2 213intergenic (+16/‑259)tse1/metHtype VI secretion system effector peptidoglycanhydrolase Tse1/methionine synthase

Marginal new junction evidence (lowest skew 10 of 261 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_008463 = 473590658 (1.290)23 (0.050) 21/260 28.4 4.1% coding (4727/7431 nt) chpA chemotaxis signal transduction system protein ChpA
?NC_008463 474032 = 512 (1.160)coding (5169/7431 nt) chpA chemotaxis signal transduction system protein ChpA
* ? NC_008463 = 1530568614 (1.200)14 (0.030) 12/272 34.1 2.7% coding (1099/1728 nt) PA14_RS07160 acyl‑CoA dehydrogenase family protein
?NC_008463 1531091 = 460 (1.000)coding (1622/1728 nt) PA14_RS07160 acyl‑CoA dehydrogenase family protein
* ? NC_008463 = 229722550 (1.080)11 (0.020) 11/270 34.6 2.2% coding (676/882 nt) PA14_RS01030 triphosphoribosyl‑dephospho‑CoA synthase
?NC_008463 230146 = 499 (1.090)coding (217/300 nt) PA14_RS01035 malonate decarboxylase subunit delta
* ? NC_008463 = 5884344529 (1.040)11 (0.020) 11/274 34.9 2.0% coding (746/1425 nt) hldE bifunctional D‑glycero‑beta‑D‑manno‑heptose‑7‑phosphate kinase/D‑glycero‑beta‑D‑manno‑heptose 1‑phosphate adenylyltransferase HldE
?NC_008463 5884597 = 610 (1.310)coding (493/1425 nt) hldE bifunctional D‑glycero‑beta‑D‑manno‑heptose‑7‑phosphate kinase/D‑glycero‑beta‑D‑manno‑heptose 1‑phosphate adenylyltransferase HldE
* ? NC_008463 = 5935541504 (0.990)10 (0.020) 9/258 35.3 2.2% coding (251/981 nt) PA14_RS27150 hypothetical protein
?NC_008463 5935938 = 465 (1.060)coding (648/981 nt) PA14_RS27150 hypothetical protein
* ? NC_008463 = 270317504 (0.990)9 (0.020) 9/272 36.2 1.8% coding (244/1185 nt) pobA 4‑hydroxybenzoate 3‑monooxygenase
?NC_008463 270615 = 520 (1.130)intergenic (‑55/‑123) pobA/PA14_RS01225 4‑hydroxybenzoate 3‑monooxygenase/helix‑turn‑helix domain‑containing protein
* ? NC_008463 = 5422623548 (1.080)10 (0.020) 9/278 36.5 2.0% coding (919/1440 nt) oprJ multidrug efflux transporter outer membrane subunit OprJ
?NC_008463 5423053 = 480 (1.020)coding (489/1440 nt) oprJ multidrug efflux transporter outer membrane subunit OprJ
* ? NC_008463 = 992802565 (1.110)9 (0.020) 9/278 36.5 1.6% coding (445/480 nt) nusB transcription antitermination factor NusB
?NC_008463 993104 = 615 (1.300)coding (249/969 nt) thiL thiamine‑phosphate kinase
* ? NC_008463 = 3304506467 (0.920)8 (0.020) 8/268 36.7 1.9% coding (1023/1449 nt) PA14_RS15105 aldehyde dehydrogenase family protein
?NC_008463 3304825 = 405 (0.890)coding (704/1449 nt) PA14_RS15105 aldehyde dehydrogenase family protein
* ? NC_008463 = 44749597 (1.170)8 (0.020) 8/268 36.7 1.8% coding (1816/10332 nt) tpsA2 two‑partner secretion system putative hemagglutinin TpsA2
?NC_008463 2855205 = 347 (0.760)coding (6638/15639 nt) PA14_RS13370 hemagglutinin repeat‑containing protein