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 I211 R1 223 18.0 1007644 96.3% 970361 85.7

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% 68.2 / NA 22E162G (GAA→GGA) mltBlytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/12);  total (10/12)

ACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATT  >  NZ_CP009273/2818365‑2818533
                                                                                      |                                                                                  
aCTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCaa                                                                                 <  2:197495/90‑1 (MQ=255)
aCTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCaa                                                                                 >  2:169766/1‑90 (MQ=255)
     gcgcgGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCg                                                                            >  1:277900/1‑90 (MQ=255)
                          gACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTcc                                                       <  1:210762/90‑1 (MQ=255)
                             aGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggc                                                    >  1:257665/1‑90 (MQ=255)
                             aGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggc                                                    >  1:484763/1‑90 (MQ=255)
                               cGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCcggcgg                                                  <  2:319536/90‑1 (MQ=255)
                                 ttGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTa                                                <  1:43125/90‑1 (MQ=255)
                                  tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               <  1:414908/90‑1 (MQ=255)
                                  tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               <  2:339301/90‑1 (MQ=255)
                                  tGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               <  2:348388/90‑1 (MQ=255)
                                       gcgcATCGAGGATGCGAGTTTTCCCCATCCCGTGCCCCCAGCGGGTTCCACCGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCAt                                          <  1:265933/90‑1 (MQ=255)
                                                    gCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               <  1:318095/72‑1 (MQ=255)
                                                    gCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               >  2:318095/1‑72 (MQ=255)
                                                      gAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCAcgcgc                           >  1:18836/1‑90 (MQ=255)
                                                               ccATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACg                  >  2:426020/1‑90 (MQ=255)
                                                                cATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGc                 <  2:484763/90‑1 (MQ=255)
                                                                     cgcgCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCtt            <  1:426020/90‑1 (MQ=255)
                                                                             cAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGa    >  1:319379/1‑90 (MQ=255)
                                                                               gCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGAtt  >  1:352732/1‑90 (MQ=255)
                                                                                   gTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               >  1:339366/1‑41 (MQ=255)
                                                                                   gTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTAc                                               <  2:339366/41‑1 (MQ=255)
                                                                                      |                                                                                  
ACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATT  >  NZ_CP009273/2818365‑2818533

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAAATACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAG  >  NZ_CP009273/2818360‑2818537
                                                                                           |                                                                                      
GAAATACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCG                                                                                >  SRR3722094.282071/1‑100 (MQ=60)
                        GTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGC                                                        >  SRR3722094.261474/1‑100 (MQ=60)
                        GTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGC                                                        >  SRR3722094.493316/1‑100 (MQ=60)
                               GACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTC                                                 <  SRR3722094.213605/100‑1 (MQ=60)
                                      TTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACAC                                          <  SRR3722094.43731/100‑1 (MQ=60)
                                       agagacaGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACC                                         <  SRR3722094.323165/93‑1 (MQ=60)
                                       TGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACC                                         <  SRR3722094.422130/100‑1 (MQ=60)
                                            GCGCATCGAGGATGCGAGTTTTCCCCATCCCGTGCCCCCAGCGGGTTCCACCGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCA                                    <  SRR3722094.269897/100‑1 (MQ=60)
                                                 TCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGC                               >  SRR3722094.19115/1‑100 (MQ=60)
                                                                        CACGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGA        >  SRR3722094.324483/1‑100 (MQ=60)
                                                                          CGCGCCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATT      >  SRR3722094.358623/1‑100 (MQ=60)
                                                                          CGCGCCCCCAGCGGGTTCCAACGCCGATAACCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATT      <  SRR3722094.433459/100‑1 (MQ=60)
                                                                              CCCCCAGCGGGTTCCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCctggctcttatacacatctgacgctgccgacgaccatct  >  SRR3722094.344953/1‑61 (MQ=60)
                                                                                           |                                                                                      
GAAATACTCCGCGCGGCGTGGGTAGTTAAATGACAGCGTTGCCAGCGCATCGAGGATGCGAGTTTTCCCCATCACGCGCCCCCAGCGGGTTTCAACGCCGATAATCCCGACGATAATTTCCGGCGGTACTCCATACACCTGCCACGCGCGATTCAACGCATCTTCATACTGATTCCAG  >  NZ_CP009273/2818360‑2818537

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