cn.MOPS: Error in getReadCountsFromBAM
0
0
Entering edit mode
@laurabuggiotti-11728
Last seen 5.6 years ago

Hi,

Im trying to run cn.MOPS using multiple bam files of different breeds and while I can successfully run it by breed, when attempting to use all of them i get the following error:

Error in getReadCountsFromBAM(BAMFiles, refSeqName = readLines("/users/lbuggiotti/WGresRussia/cn.MOPS/Chr_list.txt"),  : 
  Targets in Header file of  /users/lbuggiotti/bamKorea/bam_NCBI/SRR1262755.bam  are not identical to the header of the file /users/lbuggiotti/bamKorea/A01.realigned.rmdup.bam !
In addition: Warning messages:
1: In sn == names(targets[[i]]) :
  longer object length is not a multiple of shorter object length
2: In sl == as.integer(targets[[i]]) :
  longer object length is not a multiple of shorter object length
Execution halted

I have checked the header of the above mentioned bam files and they are different, however i dont know how to change them in order to look alike and run all of them together. Some of the bam files have been downloaded from ncbi, however the reference genome used is the same in all of them.

I hope you can help me with the issue.

Thanks a lot for your time and support,

Laura

cn.mops • 959 views
ADD COMMENT
0
Entering edit mode

I also find out about this, hoping to find the best answer! hole io

ADD REPLY

Login before adding your answer.

Traffic: 409 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6