|
From: Juan D. M. C. <jdm...@gm...> - 2016-10-10 21:15:29
|
Hi Rebeca, I think the merge command expects a coordinate sorted bam and you are sorting it by query name. You should remove the "-n" flag in your sort command and then merge the sorted files. Regards, Juan Montenegro On 11 Oct 2016 6:52 AM, "Rebecca Harris" <re...@gm...> wrote: > Hi, > > I have a couple of genomes, aligned to the same reference, that I am > trying to merge using the command: > > java -jar picard.jar MergeSamFiles MERGE_SEQUENCE_DICTIONARIES=true > I=NA18867.bam I=Denisovan.bam O=out.bam > > Upon merging, I get an error that my BAM files are in different orders - > in particular 2 and 19. So, I attempted to use samtools sort -n to sort all > my BAM files. However, this is not working as expected. After sorting, I > used samtools view -H to look at the header of my sorted BAM files (results > shown below) and find that these genomes continue to be in different orders. > > I am using samtools version 1.2 and the following command for sorting: > samtools sort -n -T human -o humans/NA18867.sorted.bam -O bam -@ 12 > humans/NA18867.bam > > Any help would be greatly appreciated. Thanks! > Rebecca > > > For one genome I get: > > @HD VN:1.0 SO:queryname > > @SQ SN:1 LN:249250621 > > @SQ SN:10 LN:135534747 > > @SQ SN:11 LN:135006516 > > @SQ SN:12 LN:133851895 > > @SQ SN:13 LN:115169878 > > @SQ SN:14 LN:107349540 > > @SQ SN:15 LN:102531392 > > @SQ SN:16 LN:90354753 > > @SQ SN:17 LN:81195210 > > @SQ SN:18 LN:78077248 > > @SQ SN:19 LN:59128983 > @SQ SN:2 LN:243199373 > > and for the other genome I get: > > @HD VN:1.0 GO:none SO:queryname > > @SQ SN:1 LN:249250621 AS:NCBI37 UR:ftp://ftp.1000genomes.ebi. > ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz M5: > 1b22b98cdeb4a9304cb5d48026a85128 SP:Human > > @SQ SN:2 LN:243199373 AS:NCBI37 UR:ftp://ftp.1000genomes.ebi. > ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz M5: > a0d9851da00400dec1098a9255ac712e SP:Human > > @SQ SN:3 LN:198022430 AS:NCBI37 UR:ftp://ftp.1000genomes.ebi. > ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz M5: > fdfd811849cc2fadebc929bb925902e5 SP:Human > > @SQ SN:4 LN:191154276 AS:NCBI37 UR:ftp://ftp.1000genomes.ebi. > ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz M5: > 23dccd106897542ad87d2765d28a19a1 SP:Human > > > ------------------------------------------------------------ > ------------------ > Check out the vibrant tech community on one of the world's most > engaging tech sites, SlashDot.org! http://sdm.link/slashdot > _______________________________________________ > Samtools-help mailing list > Sam...@li... > https://lists.sourceforge.net/lists/listinfo/samtools-help > > |