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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,061,342 T→A intergenic (‑576/‑774) BW25113_RS10415 ← / → BW25113_RS24005 IS5‑like element IS5 family transposase/TonB‑dependent receptor plug domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,061,3420TA100.0% 104.1 / NA 28intergenic (‑576/‑774)BW25113_RS10415/BW25113_RS24005IS5‑like element IS5 family transposase/TonB‑dependent receptor plug domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (8/20);  total (8/20)

AAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTT  >  NZ_CP009273/2061258‑2061428
                                                                                    |                                                                                      
aaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTa                                                                                   >  2:570866/1‑90 (MQ=255)
  gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 >  2:548496/1‑90 (MQ=255)
  gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 <  1:274021/90‑1 (MQ=255)
  gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 <  1:548496/90‑1 (MQ=255)
              ctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAg                                                                     >  2:538527/1‑90 (MQ=255)
                    ttCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAAt                                                               <  2:450758/90‑1 (MQ=255)
                    ttCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAAt                                                               <  1:538527/90‑1 (MQ=255)
                         cATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCa                                                          <  2:105386/90‑1 (MQ=255)
                            tctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                       <  1:20877/90‑1 (MQ=255)
                            tctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                       <  1:12061/90‑1 (MQ=255)
                              tcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCgg                                                     <  1:532712/90‑1 (MQ=255)
                                     caGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 <  1:462589/55‑1 (MQ=255)
                                     caGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 >  2:462589/1‑55 (MQ=255)
                                     caGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCa                                              >  1:317858/1‑90 (MQ=255)
                                       gCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 >  2:411133/1‑53 (MQ=255)
                                       gCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 <  1:411133/53‑1 (MQ=255)
                                                tCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCAc                                   <  2:470476/90‑1 (MQ=255)
                                                tCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCAc                                   <  2:521622/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                                 <  1:82907/90‑1 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTa                                >  1:572593/1‑90 (MQ=255)
                                                    tGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTAt                               >  2:500211/1‑90 (MQ=255)
                                                           ggCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata                        <  2:275217/90‑1 (MQ=255)
                                                                aaaTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg                   <  1:500211/90‑1 (MQ=255)
                                                                aaaTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg                   <  2:541645/90‑1 (MQ=255)
                                                                     cGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCt              <  1:223429/90‑1 (MQ=255)
                                                                               cAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTgtgt    <  2:100667/90‑1 (MQ=255)
                                                                               cAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTgtgt    <  2:75265/90‑1 (MQ=255)
                                                                                 ggAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCGCTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGttt  <  2:262771/90‑1 (MQ=255)
                                                                                    |                                                                                      
AAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTT  >  NZ_CP009273/2061258‑2061428

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  >  NZ_CP009273/2061245‑2061426
                                                                                                 |                                                                                    
AAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACA                                                                                    <  SRR3722109.219219/100‑1 (MQ=60)
 AGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACAT                                                                                   <  SRR3722109.240178/100‑1 (MQ=60)
    TCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAG                                                                                <  SRR3722109.301534/100‑1 (MQ=60)
               gcagcgtcagatgtgtataagagacagCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGA                                                                     <  SRR3722109.417124/73‑1 (MQ=60)
               agcgtcagatgtgtataagagacagTTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGA                                                                     <  SRR3722109.469532/75‑1 (MQ=60)
               GCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGA                                                                     <  SRR3722109.277458/100‑1 (MQ=60)
               GCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGA                                                                     <  SRR3722109.556825/100‑1 (MQ=60)
                                 TTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGG                                                   <  SRR3722109.546713/100‑1 (MQ=60)
                                        TTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCA                                            >  SRR3722109.322049/1‑100 (MQ=60)
                                         TCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAC                                           <  SRR3722109.12233/100‑1 (MQ=60)
                                         TCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAC                                           <  SRR3722109.21151/100‑1 (MQ=60)
                                           TCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTC                                         <  SRR3722109.540826/100‑1 (MQ=60)
                                                      TTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTA                              >  SRR3722109.581282/1‑100 (MQ=60)
                                                               CATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATAT                     <  SRR3722109.83929/100‑1 (MQ=60)
                                                                             AAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGT       <  SRR3722109.507765/100‑1 (MQ=60)
                                                                                  CGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  <  SRR3722109.226094/100‑1 (MQ=60)
                                                                                                 |                                                                                    
AAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  >  NZ_CP009273/2061245‑2061426

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