I've got another data from where hopefully will be easier you guys to calculate intron length. It's "gtf" extension. Shortly it look like this.

I Cufflinks transcript 218549 219145 1000 + . +gene_id "CUFF.54"; transcript_id "YAR062W"; FPKM "18.7631305676"; fra +c "0.697479"; conf_lo "12.494853"; conf_hi "25.031408"; cov "15.86482 +9"; full_read_support "yes"; I Cufflinks exon 218549 219145 1000 + . gene_i +d "CUFF.54"; transcript_id "YAR062W"; exon_number "1"; FPKM "18.76313 +05676"; frac "0.697479"; conf_lo "12.494853"; conf_hi "25.031408"; co +v "15.864829"; I Cufflinks transcript 217157 217492 1000 - . +gene_id "CUFF.55"; transcript_id "YAR060C"; FPKM "23.1145085280"; fra +c "0.302521"; conf_lo "11.245570"; conf_hi "34.983447"; cov "19.54405 +9"; full_read_support "yes"; I Cufflinks exon 217157 217492 1000 - . gene_i +d "CUFF.55"; transcript_id "YAR060C"; exon_number "1"; FPKM "23.11450 +85280"; frac "0.302521"; conf_lo "11.245570"; conf_hi "34.983447"; co +v "19.544059"; II Cufflinks transcript 7605 7733 1000 - . gen +e_id "YBL108C-A"; transcript_id "YBL108C-A"; FPKM "18.0747134240"; fr +ac "1.000000"; conf_lo "0.000000"; conf_hi "90.373567"; cov "8.950119 +"; full_read_support "yes"; II Cufflinks exon 7605 7733 1000 - . gene_id " +YBL108C-A"; transcript_id "YBL108C-A"; exon_number "1"; FPKM "18.0747 +134240"; frac "1.000000"; conf_lo "0.000000"; conf_hi "90.373567"; co +v "8.950119"; II Cufflinks transcript 8177 8482 1 + . gene_i +d "YBL108W"; transcript_id "YBL108W"; FPKM "0.0000000000"; frac "0.00 +0000"; conf_lo "0.000000"; conf_hi "0.000000"; cov "0.000000"; full_r +ead_support "no"; II Cufflinks exon 8177 8482 1 + . gene_id "YBL +108W"; transcript_id "YBL108W"; exon_number "1"; FPKM "0.0000000000"; + frac "0.000000"; conf_lo "0.000000"; conf_hi "0.000000"; cov "0.0000 +00"; II Cufflinks transcript 9268 9372 1 + . gene_i +d "YBL107W-A"; transcript_id "YBL107W-A"; FPKM "0.0000000000"; frac " +0.000000"; conf_lo "0.000000"; conf_hi "0.000000"; cov "0.000000"; fu +ll_read_support "no"; II Cufflinks exon 9268 9372 1 + . gene_id "YBL +107W-A"; transcript_id "YBL107W-A"; exon_number "1"; FPKM "0.00000000 +00"; frac "0.000000"; conf_lo "0.000000"; conf_hi "0.000000"; cov "0. +000000"; II Cufflinks transcript 9583 9666 1 + . gene_i +d "tL(UAA)B1"; transcript_id "tL(UAA)B1"; FPKM "0.0000000000"; frac " +0.000000"; conf_lo "0.000000"; conf_hi "0.000000"; cov "0.000000"; fu +ll_read_support "no"; II Cufflinks exon 9583 9666 1 + . gene_id "tL( +UAA)B1"; transcript_id "tL(UAA)B1"; exon_number "1"; FPKM "0.00000000 +00"; frac "0.000000"; conf_lo "0.000000"; conf_hi "0.000000"; cov "0. +000000"; I Cufflinks transcript 227742 228953 1000 + . +gene_id "CUFF.57"; transcript_id "YAR073W"; FPKM "5.8852029015"; frac + "0.892308"; conf_lo "3.699493"; conf_hi "8.070912"; cov "3.517134"; +full_read_support "yes"; I Cufflinks exon 227742 228953 1000 + . gene_i +d "CUFF.57"; transcript_id "YAR073W"; exon_number "1"; FPKM "5.885202 +9015"; frac "0.892308"; conf_lo "3.699493"; conf_hi "8.070912"; cov " +3.517134"; II Cufflinks transcript 142112 142868 1000 - . + gene_id "YBL040C"; transcript_id "YBL040C"; FPKM "246.2287249186"; f +rac "0.288110"; conf_lo "225.158166"; conf_hi "267.299284"; cov "208. +319539"; full_read_support "yes"; II Cufflinks exon 142112 142749 1000 - . gene_ +id "YBL040C"; transcript_id "YBL040C"; exon_number "1"; FPKM "246.228 +7249186"; frac "0.288110"; conf_lo "225.158166"; conf_hi "267.299284" +; cov "208.319539"; II Cufflinks exon 142847 142868 1000 - . gene_ +id "YBL040C"; transcript_id "YBL040C"; exon_number "2"; FPKM "246.228 +7249186"; frac "0.288110"; conf_lo "225.158166"; conf_hi "267.299284" +; cov "208.319539";
Each lines starts with "I", "II" and etc. Here they may seen written in several lines. After each transcript there comes exon. I need to minus total exon length from transcription length to get intron length in most cases it's equal to zero but some how a few transcripts include couple of exons between which can be found gaps, which will be intron length.


In reply to intron length by MBobur

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":



  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.