very low number of reads aligned in bowtie (colourspace)
0
0
Entering edit mode
2.9 years ago
Lila M ★ 1.2k

Hi there, I would like to map a file with the human genome GRCh38. They used Solid technology to generate the file, so what I did was to generate the index as follow:

bowtie-build -C GRCh38.p10.genome.fa human_index_bowtie1/colour_index

and then run bowtie1

bowtie -p 8 -t -C human_index_bowtie1/colour_index SRR11557616.fastq.gz -S  SRR11557616.sam

however, I've got a very low number of reads mapped:

# reads processed: 37684092
# reads with at least one reported alignment: 2864697 (7.60%)
# reads that failed to align: 34819395 (92.40%)
Reported 2864697 alignments

The fastqc analysis was not terrible, Per base sequence quality and kmer have some issues but I don't expect a better result if I do a trimming.

Any suggestion about how to improve the quality of the mapping? Thanks!!!

map colourspace bowtie • 676 views
ADD COMMENT
0
Entering edit mode

You should perhaps use BFAST (LINK) like the SRA entry notes. Not many aligners left that can deal with colorspace.

ADD REPLY
0
Entering edit mode

Thank you! Is there any manual or more information about how to use it?

ADD REPLY
0
Entering edit mode

Program presumably should have in-line help or a doc file. Get the code and see what you find. Since SRA entry says they used this it should work (in theory).

ADD REPLY

Login before adding your answer.

Traffic: 2051 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