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 3,103,871 (A)7→6 intergenic (+71/+78) BW25113_RS15410 → / ← yghD tRNA‑Phe/GspM family type II secretion system protein YghD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,103,8650A.100.0% 107.9 / NA 26intergenic (+65/+84)BW25113_RS15410/yghDtRNA‑Phe/GspM family type II secretion system protein YghD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (9/17);  total (9/17)

GATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103781‑3103949
                                                                                    |                                                                                    
gATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaa                                                                                 >  2:683389/1‑90 (MQ=255)
           cGGGCACCACTAATTCTTAAGAACCCTCCCACAAGGCGGGTTTTTGCTTTTGGTTCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                                     <  1:246692/90‑1 (MQ=255)
            gggCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaaa                                                                    >  2:407251/1‑90 (MQ=255)
             ggCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAg                                                                   <  2:365977/90‑1 (MQ=255)
                 ccacTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTg                                                               <  1:407251/90‑1 (MQ=255)
                             aaGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTActc                                                   <  1:250971/90‑1 (MQ=255)
                                 aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  1:428489/90‑1 (MQ=255)
                                 aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  1:772635/90‑1 (MQ=255)
                                  cccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGcc                                              >  2:65371/1‑90 (MQ=255)
                                   ccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  2:38909/90‑1 (MQ=255)
                                           aaGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                                     <  1:683389/90‑1 (MQ=255)
                                             ggCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACgag                                   <  2:193351/90‑1 (MQ=255)
                                                gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  1:99383/90‑1 (MQ=255)
                                                gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  1:95656/90‑1 (MQ=255)
                                                gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCACGCCACGTTAACGAGAGt                                >  1:358847/1‑90 (MQ=255)
                                                            ttGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGat                    <  2:492940/90‑1 (MQ=255)
                                                                aTCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTa                <  1:489655/90‑1 (MQ=255)
                                                                  cTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  1:159924/90‑1 (MQ=255)
                                                                  cTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  2:232011/90‑1 (MQ=255)
                                                                  cTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  2:79434/90‑1 (MQ=255)
                                                                        aTAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTg        >  2:95921/1‑90 (MQ=255)
                                                                          aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:81859/90‑1 (MQ=255)
                                                                              ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:336064/1‑90 (MQ=255)
                                                                              ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:346343/1‑90 (MQ=255)
                                                                              ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  1:758679/1‑90 (MQ=255)
                                                                              ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  1:509084/1‑90 (MQ=255)
                                                                              ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:600806/1‑90 (MQ=255)
                                                                                    |                                                                                    
GATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103781‑3103949

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ 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
CGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCCGTCCAAACTC  >  NZ_CP009273/3103792‑3103966
                                                                               |                                                                                               
CGGGCACCACTAATTCTTAAGAACCCTCCCACAAGGCGGGTTTTTGCTTTTGGTTCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGG                                                                            <  SRR3722073.251869/100‑1 (MQ=60)
      CCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAC                                                                      <  SRR3722073.415947/100‑1 (MQ=60)
                  AAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTT                                                          <  SRR3722073.256249/100‑1 (MQ=60)
                      ACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACG                                                      <  SRR3722073.437662/100‑1 (MQ=60)
                      ACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACG                                                      <  SRR3722073.788217/100‑1 (MQ=60)
                           CCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCACGCCACGTTAACGAGAGT                                                 >  SRR3722073.366442/1‑100 (MQ=60)
                                AAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTT                                            <  SRR3722073.697326/100‑1 (MQ=60)
                                     GGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTG                                       <  SRR3722073.101544/100‑1 (MQ=60)
                                     GGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTG                                       <  SRR3722073.97767/100‑1 (MQ=60)
                                                     ATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAT                       <  SRR3722073.500040/100‑1 (MQ=60)
                                                       CTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATAC                     <  SRR3722073.163149/100‑1 (MQ=60)
                                                         GACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT                   >  SRR3722073.519818/1‑100 (MQ=60)
                                                         GACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT                   >  SRR3722073.774011/1‑100 (MQ=60)
                                                               AACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCC             <  SRR3722073.83710/100‑1 (MQ=60)
                                                                           AAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCCGTCCAAACTC  <  SRR3722073.66902/100‑1 (MQ=60)
                                                                               |                                                                                               
CGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCCGTCCAAACTC  >  NZ_CP009273/3103792‑3103966

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