SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Updated How to convert .txt file to .bed .GFF or .BAR file format, forevermark4 Bioinformatics 2 06-30-2014 05:02 AM
How to convert .txt file to .bed or .gff, How can we use chip seq data in R software forevermark4 Bioinformatics 57 06-30-2014 05:01 AM
Is there a BED file format validator? Does a BED file have to be sorted position? LauraSmith Bioinformatics 3 05-21-2013 11:54 AM
from excel file to genomic features ohsu Bioinformatics 5 01-31-2012 10:33 PM
extract alignment from SAM with a GFF file NicoBxl Bioinformatics 4 08-02-2011 01:45 PM

Reply
 
Thread Tools
Old 06-26-2013, 06:26 AM   #1
jazz
Member
 
Location: us

Join Date: Nov 2008
Posts: 28
Default compare genomic features from gff file to a SAM/BED file

I have an alignment from Bowtie2 for my RNA-IP and I want to compile a pie chart of % of reads mapping to different genomic features like exons, introns, UTR's etc. I am planning to compare my BAM/BED file to a genome annotation gff file. I have heard of intersectbed in bedtools, however the gff files and bed files have genomic coordinates in different columns. can anyone suggest a simple fix for that? I have also heard of people using mpileup option in samtools, although I have no idea how that would work here. I am a bioinformatics newbie, so writing custom scripts is out of question at this point. Given that this is a very commonly sought after problem, I am hopeful that one of the great minds would have found an easy solution for this.
Thanks in advance!
jazz is offline   Reply With Quote
Old 06-26-2013, 02:06 PM   #2
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 311
Default

Quote:
Originally Posted by jazz View Post
I have heard of intersectbed in bedtools, however the gff files and bed files have genomic coordinates in different columns. can anyone suggest a simple fix for that?
No need to fix anything! Bedtools recognizes the file format (GFF/bed/bam) and handles it accordingly. Have you tried it at all with your files?

Dario
dariober is offline   Reply With Quote
Old 07-01-2013, 07:32 AM   #3
jazz
Member
 
Location: us

Join Date: Nov 2008
Posts: 28
Default

I am getting a bizarre result from using the intersectbed from bedtools. The output bed file contains strange characters and not the chromosomal coordinates as I expected.
Any suggestions?
Can anyone explain if the samtools mpileup would help?

Thanks
jazz is offline   Reply With Quote
Old 07-01-2013, 07:53 AM   #4
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 311
Default

Quote:
Originally Posted by jazz View Post
I am getting a bizarre result from using the intersectbed from bedtools. The output bed file contains strange characters and not the chromosomal coordinates as I expected.
Any suggestions?
Can anyone explain if the samtools mpileup would help?

Thanks
Hi there,

It would help other people answering you if you post
- The command line you executed
- A sample (a few lines) of the output including errors and messages
- A sample of the input files

Dario
dariober is offline   Reply With Quote
Old 07-01-2013, 07:59 AM   #5
jazz
Member
 
Location: us

Join Date: Nov 2008
Posts: 28
Default

Sorry!

Here is the command I tried:

intersectbed -a Inputfile.bed -b compare.gff > output.bed

The first few lines of the resulting bed file look like this:

BAM◊‹%@HD VN:1.0 SO:unsorted
@SQ SN:scaffold16 LN:1343
@SQ SN:scaffold6 LN:1021
@SQ SN:scaffold4 LN:958
@SQ SN:scaffold22 LN:923

Below, there is something like this:

$&%('%# ##ASCXScXNCXMCXOCXGCNMCMDZ20YTZUUM("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1202:3055:549570((A"))))))'''''%%#%%"""ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUUP)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2106:13176:49083 !$!H))((&'''''%%%%%ASCXSc˚XNCXMCXOCXGCNMCMDZ18YTZUUp(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1216:3600:66565@B(()()))))&'%&&%%%##"ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUUs(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1213:3015:62620BH!(B)))(()())))'&'%'##%%#"!ASCXSCXNCXMCXOCXGCNMCMDZ24YTZUUu(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1201:4444:11571PH!(B))))))))'''''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUUv(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1204:3545:82135P(H$!))(&))()&&'''%%%%%"ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUUv)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2207:17932:14595(H$!!))))))))))(''&''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ24YTZUUy)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2212:15502:16113P!(BB% $" ''%###ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU|)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2202:15199:42536PH$!!@"""%%%%%&'''')))))('(ASCXSc˙XNCXMCXOCXGCNMCMDZ21YTZUU(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1209:14989:7265P(BB %!&$'$'%#!%####ASCXNCXMCXOCXGCNMCMDZ21YTZUU(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:3945:10998`H$!!HB"""%%%%%'''''()'()))))ASCXNCXMCXOCXGCNMCMDZ22YTZUU)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2106:14259:54678`!B(HBDD)((('('((''''&%%%%%"!ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUU™)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2103:21161:300970D(H!(((&'&#%'&'%####ASc˚XSc˚XNCXMCXOCXGCNMCMDZ10G8YTZUU≤(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2109:8591:70631@DH!$H&&$''&&%'''%%%%%%ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU√(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2209:1385:21409PAB"(!(@()())))('''''%%%#%""ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU∆'IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2107:3494:9530`(!"))))())))'''''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUU∆(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2213:6747:92868`(!"'()(&&(('% '%#! %##ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUU*)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2109:16654:84992`!"$"'$&$'$!#% %!###ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUU7)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2105:14962:85910!"ABH"&'))&''''%%%%%!"!ASCXSCXNCXMCXOCXGCNMCMDZ17YTZUUM(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1115:3600:974280("!""@())()('''#'%%%%%"""ASc˚XSc˚XNCXMCXOCXGCNMCMDZ6A12YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2205:5141:38299P!B$BHD ))))))(''''''%%%%#!!ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU)#IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1103:14890:79414@D$A!B)))))))'''''#%%%%"""ASCXSc˙XNCXMCXOCXGCNMCMDZ20YTZUU)#IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1108:19803:69818@D$A!B)))))))'''''#%%%%"""ASCXSc˙XNCXMCXOCXGCNMCMDZ20YTZUU(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1105:9529:99475@D!"!"""%%%%%'''%&))(((((ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1110:4076:19142p!!(BD))))))))))'''''%%%%%"""ASCXSc˚XNCXMCXOCXGCNMCMDZ23YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2215:1801:25480`!!(BD$')&&&')(%'%%'%%%%#ASCXSc˚XNCXMCXOCXGCNMCMDZ22YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1209:5234:51548PHAB(D@)))(())))'''''%%%%%!ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2204:6685:91033PHAB(D@)))(&)))'''''%%%%%"""ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU™)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2216:18900:37651@!(BD)()))))''''&%%%%%"""ASCXSc˚XNCXMCXOCXGCNMCMDZ20YTZUU("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1104:5341:951690HAB(D@(((%(&'%#' ## %ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUU∂(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:3109:87706p(BDH"""%%%%%%''''(()))$())(ASCXNCXMCXOCXGCNMCMDZ23YTZUU√(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:9089:16084pAAHBA@)))((((((('''''%%%#%"""ASc˚XNCXMCXOCXGCNMCMDZ21T1YTZUU”) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1211:17132:48441`HDAH(HA)))))))))'''''%%%%%"""ASCXSc˚XNCXMCXOCXGCNMCMDZ22YTZUU’)#IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2108:10056:99686@HDAH(HA%%%###'''#')(()))ASCXSc˙XNCXMCXOCXGCNMCMDZ20YTZUU’) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2208:15899:52213pHDAH(HAH###!###"!%(" " %ASCXSc˚XNCXMCXOCXGCNMCMDZ23YTZUU(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1203:2989:39462`HDH)))))))(&'''''%%%%%"""ASc˙XNCXMCXOCXGCNMCMDZ0T21YTZUU)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2209:11941:93467@HD)(()((('&'&'%%%%%!ASc˙XNCXMCXOCXGCNMCMDZ0T19YTZUUfi("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1205:4676:31348@DHA!#'%"% %#" !#!ASCXSc˚XNCXMCXOCXGCNMCMDZ20YTZUU‡'"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1210:7109:8286@DHA!'#% %''#!'!%###ASCXSc˚XNCXMCXOCXGCNMCMDZ20YTZUU‰) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1208:20748:42297PDHA!@"""%%%%%'''''((((())(ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU‰) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2201:19649:94346`DHA!H"""%%%%%'''&'))))))))(ASCXSc˚XNCXMCXOCXGCNMCMDZ22YTZUU) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2111:17700:75673PAHA%#$ #'%%'')&(())))(ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2208:18875:78749`AHAA))))())()'''&&%%%%%""ASCXSc˙XNCXMCXOCXGCNMCMDZ22YTZUU(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:5271:37504PHA!H!##%#'''''&&(((('&ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUUˇ'*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1204:8193:9069pADHA """%%%%%''%''))(&')))))ASCXNCXMCXOCXGCNMCMDZ23YTZUU(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1111:5151:28312PADH(D!%%%%#%''''())'))))ASCXNCXMCXOCXGCNMCMDZ21YTZUU(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2113:4928:365430ADH(D !!#%%% "&(((&'ASCXNCXMCXOCXGCNMCMDZ19YTZUU)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1114:15413:28126PDHA(('"(((("%%''####%"ASCXNCXMCXOCXGCNMCMDZ21YTZUU)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2101:20368:77815DHA((@)))))))))()&'''''%%%%%"!ASCXNCXMCXOCXGCNMCMDZ25YTZUU(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2111:5703:20622PDHA()('$&())'% % $###ASCXNCXMCXOCXGCNMCMDZ21YTZUU(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1110:9920:44661H(D(&$)((('(&&'''''%%%%%""ASCXNCXMCXOCXGCNMCMDZ24YTZUU )*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2215:15395:89845PD "%%####'%''$'(()''%ASCXNCXMCXOCXGCNMCMDZ21YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1202:3721:24291!B"""%%%%%'''''))))))))))))ASCXSc˚XNCXMCXOCXGCNMCMDZ25YTZUU((IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1211:7680:26402pDD(ABHB))))))))))'''''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ23YTZUU))IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1102:14775:70644pHDH!(%%%%%'''')))))''())ASc˙XSc˙XNCXMCXOCXGCNMCMDZ22G0YTZUUF)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2202:16997:641860A(B("))('$$''%%"%%%%%"""ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUUP)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1102:13792:73293@("D$(&))))))'''''%%%%%""!ASCXNCXMCXOCXGCNMCMDZ20YTZUUP("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1116:6671:74058 ("D$(&$"&!$%#%%###ASCXSc˚XNCXMCXOCXGCNMCMDZ18YTZUUS)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1214:12261:89235@!ABA!""%##%%'&''')))()()ASCXNCXMCXOCXGCNMCMDZ20YTZUUS(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1215:18654:3307@!ABA!"#%%%%'''''&())(((ASCXNCXMCXOCXGCNMCMDZ20YTZUUS)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2113:13606:52861@!ABA!"#%#%%&#%%'(((())&ASCXNCXMCXOCXGCNMCMDZ20YTZUUS(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2202:5262:60786@!ABA"""%#%%%'''''))))())ASc˙XNCXMCXOCXGCNMCMDZ19C0YTZUUS)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2203:13218:67864@!ABA"%%#%%&&%''())))))ASCXNCXMCXOCXGCNMCMDZ20YTZUUT)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2105:15872:468210"D$(  #  "!%! ""ASCXNCXMCXOCXGCNMCMDZ19YTZUUw(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2203:11181:5745PH!B"$ ))))))))'''''%%%%#"!ASc˙XSc˙XNCXMCXOCXGCNMCMDZ5C15YTZUU}(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2214:7771:596390B"$!"&'%$('%'%###!!ASCXSCXNCXMCXOCXGCNMCMDZ19YTZUU)#IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2210:19479:23401 BAAB(D)))))''''&%%%%%"""ASCXSc˙XNCXMCXOCXGCNMCMDZ18YTZUU) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1107:19778:32622p!"A( """%%%%%''''')))))(()'(ASCXSc˚XNCXMCXOCXGCNMCMDZ23YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2207:2421:69609P!A!"""%%%%%'''''))))()))ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1105:6587:40443!"A($))(())))())'''''$#%#%"""AScXNCXMCXOCXGCNMCMDZ0A0A22YTZUU∆)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1209:11138:58479PA($((('%(()'''''$##%#""ASc˙XNCXMCXOCXGCNMCMDZ0C20YTZUU("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2208:7830:655820!AB(('$)(%''''%%%%#ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUU–(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2113:2124:86226pA($!$(D
"  !#ASCXNCXMCXOCXGCNMCMDZ23YTZUU—)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:10490:13666pABAA('%)))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ23YTZUU“(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1112:6679:60130`($!$(D"%%%%%''''')'(&((('(ASCXNCXMCXOCXGCNMCMDZ22YTZUU“)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2102:10863:63575`($!$(D"""%%%%%'''''))))))''(ASCXNCXMCXOCXGCNMCMDZ22YTZUU“(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1207:1741:35309P($!$(D''(&(&)(''&''#%$#%ASCXNCXMCXOCXGCNMCMDZ21YTZUU“(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2103:7127:52009P($!$(D%%'%'&(''''''%%%%%"ASCXNCXMCXOCXGCNMCMDZ21YTZUU“(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2204:5252:65081p($!$(D))())())()'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ23YTZUU“)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2212:10510:35497`($!$(D"
% $!!!!%%%!###ASCXNCXMCXOCXGCNMCMDZ22YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1102:2867:20247@ABAA')()())'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ20YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1102:8425:48665`ABAA())))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ22YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1103:13170:30492`ABAA'))))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ22YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1104:6781:65853PABAA&('$()()'%%%'#%%%%"""ASCXNCXMCXOCXGCNMCMDZ21YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1112:7147:56091@ABAA'&&))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ20YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1112:8098:65458PABAA())(()(('''''%%%%%"ASCXNCXMCXOCXGCNMCMDZ21YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:16505:73260`ABAA'&(()))''%%'&'%%%%%!ASCXNCXMCXOCXGCNMCMDZ22YTZUU”)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1211:10351:13410pBABAA(&$()()))'%'''&%%%%%"""ASc˙XNCXMCXOCXGCNMCMDZ0A22YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2104:20388:8832PABAA))))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ21YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2105:15570:76700`ABAA'))))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ22YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2108:1907:87880PABAA''&)))))'''''%%%%%ASCXNCXMCXOCXGCNMCMDZ21YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2110:7114:11279pABAA%)()))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ23YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2114:8347:67654PABAA'#(&((((&&''!###%#ASCXNCXMCXOCXGCNMCMDZ21YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2115:19327:11357`ABAA!(%)(())('%'''%%%%%""ASCXNCXMCXOCXGCNMCMDZ22YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2201:11110:50153ABAA'#(&(&#%$$ #%!%&'%%#%%ASCXNCXMCXOCXGCNMCMDZ25YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2201:12566:71306ABAA%&$)())('% "%%#%'%%%%%"!ASCXNCXMCXOCXGCNMCMDZ25YTZUU”)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2202:18198:47631ABAA(()(&()()'&'&''%$%###!!ASCXNCXMCXOCXGCNMCMDZ24YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2202:2362:66940ABAA())))))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ24YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2203:2172:23622@ABAA%'&()))'&&''%%%%%""ASCXNCXMCXOCXGCNMCMDZ20YTZUU”(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2205:4395:36383pABAA()))))))()'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ23YTZUU”)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2216:10710:74034pABAA&&&$'(&'%" #'#'$#%%#ASc˚XNCXMCXOCXGCNMCMDZ0A22YTZUU‘(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1105:3277:44016($!$(D%"'$$''$"$' ''%%$%##ASCXNCXMCXOCXGCNMCMDZ24YTZUU‘(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1208:7802:34643`($!$(D(&$$)($&&''&''%%%%%ASCXNCXMCXOCXGCNMCMDZ22YTZUU‘(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2116:7778:29884`($!$(D))))))))('''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ22YTZUU‘(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2201:2018:63939P($!$(D)'&&'((#'%#%'#%%%%"""ASCXNCXMCXOCXGCNMCMDZ21YTZUU‘(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2205:6440:23480@($!$(D)))))))'''''%%%%#"!!ASCXNCXMCXOCXGCNMCMDZ20YTZUU’)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1104:10597:245130ABAA))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ19YTZUU’)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1107:11301:46632@ABAA))))))('%''&%%%%%""!ASCXNCXMCXOCXGCNMCMDZ20YTZUU’(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1116:7714:37811PABAA!&(!((&" ''%'%#%%%"""ASCXNCXMCXOCXGCNMCMDZ21YTZUU’)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1201:14408:333510ABAA))))))'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ19YTZUU’)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1212:15336:85329PABAA%$(!$&&!%##%%%%##ASCXNCXMCXOCXGCNMCMDZ21YTZUU’(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1213:8794:577530ABAA())(()%''!##%%%#ASCXNCXMCXOCXGCNMCMDZ19YTZUU◊(*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2208:8325:13729PABAA)))))(()'''''%%%%%"""ASCXNCXMCXOCXGCNMCMDZ21YTZUU€) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2103:3275:100501PBAA"%&)))(&''&''%%%#%"ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU›) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2104:18900:88114PAA ""%%%%%'''''))())))(ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU›) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1112:15391:75683AA(B))))))))))))'''''%%%%%"""ASCXSc˙XNCXMCXOCXGCNMCMDZ25YTZUU›)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1211:21140:394260AA')()))'''''!%%%%""ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUUfi( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1209:7509:61608P$(DAB""%%%%%'&&'')))))())ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUUfi)#IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1101:19285:54092@$(DAB(&())((''''%%%%%%"ASCXSc˙XNCXMCXOCXGCNMCMDZ20YTZUU) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1106:17172:28460PA(B%#$#%''''&)))()))&ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU)"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1214:13003:11554@A(B"""%%%%%'''''))()))(ASCXSc˚XNCXMCXOCXGCNMCMDZ20YTZUU( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2213:6335:24825PA(B"""%%%%%'''&'))))))))ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUU (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1212:7719:35597@((B #$%%' %'&'(&!!&&(ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2105:12730:197780((B""%%%%%%''''&((&&'ASCXSCXNCXMCXOCXGCNMCMDZ19YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1104:10693:56496P!(AD( """%%%%%''''&))()))))ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1209:14997:555950!(AD("""%%%%%'''''(())))ASCXSCXNCXMCXOCXGCNMCMDZ19YTZUU (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2110:11871:3278P!(AD( """%%%%#'''''()')))))ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2114:18299:89262P!(AD( #####'''!%$&(&&ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2206:13913:17841P!(AD( """%%%%#'''''))))))))ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2206:20333:22711@!(AD("""%%%%#'%'''()(()))ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2106:4675:75715@"!!)))()))'''''%%#%%"!"ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU< )"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2101:16869:10135 BD(H"B)(((&'!'!&%%#%%ASCXSc˚XNCXMCXOCXGCNMCMDZ18YTZUU@ ) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2213:17512:13789PD(H"B!B ""%%#%%'''''()))))))ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUUG )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1104:16569:86225@$"$!B#"'#'&"#%##!#ASc˚XSc˚XNCXMCXOCXGCNMCMDZ11A8YTZUUG (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1116:5607:77564@$"!!B&&($(((%''###%###ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUUG )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1203:10461:40370`$"!!B&%#$!&&%%''!##%%%%ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUUG (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2102:8390:14238@$"!!B(&&%(#'#%##'#%%%%ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUUO )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1206:18423:84555`!!B!B('$'%('"$&&##' %###ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUUT (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1103:3897:51516@!B("("H)))))))'''''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUUT )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2105:19766:933400!B("("@((%))(%%%'&###%%"ASCXSCXNCXMCXOCXGCNMCMDZ19YTZUUi ) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2115:20101:40407pH(B(DB""%%%%%'&%&'%((&'&((((ASCXSc˚XNCXMCXOCXGCNMCMDZ23YTZUUr *IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2111:10371:100573P!BH(H(@&&&$!$(&'&'''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ21YTZUU ("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2211:1274:32249@!BB!((&%()&%''&'%%%##!!ASCXSc˚XNCXMCXOCXGCNMCMDZ20YTZUU≠ ) IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1103:20055:27875!H"!((''())))(('''''%%%%%!!ASCXSc˚XNCXMCXOCXGCNMCMDZ24YTZUU )"IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2206:11020:48210 "!!)))))'''''%%%%%!!!ASCXSc˚XNCXMCXOCXGCNMCMDZ18YTZUU (IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2115:12981:3475`"(((B&&&'&'$ %%!'%#%##ASCXSCXNCXMCXOCXGCNMCMDZ22YTZUU )IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2116:11206:58814@"A&$(&&&)'''%'%%%%%"""ASc˚XSc˙XNCXMCXOCXGCNMCMDZ0T19YTZUU
)*IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1206:19932:98308p$H!!"""%%%%%''''')())))))))ASCXNCXMCXOCXGCNMCMDZ23YTZUU
("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1109:1943:983830!!""""%%%%%'''''))))()ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUU
("IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2108:7758:743890!!""""%%%%%'''''))))))ASCXSc˚XNCXMCXOCXGCNMCMDZ19YTZUU
(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2111:2705:390750B!!"## ""#%$ &ASc˚XNCXMCXOCXGCNMCMDZ0T18YTZUU"
)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2213:20007:56613"!$###!%##%&'$%$&&ASCXSCXNCXMCXOCXGCNMCMDZ24YTZUU.
(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1204:5967:66521@!$A))(('))'''''%%%%#!!ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU?
(IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2111:6446:99235@""AH(())))(''#&'%%%%%"""ASc˙XSc˙XNCXMCXOCXGCNMCMDZ0C19YTZUU^
)IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2106:11632:21600@($")((((()'&'''%%%%%"""ASCXSCXNCXMCXOCXGCNMCMDZ20YTZUU_
( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:2212:8365:29810†HH(!"######%''%((&'&$'&(&'$&ASCXSc˚XNCXMCXOCXGCNMCMDZ26YTZUUs
( IˇˇˇˇˇˇˇˇJLK5VL1:383:H0TDTADXX:1:1113:5678:26750P"B(%#%%%'&%%'($ $!#'ASCXSc˚XNCXMCXOCXGCNMCMDZ21YTZUUz

These are just snippets, since the whole file is pretty big
jazz is offline   Reply With Quote
Old 07-01-2013, 08:04 AM   #6
jazz
Member
 
Location: us

Join Date: Nov 2008
Posts: 28
Default

I tried it again and this time it gave a sane output like this

scaffold68 3144 3165 JLK5VL1:383:H0TDTADXX:1:1116:17468:54631 42 +
scaffold68 3162 3184 JLK5VL1:383:H0TDTADXX:1:1205:10063:18561 42 -
scaffold68 3213 3232 JLK5VL1:383:H0TDTADXX:1:2211:8864:24912 42 -
scaffold68 3229 3250 JLK5VL1:383:H0TDTADXX:1:2103:11115:20978 42 -
scaffold68 3230 3250 JLK5VL1:383:H0TDTADXX:1:2210:10068:78079 23 -
scaffold68 3325 3346 JLK5VL1:383:H0TDTADXX:1:2204:21149:40214 24 +
scaffold68 3352 3373 JLK5VL1:383:H0TDTADXX:1:2114:18637:46559 24 -
scaffold68 3352 3373 JLK5VL1:383:H0TDTADXX:1:2216:4275:35276 42 -
scaffold68 3367 3389 JLK5VL1:383:H0TDTADXX:1:1102:12483:97822 24 +
scaffold68 3405 3427 JLK5VL1:383:H0TDTADXX:1:1215:17603:90243 42 -
jazz is offline   Reply With Quote
Old 11-04-2013, 07:38 AM   #7
lran2008
Member
 
Location: quebec

Join Date: Apr 2013
Posts: 35
Default

Hi, have you found out how to do it? I am also having this problem.
lran2008 is offline   Reply With Quote
Old 11-04-2013, 12:31 PM   #8
shi
Wei Shi
 
Location: Australia

Join Date: Feb 2010
Posts: 235
Default

You may try the featureCounts program included in the Subread package.
shi is offline   Reply With Quote
Old 11-04-2013, 11:04 PM   #9
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 311
Default

Quote:
Originally Posted by jazz View Post
Here is the command I tried:

intersectbed -a Inputfile.bed -b compare.gff > output.bed

The first few lines of the resulting bed file look like this:

BAM◊%@HD VN:1.0 SO:unsorted
@SQ SN:scaffold16 LN:1343
@SQ SN:scaffold6 LN:1021
@SQ SN:scaffold4 LN:958
@SQ SN:scaffold22 LN:923
Hi- I suspect your "Inputfile.bed" is not a bed file as the name suggests but a bam file. Could you post the output of:

Code:
head Inputfile.bed 
head compare.gff
If Inputfile.bed is indeed a bam file you should use something like (depending on exactly what output you want):
Code:
intersectBed -abam Inputfile.bam -b compare.gff -bed -wb > results.bed
dariober is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:28 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO