Unconfigured Ad

Collapse
X
Collapse
 
  • Page of 1
  • Time
  • Show
Clear All
new posts
Previous template Next
  • mapardo
    Junior Member
    • Jun 2011
    • 8
    #1

    incoherent result cufflinks

    I have done a easy test with cufflinks and the result is not coherent with the input files.

    I use strand specific data (input.sam).

    ILLUMINA+GA_0000:1:1:1079:18161#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA YccacSccc\ccccXcYScc]a]a\caScc^acccK NM:i:0 NH:i:1 XS:A:-
    ILLUMINA+GA_0000:1:1:1079:18161#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:4 NH:i:1 XS:A:-
    ILLUMINA+GA_0000:1:1:1079:18161#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA YccacSccc\ccccXcYScc]a]a\caScc^acccK NM:i:0 NH:i:1 XS:A:-
    ILLUMINA+GA_0000:1:1:1079:18161#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:4 NH:i:1 XS:A:-
    ILLUMINA+GA_0000:1:1:1079:18161#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA YccacSccc\ccccXcYScc]a]a\caScc^acccK NM:i:0 NH:i:1 XS:A:-
    ILLUMINA+GA_0000:1:1:1079:18161#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:4 NH:i:1 XS:A:-
    ILLUMINA-GA_0000:1:1:1080:14412#0 83 chr3 9861 255 36M = 9828 -69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA c]\b]JccYcccccPcaRb]]Jb]aZYYaccYaccc NM:i:0 NH:i:1 XS:A:+
    ILLUMINA-GA_0000:1:1:1080:14412#0 163 chr3 9828 255 36M = 9861 69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:3 NH:i:1 XS:A:+
    ILLUMINA-GA_0000:1:1:1080:14412#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA c]\b]JccYcccccPcaRb]]Jb]aZYYaccYaccc NM:i:0 NH:i:1 XS:A:-
    ILLUMINA-GA_0000:1:1:1080:14412#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:3 NH:i:1 XS:A:-
    ILLUMINA-GA_0000:1:1:1080:14412#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA c]\b]JccYcccccPcaRb]]Jb]aZYYaccYaccc NM:i:0 NH:i:1 XS:A:-
    ILLUMINA-GA_0000:1:1:1080:14412#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:3 NH:i:1 XS:A:-
    ILLUMINA-GA_0000:1:1:1080:14412#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA c]\b]JccYcccccPcaRb]]Jb]aZYYaccYaccc NM:i:0 NH:i:1 XS:A:-
    ILLUMINA-GA_0000:1:1:1080:14412#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB NM:i:3 NH:i:1 XS:A:-


    I sort the file with samtools (input.sort.bam) and use this command for Cufflinks:

    cufflinks --library-type fr-firststrand --min-frags-per-transfrag 3 input.sort.bam

    The result is the file transcript3.gtf, where transcripts have no sense.

    chr3 Cufflinks transcript 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
    chr3 Cufflinks exon 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; exon_number "1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
    chr3 Cufflinks transcript 9828 9896 1000 + . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; FPKM "1157774097.4185056686"; frac "1.000000"; conf_lo "0.000000"; conf_hi "2315548194.837011"; cov "583.518145";
    chr3 Cufflinks exon 9828 9896 1000 + . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; exon_number "1"; FPKM "1157774097.4185056686"; frac "1.000000"; conf_lo "0.000000"; conf_hi "2315548194.837011"; cov "583.518145";


    If I use a different value for --min-frags-per-transfrag parameter, results change to the correct one

    cufflinks --library-type fr-firststrand --min-frags-per-transfrag 1 input.sort.bam

    File transcript1.gtf.


    chr3 Cufflinks transcript 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
    chr3 Cufflinks exon 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; exon_number "1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
    chr3 Cufflinks transcript 9828 9896 1000 + . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; FPKM "289443524.3546264172"; frac "0.250000"; conf_lo "0.000000"; conf_hi "868330573.063879"; cov "145.879536";
    chr3 Cufflinks exon 9828 9896 1000 + . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; exon_number "1"; FPKM "289443524.3546264172"; frac "0.250000"; conf_lo "0.000000"; conf_hi "868330573.063879"; cov "145.879536";
    chr3 Cufflinks transcript 9828 9896 1000 - . gene_id "CUFF.3"; transcript_id "CUFF.3.1"; FPKM "868330573.0638793707"; frac "0.750000"; conf_lo "0.000000"; conf_hi "1870992353.271905"; cov "437.638609";
    chr3 Cufflinks exon 9828 9896 1000 - . gene_id "CUFF.3"; transcript_id "CUFF.3.1"; exon_number "1"; FPKM "868330573.0638793707"; frac "0.750000"; conf_lo "0.000000"; conf_hi "1870992353.271905"; cov "437.638609";


    There is something I have missed?
    Tags: cufflinks, error, result
    • swbarnes2
      Senior Member
      • May 2008
      • 910
      #2
      First thing anyone thinks when seeing a string of Bs for a quality string is that you are using an oldish dataset, and those reads are the worst possible quailty. You appear to be using phred 64 offset quality scaling. So a 'B' is a 66 in ASCII, which makes it a 2 in quality. An 'h' on the other hand, is a 104 in ASCII, which is a quality of 40.

      So try your test on higher quality reads for starters. If the quality string has a lot of h's that's a good quality read.

      Newer runs have phred offsets of 33, so that '#' is the lowest quality, and 'I' is the highest.

        Comment

        • mapardo
          Junior Member
          • Jun 2011
          • 8
          #3
          Thanks for the quick reply.

          I have changed the quality of the reads to 'h' for all the nt but the result is the same.

          INPUT FILE 1
          ILLUMINA+GA_0000:1:1:1079:18161#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18161#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18162#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18162#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18163#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18163#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14414#0 83 chr3 9861 255 36M = 9828 -69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:+
          ILLUMINA+GA_0000:1:1:1080:14414#0 163 chr3 9828 255 36M = 9861 69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:+
          ILLUMINA+GA_0000:1:1:1080:14415#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14415#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14416#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14416#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14417#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14417#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-


          RESULT 1
          chr3 Cufflinks transcript 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
          chr3 Cufflinks exon 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; exon_number "1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
          chr3 Cufflinks transcript 9828 9896 1000 + . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; FPKM "1157774097.4185056686"; frac "1.000000"; conf_lo "0.000000"; conf_hi "2315548194.837011"; cov "583.518145";
          chr3 Cufflinks exon 9828 9896 1000 + . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; exon_number "1"; FPKM "1157774097.4185056686"; frac "1.000000"; conf_lo "0.000000"; conf_hi "2315548194.837011"; cov "583.518145";




          If I change the position of the less abundant fragment, the result is the correct.

          INPUT FILE 2
          ILLUMINA+GA_0000:1:1:1079:18161#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18161#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18162#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18162#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18163#0 99 chr3 1355 255 36M = 1366 47 TTCAGGAGGCCACCCGACATAGTGCAGAGAAGCGGA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1079:18163#0 147 chr3 1366 255 36M = 1355 +47 ACCCGACATANTNCAGAGAAGCGGCCCACAATGCGN hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14414#0 83 chr3 2861 255 36M = 9828 -69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:+
          ILLUMINA+GA_0000:1:1:1080:14414#0 163 chr3 2828 255 36M = 9861 69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:+
          ILLUMINA+GA_0000:1:1:1080:14415#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14415#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14416#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14416#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14417#0 99 chr3 9828 255 36M = 9861 69 CTTCATCCGCCAACTAATATTTCACTTTACATCCAA hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-
          ILLUMINA+GA_0000:1:1:1080:14417#0 147 chr3 9861 255 36M = 9828 +69 NGTCATTATTGGCTCAACTTTCCNCNCTATCTGCTT hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh NM:i:0 NH:i:1 XS:A:-

          RESULT 2
          chr3 Cufflinks transcript 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
          chr3 Cufflinks exon 1355 1401 1000 - . gene_id "CUFF.1"; transcript_id "CUFF.1.1"; exon_number "1"; FPKM "8019117703.6926527023"; frac "1.000000"; conf_lo "0.000000"; conf_hi "17278797233.473148"; cov "4041.635323";
          chr3 Cufflinks transcript 9828 9896 1000 - . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; FPKM "868330573.0638793707"; frac "1.000000"; conf_lo "0.000000"; conf_hi "1870992353.271905"; cov "437.638609";
          chr3 Cufflinks exon 9828 9896 1000 - . gene_id "CUFF.2"; transcript_id "CUFF.2.1"; exon_number "1"; FPKM "868330573.0638793707"; frac "1.000000"; conf_lo "0.000000"; conf_hi "1870992353.271905"; cov "437.638609";

          In both runs I have been using the same parameters for cufflinks:

          cufflinks --library-type fr-firststrand -F 0.01 --min-frags-per-transfrag 2 input.sort.bam

            Comment

            Previous template Next

            Latest Articles

            Collapse
            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating              ...
              10-14-2026, 09:55 AM
            • SEQadmin2
              New Genomics Technologies Take Aim at Long-Standing Limits
              by SEQadmin2


              Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

              We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing              ...
              09-28-2026, 10:25 AM
            View All

            ad_right_rmr

            Collapse

            News

            Collapse
            Topics Statistics Last Post
            Exeter-Led Study Maps What an Average Human Genome Looks Like by SEQadmin2
            Started by SEQadmin2, Yesterday, 11:13 AM
            0 responses
            19 views
            0 reactions
            		 Last Post SEQadmin2
            by SEQadmin2
            Yesterday, 11:13 AM  
            Long-Read Sequencing Uncovers Cancer Mutations Conventional Testing Missed by SEQadmin2
            Started by SEQadmin2, Yesterday, 10:48 AM
            0 responses
            18 views
            0 reactions
            		 Last Post SEQadmin2
            by SEQadmin2
            Yesterday, 10:48 AM  
            Mayo Clinic-Led Study Detects DNA Mutations Across Noncancerous Tissues by SEQadmin2
            Started by SEQadmin2, 10-15-2026, 10:52 AM
            0 responses
            36 views
            0 reactions
            		 Last Post SEQadmin2
            by SEQadmin2
            10-15-2026, 10:52 AM  
            Study Finds Chromatin Motion Holds Steady Until Cells Face Stress by SEQadmin2
            Started by SEQadmin2, 10-14-2026, 05:39 AM
            0 responses
            42 views
            0 reactions
            		 Last Post SEQadmin2
            by SEQadmin2
            10-14-2026, 05:39 AM  
            View All
            Working...