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,038,183 (T)5→4 intergenic (+78/‑236) BW25113_RS10330 → / → yeeJ tRNA‑Asn/inverse autotransporter adhesin YeeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,038,1790T.100.0% 114.6 / NA 27intergenic (+74/‑240)BW25113_RS10330/yeeJtRNA‑Asn/inverse autotransporter adhesin YeeJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (10/17);  total (10/17)

GAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAA  >  NZ_CP009273/2038097‑2038264
                                                                                  |                                                                                     
gaggagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaaca                                                                                >  2:607875/1‑86 (MQ=255)
gaggagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaaca                                                                                <  1:613333/90‑5 (MQ=255)
gaggagCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGttttaaca                                                                                >  2:846797/1‑86 (MQ=255)
        aaaTTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAgc                                                                       <  2:435074/90‑1 (MQ=255)
              tAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTc                                                                 <  2:87285/90‑1 (MQ=255)
               aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                                <  2:761172/90‑1 (MQ=255)
               aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                                <  2:953243/90‑1 (MQ=255)
               aaaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCg                                                                <  1:607875/90‑1 (MQ=255)
                aaaaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGa                                                               <  2:395413/90‑1 (MQ=255)
                  aaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGAcc                                                             <  2:52488/90‑1 (MQ=255)
                  aaTTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGAcc                                                             <  2:282925/90‑1 (MQ=255)
                       gCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTAt                                                        >  2:322488/1‑90 (MQ=255)
                                     aTGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                          <  1:1062443/90‑1 (MQ=255)
                                     aTGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                          <  1:640139/90‑1 (MQ=255)
                                             gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTc                                                                 <  1:990019/59‑1 (MQ=255)
                                             gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTc                                                                 >  2:990019/1‑59 (MQ=255)
                                             gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                          <  2:1179152/82‑1 (MQ=255)
                                             gACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGttt                                          >  1:1179152/1‑82 (MQ=255)
                                               ccTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAc                                <  1:322507/90‑1 (MQ=255)
                                                     ttGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATc                          <  2:905236/90‑1 (MQ=255)
                                                       gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc                        >  1:676476/1‑90 (MQ=255)
                                                       gAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAc                        >  1:199747/1‑90 (MQ=255)
                                                       gAACATTGTTTTTTAACGGATAGCGGGGTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTTAGGGCGGTTTGCTATGAACGCCATCAc                        >  1:1075019/1‑90 (MQ=255)
                                                         aCATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATc                          >  2:454781/1‑86 (MQ=255)
                                                         aCATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATc                          <  1:454781/86‑1 (MQ=255)
                                                            ttGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCAttt                   >  2:821268/1‑90 (MQ=255)
                                                                 ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAggg                                               <  2:224206/57‑1 (MQ=255)
                                                                 ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAggg                                               >  1:224206/1‑57 (MQ=255)
                                                                 ttttAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc              <  1:1165205/90‑1 (MQ=255)
                                                                             gCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATaaaa  >  1:903207/1‑90 (MQ=255)
                                                                                  |                                                                                     
GAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAA  >  NZ_CP009273/2038097‑2038264

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  >  NZ_CP009273/2038086‑2038281
                                                                                                 |                                                                                                  
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAAC                                                                                                 <  SRR3722089.1032157/100‑1 (MQ=60)
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAAC                                                                                                 <  SRR3722089.82042/100‑1 (MQ=60)
           GAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGC                                                                                      <  SRR3722089.620990/100‑1 (MQ=60)
                         tcagatgtgtataagagacaGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATG                                                                        <  SRR3722089.1004712/80‑1 (MQ=60)
                          AAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGG                                                                       <  SRR3722089.615425/100‑1 (MQ=60)
                                              CAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGA                                                   >  SRR3722089.1197115/1‑100 (MQ=60)
                                                ATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAC                                                 <  SRR3722089.1078394/100‑1 (MQ=60)
                                                ATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAAC                                                 <  SRR3722089.648314/100‑1 (MQ=60)
                                                        GACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAC                                         >  SRR3722089.685244/1‑100 (MQ=60)
                                                        GACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCAC                                         >  SRR3722089.202268/1‑100 (MQ=60)
                                                        GACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGGTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTTAGGGCGGTTTGCTATGAACGCCATCAC                                         >  SRR3722089.1091185/1‑100 (MQ=60)
                                                          CCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCA                                       <  SRR3722089.326571/100‑1 (MQ=60)
                                                                TTGAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCC                                 <  SRR3722089.460226/100‑1 (MQ=60)
                                                                  GAACATTGTTTTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTActgtctcttatacacatctgacg                               >  SRR3722089.227051/1‑77 (MQ=60)
                                                                            TTTTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAA                     <  SRR3722089.1182923/100‑1 (MQ=60)
                                                                              TTAACGGATAGCGGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAA                   >  SRR3722089.916320/1‑100 (MQ=60)
                                                                                           GGGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTC       <  SRR3722089.726905/100‑1 (MQ=60)
                                                                                            GGTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCA      >  SRR3722089.370406/1‑100 (MQ=60)
                                                                                                TTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  <  SRR3722089.209834/100‑1 (MQ=60)
                                                                                                TTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  <  SRR3722089.1145140/100‑1 (MQ=60)
                                                                                                TTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  <  SRR3722089.832812/100‑1 (MQ=60)
                                                                                                 |                                                                                                  
GAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAG  >  NZ_CP009273/2038086‑2038281

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