Export reads from an alignment file
Hi! I have annotated some clustered genes in a very fragmented genome assembly based on Illumina reads. I want to study their synteny (impossible in the assembly I have) and I have a lot of nanopore reads as well. I had aligned the nanopore reads to the Illumina annotated-genes using BWA and I have now an alignment file. I now want to get the COMPLETE SEQUENCE of the reads aligning to my reference (so I can make a local assembly of those reads and hopefully get the synteny) and NOT ONLY THE PART OF THE READS THAT ALIGN TO THE REFERENCE. Is it possible to do that? How can I do it? Thanks a lot!
• 2,038 views
•
link
written
by
elena.syp •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Genome assembly reference .fasta + Illumina reads
written by A_heath •
-
Tools for Reference-Guided Assembly Using Nanopore Sequencing Reads
written by pratikshaadagale39722 •
-
Isoforms and Transcriptome assembly with direct RNA Nanopore reads
written by mariana.villegas •
-
Plasmid assembly of Oxford nanopore sequencing reads
written by MG_19 •
-
MUMmer alignments: Trying to understand the average identity result
written by caro-ca •
-
Improvement of genome assembly using illumina contigs and nanopore reads
written by KG •
-
Can any long read assemblers incorporate both nanopore and pacbio reads in the same run?
written by O.rka
-
Replacing reference with aligned reads
written by elena.syp •
-
de novo transcriptome assembly of Oxford Nanopore Technologies long reads
written by Alicia •
-
Contamination in assembly
written by elena.syp •
samtools view | cut?? #oramimissingsomething ?As per Pierre, you can just use
samtools view -b -f 0x2 MyBAM.bamin order to extract just the reads that have aligned successfully to your genome (assuming paired-end). You can then use thecutcommand in combination in order to pull out information relating to read positions if you want to group them into specific loci prior to performing local assembly.Another way to extract aligned reads over a defined genomic region is, for example: