Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • johnsequence
    Member
    • Jan 2010
    • 22

    percentage coverage after alignment

    How to calculate the percentage of alignment match after mapping in FASTQ format? Thanks.
  • imilne
    Member
    • Jan 2010
    • 68

    #2
    Not quite sure if it's what you're after but Tablet (http://bioinf.scri.ac.uk/tablet) will give you the percentage mismatch for a contig (read bases vs consensus bases, averaged over every read).
    Our software: Tablet | Flapjack | Strudel | CurlyWhirly | TOPALi

    Comment

    • drio
      Senior Member
      • Oct 2008
      • 323

      #3
      Originally posted by johnsequence View Post
      How to calculate the percentage of alignment match after mapping in FASTQ format? Thanks.
      your mapping out is in fastq format? Define % of alignment match.
      You should dump your alignments in [S|B]AM and then write your own
      tool to report stats. samtools has flagstat (too cryptic for my taste).
      dnaa (dnaa.sourceforge.net -- dnaa/dbamstats).

      Ultimately I suggest you use any of the samtools APis and write your
      own stats tool.
      -drd

      Comment

      • KevinLam
        Senior Member
        • Nov 2009
        • 204

        #4
        Originally posted by imilne View Post
        Not quite sure if it's what you're after but Tablet (http://bioinf.scri.ac.uk/tablet) will give you the percentage mismatch for a contig (read bases vs consensus bases, averaged over every read).
        Neat.
        is BAM - A BAM file is a highly compressed, binary version of SAM. supported by tablet? it doesn't say explicitly
        http://kevin-gattaca.blogspot.com/

        Comment

        • imilne
          Member
          • Jan 2010
          • 68

          #5
          Originally posted by KevinLam View Post
          Neat.
          is BAM - A BAM file is a highly compressed, binary version of SAM. supported by tablet? it doesn't say explicitly
          Yep, it's supported (since the last version). I'll update that page to say so too.
          Our software: Tablet | Flapjack | Strudel | CurlyWhirly | TOPALi

          Comment

          • bioinfosm
            Senior Member
            • Jan 2008
            • 483

            #6
            Originally posted by imilne View Post
            Yep, it's supported (since the last version). I'll update that page to say so too.
            I am very excited to use tablet, just cautious of its performance of human genome mapped reads!
            --
            bioinfosm

            Comment

            • imilne
              Member
              • Jan 2010
              • 68

              #7
              Originally posted by bioinfosm View Post
              I am very excited to use tablet, just cautious of its performance of human genome mapped reads!
              Heh, cautious in what way?

              We've had 100GB+ bam files loaded into Tablet ok on a 1GB Eee netbook. For now, bam isn't the problem, it's the reference data, which still isn't cached. We do have very efficient in-memory storage (approx 4x compression) but obviously with 100s of megabases of reference sequence that's still a *lot* of data.

              Iain
              Our software: Tablet | Flapjack | Strudel | CurlyWhirly | TOPALi

              Comment

              • KevinLam
                Senior Member
                • Nov 2009
                • 204

                #8
                Originally posted by imilne View Post
                Heh, cautious in what way?

                We've had 100GB+ bam files loaded into Tablet ok on a 1GB Eee netbook. For now, bam isn't the problem, it's the reference data, which still isn't cached. We do have very efficient in-memory storage (approx 4x compression) but obviously with 100s of megabases of reference sequence that's still a *lot* of data.

                Iain
                Wow, netbook that's pretty impressive!
                Well we just have to wait for someone to write a 'fam' for indexed reference genomes then!
                http://kevin-gattaca.blogspot.com/

                Comment

                • imilne
                  Member
                  • Jan 2010
                  • 68

                  #9
                  Originally posted by KevinLam View Post
                  Wow, netbook that's pretty impressive!
                  Well we just have to wait for someone to write a 'fam' for indexed reference genomes then!
                  I believe there is .fai for indexed fasta. We just haven't had a chance to look at it, or other options.

                  Iain
                  Our software: Tablet | Flapjack | Strudel | CurlyWhirly | TOPALi

                  Comment

                  • bioinfosm
                    Senior Member
                    • Jan 2008
                    • 483

                    #10
                    Originally posted by imilne View Post
                    Heh, cautious in what way?

                    We've had 100GB+ bam files loaded into Tablet ok on a 1GB Eee netbook. For now, bam isn't the problem, it's the reference data, which still isn't cached. We do have very efficient in-memory storage (approx 4x compression) but obviously with 100s of megabases of reference sequence that's still a *lot* of data.

                    Iain
                    some visualization tools tend to freeze the system!
                    I liked the way tablet worked..
                    --
                    bioinfosm

                    Comment

                    Latest Articles

                    Collapse

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by SEQadmin2, 06-09-2026, 11:58 AM
                    0 responses
                    25 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 06-05-2026, 10:09 AM
                    0 responses
                    31 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 06-04-2026, 08:59 AM
                    0 responses
                    39 views
                    0 reactions
                    Last Post SEQadmin2  
                    Started by SEQadmin2, 06-02-2026, 12:03 PM
                    0 responses
                    62 views
                    0 reactions
                    Last Post SEQadmin2  
                    Working...