Hi,
I have been trying to fetch the reads which are being partially mapped to the genome. Like if the total length of a sequence is 126 bases and only 66 bases out of 126 are mapping to the genome. Is it possible to fetch those sequences whose half of the bases align to the genome? I created a fake fastq file having just 4 sequences to see how a mapper is aligning the reads. I have pasted those 4 sequences below. So far I have used Bowtie2, Hisat2 and STAR. None of them can give me hits if half of the sequence is mapping to the genome.
The fake fastq file I created has following sequences (below) and I am trying to map them to Giardia lamblia genome.
@seq-1_Whole_sequence_is_Copied_from_Genome
CTTGCGGCAAATTGTGAAACAGAAGGCAATAAAAGCAATTGTGACAGCGGTCGCTGTGAAATGGTTGGCATTGTTGAAGTCTGTACTCAGTGCAAGACAAATGGACACGTCCCTATTGACGGAGTA
+
HHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJDD9?BDDDDADCDDDDBCCCFFFF
@seq-2_Left_77_bases_are_copied_from_the_genome
TCTTGCGGCAAATTGTGAAACAGAAGGCAATAAAAGCAATTGTGACAGCGGTCGCTGTGAAATGGTTGGCATTGTTATGACGATGACGATGACGATGACGTAGACGATGACGTAGCAGTAGACGAT
+
CCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDD
@seq-3_Right_77_bases_are_from_genome
GAGTGACGTAGCGATGCAGTGACGTAGCGGAGTACCAGATGATCACAAAGTCGCTGTGAAATGGTTGGCATTGTTGAAGTCTGTACTCAGTGCAAGACAAATGGACACGTCCCTATTGACGGAGTA
+
CCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDD
@seq-4_Two_different_parts_of_genome_spliced_together
TCTTGCGGCAAATTGTGAAACAGAAGGCAATAAAAGCAATTGTGACAGCGGTCGCTGTGAAATGTCCGGGCCAGCTAGGCCCATCGACATACCTCGTTGACCGGGGAAACAGCATGCCAGGACAGG
+
CCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDDBDBDDD9?BDDDDADCDDDDBCCCFFFFFHHHHHJJJIJIJHFFDCEDD
Do you know any parameter in Bowtie2, HSAT2, STAR etc which can give me hits for these type of sequences?
Thanks in Advance, Swadha
0 answers
No answers yet.
Log in to answer this question.
Can you post the command line options you used for those alignment programs?
I also suggest that you give
bbmap.shfrom BBMap suite a try. Lots of options and flexibility. Guide available here.I used default parameters. I was not sure which parameters can give hits for such kind of sequences. Like for Hisat I used:
The output I got is:
Which problem are you working on? What is the origin of these reads which align only half?
I am trying to find mRNA reads which are trans-spliced (this results when protein-coding regions from multiple RNAs transcribed from different loci can be joined). I have created a fake dataset so see if Bowtie2, HISAT2 or STAR can map these type of trans-spliced reads. I can move forward with the partial mapping as well.
Information like that is quite important, so please be as complete as possible when asking questions. So maybe something like https://genomebiology.biomedcentral.com/articles/10.1186/gb-2014-15-2-r34? Perhaps it is not necessary to reinvent the wheel.