Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 802,549 T→G T63T (ACA→ACC ybhC ← putative acyl‑CoA thioester hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273802,5490TG100.0% 67.0 / NA 20T63T (ACA→ACCybhCputative acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/11);  total (9/11)

GAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT  >  NZ_CP009273/802470‑802635
                                                                               |                                                                                      
gAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATg                                                                              <  1:17534/90‑1 (MQ=255)
 aGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGc                                                                             <  1:130054/90‑1 (MQ=255)
 aGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGc                                                                             <  2:301812/90‑1 (MQ=255)
   tGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCaa                                                                           <  1:122434/90‑1 (MQ=255)
    gCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCaaa                                                                          >  2:694423/1‑90 (MQ=255)
    gCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCaaa                                                                          >  1:453057/1‑90 (MQ=255)
      ccgccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                                        >  1:587336/1‑90 (MQ=255)
      ccgccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                                        <  1:617201/90‑1 (MQ=255)
      ccgccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                                        <  2:17994/90‑1 (MQ=255)
      ccgccgGGCCAACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                                        >  1:772787/1‑90 (MQ=255)
        gccgGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAg                                                                      <  1:393135/90‑1 (MQ=255)
                gACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGc                                                              <  1:138694/90‑1 (MQ=255)
                                       gcaggGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTg                                       >  1:360410/1‑90 (MQ=255)
                                                tCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGaa                              >  2:239394/1‑90 (MQ=255)
                                                  ggggAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATc                            <  1:133541/90‑1 (MQ=255)
                                                        ggATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                                        >  1:389393/1‑40 (MQ=38)
                                                        ggATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAt                                                                        <  2:389393/40‑1 (MQ=38)
                                                             ccACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGGGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAAGCGGgagagaggg                 >  1:196603/1‑87 (MQ=255)
                                                                        caCCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGgcg      <  1:678547/90‑1 (MQ=255)
                                                                            tGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGt  >  1:103253/1‑90 (MQ=255)
                                                                               |                                                                                      
GAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT  >  NZ_CP009273/802470‑802635

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAG  >  NZ_CP009273/802452‑802646
                                                                                                 |                                                                                                 
TGGTATGCGTTCCACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGT                                                                                                 <  SRR3722073.503996/100‑1 (MQ=60)
     TGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGG                                                                                            <  SRR3722073.532007/100‑1 (MQ=60)
            CACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAA                                                                                     >  SRR3722073.462740/1‑100 (MQ=60)
              CCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAT                                                                                   >  SRR3722073.599609/1‑100 (MQ=60)
              CCTTGAGTGCCCGCCGGGCCAACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAAT                                                                                   >  SRR3722073.788370/1‑100 (MQ=60)
                  GAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTG                                                                               <  SRR3722073.17980/100‑1 (MQ=60)
                   AGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGT                                                                              <  SRR3722073.132779/100‑1 (MQ=60)
                     TGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTG                                                                            <  SRR3722073.125027/100‑1 (MQ=60)
                        CCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGC                                                                         <  SRR3722073.630030/100‑1 (MQ=60)
                          GCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCAT                                                                       <  SRR3722073.401511/100‑1 (MQ=60)
                                  GACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTC                                                               <  SRR3722073.141577/100‑1 (MQ=60)
                                               GGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTG                                                  >  SRR3722073.368042/1‑100 (MQ=60)
                                                                AATCGGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCctgtctcttatacacatctgacgctgccgacgatatgcag                                 >  SRR3722073.397683/1‑60 (MQ=60)
                                                                    GGGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAG                             <  SRR3722073.136328/100‑1 (MQ=60)
                                                                     GGGAAGGATTCCACGCTGCAGCACCTGGGGTCAGGGATGCAAAATAGGGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAAgcgggagagaggg                            >  SRR3722073.200539/1‑87 (MQ=60)
                                                                                    CTGCAGCACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGT             >  SRR3722073.105485/1‑100 (MQ=60)
                                                                                          CACCTGGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATC       <  SRR3722073.692412/100‑1 (MQ=60)
                                                                                               GGGGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAG  <  SRR3722073.632819/100‑1 (MQ=60)
                                                                                                 |                                                                                                 
TGGTATGCGTTACACCTTGAGTGCCCGCCGGGCCGACAACAAAGTCAGGTTGCGCAGGCAGGGTAATCGGGGAAGGATTCCACGCTGCAGCACCTGGTGTCAGGGATGCAAAATAGTGTTGAGCATCGAAATTCTGCGCTTCTTTTGCCGACAGAATCGGGCGAGAAGAGGTACCAGGCGCGGTTTGATCAGAAG  >  NZ_CP009273/802452‑802646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: