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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

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% 177.4 / NA 41intergenic (+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 . (11/30);  total (11/30)

GGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103776‑3103949
                                                                                         |                                                                                    
ggTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGa                                                                                      <  2:1150777/90‑1 (MQ=255)
     gATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaa                                                                                 >  1:372270/1‑90 (MQ=255)
      aTTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                                <  2:600769/90‑2 (MQ=255)
          cGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGc                                                                           <  1:700348/90‑1 (MQ=255)
          cGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGc                                                                           <  2:452943/90‑1 (MQ=255)
               ccGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTtata                                                                      <  2:1120030/90‑1 (MQ=255)
                cGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                                     <  2:912379/90‑1 (MQ=255)
                cGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                                     <  1:881652/90‑1 (MQ=255)
                 gggCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTTACAATAACCTTCACGAAAAAATTAGCTTATaaa                                                                    >  1:1017093/1‑90 (MQ=255)
                      ccacTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTg                                                               <  1:657570/90‑1 (MQ=255)
                        acTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTggg                                                             <  2:90042/90‑1 (MQ=255)
                        acTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTggg                                                             <  2:133722/90‑1 (MQ=255)
                             ttCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGTCAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAAt                                                        <  2:593893/90‑1 (MQ=255)
                             ttCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAAt                                                        <  2:763151/90‑1 (MQ=255)
                              tCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAAtt                                                       <  1:164706/90‑1 (MQ=255)
                                ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc                                                     >  2:237863/1‑90 (MQ=255)
                                ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc                                                     >  1:891083/1‑90 (MQ=255)
                                      aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  2:186541/90‑1 (MQ=255)
                                      aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               >  2:963949/1‑90 (MQ=255)
                                       cccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               >  1:1109138/1‑89 (MQ=255)
                                       cccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  2:1109138/89‑1 (MQ=255)
                                        ccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa                                             <  2:145620/90‑1 (MQ=255)
                                         cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                                <  2:337014/55‑2 (MQ=255)
                                         cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                                >  1:337014/1‑54 (MQ=255)
                                         cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc                                            >  1:1001848/1‑90 (MQ=255)
                                         cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc                                            >  2:209947/1‑90 (MQ=255)
                                              acaAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTaa                                       >  2:553740/1‑90 (MQ=255)
                                                aaGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                                     <  1:89004/90‑1 (MQ=255)
                                                aaGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                                     <  2:480348/90‑1 (MQ=255)
                                                   gCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACgaga                                  <  1:209947/90‑1 (MQ=255)
                                                     gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  2:399091/90‑1 (MQ=255)
                                                                 ttGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGat                    <  1:1136734/90‑1 (MQ=255)
                                                                       cTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  2:22097/90‑1 (MQ=255)
                                                                       cTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt              <  2:372270/90‑1 (MQ=255)
                                                                        tGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTCCTCTCGCCCCGTTAACGAGAGTAATTTTATTGATATTAAtc             <  1:733770/90‑1 (MQ=255)
                                                                        tGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtc             <  1:362120/90‑1 (MQ=255)
                                                                        tGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtc             <  1:627006/90‑1 (MQ=255)
                                                                        tGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtc             <  1:910320/90‑1 (MQ=255)
                                                                         gACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtct            <  1:98566/90‑1 (MQ=255)
                                                                             aTAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTg        <  1:905224/90‑1 (MQ=255)
                                                                              tAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGa       >  1:937439/1‑90 (MQ=255)
                                                                               aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:1191944/90‑1 (MQ=255)
                                                                               aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:976487/90‑1 (MQ=255)
                                                                               aaCCTTCACGAAAAAATTAACTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  1:721256/90‑1 (MQ=255)
                                                                                   ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:329215/1‑90 (MQ=255)
                                                                                   ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:101084/1‑90 (MQ=255)
                                                                                         |                                                                                    
GGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103776‑3103949

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
TTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCCGTCCAAACT  >  NZ_CP009273/3103783‑3103965
                                                                                        |                                                                                              
TTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTTACAATAACCTTCACGAAAAAATTAGCTTATAAA                                                                                    >  SRR3722089.1032250/1‑100 (MQ=60)
   CGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTC                                                                                 <  SRR3722089.709477/100‑1 (MQ=60)
         CGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGG                                                                           <  SRR3722089.894350/100‑1 (MQ=60)
               CCACTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAC                                                                     <  SRR3722089.666041/100‑1 (MQ=60)
               CCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAC                                                                     >  SRR3722089.903959/1‑100 (MQ=60)
                      TTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCC                                                              >  SRR3722089.1125918/1‑100 (MQ=60)
                       TCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCA                                                             <  SRR3722089.166766/100‑1 (MQ=60)
                        CTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAActgtctcttatacacatctgacgct                                                            >  SRR3722089.341216/1‑75 (MQ=60)
                        CTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAC                                                            >  SRR3722089.1016743/1‑100 (MQ=60)
                                         AAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTT                                           <  SRR3722089.90105/100‑1 (MQ=60)
                                            GCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTAT                                        <  SRR3722089.212590/100‑1 (MQ=60)
                                                          TTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCC                          <  SRR3722089.1154027/100‑1 (MQ=60)
                                                             GATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGA                       >  SRR3722089.951160/1‑100 (MQ=60)
                                                                 TGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTCCTCTCGCCCCGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACT                   <  SRR3722089.743337/100‑1 (MQ=60)
                                                                 TGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACT                   <  SRR3722089.366581/100‑1 (MQ=60)
                                                                 TGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACT                   <  SRR3722089.634943/100‑1 (MQ=60)
                                                                 TGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACT                   <  SRR3722089.923562/100‑1 (MQ=60)
                                                                  GACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT                  <  SRR3722089.99793/100‑1 (MQ=60)
                                                                      ATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACC              <  SRR3722089.918371/100‑1 (MQ=60)
                                                                        AACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCC            <  SRR3722089.990916/100‑1 (MQ=60)
                                                                        AACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCC            <  SRR3722089.1210067/100‑1 (MQ=60)
                                                                        AACCTTCACGAAAAAATTAACTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCC            <  SRR3722089.730642/100‑1 (MQ=60)
                                                                                   AAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCCGTCCAAACT  <  SRR3722089.560169/100‑1 (MQ=60)
                                                                                        |                                                                                              
TTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTTTACCCCCGTCCAAACT  >  NZ_CP009273/3103783‑3103965

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