This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Mapping of tag sequences to the reference genome in Tassel 5 GBS v2 Pipeline

I have a question regarding appropriate bowtie2 settings for mapping of tag sequences to the reference genome in Tassel 5 GBS v2 Pipeline. According to the Tassel 5 GBS v2 Pipeline manual (enter link description here) i should run bowtie2 with the following settings:

bowtie2 -p 15 --very-sensitive -x ../AGPv3/ZMAGP3_9103 -U tagsForAlign.fq -S tagsForAlignFullvs.sam

However, i found that tag sequences in FASTQ file have old Phred +64 quality encoding. In particular, all the bases in the FASTQ file have the same f letter in the quality string (the 4-th line). My question is : Should i add --phred64 option to the my bowtie2 command? So the correct command should look like:

bowtie2 -p 15 --phred64 --very-sensitive -x ../AGPv3/ZMAGP3_9103 -U tagsForAlign.fq -S tagsForAlignFullvs.sam
genotyping genotyping-by-sequencing tassel tassel_gbs_pipeline bowtie2
• 930 views • link
written by Denis Denis

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Could not locate a Bowtie index corresponding to basename
    written by Rahul • Rahul
    Dear all, Can some one suggest me how to correct this problem. I am running wsl based Oligominer tool for probe designing. Th ecommand after …
  • rsem-calculate-expression error
    written by Bhoomika • Bhoomika
    **I have prepared reference using the command:** rsem-prepare-reference -p 5 --bowtie2 cd_hit_est.fasta ref_output_12 **And the files i obtained are:** ref_12.1.bt2 ref_12.2.bt2 ref_12.3.bt2 ref_12.4.bt2 ref_12.grp ref_12.idx.fa …
  • error with bowtie2 alignment
    written by pt.taklifi pt.taklifi
    Hi I am trying to align a fastq file to mitochondrial genome with bowtie 2 with the following command and store unmapped read to folder …
  • bowtie2 aligment problems
    written by amitpande74 • amitpande74
    Hi, I have reads from DNA sequencing,5 pools. The reads look like this: @NS500455:80:HG7TNBGXB:1:11101:17723:1055 1:N:0:ATCACG ACTTANGTGTATGTAAACTTCCGACTTCAACTGTATAGGGATCCNAGCTCCAATTCGCCCTATAGTGAGTCGTAT + /AAAA#EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE#EEEEEEEEEEEEEEEEEEEEEEEEEEEEEE It's illumina primer followed by the construct …
  • Unexpected Bowtie2 behavior with --phred64 and --phred33 flags
    written by Denis Denis
    I'm stiil not able to figure out how `Bowtie2` deal with different `Phred quality` encoding. In my recent post: [enter link description here][1] [1]: https://www.biostars.org/p/381067/ …
  • Analysis of paired-end GBS data
    written by Denis Denis
    Hi there! I'm working on GBS data analysis of plants with available reference genome. I already tried Tassel 5 GBS v2 pipeline (https://bitbucket.org/tasseladmin/tassel-5-source/wiki/Tassel5GBSv2Pipeline) and it …
  • RSEM ignores quality scores?
    written by AW AW
    Dear all, I want to quantify gene expression for my Trinity transcriptome assembly using RSEM. I am getting confusing results depending on how I specify …
  • Bowtie2 error with single end mapping
    written by H.Hasani H.Hasani
    Dear all, I'm facing a tedious error and not sure what causing it. I have single-end reads and simply trying to align it using Bowtie2 …
  • Are the BAM files from tophat and bowtie-samtools the same?
    written by jolin0701-dy jolin0701-dy
    I know we can get accepted_hits.bam from tophat like this: tophat2 --max-multihits 1 --b2-very-sensitive --GTF mm10.gtf --output-dir CZ-1_thout genome CZ1.fastq And I know we can …
  • miRNA mapping rate is very low.. (less than 0.03%)
    written by illinois.ks illinois.ks
    Hello, I haven't worked with miRNA mapping..( the specie is mouse) It is the first time to do it. I have read a lot of …