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 F3 I122 R1 193 20.9 1228102 94.7% 1163012 84.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,649,376 (T)5→4 coding (277/327 nt) hdeB ← acid‑activated periplasmic chaperone HdeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,649,3720T.100.0% 110.1 / NA 26coding (281/327 nt)hdeBacid‑activated periplasmic chaperone HdeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/14);  total (12/14)

CAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCAC  >  NZ_CP009273/3649295‑3649456
                                                                             |                                                                                    
cAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGtt                                                                         <  1:449074/90‑1 (MQ=255)
cAACTTGATAGTCCTGTTTACTTGCACCTCATAAATTCGGCAAGTCATTAGATGATTGATTTTTAAAGGTATACAAATTTTCTGCGGGtt                                                                         <  1:398719/90‑1 (MQ=255)
       cTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTa                                                                  >  2:372777/1‑90 (MQ=255)
        tACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTAc                                                                 <  2:226263/90‑1 (MQ=255)
               ttttACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTc                                                          >  2:151605/1‑90 (MQ=255)
               ttttACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTc                                                          >  1:547534/1‑90 (MQ=255)
                   aCTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATc                                                      <  1:506841/90‑1 (MQ=255)
                      tGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                                                   >  2:326371/1‑90 (MQ=255)
                      tGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                                                   >  2:188809/1‑90 (MQ=255)
                      tGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                                                   >  1:342484/1‑90 (MQ=255)
                      tGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                                                   >  1:350972/1‑90 (MQ=255)
                           cTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAgg                                              <  1:298759/90‑1 (MQ=255)
                                       gCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtgag                                  >  2:613835/1‑90 (MQ=255)
                                       gCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtgag                                  >  1:253218/1‑90 (MQ=255)
                                              aTTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt                           <  1:77978/90‑1 (MQ=255)
                                              aTTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt                           <  2:454231/90‑1 (MQ=255)
                                               ttAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATc                                                      >  2:204838/1‑62 (MQ=255)
                                               ttAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATc                                                      <  1:204838/62‑1 (MQ=255)
                                                      cTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaa                   <  1:293168/90‑1 (MQ=255)
                                                       ttGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaaa                  >  2:394809/1‑90 (MQ=255)
                                                                    gTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGc     >  1:472202/1‑90 (MQ=255)
                                                                     tataCAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGcc    <  1:394809/90‑1 (MQ=255)
                                                                     tataCAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGcc    <  1:326371/90‑1 (MQ=255)
                                                                       taCAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCAc  <  2:122002/90‑1 (MQ=255)
                                                                       taCAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCAc  <  1:222044/90‑1 (MQ=255)
                                                                       taCAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCAc  <  2:603237/90‑1 (MQ=255)
                                                                             |                                                                                    
CAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCAC  >  NZ_CP009273/3649295‑3649456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATAC  >  NZ_CP009273/3649295‑3649466
                                                                                 |                                                                                          
CAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAG                                                                         <  SRR3721989.458445/100‑1 (MQ=60)
CAACTTGATAGTCCTGTTTACTTGCACCTCATAAATTCGGCAAGTCATTAGATGATTGATTTTTAAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAG                                                                         <  SRR3721989.406834/100‑1 (MQ=60)
     TGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTC                                                                    >  SRR3721989.559361/1‑100 (MQ=60)
            CCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACT                                                             >  SRR3721989.349034/1‑100 (MQ=60)
            CCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACT                                                             >  SRR3721989.357709/1‑100 (MQ=60)
                   gatgtgtataagagacagCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAA                                                      <  SRR3721989.208552/82‑1 (MQ=60)
                   ACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAA                                                      <  SRR3721989.517650/100‑1 (MQ=60)
                           CTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTG                                              <  SRR3721989.304358/100‑1 (MQ=60)
                             CATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAG                                            >  SRR3721989.257724/1‑100 (MQ=60)
                                              ATTAGATGCTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAG                           <  SRR3721989.79375/100‑1 (MQ=60)
                                                      CTTGATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTA                   <  SRR3721989.298599/100‑1 (MQ=60)
                                                          ATTTTTGAAGGTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGC               >  SRR3721989.482129/1‑100 (MQ=60)
                                                                    GTATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATA     >  SRR3721989.591676/1‑100 (MQ=60)
                                                                     TATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATAT    <  SRR3721989.332585/100‑1 (MQ=60)
                                                                     TATACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATAT    <  SRR3721989.402812/100‑1 (MQ=60)
                                                                       TACAAATTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATAC  <  SRR3721989.226060/100‑1 (MQ=60)
                                                                                 |                                                                                          
CAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATAC  >  NZ_CP009273/3649295‑3649466

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