SAM file generation problem
my SAM files generated are of 1.1kb each.. also the alignmnt is not done.. what to do??
• 2,691 views
•
link
written
by
kuhelikadas1995 •
1 answer
Change --dta -cufflinks to --dta-cufflinks.
• 0 views
•
link
written
by
mbelmadani
Log in to answer this question.
More posts like this
-
bamdst gives error "EOF marker is absent. The input is probably truncated."
written by shpak.max •
-
Split a sam file by chromosome
written by ICfc97 •
-
BWA MEM output directory
written by sarastrafella.ss •
-
how to check the integrity of a sam file derived from bwa mem?
written by dongguiying2017 •
-
transcript assembling problem by cufflinks
written by kuhelikadas1995 •
-
Merging two fastq.gz files
written by tcf.hcdg
-
How to compare compare two Sam files to check mapping
written by madhu.9124 •
-
bowtie2-align died with signal 6 (ABRT) (core dumped)
written by micro32uvas •
-
File STARtmp remains after mapping
written by katerina192 •
-
Sorted vs unsorted SAM/BAM
written by ofonov •
Please put some more effort into your questions to get proper help.
this was the command used for aligning the adapter trimmed fastq files against human refernce genome. the SRA data here is paired end rna seq data. Command used for adapter trimming is:
In the future when you paste code examples, highlight your code and click the "101 010" button to make it a code block, makes it easier to read.
done. but it didnot work. same problem. no alignmnt done.
If it did not work then why did you accept the answer below?
I accepted the answer nd tried.. then saw that i didnt work..
Can show us the file sizes for your original and trimmed data, as well as your genome? One thing I notice is that in your cutadapt command, the quality filter (
-q 30) is after the input files whereas in the documentation it's typically before.If your trimmed files are suspiciously small I would guess that's the error.
raw data for r1 read is of size 502.3MB and that of r2 is of 498.8MB.. trimmed data for r1 is of 489.6MB and that of r2 is 486.6MB.. moreover, applying
-q 30before or after made no difference. Size of reference genome is 1.1GB(zipped file) and 4.5GB(when unzipped)As I recall HISAT2 is one of those programs that is sensitive to order of program options and input files. So you may want to pay some attention to that order. Did you generate the indexes yourself and are they known to work?
yeah.. and it did generate the indexes for genome by myself.. nothing was mentioned about the adapteer sequence.. it used AGATCGGAAGA.. for both forward and reverse end reads. should i use AGATCGGAAGAGC as the adapter sequence fr both the end reads??
Specifying an incorrect adapter sequence should have no impact on being able to produce a good SAM file from HISAT2. It sounds like that is where your problem lies. Have you looked at the log file for HISAT2 to see if there are any errors that can be addressed?
initially i had used the command