Extract read positions from a bas file and use as a list in Python
I would like to extract all the positions of the reads from a specific contig in a bam file and import them as a list in python for downstream analysis. Can someone recommend how to do that?
• 1,228 views
•
link
updated
by
Ram
•
written
by
bluebubbl19 •
•
written
by
bluebubbl19 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Extract aligned positions from reads
written by María José •
-
Extracting only reads mapped to specific genomic region
written by agata •
-
Tool: Filtering BAM files efficiently
written by ddeemer •
-
Extract quality scores from aligned BAM files
written by elb
-
Pysam VariantFile.fetch() without specifying contig
written by ricardoguerreiro2121
-
Extracting Reads from paired end data in SAM file that map to a list of contig IDs
written by yp19
-
tool for a pattern of repeats
written by 2nelly
-
aligning scaffolds to the reference genome
written by arta
-
Samtools sam-> Bam
written by V
-
Split BAM into 1 BAM per contig
written by Marvin
What have you tried?
Used samtools to extract positions of the reads of the contigs I need, directed that into a csv file and now I am stuck. I tried to read that file using read_csv in python, but it does not work.
I guess the reason for that is this is the first time I am working in Azure, and I am not sure I am indicating the path/link to the file on the drive correctly. I would really appreciate some help in this regard as well.
I'm not familiar with Azure in command line context, you may want to try and fix the problem on a local machine before porting working code to the cloud platform.
As in reads that are aligned to a specific contig reference? Otherwise this does not make sense since there will be no individual reads inside a contig.
Yes, they are aligned to that contig in a bam file