Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Paired-end reads stored in a single-end read library #58

Open
dzolier opened this issue Aug 10, 2021 · 1 comment
Open

Paired-end reads stored in a single-end read library #58

dzolier opened this issue Aug 10, 2021 · 1 comment

Comments

@dzolier
Copy link

dzolier commented Aug 10, 2021

QC'ed the reads with illumina-utils (specifically iu-filter-quality-minoche and iu-merge pairs in two different attempts), performed preliminary binning with MegaHIT; feedback at the end was as in the title. Example:

Binsanity-profile -i minoche_CONTIGS/contigs.fa -s minoche_MAPPING/ -c Binsanity_minoche/MEGAHIT_coverage --transform scale -T 1 -o Binsanity_minoche_2/

    ******************************************************
                Contigs formated to generate counts
    ******************************************************


    ==========     _____ _    _ ____  _____  ______          _____
    =====         / ____| |  | |  _ \|  __ \|  ____|   /\   |  __ \
      =====      | (___ | |  | | |_) | |__) | |__     /  \  | |  | |
        ====      \___ \| |  | |  _ <|  _  /|  __|   / /\ \ | |  | |
          ====    ____) | |__| | |_) | | \ \| |____ / ____ \| |__| |
    ==========   |_____/ \____/|____/|_|  \_\______/_/    \_\_____/
      v2.0.2

//========================== featureCounts setting ===========================\
|| ||
|| Input files : 1 BAM file ||
|| ||
|| SRR492065.bam ||
|| ||
|| Output file : SRR492065.bam.readcounts ||
|| Summary : SRR492065.bam.readcounts.summary ||
|| Paired-end : no ||
|| Count read pairs : no ||
|| Annotation : SRR492065.bam.saf (SAF) ||
|| Dir for temp files : Binsanity_minoche_2 ||
|| ||
|| Threads : 1 ||
|| Level : meta-feature level ||
|| Multimapping reads : counted ||
|| Multi-overlapping reads : counted ||
|| Min overlapping bases : 1 ||
|| ||
\============================================================================//

//================================= Running ==================================\
|| ||
|| Load annotation file SRR492065.bam.saf ... ||
|| Features : 1546 ||
|| Meta-features : 1546 ||
|| Chromosomes/contigs : 1546 ||
|| ||
|| Process BAM file SRR492065.bam... ||
ERROR: Paired-end reads were detected in single-end read library : minoche_MAPPING/SRR492065.bam

@Harpreet525
Copy link

You will have to open your Binsanity-profile script and add the parameter "-p" in line 31 to allow the script to forward the command to featureCounts to be able to read paired-end reads. It will work.
Because the script of Binsanity-profile is not prompting the featureCounts script to read paired reads.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants