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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

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% 90.2 / NA 28G594G (GGT→GGGmrdApeptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (13/15);  total (13/15)

GGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC  >  NZ_CP009273/661812‑661966
                                                                                |                                                                          
ggcagATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTcc                                                                   >  1:178403/1‑90 (MQ=255)
  cagATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAg                                                                 <  1:55712/90‑1 (MQ=255)
                tATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAg                                                   >  2:304028/1‑90 (MQ=255)
                   cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCa                                                <  2:193130/90‑1 (MQ=255)
                       cAGCATAATGTGGTCGAGGATCTGGCTCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTg                                            >  1:42426/1‑90 (MQ=255)
                        atcATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                           >  2:480792/3‑90 (MQ=255)
                            tAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                       <  1:15702/90‑1 (MQ=255)
                               tgtgGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAttgt                                    <  2:457861/90‑1 (MQ=255)
                                   gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                >  2:279124/1‑90 (MQ=255)
                                   gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                >  1:504318/1‑90 (MQ=255)
                                    tCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                               <  1:336657/90‑1 (MQ=255)
                                      gAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATgg                             >  2:442756/1‑90 (MQ=255)
                                             tGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAg                      >  2:36631/1‑90 (MQ=255)
                                                   aTCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                               >  1:167044/1‑75 (MQ=255)
                                                   aTCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                               <  2:167044/75‑1 (MQ=255)
                                                   aTCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtc                <  2:363840/90‑1 (MQ=255)
                                                      aGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc             <  1:525333/90‑1 (MQ=255)
                                                      aGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc             <  1:408848/90‑1 (MQ=255)
                                                      aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc             <  2:424638/90‑1 (MQ=255)
                                                      aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc             <  1:228477/90‑1 (MQ=255)
                                                      aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc             <  1:194561/90‑1 (MQ=255)
                                                                 cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGc                    <  2:246713/72‑1 (MQ=255)
                                                                 cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGc                    >  1:246713/1‑72 (MQ=255)
                                                                 cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc  <  1:375479/90‑1 (MQ=255)
                                                                 cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc  >  2:333765/1‑90 (MQ=255)
                                                                 cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc  >  2:50758/1‑90 (MQ=255)
                                                                   ccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAgg                     >  2:280731/1‑69 (MQ=255)
                                                                   ccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAgg                     <  1:280731/69‑1 (MQ=255)
                                                                                |                                                                          
GGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC  >  NZ_CP009273/661812‑661966

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ 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
ATTTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGT  >  NZ_CP009273/661802‑661989
                                                                                          |                                                                                                 
ATTTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCC                                                                                          >  SRR3722111.180643/1‑100 (MQ=60)
     gagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacaGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAAT                                                                                     <  SRR3722111.9726/49‑1 (MQ=60)
            CAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCG                                                                              <  SRR3722111.56377/100‑1 (MQ=60)
                       TGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCTCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTG                                                                   >  SRR3722111.42946/1‑100 (MQ=60)
                                   GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                       >  SRR3722111.512104/1‑100 (MQ=60)
                                      TAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGG                                                    <  SRR3722111.15894/100‑1 (MQ=60)
                                              TCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGG                                            <  SRR3722111.341163/100‑1 (MQ=60)
                                                   GATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCTGTC                                       >  SRR3722111.169140/1‑100 (MQ=60)
                                                        ataagagacagGTACCAACCGCCGGACCCGCCCCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAG                                  <  SRR3722111.284261/89‑1 (MQ=60)
                                                                AGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGT                          <  SRR3722111.414837/100‑1 (MQ=60)
                                                                AGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGT                          <  SRR3722111.533560/100‑1 (MQ=60)
                                                                AGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGT                          <  SRR3722111.196982/100‑1 (MQ=60)
                                                                AGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGT                          <  SRR3722111.231226/100‑1 (MQ=60)
                                                                 GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTctgtctct                         >  SRR3722111.249625/1‑92 (MQ=60)
                                                                           CGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGC               <  SRR3722111.380832/100‑1 (MQ=60)
                                                                                        CCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGT  >  SRR3722111.387043/1‑100 (MQ=60)
                                                                                          |                                                                                                 
ATTTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGT  >  NZ_CP009273/661802‑661989

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