How to convert fastq to BAM

Hi! I would like to convert fastq file to bam i in oneliner: I have tried:

bowtie2 -x {input.index}/hg19.zip -U {input.fastq} | samtools view -b -o {output}

but it does not work. Anyone know how to do it? Kindly help.

fastq bam

• 2,042 views • link

updated 4.2 years ago by cmdcolin

• written 4.2 years ago by ja4123 •

this is a useful link from the htslib docs https://www.htslib.org/workflow/fastq.html

if you wanted to use bowtie instead of minimap2, you can swap it out

• 0 views • link

written 4.2 years ago by cmdcolin

1 answer

See answer click -> : Piping bowtie2 output directly into BAM

• 0 views • link

written 4.2 years ago by GenoMax

I used: bowtie2 -x {input.index}/hg19.zip -U {input.fastq} | samtools view -bS -> {output}.bam from the site you have sent but I still have an error: Missing input files for rule fastq2bam: file.sam

• 0 views • link

written 4.2 years ago by Jimpix •

but I still have an error: Missing input files for rule fastq2bam: file.sam

it's unrelated to your question. You have a problem with snakemake, not bowtie.

• 0 views • link

written 4.2 years ago by Pierre Lindenbaum

Thanks! I know where was a problem

• 0 views • link

written 4.2 years ago by Jimpix •

More posts like this

samtools idxstats not removing ChrM
written 3.4 years ago by Jen •

I am trying to remove ChrM from my ChIP-seq data. Below is my pipeline for one sample up to where I am having the issue …
Error when runing Bowtie2: (ERR): bowtie2-align exited with value 1
written 3.5 years ago by luzglongoria •

Hi there, I'm running this command with Bowtie2: bowtie2 --threads 20 --no-unal -x /Path/to/Indexed/Genome/file.fna -q -k 1 --al file.fastq -U -b --align-paired-reads input.bam | samtools …
Filter out Bam not overlapping with Bed File (keeping a read and its mate)?
written 3.8 years ago by Eliveri •

I have a Bam file for which I want to REMOVE reads overlapping with a given Bed file (so keep all other reads/alignments). I also …
Difference between commands [bowtie2, samtools]
written 4.2 years ago by Jimpix •

Hi! I would like to ask what is the difference between commands bellow. This wrote in oneliner: bowtie2 -x {input.index} -U {input.fastq} | samtools view …
Alignment fastq files
written 4.7 years ago by Jimpix •

I have a question. I need to align and convert fastq files (unpaired) into bam file. If to beo sure I need to ast first …
Converting CRAM to BAM without reference fasta
written 5.5 years ago by Ld_60 •

Hello, I have a couple of CRAM files which I am aiming to convert to BAM files. I tried with these two commands: samtools view …
How to break script if Bowtie2 errors?
written 6.9 years ago by dfajar2

Hi, I'll try to explain my problem the best I can. I have a script that runs multiple bowtie2 alignments. The problem is that it …
Piping bowtie2 output directly into BAM
written 7.4 years ago by hersh •

I am using bowtie2, and due to the size of the SAM output, I would like to entirely avoid the SAM file, and pipe the …
How to obtain insert sequence from paired-end mapping result (SAM/BAM)
written 9.9 years ago by gundalav

I have paired-end reads R1 and R2. |-----R1-----|----------------------------|-----R2-----| A B C Total insert here is A + B + C. How can I get this …
Are the BAM files from tophat and bowtie-samtools the same?
written 10.3 years ago by jolin0701-dy

I know we can get accepted_hits.bam from tophat like this: tophat2 --max-multihits 1 --b2-very-sensitive --GTF mm10.gtf --output-dir CZ-1_thout genome CZ1.fastq And I know we can …