View Single Post
Old 01-12-2012, 01:42 AM   #3
Junior Member
Location: India

Join Date: Jan 2012
Posts: 6


I generated the SAM file using bowtie version0.12.7 but i am not able to generate a proper bcf file.BCF file is genertaed within 2 minutes and normally it will take time.The commands i gave is

time ./fastq-dump SRR068289.sra

I have falidated the fastaq using fastaqvalidator and it was successful

time ./bowtie -k 2 -v 2 hg19 -q /store/data/amruta/sratoolkit.2.1.8-centos_linux64/12JAN/SRR068289.fastq -S > /store/data/amruta/sratoolkit.2.1.8-centos_linux64/12JAN/bowtie1-amruta.sam

From sam file bam was generated and using samtools bcf and vcf were generated
./samtools view -bS bowtie1-amruta.sam > bowtie1-amruta.sam
samtools sort filename.bam sorted
samtools faidx input.reference.fa
samtools index sortedhg19.bam
samtools mpileup -d 8000 -uf input.reference.fa my.sortedhg19.bam | bcftools view -bvcg -> my.raw.bcf
bcftools view my.raw.bcf | perl varFilter -D8000 > my.flt.vcf

My fastaq file is a huge in size.So i tried with bowtie 2 beta version.But i am not able to run bowtie2

./bowtie2 -x hg19 -u file.fataq -S out.sam

The error was "Error: Encountered internal Bowtie 2 exception (#1)"
Can anyone please help is offline   Reply With Quote