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 524,350 T→C F85F (TTT→TTC BW25113_RS24960 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273524,3500TC100.0% 100.4 / NA 27F85F (TTT→TTCBW25113_RS24960hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/8);  total (19/8)

CATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAAT  >  NZ_CP009273/524265‑524436
                                                                                     |                                                                                      
cATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTa                                                                                   <  2:68401/91‑1 (MQ=255)
   ggCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACa                                                                                <  2:860862/91‑1 (MQ=255)
   ggCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACa                                                                                >  1:860862/1‑91 (MQ=255)
   ggCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACa                                                                                >  2:120751/1‑91 (MQ=255)
                 tttCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATggggg                                                                  <  2:487591/91‑1 (MQ=255)
                     cGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAgg                                                              >  2:856828/1‑91 (MQ=255)
                           tAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAg                                                        >  2:537388/1‑91 (MQ=255)
                                  ggATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACAt                                                 >  1:416204/1‑91 (MQ=255)
                                     tatatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGa                                              >  2:437720/1‑91 (MQ=255)
                                       tatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACa                                            >  1:290696/1‑91 (MQ=255)
                                       tatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACa                                            <  2:837796/91‑1 (MQ=255)
                                            atatACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTccc                                                                                      <  1:848513/44‑1 (MQ=255)
                                            atatACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAg                                       >  2:482923/1‑91 (MQ=255)
                                            atatACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAg                                       >  1:644550/1‑91 (MQ=255)
                                             tataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAgg                                      >  2:323185/1‑91 (MQ=255)
                                              ataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAACTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAggg                                     >  2:175738/1‑91 (MQ=255)
                                              ataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAggg                                     >  2:515537/1‑91 (MQ=255)
                                                 cacaAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCa                                  >  1:293710/1‑91 (MQ=255)
                                                  acaAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAg                                 <  2:167460/91‑1 (MQ=255)
                                                      aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTa                                                                                   >  2:617703/1‑37 (MQ=38)
                                                      aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTa                                                                                   <  1:617703/37‑1 (MQ=38)
                                                      aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGaaa                                                     <  1:133709/67‑1 (MQ=255)
                                                      aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGaaa                                                     >  2:133709/1‑67 (MQ=255)
                                                        gTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACa                           >  2:487656/1‑91 (MQ=255)
                                                             agaTGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTgg                      >  2:7726/1‑91 (MQ=255)
                                                                  attattTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGaaa                 >  2:425963/1‑91 (MQ=255)
                                                                     attTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATc              >  2:395187/1‑91 (MQ=255)
                                                                            tCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCaa       <  1:395187/91‑1 (MQ=255)
                                                                                 ttttcCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAAt  <  1:537388/91‑1 (MQ=255)
                                                                                     |                                                                                      
CATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAAT  >  NZ_CP009273/524265‑524436

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
GTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAA  >  NZ_CP009273/524258‑524450
                                                                                            |                                                                                                    
GTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACA                                                                                              >  SRR3722118.880012/1‑101 (MQ=60)
    tccgtctcgtgggctcggagatgagtataagagacagGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTCTCTTATTTTCCCCTAACAAATG                                                                                          <  SRR3722118.867302/64‑1 (MQ=60)
       ataggagtccgtctcgtgggctcggagatgtgtataagagacagATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATG                                                                                       <  SRR3722118.630940/57‑1 (MQ=60)
                               ACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACAT                                                               >  SRR3722118.424966/1‑101 (MQ=60)
                                    GCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACA                                                          >  SRR3722118.296654/1‑101 (MQ=60)
                                     tgtataagagacagATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAG                                                         <  SRR3722118.136241/87‑1 (MQ=60)
                                         GGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAG                                                     >  SRR3722118.658454/1‑101 (MQ=60)
                                              TATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCA                                                >  SRR3722118.299726/1‑101 (MQ=60)
                                                                                   TCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCC           <  SRR3722118.403482/101‑1 (MQ=60)
                                                                                        TTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATA      <  SRR3722118.548790/101‑1 (MQ=60)
                                                                                            CCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGctgtctcttatacacatctccgagcccacgagacgga  >  SRR3722118.734779/1‑64 (MQ=60)
                                                                                            CCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAA  >  SRR3722118.709819/1‑101 (MQ=60)
                                                                                            |                                                                                                    
GTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAA  >  NZ_CP009273/524258‑524450

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