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 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% 197.8 / NA 51intergenic (‑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 (21/30);  total (21/30)

AAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTT  >  NZ_CP009273/2061257‑2061428
                                                                                     |                                                                                      
aaaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACAtt                                                                                    >  2:727840/1‑90 (MQ=255)
 aaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTa                                                                                   >  2:76346/1‑90 (MQ=255)
   gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAgt                                                                                 >  1:102139/1‑90 (MQ=255)
              tctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGaa                                                                      <  2:785772/90‑1 (MQ=255)
                     ttCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAAt                                                               <  1:749557/90‑1 (MQ=255)
                          cATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCa                                                          <  2:84037/90‑1 (MQ=255)
                          cATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCa                                                          <  2:1041538/90‑1 (MQ=255)
                             tctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                       <  2:47808/90‑1 (MQ=255)
                             tctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                       <  1:361626/90‑1 (MQ=255)
                             tctcAGACACAGCCTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                       <  1:1173220/90‑1 (MQ=255)
                              ctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCg                                                      >  2:1186674/1‑90 (MQ=255)
                               tcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCgg                                                     <  1:110070/90‑1 (MQ=255)
                               tcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCgg                                                     <  2:1232523/90‑1 (MQ=255)
                                 aGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAg                                                   >  2:570756/1‑90 (MQ=255)
                                    cacaGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAc                                                <  2:1227915/90‑1 (MQ=255)
                                      caGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCa                                              >  2:422324/1‑90 (MQ=255)
                                        gCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACt                                            >  2:900779/1‑90 (MQ=255)
                                           ttCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGt                                         >  2:842063/1‑90 (MQ=255)
                                               aCTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcatc                                     <  2:683346/90‑1 (MQ=255)
                                               aCTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcatc                                     <  1:384660/90‑1 (MQ=255)
                                                 tCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTTGTCATCAc                                   <  1:591429/90‑1 (MQ=255)
                                                 tCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCAc                                   <  2:233100/90‑1 (MQ=255)
                                                  ccATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACg                                  <  2:1075700/90‑1 (MQ=255)
                                                  ccATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACg                                  >  2:22405/1‑90 (MQ=255)
                                                   cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                                 <  1:727840/90‑1 (MQ=255)
                                                   cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                                 <  1:116420/90‑1 (MQ=255)
                                                   cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                                 <  2:333422/90‑1 (MQ=255)
                                                    aTGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                                 >  2:486656/1‑89 (MQ=255)
                                                    aTGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                                 <  1:486656/89‑1 (MQ=255)
                                                    aTGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTa                                >  2:1220440/1‑90 (MQ=255)
                                                       gTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGGGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTAttt                             >  1:910111/1‑90 (MQ=255)
                                                        tAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTta                            <  2:1148187/90‑1 (MQ=255)
                                                         aaTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtat                           >  2:1209224/1‑90 (MQ=255)
                                                          aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCa                                                          <  2:1143940/58‑1 (MQ=255)
                                                          aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCa                                                          >  1:1143940/1‑58 (MQ=255)
                                                          aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtata                          >  1:1149655/1‑90 (MQ=255)
                                                            ggCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata                        >  1:7720/1‑90 (MQ=255)
                                                                cAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCAt                    <  1:1183990/90‑1 (MQ=255)
                                                                 aaaTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg                   >  2:855/1‑90 (MQ=255)
                                                                 aaaTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg                   >  2:80488/1‑90 (MQ=255)
                                                                 aaaTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg                   >  2:99337/1‑90 (MQ=255)
                                                                 aaaTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATg                   >  2:521976/1‑90 (MQ=255)
                                                                   aTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGaa                 <  2:147242/90‑1 (MQ=255)
                                                                   aTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGaa                 <  1:941279/90‑1 (MQ=255)
                                                                   aTACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGaa                 <  1:714038/90‑1 (MQ=255)
                                                                        aaGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAt            <  2:150526/90‑1 (MQ=255)
                                                                        aaGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAt            <  1:326645/90‑1 (MQ=255)
                                                                             gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg       >  2:24272/1‑90 (MQ=255)
                                                                                cAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTgtgt    <  1:451089/90‑1 (MQ=255)
                                                                                cAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTgtgt    <  1:76346/90‑1 (MQ=255)
                                                                                  ggAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGttt  <  1:22405/90‑1 (MQ=255)
                                                                                     |                                                                                      
AAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTT  >  NZ_CP009273/2061257‑2061428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
GAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAG  >  NZ_CP009273/2061244‑2061438
                                                                                                  |                                                                                                
GAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAAC                                                                                                 <  SRR3722089.818272/100‑1 (MQ=60)
     TCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAG                                                                                            <  SRR3722089.104207/100‑1 (MQ=60)
      CTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGT                                                                                           >  SRR3722089.103410/1‑100 (MQ=60)
                                  TTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGG                                                               <  SRR3722089.759477/100‑1 (MQ=60)
                                          TCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAC                                                       <  SRR3722089.366084/100‑1 (MQ=60)
                                          TCTCAGACACAGCCTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAC                                                       <  SRR3722089.1191083/100‑1 (MQ=60)
                                            TCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTC                                                     <  SRR3722089.111417/100‑1 (MQ=60)
                                                          CTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGGGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTT                                       >  SRR3722089.923347/1‑100 (MQ=60)
                                                            ACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTAT                                     <  SRR3722089.389352/100‑1 (MQ=60)
                                                             CTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCtgtctcttatacacatctgac                                    >  SRR3722089.1161350/1‑79 (MQ=60)
                                                             CTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATA                                    >  SRR3722089.1167157/1‑100 (MQ=60)
                                                              TCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTTGTCATCACGTATTTATAT                                   <  SRR3722089.598603/100‑1 (MQ=60)
                                                               CCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATA                                  >  SRR3722089.7806/1‑100 (MQ=60)
                                                                CATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATAT                                 <  SRR3722089.492432/100‑1 (MQ=60)
                                                                CATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATAT                                 <  SRR3722089.737303/100‑1 (MQ=60)
                                                                CATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATAT                                 <  SRR3722089.117840/100‑1 (MQ=60)
                                                                             CAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTG                    <  SRR3722089.1202016/100‑1 (MQ=60)
                                                                                ATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTG                 <  SRR3722089.723326/100‑1 (MQ=60)
                                                                                ATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTG                 <  SRR3722089.955067/100‑1 (MQ=60)
                                                                                     AAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTT            <  SRR3722089.330748/100‑1 (MQ=60)
                                                                                             CAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATC    <  SRR3722089.456501/100‑1 (MQ=60)
                                                                                             CAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATC    <  SRR3722089.77275/100‑1 (MQ=60)
                                                                                               GGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAG  <  SRR3722089.22672/100‑1 (MQ=60)
                                                                                                  |                                                                                                
GAAGCTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAG  >  NZ_CP009273/2061244‑2061438

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