BuildBamIndex Error - results from all directories
0
0
Entering edit mode
5.4 years ago
woojoy14 ▴ 10
gatk BuildBamIndex --INPUT 05_DeDuplicate/${sample}_dedup.bam && cp
/home/ubuntu/00_NA12878/00_References/chr.list ./ &&
for i in `cat chr.list`; do samtools view -bh 
05_DeDuplicate/${sample}_dedup.bam ${i} > 
05_DeDuplicate/${sample}_dedup.${i}.bam ; done  23 for i in `cat chr.list`; 
do gatk BuildBamIndex --INPUT 05_DeDuplicate/${sample}_dedup.${i}.bam; done



  >

NA12878_dedup.00_RawData.bai
NA12878_dedup.00_RawData.bam
 NA12878_dedup.00_References.bai
NA12878_dedup.00_References.bam
NA12878_dedup.01_FASTQC.bai
NA12878_dedup.01_FASTQC.bam
NA12878_dedup.02_Trimmomatic.bai
NA12878_dedup.02_Trimmomatic.bam
NA12878_dedup.03_BWA.bai
NA12878_dedup.03_BWA.bam
NA12878_dedup.04_Sort.bai
NA12878_dedup.04_Sort.bam
NA12878_dedup.05_DeDuplicate.bai
NA12878_dedup.05_DeDuplicate.bam
NA12878_dedup.bai
NA12878_dedup.bam

I run the code above and got the .bam and .bai files of all directories existing in my folder including results by chromosomes which is what I would like to get only in the first place.

I made vi file for this and run it as bash ./{file name} {sample name}

One thing I'm doubting is the name of my root directory started with the same name of sample name. But I indicate the location as INPUT 05_DeDuplicate/${sample}_dedup.bam in the first line of the code. I have no idea why this happened so far.

Do you see something wrong with the code above or is it something about the sample name?

I will really appreciate your help!

BuildBamIndex • 468 views
ADD COMMENT

Login before adding your answer.

Traffic: 1723 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