Question: Partial amplicon region missing from mapping data
0
gravatar for lakhujanivijay
10 months ago by
lakhujanivijay4.8k
India
lakhujanivijay4.8k wrote:

I have a 1 kb amplicon sequenced on Ion Torrent platform that belongs to a virus. Statistics are given below

file            virus_genome.fa     data.fastq
format          FASTA               FASTQ
type            DNA                 DNA
num_seqs        1                   60231
sum_len         9181                9416001
min_len         9181                25
avg_len         9181                156.3
max_len         9181                344

I have mapped those reads on to the virus genome using bwa with 88% mapping, but those reads only mapped to around 600-700 bps out of the entire genome. Look at below image from the IGV viewer.

Screenshot-from-2019-04-05-11-04-34

I am wondering why there is no mapping covering the entire 1 kb amplicon; i.e. I am still missing ~ 400bps from my 1kb amplicon. Is it something to do with the wetlab or is it about the mapping (soft clipped bases somewhere) ?

Any help is much appreciated!

igv bwa mapping virus • 344 views
ADD COMMENTlink modified 10 months ago by swbarnes27.5k • written 10 months ago by lakhujanivijay4.8k
0
gravatar for swbarnes2
10 months ago by
swbarnes27.5k
United States
swbarnes27.5k wrote:

Your reads diverge a lot from your reference. At first glance, I'd say you didn't sequence what you thought you sequenced. I'd assemble them and blast that against nr, and see what you really have,

ADD COMMENTlink written 10 months ago by swbarnes27.5k

That's exactly what I did and unfortunately the assembly is very much fragmented. How can I check read quality outside torrent suite? Would fastqc be suffice for that ?

ADD REPLYlink written 10 months ago by lakhujanivijay4.8k

Read quality is probably not your problem. That region is small, make the consensus manually

ADD REPLYlink written 10 months ago by swbarnes27.5k

Any idea how do I proceed about that ?

ADD REPLYlink written 10 months ago by lakhujanivijay4.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1288 users visited in the last hour