breseq  version 0.28.1  revision 52323b585d6f
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 252 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,559,2380TC56.4% 207.7 / inf 536intergenic (‑395/‑429)ykgL/ecpRhypothetical protein/helix‑turn‑helix transcriptional regulator
*DH13,184,6490AG55.1% 36.8 / inf 334noncoding (35/75 nt)ECDH1_RS15800tRNA‑TERM
*DH14,165,3220AG54.9% 191.7 / inf 576pseudogene (50/129 nt)dcuRglutamate decarboxylase
*DH1613,2670AG54.0% 101.9 / inf 423intergenic (+128/+895)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
*DH13,557,4120TC53.8% 110.2 / inf 463intergenic (‑814/‑301)ECDH1_RS17595/ykgRattachment protein/hypothetical protein
*DH14,026,7830AG53.5% 133.3 / inf 560C269R (TGT→CGT) lptFLPS export ABC transporter permease LptF
*DH11,060,4750AG53.4% 120.1 / inf 471D161G (GAC→GGC) recADNA recombination/repair protein RecA
*DH14,081,1220AG53.4% 113.1 / inf 502D189G (GAC→GGC) qorBNAD(P)‑dependent oxidoreductase
*DH14,111,5710AG53.2% 102.8 / inf 496C325R (TGT→CGT) hflKprotease modulator HflK
*DH13,542,9560TC53.2% 70.6 / inf 436S245S (TCT→TCCbetBbetaine‑aldehyde dehydrogenase
*DH14,245,0210AG53.2% 90.0 / inf 456C12C (TGT→TGCaphAacid phosphatase AphA
*DH11,831,5490AC53.0% 217.3 / inf 456V523V (GTT→GTGyeeJinverse autotransporter adhesin‑like protein YeeJ
*DH1792,7950TC52.9% 88.6 / inf 525K208E (AAG→GAG) endAdeoxyribonuclease I
*DH12,476,5850TC52.8% 63.2 / inf 414pseudogene (1363/1446 nt)abgApeptidase M20
*DH1762,7300GA52.6% 111.0 / inf 481V83I (GTT→ATT) glcAglycolate permease GlcA
*DH13,553,1670AG52.6% 91.5 / inf 494T178A (ACG→GCG) rclCprotein RclC
*DH13,874,7890TC52.6% 77.7 / inf 545L21L (TTA→TTGcreDprotein CreD
*DH12,732,5430AG52.4% 72.0 / inf 495S135P (TCG→CCG) yceDhypothetical protein
*DH13,869,4740TC52.4% 60.8 / inf 478K523E (AAA→GAA) thrAbifunctional aspartate kinase/homoserine dehydrogenase I
*DH12,999,5880AT52.4% 92.2 / inf 380F269I (TTC→ATC) dacCD‑alanyl‑D‑alanine carboxypeptidase

Marginal new junction evidence (lowest skew 10 of 67 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? DH1 1586746 =NA (NA)382 (0.800) 162/292 3.1 57.7% noncoding (1/1198 nt) repeat_region repeat region
?DH1 = 1971030 280 (0.590)coding (144/144 nt) yobF DUF2527 domain‑containing protein
* ? DH1 = 300335NA (NA)279 (0.570) 167/300 3.1 50.4% noncoding (772/773 nt) repeat_region repeat region
?DH1 490607 = 274 (0.560)coding (3102/3801 nt) yhdP DUF3971 domain‑containing protein
* ? DH1 = 3575873252 (0.520)249 (0.520) 160/292 3.3 50.6% noncoding (3/1201 nt) repeat_region repeat region
?DH1 3577073 = 241 (0.510)intergenic (+151/‑282) insH1/yagM transposase/hypothetical protein
* ? DH1 = 2201911194 (0.400)249 (0.510) 163/300 3.5 66.9% coding (1337/1533 nt) pntA NAD(P) transhydrogenase subunit alpha
?DH1 2204947 = 52 (0.110)intergenic (‑382/‑30) tqsA/mdtJ AI‑2 transporter TqsA/spermidine export protein MdtJ
* ? DH1 3203897 =317 (0.650)231 (0.500) 153/284 3.6 44.4% coding (512/708 nt) ybeR hypothetical protein
?DH1 = 3203926 279 (0.600)coding (483/708 nt) ybeR hypothetical protein
* ? DH1 = 4079176279 (0.570)228 (0.490) 152/284 3.7 46.1% noncoding (1/771 nt) repeat_region repeat region
?DH1 4079953 = 269 (0.580)intergenic (+24/+135) insB1/ytfH transposase/HxlR family transcriptional regulator
* ? DH1 = 1702500261 (0.530)236 (0.500) 157/292 3.7 48.0% noncoding (5/1205 nt) repeat_region repeat region
?DH1 1703700 = 257 (0.540)noncoding (1205/1205 nt) repeat_region repeat region
* ? DH1 = 2678609234 (0.480)246 (0.520) 156/292 3.8 49.8% noncoding (1/1259 nt) repeat_region repeat region
?DH1 2679871 = 269 (0.570)intergenic (‑69/+280) insE1/lit transposase/bacteriophage T4 late gene expression‑blocking protein
* ? DH1 = 2788747252 (0.520)221 (0.470) 153/292 4.2 49.2% noncoding (5/1206 nt) repeat_region repeat region
?DH1 2789947 = 211 (0.440)noncoding (1205/1206 nt) repeat_region repeat region
* ? DH1 = 1967355225 (0.460)208 (0.450) 148/286 4.4 48.2% noncoding (6/1206 nt) repeat_region repeat region
?DH1 1968555 = 232 (0.500)noncoding (1206/1206 nt) repeat_region repeat region