breseq  version 0.38.1  
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Read File Information

read filereadsbasespassed filtersaveragelongestmapped
errorsUnknown_CH607-002R0002_good_1.fq3,838,090572,159,134100.0%149.1 bases150 bases96.8%
total3,838,090572,159,134100.0%149.1 bases150 bases96.8%

Reference Sequence Information

seq idlengthfit meanfit relative_variance% mapped readsdescription
coveragedistributionXen14_1242,55667.42.53.1%na
coveragedistributionXen14_2185,21063.72.32.1%na
coveragedistributionXen14_3166,311133.83.04.0%na
coveragedistributionXen14_4151,52490.42.72.5%na
coveragedistributionXen14_5149,06295.52.92.6%na
coveragedistributionXen14_6146,760117.03.43.1%na
coveragedistributionXen14_7136,51471.52.61.8%na
coveragedistributionXen14_8132,518147.33.73.5%na
coveragedistributionXen14_9129,936148.63.73.5%na
coveragedistributionXen14_10107,217160.33.23.2%na
coveragedistributionXen14_11106,943113.22.72.2%na
coveragedistributionXen14_12100,89875.62.61.4%na
coveragedistributionXen14_1398,31886.42.51.5%na
coveragedistributionXen14_1489,144165.23.52.7%na
coveragedistributionXen14_1583,86958.72.50.9%na
coveragedistributionXen14_1681,941100.93.01.5%na
coveragedistributionXen14_1780,390159.24.02.4%na
coveragedistributionXen14_1878,780131.83.51.9%na
coveragedistributionXen14_1977,008137.63.51.9%na
coveragedistributionXen14_2075,56379.92.41.2%na
coveragedistributionXen14_2174,83760.62.20.8%na
coveragedistributionXen14_2274,457164.43.72.2%na
coveragedistributionXen14_2373,975110.42.91.5%na
coveragedistributionXen14_2473,66884.32.51.1%na
coveragedistributionXen14_2573,24056.72.70.8%na
coveragedistributionXen14_2673,210136.73.11.8%na
coveragedistributionXen14_2771,26084.73.31.2%na
coveragedistributionXen14_2864,05458.62.10.7%na
coveragedistributionXen14_2962,12354.82.30.6%na
coveragedistributionXen14_3061,66187.22.41.0%na
coveragedistributionXen14_3158,819176.44.31.9%na
coveragedistributionXen14_3255,39156.82.30.6%na
coveragedistributionXen14_3353,24178.72.40.8%na
coveragedistributionXen14_3453,084125.53.41.2%na
coveragedistributionXen14_3552,845154.34.01.5%na
coveragedistributionXen14_3651,934107.02.81.0%na
coveragedistributionXen14_3751,42488.12.70.8%na
coveragedistributionXen14_3848,83698.03.31.0%na
coveragedistributionXen14_3948,48569.64.30.6%na
coveragedistributionXen14_4046,79154.72.40.5%na
coveragedistributionXen14_4141,431100.43.10.7%na
coveragedistributionXen14_4241,403152.34.01.2%na
coveragedistributionXen14_4341,29689.13.40.7%na
coveragedistributionXen14_4440,957120.03.10.9%na
coveragedistributionXen14_4540,401155.83.21.1%na
coveragedistributionXen14_4640,351101.72.40.7%na
coveragedistributionXen14_4736,274141.53.20.9%na
coveragedistributionXen14_4834,806139.13.20.9%na
coveragedistributionXen14_4934,102151.72.90.9%na
coveragedistributionXen14_5033,369113.72.50.7%na
coveragedistributionXen14_5132,311159.63.90.9%na
coveragedistributionXen14_5231,365136.74.00.8%na
coveragedistributionXen14_5330,293145.23.10.8%na
coveragedistributionXen14_5428,579163.15.10.8%na
coveragedistributionXen14_5526,71587.32.80.4%na
coveragedistributionXen14_5626,428144.52.90.7%na
coveragedistributionXen14_5726,24281.12.70.4%na
coveragedistributionXen14_5824,572130.63.30.6%na
coveragedistributionXen14_5924,177105.63.20.5%na
coveragedistributionXen14_6023,18778.32.20.3%na
coveragedistributionXen14_6121,383117.83.00.5%na
coveragedistributionXen14_6220,21059.71.90.2%na
coveragedistributionXen14_6319,884100.62.40.4%na
coveragedistributionXen14_6418,26962.22.50.2%na
coveragedistributionXen14_6517,85680.12.30.3%na
coveragedistributionXen14_6617,785107.52.90.3%na
coveragedistributionXen14_6717,70491.62.30.3%na
coveragedistributionXen14_6816,382162.23.20.5%na
coveragedistributionXen14_6914,62597.84.50.3%na
coveragedistributionXen14_7014,40254.31.80.1%na
coveragedistributionXen14_7114,30484.64.80.2%na
coveragedistributionXen14_7213,687116.83.60.3%na
coveragedistributionXen14_7313,160131.74.20.3%na
coveragedistributionXen14_7413,069171.43.90.4%na
coveragedistributionXen14_7513,05294.02.60.2%na
coveragedistributionXen14_7612,93198.92.10.2%na
coveragedistributionXen14_7712,66988.23.00.2%na
coveragedistributionXen14_7812,559175.33.50.4%na
coveragedistributionXen14_7912,553134.72.40.3%na
coveragedistributionXen14_8012,035129.04.20.3%na
coveragedistributionXen14_8112,02081.23.20.2%na
coveragedistributionXen14_8211,983100.32.60.2%na
coveragedistributionXen14_8311,35696.94.20.2%na
coveragedistributionXen14_8411,35671.63.90.2%na
coveragedistributionXen14_8511,197125.22.70.3%na
coveragedistributionXen14_8610,95383.32.50.2%na
coveragedistributionXen14_8710,348227.54.00.4%na
coveragedistributionXen14_8810,25967.82.00.1%na
coveragedistributionXen14_8910,105136.33.40.3%na
coveragedistributionXen14_909,218113.82.20.2%na
coveragedistributionXen14_918,91391.43.90.1%na
coveragedistributionXen14_928,23263.22.30.1%na
coveragedistributionXen14_937,77764.82.20.1%na
coveragedistributionXen14_947,58462.12.00.1%na
coveragedistributionXen14_957,34260.52.20.1%na
coveragedistributionXen14_966,91671.12.40.1%na
coveragedistributionXen14_976,746158.64.10.2%na
coveragedistributionXen14_986,594143.12.30.2%na
coveragedistributionXen14_996,04567.11.70.1%na
coveragedistributionXen14_1005,846129.33.00.1%na
coveragedistributionXen14_1015,615910.138.00.9%na
coveragedistributionXen14_1025,55259.33.40.1%na
coveragedistributionXen14_1035,30775.82.60.1%na
coveragedistributionXen14_1045,096104.43.30.1%na
coveragedistributionXen14_1055,06193.11.00.1%na
coveragedistributionXen14_1064,99597.12.10.1%na
coveragedistributionXen14_1074,93265.71.50.1%na
coveragedistributionXen14_1084,81693.02.50.1%na
coveragedistributionXen14_1094,61689.73.60.1%na
coveragedistributionXen14_1104,57849.51.70.0%na
coveragedistributionXen14_1114,562729.419.40.6%na
coveragedistributionXen14_1124,48681.03.70.1%na
coveragedistributionXen14_1134,275160.54.00.1%na
coveragedistributionXen14_1144,03961.63.90.0%na
coveragedistributionXen14_1153,633128.83.30.1%na
coveragedistributionXen14_1163,428131.92.60.1%na
coveragedistributionXen14_1173,42779.82.60.0%na
coveragedistributionXen14_1183,35559.82.70.0%na
coveragedistributionXen14_1193,20890.62.90.0%na
coveragedistributionXen14_1203,12354.03.20.0%na
coveragedistributionXen14_1213,07686.31.70.0%na
coveragedistributionXen14_1222,908106.84.20.1%na
coveragedistributionXen14_1232,7141019.87.30.5%na
coveragedistributionXen14_1242,576127.14.90.1%na
coveragedistributionXen14_1252,519129.12.10.1%na
coveragedistributionXen14_1262,515125.73.10.1%na
coveragedistributionXen14_1272,1661058.014.00.4%na
coveragedistributionXen14_1282,14380.92.80.0%na
coveragedistributionXen14_1292,049125.63.80.0%na
coveragedistributionXen14_1302,000119.24.90.0%na
coveragedistributionXen14_1311,92251.02.40.0%na
coveragedistributionXen14_1321,806135.94.20.1%na
coveragedistributionXen14_1331,721246.75.00.1%na
coveragedistributionXen14_1341,65792.05.10.0%na
coveragedistributionXen14_1351,633150.62.30.0%na
coveragedistributionXen14_1361,567121.52.30.0%na
coveragedistributionXen14_1371,544217.82.20.1%na
coveragedistributionXen14_1381,361119.81.70.0%na
coveragedistributionXen14_1391,349153.61.70.0%na
coveragedistributionXen14_1401,298947.38.30.2%na
coveragedistributionXen14_1411,2922027.646.70.4%na
coveragedistributionXen14_1421,28157.02.80.0%na
coveragedistributionXen14_1431,1071196.43.80.1%na
coveragedistributionXen14_1441,091692.011.70.1%na
coveragedistributionXen14_1451,062696.92.30.1%na
coveragedistributionXen14_1461,046113.62.20.0%na
coveragedistributionXen14_1471,007187.72.50.0%na
coveragedistributionXen14_14898948.02.10.0%na
coveragedistributionXen14_149980112.51.50.0%na
coveragedistributionXen14_150945158.31.00.0%na
coveragedistributionXen14_151938163.91.40.0%na
coveragedistributionXen14_152934128.623.60.0%na
coveragedistributionXen14_1539292756.3134.80.4%na
coveragedistributionXen14_154920218.51.10.0%na
coveragedistributionXen14_15590073.311.00.0%na
coveragedistributionXen14_156862194.84.50.0%na
coveragedistributionXen14_157857209.17.70.0%na
coveragedistributionXen14_158799109.58.70.0%na
coveragedistributionXen14_159798137.86.90.0%na
coveragedistributionXen14_160780196.524.40.0%na
coveragedistributionXen14_161770119.81.10.0%na
coveragedistributionXen14_162765181.62.80.0%na
coveragedistributionXen14_163753263.12.40.0%na
coveragedistributionXen14_16474685.05.60.0%na
coveragedistributionXen14_165691150.85.10.0%na
coveragedistributionXen14_16668511009.87409.30.2%na
coveragedistributionXen14_167666886.132.80.1%na
coveragedistributionXen14_1686292900.8395.30.1%na
coveragedistributionXen14_169626197.92.00.0%na
coveragedistributionXen14_170619100.61.70.0%na
coveragedistributionXen14_171569133.61.00.0%na
coveragedistributionXen14_17256940.11.30.0%na
coveragedistributionXen14_173564200.38.40.0%na
coveragedistributionXen14_174564130.61.20.0%na
coveragedistributionXen14_17556099.31.30.0%na
coveragedistributionXen14_176552[163.7][28.7]0.0%na
coveragedistributionXen14_177549173.611.90.0%na
coveragedistributionXen14_178541123.79.30.0%na
coveragedistributionXen14_179536146.41.00.0%na
coveragedistributionXen14_18052699.12.20.0%na
coveragedistributionXen14_181511[140.8][23.1]0.0%na
coveragedistributionXen14_182478167.612.40.0%na
coveragedistributionXen14_183477136.73.60.0%na
coveragedistributionXen14_1844731351.056.60.1%na
coveragedistributionXen14_185471[366.7][50.9]0.0%na
coveragedistributionXen14_186465132.518.30.0%na
coveragedistributionXen14_187464372.818.40.0%na
coveragedistributionXen14_188458[235.6][25.8]0.0%na
coveragedistributionXen14_189453194.18.50.0%na
coveragedistributionXen14_19044870.013.30.0%na
coveragedistributionXen14_191447[322.6][159.4]0.0%na
coveragedistributionXen14_192445128.21.90.0%na
coveragedistributionXen14_193444343.111.90.0%na
coveragedistributionXen14_19444281.81.00.0%na
coveragedistributionXen14_19544173.62.10.0%na
coveragedistributionXen14_196436137.71.00.0%na
coveragedistributionXen14_197433101.02.30.0%na
coveragedistributionXen14_19843196.31.40.0%na
coveragedistributionXen14_199428156.12.80.0%na
coveragedistributionXen14_2004231099.029.50.1%na
coveragedistributionXen14_201421436.0137.70.0%na
coveragedistributionXen14_2024211597.4434.60.0%na
coveragedistributionXen14_203416106.98.50.0%na
coveragedistributionXen14_204415[103.4][22.8]0.0%na
coveragedistributionXen14_2054041004.235.70.0%na
coveragedistributionXen14_206403[945.6][106.5]0.1%na
total4,963,538100.0%

fit relative_variance is the ratio of the variance to the mean for the negative binomial fit. It is =1 for Poisson and >1 for over-dispersed data.

Fit failed Negative binomial fit failed for this reference sequence. It may have an unusual coverage depth distribution. JC and MC predictions may be less accurate.

New Junction Evidence

Junction Candidates Tested

optionlimitactual
Number of alignment pairs examined for constructing junction candidates≤ 10000076784
Coverage evenness (position-hash) score of junction candidates≥ 2≥ 2
Test this many junction candidates (n). May be smaller if not enough passed the coverage evenness threshold100 ≤ n ≤ 5000418
Total length of all junction candidates (factor times the reference genome length)≤ 0.10.024

Junction Skew Score Calculation

reference sequencepr(no read start)
Xen14_10.84506
Xen14_20.85381
Xen14_30.72748
Xen14_40.80065
Xen14_50.79186
Xen14_60.75465
Xen14_70.83775
Xen14_80.70884
Xen14_90.70581
Xen14_100.68437
Xen14_110.76254
Xen14_120.82849
Xen14_130.80917
Xen14_140.68168
Xen14_150.86365
Xen14_160.78360
Xen14_170.68570
Xen14_180.72897
Xen14_190.72437
Xen14_200.82038
Xen14_210.86050
Xen14_220.68134
Xen14_230.76644
Xen14_240.81378
Xen14_250.86917
Xen14_260.72113
Xen14_270.80708
Xen14_280.86216
Xen14_290.87292
Xen14_300.80879
Xen14_310.66912
Xen14_320.86714
Xen14_330.82435
Xen14_340.74312
Xen14_350.69656
Xen14_360.77526
Xen14_370.80677
Xen14_380.78478
Xen14_390.84079
Xen14_400.87248
Xen14_410.78718
Xen14_420.69460
Xen14_430.80585
Xen14_440.75082
Xen14_450.69393
Xen14_460.78043
Xen14_470.71897
Xen14_480.71814
Xen14_490.69745
Xen14_500.75969
Xen14_510.69226
Xen14_520.72705
Xen14_530.70544
Xen14_540.68647
Xen14_550.80502
Xen14_560.71184
Xen14_570.81832
Xen14_580.73028
Xen14_590.77106
Xen14_600.82211
Xen14_610.75642
Xen14_620.85997
Xen14_630.78108
Xen14_640.85108
Xen14_650.82185
Xen14_660.77231
Xen14_670.80437
Xen14_680.68044
Xen14_690.79135
Xen14_700.87340
Xen14_710.81708
Xen14_720.75597
Xen14_730.73970
Xen14_740.67652
Xen14_750.79842
Xen14_760.78629
Xen14_770.80488
Xen14_780.66462
Xen14_790.72763
Xen14_800.73639
Xen14_810.82899
Xen14_820.78503
Xen14_830.79372
Xen14_840.83471
Xen14_850.73230
Xen14_860.80510
Xen14_870.60084
Xen14_880.83439
Xen14_890.72241
Xen14_900.76508
Xen14_910.80809
Xen14_920.85836
Xen14_930.85631
Xen14_940.85951
Xen14_950.86414
Xen14_960.77345
Xen14_970.68255
Xen14_980.70890
Xen14_990.85054
Xen14_1000.74359
Xen14_1010.23393
Xen14_1020.86635
Xen14_1030.82900
Xen14_1040.79003
Xen14_1050.75598
Xen14_1060.79920
Xen14_1070.85148
Xen14_1080.79662
Xen14_1090.81161
Xen14_1100.88041
Xen14_1110.30173
Xen14_1120.81621
Xen14_1130.69626
Xen14_1140.85999
Xen14_1150.75062
Xen14_1160.74460
Xen14_1170.83149
Xen14_1180.86587
Xen14_1190.81016
Xen14_1200.87096
Xen14_1210.82152
Xen14_1220.77957
Xen14_1230.22476
Xen14_1240.75699
Xen14_1250.74931
Xen14_1260.70636
Xen14_1270.22922
Xen14_1280.83179
Xen14_1290.75939
Xen14_1300.77025
Xen14_1310.86157
Xen14_1320.68810
Xen14_1330.60895
Xen14_1340.80396
Xen14_1350.70667
Xen14_1360.77186
Xen14_1370.63245
Xen14_1380.77443
Xen14_1390.72313
Xen14_1400.28005
Xen14_1410.16602
Xen14_1420.86127
Xen14_1430.50469
Xen14_1440.37672
Xen14_1450.34934
Xen14_1460.78824
Xen14_1470.64548
Xen14_1480.88423
Xen14_1490.80000
Xen14_1500.74392
Xen14_1510.72388
Xen14_1520.81317
Xen14_1530.18030
Xen14_1540.68587
Xen14_1550.89111
Xen14_1560.70584
Xen14_1570.67678
Xen14_1580.81164
Xen14_1590.77130
Xen14_1600.75128
Xen14_1610.78442
Xen14_1620.73856
Xen14_1630.64011
Xen14_1640.84853
Xen14_1650.74891
Xen14_1660.41826
Xen14_1670.37162
Xen14_1680.29491
Xen14_1690.61994
Xen14_1700.81664
Xen14_1710.78032
Xen14_1720.91564
Xen14_1730.70922
Xen14_1740.78191
Xen14_1750.83393
Xen14_176NA
Xen14_1770.75137
Xen14_1780.80314
Xen14_1790.77892
Xen14_1800.83650
Xen14_181NA
Xen14_1820.76829
Xen14_1830.73009
Xen14_1840.39958
Xen14_185NA
Xen14_1860.83656
Xen14_1870.65409
Xen14_188NA
Xen14_1890.74172
Xen14_1900.87946
Xen14_191NA
Xen14_1920.81798
Xen14_1930.66216
Xen14_1940.86312
Xen14_1950.88869
Xen14_1960.85570
Xen14_1970.77835
Xen14_1980.84687
Xen14_1990.82092
Xen14_2000.20364
Xen14_2010.85273
Xen14_2020.55107
Xen14_2030.85457
Xen14_204NA
Xen14_2050.42946
Xen14_206NA

pr(no read start) is the probability that there will not be an aligned read whose first base matches a given position on a given strand.

Final Junction Predictions

optionvalue
Coverage evenness (position-hash) score of predicted junctions must be≥ 3
Minimum probablilty assigned that no mapped read will start at a given position and strand for junction prediction0.1
Junction allow suboptimal matchesFALSE
Skew score of predicted junction (−log10 probability of unusual coverage evenness) must be≤ 3
Number of bases that at least one read must overlap each uniquely aligned side of a predicted junction≥ 1

Read Alignment Evidence

optionvalue
ModeConsensus/Mixed Base
Ploidy1 (haploid)
Consensus mutation E-value cutoff10
Consensus frequency cutoff0.8
Consensus minimum variant coverage each strandOFF
Consensus minimum total coverage each strandOFF
Consensus minimum variant coverageOFF
Consensus minimum total coverageOFF
Polymorphism E-value cutoff10
Polymorphism frequency cutoff0.2
Polymorphism minimum variant coverage each strandOFF
Polymorphism minimum total coverage each strandOFF
Polymorphism minimum variant coverageOFF
Polymorphism minimum total coverageOFF
Polymorphism bias cutoffOFF
Predict indel polymorphismsYES
Skip indel polymorphisms in homopolymers runs ofOFF
Skip base substitutions when they create a homopolymer flanked on each side byOFF

Software Versions

programversion
bowtie2s.2023.01/rocky8/x86_64/amd/zen3/software/Bowtie2/2.4.5-GCC-11.3.0/bin/bowtie2-align-s
R4.2.1

Execution Times

stepstartendelapsed
Read and reference sequence file input21:16:59 30 Oct 202421:18:59 30 Oct 20242 minutes 0 seconds
Read alignment to reference genome21:19:01 30 Oct 202421:20:54 30 Oct 20241 minute 53 seconds
Preprocessing alignments for candidate junction identification21:20:54 30 Oct 202421:21:37 30 Oct 202443 seconds
Preliminary analysis of coverage distribution21:21:38 30 Oct 202421:32:30 30 Oct 202410 minutes 52 seconds
Identifying junction candidates21:32:30 30 Oct 202421:33:16 30 Oct 202446 seconds
Re-alignment to junction candidates21:33:17 30 Oct 202421:33:52 30 Oct 202435 seconds
Resolving best read alignments21:33:52 30 Oct 202421:35:13 30 Oct 20241 minute 21 seconds
Creating BAM files21:35:13 30 Oct 202421:36:07 30 Oct 202454 seconds
Tabulating error counts21:36:07 30 Oct 202421:38:17 30 Oct 20242 minutes 10 seconds
Re-calibrating base error rates21:38:17 30 Oct 202421:43:08 30 Oct 20244 minutes 51 seconds
Examining read alignment evidence21:43:08 30 Oct 202421:51:14 30 Oct 20248 minutes 6 seconds
Polymorphism statistics21:51:14 30 Oct 202421:51:17 30 Oct 20243 seconds
Output21:51:18 30 Oct 202421:51:45 30 Oct 202427 seconds
Output :: Mutation Prediction21:51:18 30 Oct 202421:51:20 30 Oct 20242 seconds
Output :: Mutation Annotation21:51:20 30 Oct 202421:51:20 30 Oct 20240 seconds
Total 34 minutes 43 seconds