Deleted:How to convert BAM to paired FASTQ from Nanopore data?
2
0
Entering edit mode
3.0 years ago
ej6474 ▴ 10

I have a set of coordinate sorted, mapped BAM files for SARS-CoV-2 sequencing data through Nanopore MinION. that I need to generate paired fastq files for.

When I traditionally name sort the BAM files, run samtools fixmate and then bedtools bamtofastq, my paired fastq files are empty. When I run samtools flagstat on my bam files, I see that read1 and read 2 are 0:

52244 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 secondary
737 + 0 supplementary
0 + 0 duplicates
52244 + 0 mapped (100.00% : N/A)
0 + 0 paired in sequencing
0 + 0 read1
0 + 0 read2
0 + 0 properly paired (N/A : N/A)
0 + 0 with itself and mate mapped
0 + 0 singletons (N/A : N/A)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

I am not fully familiar with Nanopore data so what is the best way to extract paired reads from my BAM files or are there any intermediate steps I am missing?

nanopore • 2.1k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2150 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6