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,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% 113.4 / NA 27intergenic (+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/18);  total (9/18)

ATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103782‑3103949
                                                                                   |                                                                                    
aTTCCGCGTCCGGGCCCCACTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                                <  2:70547/90‑2 (MQ=255)
aTTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                                >  2:118741/1‑89 (MQ=255)
     gAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCt                                                                          >  1:443120/1‑90 (MQ=255)
          cGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                                     <  2:38421/90‑1 (MQ=255)
                 cacTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTgg                                                              >  1:708850/1‑90 (MQ=255)
                  acTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTggg                                                             <  1:432469/90‑1 (MQ=255)
                  acTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTggg                                                             <  2:422011/90‑1 (MQ=255)
                     aaTTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGa                                                          >  2:113893/1‑90 (MQ=255)
                          ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc                                                     >  2:4707/1‑90 (MQ=255)
                                aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  1:605000/90‑1 (MQ=255)
                                aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  2:401209/90‑1 (MQ=255)
                                aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  2:575963/90‑1 (MQ=255)
                                  ccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  1:268943/90‑1 (MQ=255)
                                  ccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  2:125606/90‑1 (MQ=255)
                                  ccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  2:33989/90‑1 (MQ=255)
                                   cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc                                            >  2:87888/1‑90 (MQ=255)
                                   cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc                                            >  2:56113/1‑90 (MQ=255)
                                          aaGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                                     <  1:130159/90‑1 (MQ=255)
                                                       gCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTAt                        <  1:4707/90‑1 (MQ=255)
                                                        cTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTAtt                       <  2:692631/90‑1 (MQ=255)
                                                             ggATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGatat                  <  2:708850/90‑1 (MQ=255)
                                                                  tGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtc             <  1:56113/90‑1 (MQ=255)
                                                                    aCAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtctc           >  2:342730/1‑90 (MQ=255)
                                                                         aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:87888/90‑1 (MQ=255)
                                                                         aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:41155/90‑1 (MQ=255)
                                                                          aCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATa     <  1:629190/90‑1 (MQ=255)
                                                                             ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             >  1:716016/1‑47 (MQ=255)
                                                                             ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  2:716016/47‑1 (MQ=255)
                                                                             ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  1:199638/1‑90 (MQ=255)
                                                                                   |                                                                                    
ATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103782‑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
GTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCC  >  NZ_CP009273/3103777‑3103956
                                                                                              |                                                                                     
GTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCT                                                                                 >  SRR3722112.449264/1‑100 (MQ=60)
            GTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGG                                                                     >  SRR3722112.719928/1‑100 (MQ=60)
                       ACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTC                                                          <  SRR3722112.438361/100‑1 (MQ=60)
                                     ACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACG                                            <  SRR3722112.614274/100‑1 (MQ=60)
                                       CCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAG                                          <  SRR3722112.272096/100‑1 (MQ=60)
                                               AAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTT                                  <  SRR3722112.131775/100‑1 (MQ=60)
                                                            GCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAT                     <  SRR3722112.4758/100‑1 (MQ=60)
                                                                       TGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACT          <  SRR3722112.56800/100‑1 (MQ=60)
                                                                        GACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGctgtctcctatacacatctgacgctgccgacga         >  SRR3722112.727198/1‑67 (MQ=60)
                                                                        GACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT         >  SRR3722112.202068/1‑100 (MQ=60)
                                                                              AACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCC   <  SRR3722112.41649/100‑1 (MQ=60)
                                                                              AACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCC   <  SRR3722112.88988/100‑1 (MQ=60)
                                                                               ACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCC  <  SRR3722112.638919/100‑1 (MQ=60)
                                                                                              |                                                                                     
GTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCC  >  NZ_CP009273/3103777‑3103956

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