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 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 611,351 Δ1 bp coding (1739/3882 nt) entF → enterobactin non‑ribosomal peptide synthetase EntF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273611,3480C.100.0% 101.2 / NA 23coding (1736/3882 nt)entFenterobactin non‑ribosomal peptide synthetase EntF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (11/12);  total (11/12)

AAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCT  >  NZ_CP009273/611258‑611428
                                                                                          |                                                                                
aaGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGc                                                                                  <  2:560019/91‑1 (MQ=255)
 aGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGcc                                                                                 <  2:786426/91‑1 (MQ=255)
        cgtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTac                                                                         >  1:13052/1‑91 (MQ=255)
        cgtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTac                                                                         <  2:389684/91‑1 (MQ=255)
                           acagcgacgaTTAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGcgccg                                                      >  1:37151/5‑91 (MQ=255)
                             accGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCt                                                    >  2:50703/1‑91 (MQ=255)
                             accGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCt                                                    >  1:247266/1‑91 (MQ=255)
                                 acgaTCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCaa                                                <  1:38211/91‑1 (MQ=255)
                                           gCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAAc                                      <  2:22721/91‑1 (MQ=255)
                                              gcgcTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGc                                   >  1:634907/1‑91 (MQ=255)
                                              gcgcTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGc                                   >  2:478779/1‑91 (MQ=255)
                                                     aGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCacac                            <  1:419704/91‑1 (MQ=255)
                                                        gATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGc                         <  1:834051/91‑1 (MQ=255)
                                                         aTGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCt                        <  1:804133/91‑1 (MQ=255)
                                                         aTGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCt                        <  1:193200/91‑1 (MQ=255)
                                                           gTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTta                      >  1:101445/1‑91 (MQ=255)
                                                              cccAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat                   <  2:368041/91‑1 (MQ=255)
                                                              cccAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat                   <  1:20731/91‑1 (MQ=255)
                                                                     tAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCttt            <  1:289818/91‑1 (MQ=255)
                                                                     tAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCttt            >  2:685997/1‑91 (MQ=255)
                                                                          aGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGc                                   <  2:250654/63‑1 (MQ=255)
                                                                          aGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGc                                   >  1:250654/1‑63 (MQ=255)
                                                                          aGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCacac                            >  2:138581/1‑70 (MQ=255)
                                                                          aGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCacac                            <  1:138581/70‑1 (MQ=255)
                                                                               ttGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCt  >  1:118840/1‑91 (MQ=255)
                                                                                          |                                                                                
AAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCT  >  NZ_CP009273/611258‑611428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
GGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGG  >  NZ_CP009273/611256‑611449
                                                                                               |                                                                                                  
GGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTAC                                                                                              >  SRR3722118.13303/1‑101 (MQ=60)
                   gtggtaattaacagCGACGATTAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCG                                                                           >  SRR3722118.37828/15‑101 (MQ=60)
                     GTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCT                                                                         >  SRR3722118.252276/1‑101 (MQ=60)
                                   ACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAAC                                                           <  SRR3722118.38911/101‑1 (MQ=60)
                                      ATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGC                                                        >  SRR3722118.648597/1‑101 (MQ=60)
                                                   CTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTA                                           >  SRR3722118.103308/1‑101 (MQ=60)
                                                       ataagagacagCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATC                                       <  SRR3722118.141217/90‑1 (MQ=60)
                                                       AGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATC                                       <  SRR3722118.428564/101‑1 (MQ=60)
                                                          GATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCCTATATCATC                                    <  SRR3722118.852434/101‑1 (MQ=60)
                                                           ATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCT                                   <  SRR3722118.197028/101‑1 (MQ=60)
                                                           ATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCT                                   <  SRR3722118.821790/101‑1 (MQ=60)
                                                                CCCAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACC                              <  SRR3722118.21117/101‑1 (MQ=60)
                                                                  CAATTTAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCctgtctcttatacacatc                            >  SRR3722118.255754/1‑83 (MQ=60)
                                                                       TAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCT                       <  SRR3722118.295763/101‑1 (MQ=60)
                                                                       TAACAAGCCTTTGCTATAACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCT                       >  SRR3722118.121068/1‑101 (MQ=60)
                                                                                         ACGCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAG     <  SRR3722118.700922/101‑1 (MQ=60)
                                                                                           cgCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGG    <  SRR3722118.421933/99‑1 (MQ=60)
                                                                                             CCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGG  >  SRR3722118.568899/1‑101 (MQ=60)
                                                                                               |                                                                                                  
GGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGG  >  NZ_CP009273/611256‑611449

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