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

Read File Information

read filereadsbasespassed filtersaveragelongestmapped
errorsUnknown_CH607-002R0003_good_1.fq4,487,615668,794,458100.0%149.0 bases150 bases92.1%
total4,487,615668,794,458100.0%149.0 bases150 bases92.1%

Reference Sequence Information

seq idlengthfit meanfit relative_variance% mapped readsdescription
coveragedistributionXen14_1242,55688.72.63.6%na
coveragedistributionXen14_2185,21063.92.31.9%na
coveragedistributionXen14_3166,311153.43.44.1%na
coveragedistributionXen14_4151,524105.83.12.6%na
coveragedistributionXen14_5149,062108.92.72.6%na
coveragedistributionXen14_6146,760133.63.13.2%na
coveragedistributionXen14_7136,51469.42.31.5%na
coveragedistributionXen14_8132,518164.33.83.5%na
coveragedistributionXen14_9129,936164.34.23.5%na
coveragedistributionXen14_10107,217176.33.53.1%na
coveragedistributionXen14_11106,943127.02.82.2%na
coveragedistributionXen14_12100,89888.22.71.5%na
coveragedistributionXen14_1398,31898.72.71.6%na
coveragedistributionXen14_1489,144175.74.12.5%na
coveragedistributionXen14_1583,86975.92.31.0%na
coveragedistributionXen14_1681,941115.42.81.5%na
coveragedistributionXen14_1780,390174.13.52.4%na
coveragedistributionXen14_1878,780146.83.11.9%na
coveragedistributionXen14_1977,008143.13.21.8%na
coveragedistributionXen14_2075,56392.92.31.3%na
coveragedistributionXen14_2174,83760.81.90.7%na
coveragedistributionXen14_2274,457180.04.22.2%na
coveragedistributionXen14_2373,975126.42.71.5%na
coveragedistributionXen14_2473,66881.32.41.0%na
coveragedistributionXen14_2573,24068.92.10.8%na
coveragedistributionXen14_2673,210152.63.51.8%na
coveragedistributionXen14_2771,260126.84.21.5%na
coveragedistributionXen14_2864,05463.92.10.7%na
coveragedistributionXen14_2962,12365.22.10.7%na
coveragedistributionXen14_3061,66199.72.61.0%na
coveragedistributionXen14_3158,819159.86.11.5%na
coveragedistributionXen14_3255,39165.52.00.6%na
coveragedistributionXen14_3353,24193.32.30.8%na
coveragedistributionXen14_3453,084138.53.11.2%na
coveragedistributionXen14_3552,845168.23.41.4%na
coveragedistributionXen14_3651,934121.02.81.0%na
coveragedistributionXen14_3751,424101.72.80.8%na
coveragedistributionXen14_3848,836111.63.51.0%na
coveragedistributionXen14_3948,48568.43.40.5%na
coveragedistributionXen14_4046,79169.42.10.6%na
coveragedistributionXen14_4141,431111.34.00.7%na
coveragedistributionXen14_4241,403168.74.51.1%na
coveragedistributionXen14_4341,296124.44.00.8%na
coveragedistributionXen14_4440,957134.82.70.9%na
coveragedistributionXen14_4540,401168.13.81.1%na
coveragedistributionXen14_4640,351116.42.70.8%na
coveragedistributionXen14_4736,274159.03.60.9%na
coveragedistributionXen14_4834,806153.83.50.9%na
coveragedistributionXen14_4934,102168.94.40.9%na
coveragedistributionXen14_5033,369129.53.10.7%na
coveragedistributionXen14_5132,311172.93.80.9%na
coveragedistributionXen14_5231,365150.54.80.7%na
coveragedistributionXen14_5330,293163.43.80.8%na
coveragedistributionXen14_5428,579175.73.30.8%na
coveragedistributionXen14_5526,71599.32.70.4%na
coveragedistributionXen14_5626,428160.32.80.7%na
coveragedistributionXen14_5726,24280.72.20.3%na
coveragedistributionXen14_5824,572141.63.00.6%na
coveragedistributionXen14_5924,177124.62.70.5%na
coveragedistributionXen14_6023,18776.42.30.3%na
coveragedistributionXen14_6121,383124.42.60.4%na
coveragedistributionXen14_6220,21075.32.20.2%na
coveragedistributionXen14_6319,884113.22.20.4%na
coveragedistributionXen14_6418,26984.21.70.3%na
coveragedistributionXen14_6517,85676.92.10.2%na
coveragedistributionXen14_6617,785122.73.30.3%na
coveragedistributionXen14_6717,704102.62.90.3%na
coveragedistributionXen14_6816,382175.53.20.5%na
coveragedistributionXen14_6914,625137.14.30.3%na
coveragedistributionXen14_7014,40270.22.10.2%na
coveragedistributionXen14_7114,304122.65.00.3%na
coveragedistributionXen14_7213,687132.22.20.3%na
coveragedistributionXen14_7313,160134.34.00.3%na
coveragedistributionXen14_7413,069181.62.20.4%na
coveragedistributionXen14_7513,052103.53.70.2%na
coveragedistributionXen14_7612,931115.83.80.2%na
coveragedistributionXen14_7712,669100.83.70.2%na
coveragedistributionXen14_7812,559188.64.00.4%na
coveragedistributionXen14_7912,553150.12.40.3%na
coveragedistributionXen14_8012,035138.94.30.3%na
coveragedistributionXen14_8112,02072.23.60.1%na
coveragedistributionXen14_8211,983120.73.90.2%na
coveragedistributionXen14_8311,356133.23.10.2%na
coveragedistributionXen14_8411,356122.32.40.2%na
coveragedistributionXen14_8511,197137.22.30.3%na
coveragedistributionXen14_8610,95375.02.30.2%na
coveragedistributionXen14_8710,348244.16.00.4%na
coveragedistributionXen14_8810,25966.12.00.1%na
coveragedistributionXen14_8910,105134.86.30.2%na
coveragedistributionXen14_909,218126.73.30.2%na
coveragedistributionXen14_918,913138.34.90.2%na
coveragedistributionXen14_928,23264.62.10.1%na
coveragedistributionXen14_937,77775.72.30.1%na
coveragedistributionXen14_947,58470.34.30.1%na
coveragedistributionXen14_957,34264.21.60.1%na
coveragedistributionXen14_966,91679.56.00.1%na
coveragedistributionXen14_976,746166.65.10.2%na
coveragedistributionXen14_986,594154.73.60.2%na
coveragedistributionXen14_996,04576.51.60.1%na
coveragedistributionXen14_1005,846146.83.40.1%na
coveragedistributionXen14_1015,615865.030.20.8%na
coveragedistributionXen14_1025,55264.92.60.1%na
coveragedistributionXen14_1035,30773.52.10.1%na
coveragedistributionXen14_1045,09697.52.20.1%na
coveragedistributionXen14_1055,061112.91.60.1%na
coveragedistributionXen14_1064,995130.34.40.1%na
coveragedistributionXen14_1074,93285.72.70.1%na
coveragedistributionXen14_1084,816128.23.00.1%na
coveragedistributionXen14_1094,61699.32.80.1%na
coveragedistributionXen14_1104,57859.91.70.0%na
coveragedistributionXen14_1114,562695.522.40.5%na
coveragedistributionXen14_1124,48677.21.80.1%na
coveragedistributionXen14_1134,275189.22.70.1%na
coveragedistributionXen14_1144,03975.82.20.0%na
coveragedistributionXen14_1153,633131.64.10.1%na
coveragedistributionXen14_1163,428145.11.60.1%na
coveragedistributionXen14_1173,42771.41.80.0%na
coveragedistributionXen14_1183,35577.04.60.0%na
coveragedistributionXen14_1193,208101.73.70.1%na
coveragedistributionXen14_1203,12367.92.60.0%na
coveragedistributionXen14_1213,076100.92.10.0%na
coveragedistributionXen14_1222,908124.21.30.1%na
coveragedistributionXen14_1232,7141116.84.30.5%na
coveragedistributionXen14_1242,576126.02.50.1%na
coveragedistributionXen14_1252,519124.02.90.0%na
coveragedistributionXen14_1262,515130.11.40.1%na
coveragedistributionXen14_1272,166692.619.90.2%na
coveragedistributionXen14_1282,14389.83.00.0%na
coveragedistributionXen14_1292,049135.31.20.0%na
coveragedistributionXen14_1302,000136.62.10.0%na
coveragedistributionXen14_1311,92259.92.70.0%na
coveragedistributionXen14_1321,806166.92.00.1%na
coveragedistributionXen14_1331,721296.02.30.1%na
coveragedistributionXen14_1341,65773.25.90.0%na
coveragedistributionXen14_1351,633145.01.20.0%na
coveragedistributionXen14_1361,567125.82.00.0%na
coveragedistributionXen14_1371,544240.64.20.1%na
coveragedistributionXen14_1381,361138.53.60.0%na
coveragedistributionXen14_1391,349143.61.60.0%na
coveragedistributionXen14_1401,2981042.43.20.2%na
coveragedistributionXen14_1411,2922411.079.00.4%na
coveragedistributionXen14_1421,28159.71.20.0%na
coveragedistributionXen14_1431,1071408.310.90.1%na
coveragedistributionXen14_1441,091754.721.20.1%na
coveragedistributionXen14_1451,062802.34.50.1%na
coveragedistributionXen14_1461,046126.41.80.0%na
coveragedistributionXen14_1471,007204.23.40.0%na
coveragedistributionXen14_14898976.73.80.0%na
coveragedistributionXen14_149980139.73.40.0%na
coveragedistributionXen14_150945161.81.00.0%na
coveragedistributionXen14_151938198.41.70.0%na
coveragedistributionXen14_152934143.53.60.0%na
coveragedistributionXen14_1539292848.485.20.4%na
coveragedistributionXen14_154920570.6116.50.0%na
coveragedistributionXen14_15590061.01.80.0%na
coveragedistributionXen14_156862239.42.10.0%na
coveragedistributionXen14_157857362.733.30.0%na
coveragedistributionXen14_158799153.76.10.0%na
coveragedistributionXen14_159798139.82.00.0%na
coveragedistributionXen14_160780177.21.00.0%na
coveragedistributionXen14_161770139.77.00.0%na
coveragedistributionXen14_162765161.21.00.0%na
coveragedistributionXen14_163753297.28.70.0%na
coveragedistributionXen14_16474678.42.40.0%na
coveragedistributionXen14_165691161.68.80.0%na
coveragedistributionXen14_1666852463.5335.00.2%na
coveragedistributionXen14_1676661083.818.10.1%na
coveragedistributionXen14_1686293509.7557.40.1%na
coveragedistributionXen14_169626243.31.40.0%na
coveragedistributionXen14_170619160.75.70.0%na
coveragedistributionXen14_171569160.01.80.0%na
coveragedistributionXen14_17256961.91.00.0%na
coveragedistributionXen14_173564[235.8][36.4]0.0%na
coveragedistributionXen14_174564159.31.00.0%na
coveragedistributionXen14_175560151.35.90.0%na
coveragedistributionXen14_176552248.41.00.0%na
coveragedistributionXen14_177549162.66.80.0%na
coveragedistributionXen14_178541129.75.60.0%na
coveragedistributionXen14_179536141.61.00.0%na
coveragedistributionXen14_180526165.71.30.0%na
coveragedistributionXen14_181511196.226.90.0%na
coveragedistributionXen14_182478210.210.70.0%na
coveragedistributionXen14_183477132.81.20.0%na
coveragedistributionXen14_1844731214.12.00.1%na
coveragedistributionXen14_185471[427.6][67.1]0.0%na
coveragedistributionXen14_186465153.334.30.0%na
coveragedistributionXen14_187464[270.1][26.9]0.0%na
coveragedistributionXen14_188458565.667.30.0%na
coveragedistributionXen14_189453278.343.60.0%na
coveragedistributionXen14_19044872.85.10.0%na
coveragedistributionXen14_191447[428.3][247.4]0.0%na
coveragedistributionXen14_192445[99.7][19.5]0.0%na
coveragedistributionXen14_193444379.34.90.0%na
coveragedistributionXen14_19444223.01.20.0%na
coveragedistributionXen14_195441[46.9][14.9]0.0%na
coveragedistributionXen14_196436218.521.60.0%na
coveragedistributionXen14_197433134.61.20.0%na
coveragedistributionXen14_198431[73.3][15.4]0.0%na
coveragedistributionXen14_199428205.31.00.0%na
coveragedistributionXen14_200423[1023.0][122.2]0.1%na
coveragedistributionXen14_201421141.61.00.0%na
coveragedistributionXen14_2024211399.8176.10.0%na
coveragedistributionXen14_203416[113.7][21.1]0.0%na
coveragedistributionXen14_204415166.71.70.0%na
coveragedistributionXen14_205404[856.0][83.6]0.0%na
coveragedistributionXen14_2064031393.24.40.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≤ 100000100000
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 ≤ 5000532
Total length of all junction candidates (factor times the reference genome length)≤ 0.10.030

Junction Skew Score Calculation

reference sequencepr(no read start)
Xen14_10.79116
Xen14_20.84258
Xen14_30.67839
Xen14_40.75841
Xen14_50.75192
Xen14_60.70762
Xen14_70.82979
Xen14_80.66137
Xen14_90.66105
Xen14_100.64311
Xen14_110.71971
Xen14_120.79252
Xen14_130.77088
Xen14_140.64686
Xen14_150.81469
Xen14_160.73897
Xen14_170.64292
Xen14_180.68789
Xen14_190.69623
Xen14_200.78145
Xen14_210.84816
Xen14_220.63669
Xen14_230.72149
Xen14_240.80643
Xen14_250.83009
Xen14_260.67735
Xen14_270.71835
Xen14_280.84091
Xen14_290.83825
Xen14_300.77143
Xen14_310.67421
Xen14_320.83804
Xen14_330.78359
Xen14_340.70132
Xen14_350.65694
Xen14_360.73204
Xen14_370.76802
Xen14_380.74270
Xen14_390.83468
Xen14_400.83074
Xen14_410.75118
Xen14_420.65151
Xen14_430.72569
Xen14_440.70979
Xen14_450.65512
Xen14_460.73900
Xen14_470.67164
Xen14_480.67713
Xen14_490.65294
Xen14_500.71710
Xen14_510.65201
Xen14_520.69005
Xen14_530.65933
Xen14_540.65172
Xen14_550.76946
Xen14_560.66708
Xen14_570.80512
Xen14_580.69435
Xen14_590.72761
Xen14_600.81432
Xen14_610.72682
Xen14_620.81608
Xen14_630.74407
Xen14_640.79890
Xen14_650.81575
Xen14_660.73092
Xen14_670.76528
Xen14_680.64553
Xen14_690.70714
Xen14_700.83014
Xen14_710.73147
Xen14_720.71064
Xen14_730.72033
Xen14_740.64511
Xen14_750.76061
Xen14_760.74000
Xen14_770.77149
Xen14_780.62704
Xen14_790.68275
Xen14_800.70274
Xen14_810.83648
Xen14_820.73454
Xen14_830.71786
Xen14_840.74949
Xen14_850.68786
Xen14_860.80171
Xen14_870.55861
Xen14_880.82521
Xen14_890.72014
Xen14_900.72559
Xen14_910.70885
Xen14_920.84603
Xen14_930.81953
Xen14_940.83742
Xen14_950.84241
Xen14_960.73936
Xen14_970.67099
Xen14_980.67497
Xen14_990.82035
Xen14_1000.68645
Xen14_1010.21380
Xen14_1020.84609
Xen14_1030.81807
Xen14_1040.77698
Xen14_1050.70846
Xen14_1060.72703
Xen14_1070.80647
Xen14_1080.71782
Xen14_1090.78189
Xen14_1100.84851
Xen14_1110.29614
Xen14_1120.80915
Xen14_1130.64877
Xen14_1140.82310
Xen14_1150.72695
Xen14_1160.70405
Xen14_1170.83309
Xen14_1180.82876
Xen14_1190.76683
Xen14_1200.83333
Xen14_1210.77259
Xen14_1220.73986
Xen14_1230.17668
Xen14_1240.73447
Xen14_1250.74137
Xen14_1260.67674
Xen14_1270.29963
Xen14_1280.80339
Xen14_1290.71547
Xen14_1300.71425
Xen14_1310.81556
Xen14_1320.62138
Xen14_1330.53312
Xen14_1340.78600
Xen14_1350.67425
Xen14_1360.74505
Xen14_1370.58290
Xen14_1380.69583
Xen14_1390.71868
Xen14_1400.22612
Xen14_1410.14745
Xen14_1420.82365
Xen14_1430.41328
Xen14_1440.31072
Xen14_1450.31733
Xen14_1460.76769
Xen14_1470.62399
Xen14_1480.83670
Xen14_1490.73827
Xen14_1500.73016
Xen14_1510.67005
Xen14_1520.72270
Xen14_1530.15985
Xen14_1540.63424
Xen14_1550.86389
Xen14_1560.68047
Xen14_1570.58985
Xen14_1580.73529
Xen14_1590.76504
Xen14_1600.71154
Xen14_1610.76558
Xen14_1620.71503
Xen14_1630.59429
Xen14_1640.84406
Xen14_1650.73227
Xen14_1660.39557
Xen14_1670.31907
Xen14_1680.27027
Xen14_1690.70573
Xen14_1700.73667
Xen14_1710.73989
Xen14_1720.88840
Xen14_173NA
Xen14_1740.73138
Xen14_1750.75446
Xen14_1760.64855
Xen14_1770.74681
Xen14_1780.77542
Xen14_1790.76866
Xen14_1800.74905
Xen14_1810.71918
Xen14_1820.67573
Xen14_1830.69174
Xen14_1840.36469
Xen14_185NA
Xen14_1860.80323
Xen14_187NA
Xen14_1880.60153
Xen14_1890.71854
Xen14_1900.86942
Xen14_191NA
Xen14_192NA
Xen14_1930.59685
Xen14_1940.87783
Xen14_195NA
Xen14_1960.75573
Xen14_1970.69759
Xen14_198NA
Xen14_1990.79433
Xen14_200NA
Xen14_2010.79572
Xen14_2020.51306
Xen14_203NA
Xen14_2040.75301
Xen14_205NA
Xen14_2060.40074

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:51:46 30 Oct 202421:53:09 30 Oct 20241 minute 23 seconds
Read alignment to reference genome21:53:10 30 Oct 202421:55:22 30 Oct 20242 minutes 12 seconds
Preprocessing alignments for candidate junction identification21:55:22 30 Oct 202421:56:09 30 Oct 202447 seconds
Preliminary analysis of coverage distribution21:56:09 30 Oct 202422:03:27 30 Oct 20247 minutes 18 seconds
Identifying junction candidates22:03:27 30 Oct 202422:04:37 30 Oct 20241 minute 10 seconds
Re-alignment to junction candidates22:04:37 30 Oct 202422:05:23 30 Oct 202446 seconds
Resolving best read alignments22:05:23 30 Oct 202422:06:42 30 Oct 20241 minute 19 seconds
Creating BAM files22:06:42 30 Oct 202422:07:43 30 Oct 20241 minute 1 second
Tabulating error counts22:07:43 30 Oct 202422:09:33 30 Oct 20241 minute 50 seconds
Re-calibrating base error rates22:09:33 30 Oct 202422:13:56 30 Oct 20244 minutes 23 seconds
Examining read alignment evidence22:13:56 30 Oct 202422:23:21 30 Oct 20249 minutes 25 seconds
Polymorphism statistics22:23:21 30 Oct 202422:23:26 30 Oct 20245 seconds
Output22:23:26 30 Oct 202422:23:59 30 Oct 202433 seconds
Output :: Mutation Prediction22:23:26 30 Oct 202422:23:31 30 Oct 20245 seconds
Output :: Mutation Annotation22:23:31 30 Oct 202422:23:31 30 Oct 20240 seconds
Total 32 minutes 17 seconds