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 F3 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,543,103 A→G E129G (GAA→GGA)  uxuB → fructuronate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,543,1030AG100.0% 117.5 / NA 35E129G (GAA→GGA) uxuBfructuronate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/14);  total (21/14)

CGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543016‑4543187
                                                                                       |                                                                                    
cgcTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGgaa                                                                                    <  2:809175/90‑1 (MQ=255)
        ccGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            <  2:91941/90‑1 (MQ=255)
         cGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACt                                                                           >  2:311079/1‑90 (MQ=255)
           gAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTgc                                                                         >  1:141110/1‑90 (MQ=255)
            aaaTCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTgcg                                                                        >  2:111526/1‑90 (MQ=255)
             aaTCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTgcgc                                                                       <  2:93623/90‑1 (MQ=255)
               tCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTg                                                                     <  1:697024/90‑1 (MQ=255)
                     gCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCgg                                                               >  1:456638/1‑90 (MQ=255)
                        gCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCaa                                                            >  1:732989/1‑90 (MQ=255)
                        gCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCaa                                                            >  1:856572/1‑90 (MQ=255)
                            aGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCgg                                                        <  2:406792/90‑1 (MQ=255)
                                aTTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCaa                                                    >  1:263680/1‑90 (MQ=255)
                                 ttCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAAc                                                   <  2:141110/90‑1 (MQ=255)
                                       aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:490854/1‑90 (MQ=255)
                                       aCGCGATGGAGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  1:905012/1‑90 (MQ=255)
                                          cgATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCa                                          >  1:152411/1‑90 (MQ=255)
                                          cgATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCa                                          >  1:235506/1‑90 (MQ=255)
                                            aTGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGT‑‑CGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATaa                                      >  1:963404/1‑90 (MQ=255)
                                               gcgcGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAAc                                     >  2:188484/1‑90 (MQ=255)
                                                 gcgTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAAccc                                   >  1:483114/1‑90 (MQ=255)
                                                   gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            >  2:470781/1‑47 (MQ=255)
                                                   gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            <  1:470781/47‑1 (MQ=255)
                                                   gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 >  2:121038/1‑90 (MQ=255)
                                                   gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  2:101949/90‑1 (MQ=255)
                                                            cGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGc                        <  2:856572/90‑1 (MQ=255)
                                                            cGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGc                        <  1:389684/90‑1 (MQ=255)
                                                             ggCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCa                       <  2:483114/90‑1 (MQ=255)
                                                             ggCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCa                       <  1:674784/90‑1 (MQ=255)
                                                              gCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  2:122275/79‑1 (MQ=255)
                                                              gCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 >  1:122275/1‑79 (MQ=255)
                                                                 aTTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGAt                   <  2:31517/90‑1 (MQ=255)
                                                                        ctctAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGaaa            >  2:543470/1‑90 (MQ=255)
                                                                               ggTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAAccc                                   >  2:510422/1‑60 (MQ=255)
                                                                               ggTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAAccc                                   <  1:510422/60‑1 (MQ=255)
                                                                                  cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  2:231255/1‑90 (MQ=255)
                                                                                  cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:921255/1‑90 (MQ=255)
                                                                                       |                                                                                    
CGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543016‑4543187

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
AGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAGT  >  NZ_CP009273/4543011‑4543195
                                                                                            |                                                                                            
ttcgtcggcagcgtcagatgtgtataagagacagGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGC                                                                                       <  SRR3722236.313717/66‑1 (MQ=60)
      GCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGC                                                                                 >  SRR3722236.143150/1‑100 (MQ=60)
             tcgtcggcagcgtcagatgtgtataagagacagGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATG                                                                          <  SRR3722236.477759/67‑1 (MQ=60)
                GAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGG                                                                       >  SRR3722236.463281/1‑100 (MQ=60)
                   ATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAA                                                                    >  SRR3722236.871575/1‑100 (MQ=60)
                   ATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAA                                                                    >  SRR3722236.745474/1‑100 (MQ=60)
                    TCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAG                                                                   <  SRR3722236.708774/100‑1 (MQ=60)
                           CTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAA                                                            >  SRR3722236.267543/1‑100 (MQ=60)
                                  GGTATTCTCAACGCGATGGAGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATC                                                     >  SRR3722236.920969/1‑100 (MQ=60)
                                     ATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCA                                                  >  SRR3722236.239012/1‑100 (MQ=60)
                                     ATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCA                                                  >  SRR3722236.154634/1‑100 (MQ=60)
                                       TCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGT‑‑CGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAA                                              >  SRR3722236.980398/1‑100 (MQ=60)
                                            ACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCC                                           >  SRR3722236.490365/1‑100 (MQ=60)
                                                      cagatgtgtataagagacaGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAG                                 <  SRR3722236.518172/81‑1 (MQ=60)
                                                         TCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCAC                              >  SRR3722236.123996/1‑100 (MQ=60)
                                                                 CGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGA                      <  SRR3722236.395142/100‑1 (MQ=60)
                                                                  GGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAA                     <  SRR3722236.686025/100‑1 (MQ=60)
                                                                             CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC          >  SRR3722236.937527/1‑100 (MQ=60)
                                                                               CTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGC        >  SRR3722236.161528/1‑100 (MQ=60)
                                                                                     GTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAGT  >  SRR3722236.310812/1‑100 (MQ=60)
                                                                                            |                                                                                            
AGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAGT  >  NZ_CP009273/4543011‑4543195

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: