|
FSplice |
Program provides the possibility to search for both donor and acceptor sites, and to define thresholds for them independently. Program allows to search minor variants of splicing donor site (GC-site) as well.
Output example
FSplice 1.0. Prediction of potential splice sites in Homo_sapiens genomic DNA Seq name: NM_000449 chr 1 - 148089557 148094091 4535 Length of sequence: 4535 Direct chain. Acceptor(AG) sites. Treshold 4.175 (90%). 1 P: 187 W: 7.47 Seq: attctAGccctc 2 P: 296 W: 6.42 Seq: tcttcAGaggct 3 P: 495 W: 7.30 Seq: tccctAGcagtc 4 P: 498 W: 5.72 Seq: ctagcAGtcaga 5 P: 559 W: 14.18 Seq: cccacAGcaagg 6 P: 847 W: 6.42 Seq: atggtAGcctat 7 P: 1332 W: 9.70 Seq: acctcAGcaaga 8 P: 1383 W: 9.25 Seq: ccttcAGctccc 9 P: 1393 W: 5.38 Seq: ccctcAGgaccc 10 P: 1673 W: 9.95 Seq: tctgtAGctcag 11 P: 1721 W: 4.72 Seq: cctatAGgtgga 12 P: 1916 W: 6.72 Seq: tccctAGggact 13 P: 1984 W: 9.70 Seq: cactcAGgaagt 14 P: 2366 W: 12.18 Seq: ctcccAGgtaaa 15 P: 2467 W: 7.12 Seq: cctgtAGctgag 16 P: 2638 W: 7.42 Seq: acttcAGccaga 17 P: 2779 W: 6.42 Seq: gctacAGcagca 18 P: 2867 W: 6.42 Seq: gtctcAGcaacc 19 P: 2995 W: 5.03 Seq: ctaccAGtcagt 20 P: 3033 W: 5.85 Seq: tcctcAGtttcc 21 P: 3078 W: 9.68 Seq: tctgcAGaagag 22 P: 3342 W: 9.88 Seq: tttttAGcctcc 23 P: 3545 W: 8.12 Seq: cccccAGgcttt 24 P: 4435 W: 6.70 Seq: tcctaAGgaagt 25 P: 4458 W: 6.65 Seq: tgtacAGacagc 26 P: 4513 W: 5.65 Seq: ttttcAGcttga 27 P: 4533 W: 4.58 Seq: gctttAGtg--- Donor(GT) sites. Treshold 6.099 (90%). 1 P: 40 W: 8.20 Seq: aagtgGTgagaa 2 P: 150 W: 7.50 Seq: ccagtGTgagtt 3 P: 307 W: 7.64 Seq: ccgagGTaccat 4 P: 317 W: 9.32 Seq: atttcGTaagta 5 P: 594 W: 15.48 Seq: tcctgGTaagtg 6 P: 691 W: 9.60 Seq: gagagGTagggt 7 P: 1416 W: 13.38 Seq: aaaagGTaggtt 8 P: 1794 W: 7.36 Seq: tatcgGTgggtg 9 P: 2325 W: 10.44 Seq: agagtGTaagta 10 P: 2367 W: 13.10 Seq: cccagGTaaaag 11 P: 2438 W: 8.06 Seq: tctagGTatgat 12 P: 2841 W: 7.36 Seq: cgctgGTgtgtt 13 P: 3180 W: 14.08 Seq: cccagGTaagga 14 P: 3733 W: 10.16 Seq: gagagGTaggca 15 P: 3796 W: 8.62 Seq: tacctGTgagtg 16 P: 4177 W: 11.56 Seq: caaaaGTgagtg 17 P: 4237 W: 6.38 Seq: gagagGTagaca 18 P: 4341 W: 8.06 Seq: tacagGTctgtg Reverse chain. Acceptor(AG) sites. Treshold 4.175 (90%). 1 P: 193 W: 6.42 Seq: cccacAGacctg 2 P: 292 W: 5.40 Seq: ggtgcAGtgtct 3 P: 316 W: 4.58 Seq: gccaaAGgaaaa 4 P: 481 W: 8.07 Seq: ttttcAGcctct 5 P: 517 W: 10.38 Seq: cctccAGctgag 6 P: 646 W: 4.17 Seq: tttcgAGggcgc 7 P: 709 W: 7.05 Seq: gctttAGctggt 8 P: 742 W: 6.70 Seq: ctcacAGgtact 9 P: 1424 W: 5.67 Seq: ggtttAGatgac 10 P: 1463 W: 6.97 Seq: tctgcAGaggta 11 P: 1964 W: 7.45 Seq: ttgtcAGagatc 12 P: 2035 W: 6.78 Seq: attgcAGaagcc 13 P: 2068 W: 7.25 Seq: gcctcAGctaca 14 P: 2287 W: 4.72 Seq: actgtAGcaata 15 P: 2397 W: 9.20 Seq: ctcccAGgtcct 16 P: 2421 W: 4.40 Seq: tctctAGtcaag 17 P: 2748 W: 5.08 Seq: ccgatAGgcatc 18 P: 2798 W: 5.47 Seq: cttccAGgtggt 19 P: 3064 W: 6.58 Seq: ttcccAGtgaac 20 P: 3133 W: 10.05 Seq: tctccAGtggtg 21 P: 3901 W: 9.50 Seq: ccctcAGcattt 22 P: 3945 W: 6.03 Seq: ttaccAGgatcc 23 P: 4298 W: 4.72 Seq: cccccAGtcttg 24 P: 4406 W: 11.57 Seq: tccccAGaaggc 25 P: 4440 W: 9.12 Seq: tacccAGaaagg Donor(GT) sites. Treshold 6.099 (90%). 1 P: 31 W: 8.48 Seq: aaaagGTcagag 2 P: 49 W: 10.02 Seq: accagGTactaa 3 P: 400 W: 7.08 Seq: ctttgGTatgct 4 P: 743 W: 10.02 Seq: cacagGTacttc 5 P: 832 W: 6.80 Seq: gctgaGTgagtc 6 P: 896 W: 12.40 Seq: agttgGTaagat 7 P: 1218 W: 7.64 Seq: acacaGTaaggt 8 P: 1223 W: 8.90 Seq: gtaagGTgtgaa 9 P: 1466 W: 7.64 Seq: cagagGTaccaa 10 P: 1477 W: 12.26 Seq: aaaagGTaatag 11 P: 1491 W: 11.84 Seq: tgaagGTgagga 12 P: 1830 W: 7.64 Seq: cacagGTcaggg 13 P: 2196 W: 6.94 Seq: ggaagGTgattt 14 P: 2686 W: 6.80 Seq: catggGTgaggg 15 P: 2982 W: 7.22 Seq: ccctgGTaaacc 16 P: 3159 W: 9.32 Seq: tgaagGTagaga 17 P: 3209 W: 10.16 Seq: ctgagGTaggag 18 P: 3773 W: 6.80 Seq: atcaaGTgagag 19 P: 4253 W: 8.34 Seq: gggtgGTaggtt
Acceptor(AG) sites. - the type of splicing sites. For the current case "Acceptor(AG)" means the U2-type acceptor site. Possible variants: Donor(GT) sites. means U2-type donor GT-site (Major variant). Donor(GC) sites. means U2-type donor GC- site (Minor variant).
Treshold 4.175 (90%) - means that for the current threshold value (4.175) 90% of true splicing sites are being classified as true.
P: 187 - position of splicing site
W: - weight of site.