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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,439,101 T→C K255E (AAG→GAG)  mepA ← penicillin‑insensitive murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,439,1010TC100.0% 93.0 / NA 27K255E (AAG→GAG) mepApenicillin‑insensitive murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/11);  total (16/11)

ATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGT  >  NZ_CP009273/2439014‑2439186
                                                                                       |                                                                                     
aTAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGtcttctcct                                                                                     >  2:637309/1‑90 (MQ=255)
  aaaCAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCa                                                                                   <  1:648665/90‑1 (MQ=255)
         cTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTttg                                                                            >  2:831843/1‑90 (MQ=255)
            ttAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTtttttg                                                                         >  2:475310/1‑90 (MQ=255)
               aaCGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTc                                                                      >  1:215204/1‑90 (MQ=255)
                aCGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTcc                                                                     <  1:861190/90‑1 (MQ=255)
                 cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                                    >  2:861926/1‑90 (MQ=255)
                 cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                                    >  2:411148/1‑90 (MQ=255)
                 cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                                    <  1:269262/90‑1 (MQ=255)
                 cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                                    >  1:537161/1‑90 (MQ=255)
                   tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                                  >  1:149788/1‑90 (MQ=255)
                   tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                                  >  2:479493/1‑90 (MQ=255)
                   tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                                  >  2:425127/1‑90 (MQ=255)
                          tAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGa                                                           <  2:215204/90‑1 (MQ=255)
                             aTCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGt                                                        <  2:25619/90‑1 (MQ=255)
                                    gCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGa                                                           >  1:350688/1‑80 (MQ=255)
                                    gCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGa                                                           <  2:350688/80‑1 (MQ=255)
                                    gCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAAcc                                                 <  1:530900/90‑1 (MQ=255)
                                    gCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGGTCAAAcc                                                 >  2:366382/1‑90 (MQ=255)
                                          ccAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCttt                                           >  1:664529/1‑90 (MQ=255)
                                                 cgcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTc                                    <  1:639315/90‑1 (MQ=255)
                                                  gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                                   >  2:410415/1‑90 (MQ=255)
                                                  gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                                   >  1:680814/1‑90 (MQ=255)
                                                  gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGGTCAAACCAGCTTTGCAGTTCt                                   >  2:545732/1‑90 (MQ=255)
                                                  gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTg                                  >  1:574435/1‑90 (MQ=255)
                                                         aGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGc                            <  2:285562/90‑1 (MQ=255)
                                                              ggcggcAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCa                       <  2:664529/90‑1 (MQ=255)
                                                                              gtgtCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTa       <  2:67468/90‑1 (MQ=255)
                                                                                   ttctcctCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGt  <  1:753796/90‑1 (MQ=255)
                                                                                   ttctcctCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGt  <  2:815399/90‑1 (MQ=255)
                                                                                       |                                                                                     
ATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGT  >  NZ_CP009273/2439014‑2439186

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ 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
AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTC  >  NZ_CP009273/2439016‑2439200
                                                                                     |                                                                                                   
AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTG                                                                                       <  SRR3722072.660251/100‑1 (MQ=60)
   CAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTC                                                                                    >  SRR3722072.218858/1‑100 (MQ=60)
     GGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCC                                                                                  >  SRR3722072.546831/1‑100 (MQ=60)
       CTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGG                                                                                >  SRR3722072.152431/1‑100 (MQ=60)
              ACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGA                                                                         <  SRR3722072.876049/100‑1 (MQ=60)
               CGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAG                                                                        <  SRR3722072.273840/100‑1 (MQ=60)
                        TAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACC                                                               >  SRR3722072.356868/1‑100 (MQ=60)
                              ACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTT                                                         >  SRR3722072.676412/1‑100 (MQ=60)
                                  GCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAG                                                     <  SRR3722072.540466/100‑1 (MQ=60)
                                      ATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTctgtctcttatacacatctccgagcccacgagacgtccgattat                                                 >  SRR3722072.510284/1‑56 (MQ=60)
                                      ATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCT                                                 >  SRR3722072.692995/1‑100 (MQ=60)
                                      ATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTG                                                >  SRR3722072.584705/1‑100 (MQ=60)
                                               CGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAA                                        <  SRR3722072.650725/100‑1 (MQ=60)
                                                                                 TTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCAC      <  SRR3722072.767109/100‑1 (MQ=60)
                                                                                     CCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTC  >  SRR3722072.104085/1‑100 (MQ=60)
                                                                                     |                                                                                                   
AAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACACTC  >  NZ_CP009273/2439016‑2439200

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