Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • #16
    Was your human genome FASTA indexed previously using a different version?

    Comment


    • #17
      Are you running it on a server or a PC, and do you have more than 3 Gb of memory available?

      Comment


      • #18
        I'm running on a PC using VMware wokstation and I allowed 4 gb.

        Comment


        • #19
          Originally posted by AJenkins View Post
          Was your human genome FASTA indexed previously using a different version?
          No. I indexed it with the same version.

          Comment


          • #20
            @FrankiB: Just want to confirm that you are using 64-bit unix of some kind?

            How much /tmp space have you allocated? Is adequate disk space allocated for the VM instance?

            Comment


            • #21
              Originally posted by GenoMax View Post
              @FrankiB: Just want to confirm that you are using 64-bit unix of some kind?

              How much /tmp space have you allocated? Is adequate disk space allocated for the VM instance?
              Yes I'm using a 64-bit ansd 250 gb allocated for VM.

              Comment


              • #22
                Can you increase the RAM allocated to VM (6 or 8 GB) if possible and try? Allocating 4 GB to the VM may mean that a single process may not be able to get ~3.2 GB, which is what you need for an alignment job: http://bio-bwa.sourceforge.net/bwa.shtml#8

                Comment


                • #23
                  Originally posted by GenoMax View Post
                  Can you increase the RAM allocated to VM (6 or 8 GB) if possible and try? Allocating 4 GB to the VM may mean that a single process may not be able to get ~3.2 GB, which is what you need for an alignment job: http://bio-bwa.sourceforge.net/bwa.shtml#8
                  Yes I tried to allocated 7 gb and I got the same result...

                  Comment


                  • #24
                    Can you try to capture the standard error output to see if there is any additional useful information? The exact procedure would differ based on the shell you are using (e.g. http://www.cyberciti.biz/faq/redirec...err-to-stdout/)

                    Comment


                    • #25
                      Originally posted by GenoMax View Post
                      Can you try to capture the standard error output to see if there is any additional useful information? The exact procedure would differ based on the shell you are using (e.g. http://www.cyberciti.biz/faq/redirec...err-to-stdout/)
                      I got about the same information:
                      @SQ SN:1 LN:249250621
                      @SQ SN:2 LN:243199373
                      @SQ SN:3 LN:198022430
                      @SQ SN:4 LN:191154276
                      @SQ SN:5 LN:180915260
                      @SQ SN:6 LN:171115067
                      @SQ SN:7 LN:159138663
                      @SQ SN:8 LN:146364022
                      @SQ SN:9 LN:141213431
                      @SQ SN:10 LN:135534747
                      @SQ SN:11 LN:135006516
                      @SQ SN:12 LN:133851895
                      @SQ SN:13 LN:115169878
                      @SQ SN:14 LN:107349540
                      @SQ SN:15 LN:102531392
                      @SQ SN:16 LN:90354753
                      @SQ SN:17 LN:81195210
                      @SQ SN:18 LN:78077248
                      @SQ SN:19 LN:59128983
                      @SQ SN:20 LN:63025520
                      @SQ SN:21 LN:48129895
                      @SQ SN:22 LN:51304566
                      @SQ SN:X LN:155270560
                      @SQ SN:Y LN:59373566
                      @SQ SN:MT LN:16569
                      @SQ SN:GL000207.1 LN:4262
                      @SQ SN:GL000226.1 LN:15008
                      @SQ SN:GL000229.1 LN:19913
                      @SQ SN:GL000231.1 LN:27386
                      @SQ SN:GL000210.1 LN:27682
                      @SQ SN:GL000239.1 LN:33824
                      @SQ SN:GL000235.1 LN:34474
                      @SQ SN:GL000201.1 LN:36148
                      @SQ SN:GL000247.1 LN:36422
                      @SQ SN:GL000245.1 LN:36651
                      @SQ SN:GL000197.1 LN:37175
                      @SQ SN:GL000203.1 LN:37498
                      @SQ SN:GL000246.1 LN:38154
                      @SQ SN:GL000249.1 LN:38502
                      @SQ SN:GL000196.1 LN:38914
                      @SQ SN:GL000248.1 LN:39786
                      @SQ SN:GL000244.1 LN:39929
                      @SQ SN:GL000238.1 LN:39939
                      @SQ SN:GL000202.1 LN:40103
                      @SQ SN:GL000234.1 LN:40531
                      @SQ SN:GL000232.1 LN:40652
                      @SQ SN:GL000206.1 LN:41001
                      @SQ SN:GL000240.1 LN:41933
                      @SQ SN:GL000236.1 LN:41934
                      @SQ SN:GL000241.1 LN:42152
                      @SQ SN:GL000243.1 LN:43341
                      @SQ SN:GL000242.1 LN:43523
                      @SQ SN:GL000230.1 LN:43691
                      @SQ SN:GL000237.1 LN:45867
                      @SQ SN:GL000233.1 LN:45941
                      @SQ SN:GL000204.1 LN:81310
                      @SQ SN:GL000198.1 LN:90085
                      @SQ SN:GL000208.1 LN:92689
                      @SQ SN:GL000191.1 LN:106433
                      @SQ SN:GL000227.1 LN:128374
                      @SQ SN:GL000228.1 LN:129120
                      @SQ SN:GL000214.1 LN:137718
                      @SQ SN:GL000221.1 LN:155397
                      @SQ SN:GL000209.1 LN:159169
                      @SQ SN:GL000218.1 LN:161147
                      @SQ SN:GL000220.1 LN:161802
                      @SQ SN:GL000213.1 LN:164239
                      @SQ SN:GL000211.1 LN:166566
                      @SQ SN:GL000199.1 LN:169874
                      @SQ SN:GL000217.1 LN:172149
                      @SQ SN:GL000216.1 LN:172294
                      @SQ SN:GL000215.1 LN:172545
                      @SQ SN:GL000205.1 LN:174588
                      @SQ SN:GL000219.1 LN:179198
                      @SQ SN:GL000224.1 LN:179693
                      @SQ SN:GL000223.1 LN:180455
                      @SQ SN:GL000195.1 LN:182896
                      @SQ SN:GL000212.1 LN:186858
                      @SQ SN:GL000222.1 LN:186861
                      @SQ SN:GL000200.1 [M::main_mem] read 82646 sequences (10000166 bp)...

                      Comment


                      • #26
                        Re-reading your previous posts I realized that in post #15 you said that index generated 2 files. That does not seem right. I have a feeling the indexing is the problem and did not work right the first time. Can you post a listing of index files for the genome? You may also want to try re-indexing.

                        There should be 5 index related files (*.amb, *.sa, *.bwt, *.ann, *.pac) besides the fasta file.
                        Last edited by GenoMax; 01-29-2014, 08:45 AM.

                        Comment


                        • #27
                          Sorry to have mistaken you, but I have generated the 5 files you are refering to.

                          Comment


                          • #28
                            Is it working now? You deleted a previous post that said that "aln" option was working but "mem" was not.

                            Comment


                            • #29
                              Yea. bwa aln is working well (it seems). Now I'm running bwa sampe to combine the *.sai files.

                              Comment


                              • #30
                                BWA.tar installation on windows

                                Hello Everyone!
                                I am looking a best software to align paired end reads from illumina MiSeq. Someone told me that I can use BWA-MEM or Subread to align paired end reads to the reference genome. I downloaded the BWA-MEM but I can't execute that because my computer has windows 32 bit version. the installation file has .tar.gz extension. Can anyone tell me how to open the execute the software on windows?
                                Many thanks!
                                MZB

                                Comment

                                Latest Articles

                                Collapse

                                • seqadmin
                                  Strategies for Sequencing Challenging Samples
                                  by seqadmin


                                  Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                                  03-22-2024, 06:39 AM
                                • seqadmin
                                  Techniques and Challenges in Conservation Genomics
                                  by seqadmin



                                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                                  Avian Conservation
                                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                                  03-08-2024, 10:41 AM

                                ad_right_rmr

                                Collapse

                                News

                                Collapse

                                Topics Statistics Last Post
                                Started by seqadmin, Yesterday, 06:37 PM
                                0 responses
                                11 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, Yesterday, 06:07 PM
                                0 responses
                                10 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-22-2024, 10:03 AM
                                0 responses
                                51 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 03-21-2024, 07:32 AM
                                0 responses
                                68 views
                                0 likes
                                Last Post seqadmin  
                                Working...
                                X