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 I230 R1 226 18.8 1048726 96.7% 1014118 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,098,336 Δ1 bp coding (545/1080 nt) mltC → membrane‑bound lytic murein transglycosylase MltC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,098,3350C.100.0% 72.2 / NA 18coding (544/1080 nt)mltCmembrane‑bound lytic murein transglycosylase MltC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/9);  total (9/9)

CTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGC  >  NZ_CP009273/3098250‑3098416
                                                                                     |                                                                                 
cTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGcgaac                                                                               <  1:128318/90‑5 (MQ=255)
  tCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAacca                                                                            <  1:113250/90‑1 (MQ=255)
     cgGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCt                                                                         >  1:487710/1‑90 (MQ=255)
                 tgAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTgc                                                             >  2:477267/1‑90 (MQ=255)
                 tgAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTgc                                                             <  1:224779/90‑1 (MQ=255)
                   aaaaaCCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTgcgc                                                           <  2:65572/90‑1 (MQ=255)
                         cGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCCCCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAat                                                     >  1:224413/1‑90 (MQ=255)
                         cGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAat                                                     <  1:477267/90‑1 (MQ=255)
                                      gcagcaACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGt                                        <  1:98824/90‑1 (MQ=255)
                                           aaCGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGcc                                   >  1:83130/1‑90 (MQ=255)
                                                 cTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGt                             >  1:128119/1‑90 (MQ=255)
                                                  tGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTc                            <  1:487140/90‑1 (MQ=255)
                                                    cGTATCATCTACAGCGTCCCCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCAc                          <  2:224413/90‑1 (MQ=255)
                                                         catcTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGaaa                     >  1:24755/1‑90 (MQ=255)
                                                            cTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAatat                  <  2:498067/90‑1 (MQ=255)
                                                             tACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATg                 >  2:114540/1‑90 (MQ=255)
                                                                 gCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGt             >  2:32297/1‑90 (MQ=255)
                                                                   gTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTtg           >  1:111227/1‑90 (MQ=255)
                                                                            aaCATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGc  <  1:184824/90‑1 (MQ=255)
                                                                                     |                                                                                 
CTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGC  >  NZ_CP009273/3098250‑3098416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTAT  >  NZ_CP009273/3098245‑3098432
                                                                                           |                                                                                                
AGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCT                                                                                         >  SRR3722114.493815/1‑100 (MQ=60)
     CTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATA                                                                                    <  SRR3722114.129653/100‑1 (MQ=60)
       TCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAA                                                                                  <  SRR3722114.114448/100‑1 (MQ=60)
                    GCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCCCCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAAT                                                                     >  SRR3722114.226544/1‑100 (MQ=60)
                      TGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATAT                                                                   <  SRR3722114.226913/100‑1 (MQ=60)
                              agatctagtacggtctcgtgggctcggagatgtgtataagagacagCCATTAACATGGTGCGAACCTCCTTGATAAACGTGCGCACAAATATCTCGGCAT                                                           <  SRR3722114.375360/54‑1 (MQ=60)
                              CGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCAT                                                           <  SRR3722114.483246/100‑1 (MQ=60)
                               gacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacaggtgCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCAT                                                           <  SRR3722114.284885/40‑1 (MQ=60)
                                      GAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCC                                                   >  SRR3722114.84006/1‑100 (MQ=60)
                                           GCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCG                                              <  SRR3722114.99865/100‑1 (MQ=60)
                                            CAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGT                                             >  SRR3722114.129453/1‑100 (MQ=60)
                                                    GGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAA                                     >  SRR3722114.24997/1‑100 (MQ=60)
                                                       TGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATAT                                  <  SRR3722114.493239/100‑1 (MQ=60)
                                                              CATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTG                           >  SRR3722114.112402/1‑100 (MQ=60)
                                                                            CCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATT             >  SRR3722114.333/1‑100 (MQ=60)
                                                                            CCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATT             >  SRR3722114.111085/1‑100 (MQ=60)
                                                                               TTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTG          >  SRR3722114.415896/1‑100 (MQ=60)
                                                                                 AACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGC        <  SRR3722114.186644/100‑1 (MQ=60)
                                                                                       GTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTAT  >  SRR3722114.166204/1‑100 (MQ=60)
                                                                                           |                                                                                                
AGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTAT  >  NZ_CP009273/3098245‑3098432

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