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

Read File Information

read filereadsbasespassed filtersaveragelongestmapped
errorsUnknown_CH607-002R0004_good_1.fq4,263,449636,038,093100.0%149.2 bases150 bases96.9%
total4,263,449636,038,093100.0%149.2 bases150 bases96.9%

Reference Sequence Information

seq idlengthfit meanfit relative_variance% mapped readsdescription
coveragedistributionXen14_1242,55673.22.73.0%na
coveragedistributionXen14_2185,21070.42.62.1%na
coveragedistributionXen14_3166,311148.73.94.0%na
coveragedistributionXen14_4151,52499.53.52.4%na
coveragedistributionXen14_5149,062103.83.72.5%na
coveragedistributionXen14_6146,760130.83.83.1%na
coveragedistributionXen14_7136,51476.72.61.7%na
coveragedistributionXen14_8132,518159.94.23.4%na
coveragedistributionXen14_9129,936161.24.93.4%na
coveragedistributionXen14_10107,217187.14.43.3%na
coveragedistributionXen14_11106,943125.63.42.2%na
coveragedistributionXen14_12100,89881.32.81.4%na
coveragedistributionXen14_1398,31893.93.21.5%na
coveragedistributionXen14_1489,144186.95.32.7%na
coveragedistributionXen14_1583,86965.22.40.9%na
coveragedistributionXen14_1681,941110.83.01.5%na
coveragedistributionXen14_1780,390185.35.12.5%na
coveragedistributionXen14_1878,780144.53.91.9%na
coveragedistributionXen14_1977,008149.34.71.9%na
coveragedistributionXen14_2075,56385.42.51.2%na
coveragedistributionXen14_2174,83767.22.30.8%na
coveragedistributionXen14_2274,457193.94.92.3%na
coveragedistributionXen14_2373,975125.74.51.5%na
coveragedistributionXen14_2473,66889.23.21.1%na
coveragedistributionXen14_2573,24062.62.10.8%na
coveragedistributionXen14_2673,210150.13.51.8%na
coveragedistributionXen14_2771,260101.74.71.2%na
coveragedistributionXen14_2864,05464.12.30.7%na
coveragedistributionXen14_2962,12358.52.30.6%na
coveragedistributionXen14_3061,66196.23.01.0%na
coveragedistributionXen14_3158,819210.25.92.0%na
coveragedistributionXen14_3255,39161.02.40.6%na
coveragedistributionXen14_3353,24185.62.80.7%na
coveragedistributionXen14_3453,084135.24.51.2%na
coveragedistributionXen14_3552,845172.15.01.5%na
coveragedistributionXen14_3651,934115.14.51.0%na
coveragedistributionXen14_3751,42496.03.30.8%na
coveragedistributionXen14_3848,836110.23.31.0%na
coveragedistributionXen14_3948,48573.36.60.6%na
coveragedistributionXen14_4046,79161.42.50.5%na
coveragedistributionXen14_4141,431110.72.80.7%na
coveragedistributionXen14_4241,403172.45.21.2%na
coveragedistributionXen14_4341,296107.14.60.7%na
coveragedistributionXen14_4440,957131.24.20.9%na
coveragedistributionXen14_4540,401173.84.71.1%na
coveragedistributionXen14_4640,351112.73.70.7%na
coveragedistributionXen14_4736,274153.94.80.9%na
coveragedistributionXen14_4834,806156.83.50.9%na
coveragedistributionXen14_4934,102172.64.81.0%na
coveragedistributionXen14_5033,369127.93.40.7%na
coveragedistributionXen14_5132,311172.14.80.9%na
coveragedistributionXen14_5231,365152.27.10.8%na
coveragedistributionXen14_5330,293161.84.80.8%na
coveragedistributionXen14_5428,579189.25.80.9%na
coveragedistributionXen14_5526,71593.02.80.4%na
coveragedistributionXen14_5626,428167.23.10.7%na
coveragedistributionXen14_5726,24284.62.80.4%na
coveragedistributionXen14_5824,572134.53.10.5%na
coveragedistributionXen14_5924,177123.14.60.5%na
coveragedistributionXen14_6023,18784.52.30.3%na
coveragedistributionXen14_6121,383136.52.50.5%na
coveragedistributionXen14_6220,21063.12.50.2%na
coveragedistributionXen14_6319,884110.32.60.4%na
coveragedistributionXen14_6418,26971.82.20.2%na
coveragedistributionXen14_6517,85683.93.00.2%na
coveragedistributionXen14_6617,785115.94.40.3%na
coveragedistributionXen14_6717,70499.92.90.3%na
coveragedistributionXen14_6816,382183.02.40.5%na
coveragedistributionXen14_6914,625118.56.30.3%na
coveragedistributionXen14_7014,40258.32.40.1%na
coveragedistributionXen14_7114,30497.83.90.2%na
coveragedistributionXen14_7213,687132.63.30.3%na
coveragedistributionXen14_7313,160157.84.30.3%na
coveragedistributionXen14_7413,069193.95.40.4%na
coveragedistributionXen14_7513,052124.92.50.3%na
coveragedistributionXen14_7612,931106.23.60.2%na
coveragedistributionXen14_7712,669101.32.70.2%na
coveragedistributionXen14_7812,559205.44.40.4%na
coveragedistributionXen14_7912,553146.62.60.3%na
coveragedistributionXen14_8012,035139.15.80.3%na
coveragedistributionXen14_8112,02091.73.10.2%na
coveragedistributionXen14_8211,983123.62.90.2%na
coveragedistributionXen14_8311,356125.45.70.2%na
coveragedistributionXen14_8411,35690.77.10.2%na
coveragedistributionXen14_8511,197136.42.50.3%na
coveragedistributionXen14_8610,95389.83.40.2%na
coveragedistributionXen14_8710,348241.35.40.4%na
coveragedistributionXen14_8810,25969.12.00.1%na
coveragedistributionXen14_8910,105168.93.90.3%na
coveragedistributionXen14_909,218132.83.80.2%na
coveragedistributionXen14_918,913110.23.20.2%na
coveragedistributionXen14_928,23269.82.00.1%na
coveragedistributionXen14_937,77772.23.50.1%na
coveragedistributionXen14_947,58474.62.30.1%na
coveragedistributionXen14_957,34268.82.10.1%na
coveragedistributionXen14_966,916165.38.40.1%na
coveragedistributionXen14_976,746187.05.60.2%na
coveragedistributionXen14_986,594162.66.00.2%na
coveragedistributionXen14_996,04576.22.10.1%na
coveragedistributionXen14_1005,846145.83.30.1%na
coveragedistributionXen14_1015,6151132.068.21.0%na
coveragedistributionXen14_1025,55272.61.90.1%na
coveragedistributionXen14_1035,30777.51.90.1%na
coveragedistributionXen14_1045,096107.13.20.1%na
coveragedistributionXen14_1055,061108.92.50.1%na
coveragedistributionXen14_1064,995121.53.50.1%na
coveragedistributionXen14_1074,93269.02.50.1%na
coveragedistributionXen14_1084,816112.23.90.1%na
coveragedistributionXen14_1094,61691.44.90.1%na
coveragedistributionXen14_1104,57856.82.20.0%na
coveragedistributionXen14_1114,562793.115.00.6%na
coveragedistributionXen14_1124,48684.42.00.1%na
coveragedistributionXen14_1134,275193.05.40.1%na
coveragedistributionXen14_1144,03968.94.00.0%na
coveragedistributionXen14_1153,633157.93.30.1%na
coveragedistributionXen14_1163,428145.63.20.1%na
coveragedistributionXen14_1173,42793.14.90.1%na
coveragedistributionXen14_1183,35560.93.40.0%na
coveragedistributionXen14_1193,20896.32.50.0%na
coveragedistributionXen14_1203,12366.93.10.0%na
coveragedistributionXen14_1213,07693.52.50.0%na
coveragedistributionXen14_1222,908117.14.50.1%na
coveragedistributionXen14_1232,7141219.010.80.5%na
coveragedistributionXen14_1242,576136.93.00.1%na
coveragedistributionXen14_1252,519152.23.60.1%na
coveragedistributionXen14_1262,515141.21.80.1%na
coveragedistributionXen14_1272,1661129.216.20.4%na
coveragedistributionXen14_1282,14387.52.90.0%na
coveragedistributionXen14_1292,049142.46.10.0%na
coveragedistributionXen14_1302,000130.82.10.0%na
coveragedistributionXen14_1311,92256.71.00.0%na
coveragedistributionXen14_1321,806161.42.90.1%na
coveragedistributionXen14_1331,721259.91.80.1%na
coveragedistributionXen14_1341,65779.45.70.0%na
coveragedistributionXen14_1351,633155.72.40.0%na
coveragedistributionXen14_1361,567123.01.00.0%na
coveragedistributionXen14_1371,544196.41.90.0%na
coveragedistributionXen14_1381,361136.52.50.0%na
coveragedistributionXen14_1391,349177.42.60.0%na
coveragedistributionXen14_1401,2981181.83.50.2%na
coveragedistributionXen14_1411,2922164.626.30.4%na
coveragedistributionXen14_1421,28157.81.00.0%na
coveragedistributionXen14_1431,1071278.37.40.1%na
coveragedistributionXen14_1441,091792.417.50.1%na
coveragedistributionXen14_1451,062823.42.80.1%na
coveragedistributionXen14_1461,046143.91.00.0%na
coveragedistributionXen14_1471,007191.62.00.0%na
coveragedistributionXen14_14898979.45.20.0%na
coveragedistributionXen14_149980133.54.50.0%na
coveragedistributionXen14_150945170.34.70.0%na
coveragedistributionXen14_151938169.93.80.0%na
coveragedistributionXen14_152934142.07.00.0%na
coveragedistributionXen14_1539293044.5147.20.4%na
coveragedistributionXen14_154920210.65.90.0%na
coveragedistributionXen14_15590078.01.00.0%na
coveragedistributionXen14_156862220.22.50.0%na
coveragedistributionXen14_157857264.46.80.0%na
coveragedistributionXen14_158799123.42.30.0%na
coveragedistributionXen14_159798169.96.50.0%na
coveragedistributionXen14_160780199.27.30.0%na
coveragedistributionXen14_161770146.83.40.0%na
coveragedistributionXen14_162765185.51.80.0%na
coveragedistributionXen14_163753272.03.60.0%na
coveragedistributionXen14_16474683.21.60.0%na
coveragedistributionXen14_165691176.44.00.0%na
coveragedistributionXen14_1666852983.7751.70.2%na
coveragedistributionXen14_1676661020.527.90.1%na
coveragedistributionXen14_1686294474.2915.80.1%na
coveragedistributionXen14_169626205.31.00.0%na
coveragedistributionXen14_170619128.11.00.0%na
coveragedistributionXen14_171569154.31.80.0%na
coveragedistributionXen14_17256952.02.90.0%na
coveragedistributionXen14_173564331.429.90.0%na
coveragedistributionXen14_174564164.41.00.0%na
coveragedistributionXen14_175560123.51.20.0%na
coveragedistributionXen14_176552274.231.40.0%na
coveragedistributionXen14_177549400.278.50.0%na
coveragedistributionXen14_178541122.92.40.0%na
coveragedistributionXen14_179536144.54.30.0%na
coveragedistributionXen14_180526130.41.00.0%na
coveragedistributionXen14_181511163.415.10.0%na
coveragedistributionXen14_182478206.66.10.0%na
coveragedistributionXen14_183477131.74.10.0%na
coveragedistributionXen14_1844731302.03.30.1%na
coveragedistributionXen14_185471515.12.20.0%na
coveragedistributionXen14_186465122.44.90.0%na
coveragedistributionXen14_187464349.12.40.0%na
coveragedistributionXen14_188458641.769.80.0%na
coveragedistributionXen14_189453210.61.00.0%na
coveragedistributionXen14_190448[85.5][19.3]0.0%na
coveragedistributionXen14_191447[423.1][219.7]0.0%na
coveragedistributionXen14_192445125.51.10.0%na
coveragedistributionXen14_193444330.11.00.0%na
coveragedistributionXen14_194442119.815.30.0%na
coveragedistributionXen14_19544182.013.10.0%na
coveragedistributionXen14_196436154.213.60.0%na
coveragedistributionXen14_197433112.51.10.0%na
coveragedistributionXen14_198431115.71.30.0%na
coveragedistributionXen14_199428162.41.00.0%na
coveragedistributionXen14_2004231893.3132.90.1%na
coveragedistributionXen14_201421[84.6][16.5]0.0%na
coveragedistributionXen14_2024211270.3147.10.0%na
coveragedistributionXen14_203416[81.2][18.1]0.0%na
coveragedistributionXen14_204415[120.1][28.3]0.0%na
coveragedistributionXen14_2054041179.17.90.1%na
coveragedistributionXen14_2064031324.713.70.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≤ 10000084324
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 ≤ 5000460
Total length of all junction candidates (factor times the reference genome length)≤ 0.10.026

Junction Skew Score Calculation

reference sequencepr(no read start)
Xen14_10.82279
Xen14_20.82954
Xen14_30.68674
Xen14_40.77225
Xen14_50.76341
Xen14_60.71529
Xen14_70.81697
Xen14_80.67011
Xen14_90.66823
Xen14_100.63189
Xen14_110.72521
Xen14_120.80671
Xen14_130.78401
Xen14_140.63069
Xen14_150.84110
Xen14_160.75003
Xen14_170.63109
Xen14_180.69318
Xen14_190.68679
Xen14_200.79618
Xen14_210.83514
Xen14_220.62216
Xen14_230.72638
Xen14_240.79197
Xen14_250.84597
Xen14_260.68351
Xen14_270.76372
Xen14_280.84174
Xen14_290.85584
Xen14_300.77833
Xen14_310.60600
Xen14_320.84724
Xen14_330.79935
Xen14_340.71093
Xen14_350.65254
Xen14_360.74118
Xen14_370.78030
Xen14_380.74905
Xen14_390.82498
Xen14_400.84974
Xen14_410.75506
Xen14_420.65128
Xen14_430.76176
Xen14_440.71894
Xen14_450.65203
Xen14_460.74987
Xen14_470.68227
Xen14_480.67500
Xen14_490.65218
Xen14_500.72458
Xen14_510.65012
Xen14_520.68937
Xen14_530.66397
Xen14_540.63380
Xen14_550.78441
Xen14_560.66068
Xen14_570.80025
Xen14_580.70853
Xen14_590.73316
Xen14_600.80088
Xen14_610.70862
Xen14_620.84359
Xen14_630.75101
Xen14_640.82504
Xen14_650.80101
Xen14_660.74487
Xen14_670.77590
Xen14_680.63881
Xen14_690.73831
Xen14_700.85205
Xen14_710.77202
Xen14_720.71696
Xen14_730.68438
Xen14_740.62744
Xen14_750.73413
Xen14_760.75679
Xen14_770.77449
Xen14_780.60658
Xen14_790.68892
Xen14_800.70843
Xen14_810.79359
Xen14_820.72774
Xen14_830.73234
Xen14_840.79933
Xen14_850.69106
Xen14_860.78133
Xen14_870.57335
Xen14_880.81681
Xen14_890.66259
Xen14_900.72169
Xen14_910.76080
Xen14_920.83157
Xen14_930.83188
Xen14_940.82641
Xen14_950.83642
Xen14_960.74971
Xen14_970.63897
Xen14_980.66822
Xen14_990.82060
Xen14_1000.69911
Xen14_1010.18415
Xen14_1020.83177
Xen14_1030.80535
Xen14_1040.76442
Xen14_1050.71784
Xen14_1060.74615
Xen14_1070.83394
Xen14_1080.75311
Xen14_1090.79455
Xen14_1100.85671
Xen14_1110.26600
Xen14_1120.79665
Xen14_1130.63953
Xen14_1140.83276
Xen14_1150.68869
Xen14_1160.69822
Xen14_1170.79268
Xen14_1180.84978
Xen14_1190.78647
Xen14_1200.84662
Xen14_1210.79096
Xen14_1220.74209
Xen14_1230.16949
Xen14_1240.72069
Xen14_1250.70013
Xen14_1260.66620
Xen14_1270.20776
Xen14_1280.81616
Xen14_1290.71645
Xen14_1300.72675
Xen14_1310.84002
Xen14_1320.64780
Xen14_1330.57379
Xen14_1340.78365
Xen14_1350.66622
Xen14_1360.74697
Xen14_1370.62565
Xen14_1380.73292
Xen14_1390.68087
Xen14_1400.21225
Xen14_1410.15441
Xen14_1420.83949
Xen14_1430.46102
Xen14_1440.32493
Xen14_1450.29614
Xen14_1460.74857
Xen14_1470.64052
Xen14_1480.83771
Xen14_1490.74388
Xen14_1500.69312
Xen14_1510.69883
Xen14_1520.74197
Xen14_1530.16416
Xen14_1540.65054
Xen14_1550.83944
Xen14_1560.68175
Xen14_1570.61669
Xen14_1580.76846
Xen14_1590.73120
Xen14_1600.69423
Xen14_1610.74610
Xen14_1620.70980
Xen14_1630.60093
Xen14_1640.84048
Xen14_1650.71274
Xen14_1660.38745
Xen14_1670.31532
Xen14_1680.26630
Xen14_1690.71694
Xen14_1700.77948
Xen14_1710.75659
Xen14_1720.89895
Xen14_1730.65869
Xen14_1740.74557
Xen14_1750.79018
Xen14_1760.71014
Xen14_1770.71038
Xen14_1780.78373
Xen14_1790.76679
Xen14_1800.73265
Xen14_1810.74238
Xen14_1820.68933
Xen14_1830.71829
Xen14_1840.35941
Xen14_1850.52548
Xen14_1860.79677
Xen14_1870.60884
Xen14_1880.62664
Xen14_1890.73400
Xen14_190NA
Xen14_191NA
Xen14_1920.79663
Xen14_1930.65878
Xen14_1940.85294
Xen14_1950.85689
Xen14_1960.76032
Xen14_1970.83141
Xen14_1980.81671
Xen14_1990.79433
Xen14_2000.15172
Xen14_201NA
Xen14_2020.52138
Xen14_203NA
Xen14_204NA
Xen14_2050.43193
Xen14_2060.40695

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 input22:24:01 30 Oct 202422:25:36 30 Oct 20241 minute 35 seconds
Read alignment to reference genome22:25:37 30 Oct 202422:27:15 30 Oct 20241 minute 38 seconds
Preprocessing alignments for candidate junction identification22:27:15 30 Oct 202422:28:01 30 Oct 202446 seconds
Preliminary analysis of coverage distribution22:28:01 30 Oct 202422:34:44 30 Oct 20246 minutes 43 seconds
Identifying junction candidates22:34:44 30 Oct 202422:35:40 30 Oct 202456 seconds
Re-alignment to junction candidates22:35:40 30 Oct 202422:36:11 30 Oct 202431 seconds
Resolving best read alignments22:36:11 30 Oct 202422:37:25 30 Oct 20241 minute 14 seconds
Creating BAM files22:37:26 30 Oct 202422:38:28 30 Oct 20241 minute 2 seconds
Tabulating error counts22:38:28 30 Oct 202422:39:46 30 Oct 20241 minute 18 seconds
Re-calibrating base error rates22:39:46 30 Oct 202422:43:58 30 Oct 20244 minutes 12 seconds
Examining read alignment evidence22:43:58 30 Oct 202422:52:12 30 Oct 20248 minutes 14 seconds
Polymorphism statistics22:52:12 30 Oct 202422:52:14 30 Oct 20242 seconds
Output22:52:15 30 Oct 202422:52:51 30 Oct 202436 seconds
Output :: Mutation Prediction22:52:15 30 Oct 202422:52:17 30 Oct 20242 seconds
Output :: Mutation Annotation22:52:17 30 Oct 202422:52:17 30 Oct 20240 seconds
Total 28 minutes 49 seconds