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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,818,451 T→C E162G (GAA→GGA)  mltB ← lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,818,4510TC100.0% 74.3 / NA 24E162G (GAA→GGA) mltBlytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/16);  total (8/16)

GCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAT  >  NZ_CP009273/2818370‑2818520
                                                                                 |                                                                     
gcgcgGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTTCCCATCACGCGCCCCCAGCGGGTTCCAACGctg                                                               >  2:32052/1‑88 (MQ=255)
gcgcgGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCg                                                               >  2:59584/1‑90 (MQ=255)
gcgcgGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCg                                                               >  2:451623/1‑90 (MQ=255)
 cgcgGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGa                                                              <  1:203265/90‑1 (MQ=255)
                 aaaTGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAAt                                              >  2:52814/1‑90 (MQ=255)
                   aTGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAAttt                                            >  2:621212/1‑90 (MQ=255)
                     gACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTcc                                          <  1:14870/90‑1 (MQ=255)
                      aCAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcg                                         >  1:408577/1‑90 (MQ=255)
                        aGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggc                                       <  1:521312/90‑1 (MQ=255)
                        aGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggc                                       <  2:555745/90‑1 (MQ=255)
                        aGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggc                                       <  1:252912/90‑1 (MQ=255)
                          cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg                                     <  2:295150/90‑1 (MQ=255)
                             tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                  <  1:516271/90‑1 (MQ=255)
                             tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                  <  2:392694/90‑1 (MQ=255)
                             tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                  <  2:474297/90‑1 (MQ=255)
                               ccAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTc                                <  1:598226/90‑1 (MQ=255)
                                  gcgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg                                     <  1:393220/82‑1 (MQ=255)
                                  gcgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg                                     >  2:393220/1‑82 (MQ=255)
                                  gcgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAt                             <  2:43704/90‑1 (MQ=255)
                                   cgcATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATa                            >  2:648285/1‑90 (MQ=255)
                                      aTCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATacac                         <  2:408577/90‑1 (MQ=255)
                                                         cccATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAAc      <  1:451623/90‑1 (MQ=255)
                                                           cATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGc    <  1:625747/90‑1 (MQ=255)
                                                             tCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAt  <  1:520603/90‑1 (MQ=255)
                                                                                 |                                                                     
GCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCAT  >  NZ_CP009273/2818370‑2818520

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
CAGAGAAATACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAAC  >  NZ_CP009273/2818356‑2818542
                                                                                               |                                                                                           
CAGAGAAATACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAAC                                                                                         <  SRR3722088.544726/100‑1 (MQ=60)
               CGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGAC                                                                          <  SRR3722088.205302/100‑1 (MQ=60)
                          TAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCG                                                               >  SRR3722088.413075/1‑100 (MQ=60)
                                   GACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTC                                                      <  SRR3722088.15038/100‑1 (MQ=60)
                                      AGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAT                                                   <  SRR3722088.255376/100‑1 (MQ=60)
                                      AGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAT                                                   <  SRR3722088.527606/100‑1 (MQ=60)
                                        CGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATAC                                                 <  SRR3722088.397440/100‑1 (MQ=60)
                                           TGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACC                                              <  SRR3722088.522500/100‑1 (MQ=60)
                                             CCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTG                                            <  SRR3722088.605591/100‑1 (MQ=60)
                                                                       CCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCAT                  <  SRR3722088.456887/100‑1 (MQ=60)
                                                                         CATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATAC                <  SRR3722088.633476/100‑1 (MQ=60)
                                                                           TCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTG              <  SRR3722088.526887/100‑1 (MQ=60)
                                                                                  CCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAG       >  SRR3722088.76071/1‑100 (MQ=60)
                                                                                       AGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAAC  >  SRR3722088.626357/1‑100 (MQ=60)
                                                                                               |                                                                                           
CAGAGAAATACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAGAAAAC  >  NZ_CP009273/2818356‑2818542

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