Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

add support for paired-ends #35

Open
JoseBlanca opened this issue Feb 28, 2011 · 0 comments
Open

add support for paired-ends #35

JoseBlanca opened this issue Feb 28, 2011 · 0 comments
Labels

Comments

@JoseBlanca
Copy link
Owner

From Giuseppe D'Auria

para que funcione el bwa con paired-ends tengo que hacer eso:

indexando la secuencia de referencia

bwa index -a bwtsw sequence.fasta

alineando los dos ficheros de forward y reverse

bwa aln sequence.fasta NG-5400_Y3_read_1.fastq> reads_y3_1.sai
bwa aln sequence.fasta NG-5400_Y3_read_2.fastq> reads_y3_2.sai

mapeando los dos ficheros

bwa sampe -f aln_y3.sam sequence.fasta reads_y3_1.sai reads_y3_2.sai
NG-5400_Y3_read_1.fastq NG-5400_Y3_read_2.fastq

creando el bam

samtools view -b -S aln_y3.sam> aln_y3.bam

sorteando el bam

samtools sort aln_y3.bam aln_y3.sorted

indexando el bam

samtools index aln_y3.sorted.bam

Tendriamos que tener en el backbone_analysis una opción para que el
pipeline funcione de esta manera.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant