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 I228 R1 214 26.5 1448192 97.1% 1406194 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% 86.2 / NA 21coding (544/1080 nt)mltCmembrane‑bound lytic murein transglycosylase MltC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (8/13);  total (8/13)

TTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCG  >  NZ_CP009273/3098251‑3098415
                                                                                    |                                                                                
ttCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAacc                                                                            <  2:138591/90‑1 (MQ=255)
     ggATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCtt                                                                       >  2:663953/1‑90 (MQ=255)
                cgAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTgc                                                            <  1:442077/89‑1 (MQ=255)
                        cGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAat                                                    >  2:196764/1‑90 (MQ=255)
                                   ccGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATg                                         <  1:395/90‑1 (MQ=255)
                                    cGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATgg                                        <  2:490716/90‑1 (MQ=255)
                                              ggCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGc                              <  2:106211/90‑1 (MQ=255)
                                                 tGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTc                           <  1:136887/90‑1 (MQ=255)
                                                  gCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCa                          >  2:364950/1‑90 (MQ=255)
                                                  gCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCa                          >  1:484192/1‑90 (MQ=255)
                                                   cGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCAc                         <  1:692972/90‑1 (MQ=255)
                                                   cGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCAc                         <  1:500980/90‑1 (MQ=255)
                                                           cTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAatat                 <  2:189584/90‑1 (MQ=255)
                                                           cTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAatat                 <  1:470095/90‑1 (MQ=255)
                                                            tACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATg                >  2:35408/1‑90 (MQ=255)
                                                                  gTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCATAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGGGTtg          >  1:551002/1‑90 (MQ=255)
                                                                   tCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTtga         >  2:723551/1‑90 (MQ=255)
                                                                     aCCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTtgatg       <  1:723551/90‑1 (MQ=255)
                                                                     aCCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTtgatg       <  1:444/90‑1 (MQ=255)
                                                                     aCCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTtgatg       <  1:378436/90‑1 (MQ=255)
                                                                          tAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCg  >  2:29468/1‑90 (MQ=255)
                                                                                    |                                                                                
TTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCG  >  NZ_CP009273/3098251‑3098415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
CGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGC  >  NZ_CP009273/3098241‑3098434
                                                                                               |                                                                                                  
CGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACC                                                                                               <  SRR3722112.522002/100‑1 (MQ=60)
CGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACC                                                                                               <  SRR3722112.485557/100‑1 (MQ=60)
             atctacacattgacattcgtcggcagcgtcagatgtgtataagagacagGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACG                                                                                  <  SRR3722112.139968/51‑1 (MQ=60)
                          CGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATAT                                                                     <  SRR3722112.448200/100‑1 (MQ=60)
                                             attcgccggcagcgtcagatgtgtataagagacagCCATTAACATGGTGCGAACCGCCTTGATAAACGTGCGCCCAAATATCTCGGCATGGTCCGCCAGG                                                  <  SRR3722112.392509/65‑1 (MQ=60)
                                             CCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGG                                                  <  SRR3722112.401/100‑1 (MQ=60)
                                                  GCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCA                                             >  SRR3722112.491152/1‑100 (MQ=60)
                                                           TGCGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATAT                                    <  SRR3722112.138590/100‑1 (MQ=60)
                                                             CGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGG                                  <  SRR3722112.703795/100‑1 (MQ=60)
                                                             CGTATCATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGG                                  <  SRR3722112.508239/100‑1 (MQ=60)
                                                                  CATCTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCATAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGGGTTG                             >  SRR3722112.559317/1‑100 (MQ=60)
                                                                     CTACAGCGTCACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATG                          <  SRR3722112.476796/100‑1 (MQ=60)
                                                                               ACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGAT                <  SRR3722112.450/100‑1 (MQ=60)
                                                                               ACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGAT                <  SRR3722112.383350/100‑1 (MQ=60)
                                                                               ACCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGAT                <  SRR3722112.734843/100‑1 (MQ=60)
                                                                                CCATTAACATGGTGCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATT               >  SRR3722112.18718/1‑100 (MQ=60)
                                                                                              GCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGC  <  SRR3722112.468981/100‑1 (MQ=60)
                                                                                              GCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGC  <  SRR3722112.628541/100‑1 (MQ=60)
                                                                                               |                                                                                                  
CGCAAGCAACTTCGCGGATTATCTGCTGAAAAACCGTCTGAAGAGCCGCAGCAACGGGCTGCGTATCATCTACAGCGTCACCATTAACATGGTGCCGAACCACCTTGATAAACGTGCGCACAAATATCTCGGCATGGTCCGCCAGGCGTCACGGAAATATGGCGTTGATGAGTCGCTGATTCTGGCAATTATGC  >  NZ_CP009273/3098241‑3098434

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