Hello,
I was given fastq files of 150 bp illumina sequenced reads after PCR was performed on the region of interest and I was given a reference sequence of the gene. I used BWA mem to align my sequences to the reference and only got 60% mapped. I was wondering if I am using the correct algorithm for PCR produced reads and some explanation as to what could cause not all of them to map? I did not adapter trim because BWA does soft-clipping and I know some reads did sequence into the adapters because they were short. Could someone else explain as to why not all would map? Could this be PCR/sequencing error or am I doing something wrong?
Thank You
1 answer
Did you check to see if there are phiX sequences in your data? Often times phiX is spiked in to normalize the base distribution when sequencing a low nucleotide diversity sample (amplicons would qualify for this). Generally phiX data is removed during demultiplexing but if your samples were not multiplexed phiX reads may still be there. If that is not the case grab a few reads (convert them to fasta format) that do not map and then blast them at NCBI to find out what they are.
Log in to answer this question.
-
clipping penalty
written by faravarkhordadpoor •
-
cutadapt: how to trim adapter dimers?
written by Ana •
-
Is the RNA- sequencing quality too bad for analysis?
written by Saran •
-
which sequence adapter trimming tool should i use for my illumina 101pb reads ?
written by v.berriosfarias •
-
Tool: Trimviz, a lightweight QC tool to assess impact of read-trimming and aligner soft-clipping
written by stuart archer •
-
Does Cutadapt change the x,y coordinates in the fastq read name?
written by drkennetz
-
Soft-clipping of reads in Amplicon-sequenced data
written by jsneaththompson
-
Extract soft-clipped reads from BWA -generated SAM file
written by ThePresident
-
How can I deal with adapter contamination in next-gen sequencing reads?
written by gbdias
-
Fastx Clipper Not Accurate In Adapter Clipping For 100Nt Reads
written by suzanne.mcgaugh
