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 661,892 A→C G594G (GGT→GGG mrdA ← peptidoglycan DD‑transpeptidase MrdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273661,8920AC100.0% 129.8 / NA 37G594G (GGT→GGGmrdApeptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/18);  total (19/18)

TTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTC  >  NZ_CP009273/661805‑661978
                                                                                       |                                                                                      
ttCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCaccccc                                                                                      <  2:458653/90‑1 (MQ=255)
  ccGCAGGCAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                    <  1:737030/90‑1 (MQ=255)
   cGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGtt                                                                                   <  2:383192/90‑1 (MQ=255)
                            cccAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACt                                                          >  1:475959/1‑90 (MQ=255)
                              cAGCATAATGTGGTCGAGGTTCTGGCGCATCAGTGTCCCACCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTg                                                        <  1:601429/90‑1 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  2:734949/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  2:280940/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  2:168195/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  2:15337/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                       >  1:604021/1‑90 (MQ=255)
                                gCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCg                                                      >  1:659011/1‑90 (MQ=255)
                                 cATAATGGGGTCGAGGATTTGGCGCATCCGTGAACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     <  1:337330/90‑1 (MQ=255)
                                   tAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   <  1:580332/90‑1 (MQ=255)
                                   tAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   <  1:424812/90‑1 (MQ=255)
                                   tAATGTGGTCGAGGATCGGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                   <  1:603484/90‑1 (MQ=255)
                                     aTGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAttg                                                 >  2:464297/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  2:158050/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  1:139578/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  2:671841/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  2:478212/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                            >  1:377329/1‑90 (MQ=255)
                                           tCGAGGATCTGGCGCATCAGTGTACCACCCGCCGGACCCGCCCCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                                           <  2:102460/90‑1 (MQ=255)
                                           tCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                                           <  1:464297/90‑1 (MQ=255)
                                           tCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                                           <  2:642727/90‑1 (MQ=255)
                                               ggATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGgcg                                       >  1:367591/1‑90 (MQ=255)
                                                    tGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAg                                  >  1:688851/1‑90 (MQ=255)
                                                     ggCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAgg                                 <  2:688851/90‑1 (MQ=255)
                                                      gcgcATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGc                                <  2:276998/90‑1 (MQ=255)
                                                       cgcATAAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCg                               <  1:671841/90‑1 (MQ=255)
                                                            cAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcat                          <  2:367591/90‑1 (MQ=255)
                                                             aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc                         <  2:241923/90‑1 (MQ=255)
                                                                     aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     >  2:333071/1‑54 (MQ=255)
                                                                     aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                     <  1:333071/54‑1 (MQ=255)
                                                                     aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTtgtg                 >  1:401067/1‑90 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              <  2:139578/90‑1 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              >  1:398926/1‑90 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              <  1:390030/90‑1 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc              >  2:132138/1‑90 (MQ=255)
                                                                                    acccccGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTc  <  1:430389/90‑1 (MQ=255)
                                                                                       |                                                                                      
TTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTC  >  NZ_CP009273/661805‑661978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ 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
CCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  >  NZ_CP009273/661807‑661988
                                                                                     |                                                                                                
CCGCAGGCAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAAT                                                                                    <  SRR3722073.751951/100‑1 (MQ=60)
CCGCAGGCAGATCGGTGTTGTTATCCCCCAGCAAAATGTGGTCGAGGTTCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAAT                                                                                    <  SRR3722073.531137/100‑1 (MQ=60)
                GTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACT                                                                    >  SRR3722073.486095/1‑100 (MQ=60)
                   GTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGC                                                                 >  SRR3722073.616600/1‑100 (MQ=60)
                    TTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCG                                                                >  SRR3722073.672513/1‑100 (MQ=60)
                            caggctaatgtgctgaaggacctgggccatcagtttcCCAACCGCCGGACCCGCACCCCCGTTCTCCAGCATCATGGCGACAGCCACTTGCGGATTGTTG                                                        <  SRR3722073.494885/63‑1 (MQ=60)
                            CAGCATAATGTGGTCGAGGTTCTGGCGCATCAGTGTCCCACCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTG                                                        <  SRR3722073.613973/100‑1 (MQ=60)
                              GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                      >  SRR3722073.385366/1‑100 (MQ=60)
                              GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                      >  SRR3722073.142481/1‑100 (MQ=60)
                               cagcgtcagatgtgtataagagacagTCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                     <  SRR3722073.340145/74‑1 (MQ=60)
                               CATAATGGGGTCGAGGATTTGGCGCATCCGTGAACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                     <  SRR3722073.344495/100‑1 (MQ=60)
                                 TAATGTGGTCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGG                                                   <  SRR3722073.592445/100‑1 (MQ=60)
                                 TAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGG                                                   <  SRR3722073.433899/100‑1 (MQ=60)
                                 TAATGTGGTCGAGGATCGGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGG                                                   <  SRR3722073.616052/100‑1 (MQ=60)
                                   ATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCG                                                 >  SRR3722073.375394/1‑100 (MQ=60)
                                        GTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAG                                            >  SRR3722073.702885/1‑100 (MQ=60)
                                         TCGAGGATCTGGCGCATCAGTGTACCACCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGG                                           <  SRR3722073.474203/100‑1 (MQ=60)
                                                     CGAATAAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTT                               <  SRR3722073.685580/100‑1 (MQ=60)
                                                         TCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTG                           >  SRR3722073.409606/1‑100 (MQ=60)
                                                            GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC                        >  SRR3722073.407418/1‑100 (MQ=60)
                                                                      CGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGC              <  SRR3722073.398341/100‑1 (MQ=60)
                                                                           GACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGC         >  SRR3722073.542645/1‑100 (MQ=60)
                                                                                  ACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  <  SRR3722073.439601/100‑1 (MQ=60)
                                                                                     |                                                                                                
CCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTG  >  NZ_CP009273/661807‑661988

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