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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

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,1030AG94.1% 94.8 / ‑0.4 34E129G (GAA→GGA) uxuBfructuronate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (16/16);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-01

TGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543019‑4543187
                                                                                    |                                                                                    
tGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAg                                                                                 <  2:40111/90‑1 (MQ=255)
tGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAg                                                                                 <  2:284526/90‑1 (MQ=255)
tGCATCCGGAAATCGATGGCTGCGAAGGTACTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAg                                                                                 <  2:677739/90‑1 (MQ=255)
 gCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAgg                                                                                >  1:635731/1‑90 (MQ=255)
 gCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAgg                                                                                >  2:314224/1‑90 (MQ=255)
     ccGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            <  1:672043/90‑1 (MQ=255)
        gAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTgc                                                                         <  2:462966/90‑1 (MQ=255)
            tCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTg                                                                     <  2:455895/90‑1 (MQ=255)
                             aTTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCaa                                                    >  1:499123/1‑90 (MQ=255)
                             aTTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCaa                                                    >  1:607520/1‑90 (MQ=255)
                              ttCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAAc                                                   <  2:378532/90‑1 (MQ=255)
                                   aaCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAt                                              <  1:517182/90‑1 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             <  1:357983/90‑1 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:441252/1‑90 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:220240/1‑90 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:141729/1‑90 (MQ=255)
                                         aTGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAAt                                        >  1:608824/1‑90 (MQ=255)
                                          tGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATa                                       <  1:709923/90‑1 (MQ=255)
                                            gcgcGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAAc                                     >  2:654401/1‑90 (MQ=255)
                                                gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  1:753596/90‑1 (MQ=255)
                                                gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  1:382898/90‑1 (MQ=255)
                                                gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  1:453566/90‑1 (MQ=255)
                                                    gCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGAt                             >  2:765619/1‑90 (MQ=255)
                                                           gCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCAc                      <  1:630107/90‑1 (MQ=255)
                                                           gCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCAc                      <  2:455418/90‑1 (MQ=255)
                                                                  tctctcTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTgg               <  2:499123/90‑1 (MQ=255)
                                                                  tctctcTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTgg               <  1:771166/90‑1 (MQ=255)
                                                                     ctctAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGaaa            >  2:474572/1‑90 (MQ=255)
                                                                           cGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGa      >  2:3749/1‑90 (MQ=255)
                                                                            ggTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGAc     <  1:513500/90‑1 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  2:447841/1‑90 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  2:352697/1‑90 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:713641/1‑90 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:283711/1‑90 (MQ=255)
                                                                                    |                                                                                    
TGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543019‑4543187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
GATGAAAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  >  NZ_CP009273/4543005‑4543194
                                                                                                  |                                                                                           
GATGAAAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGA                                                                                            >  SRR3722079.811473/1‑100 (MQ=60)
GATGAAAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGA                                                                                            >  SRR3722079.747034/1‑100 (MQ=60)
     AAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGG                                                                                       >  SRR3722079.642259/1‑100 (MQ=60)
                   CCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATG                                                                         <  SRR3722079.679074/100‑1 (MQ=60)
                                 CTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAA                                                           >  SRR3722079.503591/1‑100 (MQ=60)
                                 CTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAA                                                           >  SRR3722079.613646/1‑100 (MQ=60)
                                             TCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAAT                                               >  SRR3722079.614978/1‑100 (MQ=60)
                                                 AACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACC                                           <  SRR3722079.521934/100‑1 (MQ=60)
                                                  ACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCC                                          <  SRR3722079.360785/100‑1 (MQ=60)
                                                        TGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGAT                                    <  SRR3722079.717462/100‑1 (MQ=60)
                                                              GTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCA                              <  SRR3722079.457537/100‑1 (MQ=60)
                                                              GTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCA                              <  SRR3722079.386022/100‑1 (MQ=60)
                                                              GTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCA                              <  SRR3722079.761650/100‑1 (MQ=60)
                                                                         GCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAA                   <  SRR3722079.636545/100‑1 (MQ=60)
                                                                                TCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGAC            <  SRR3722079.779423/100‑1 (MQ=60)
                                                                                   CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722079.721240/1‑100 (MQ=60)
                                                                                   CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722079.286063/1‑100 (MQ=60)
                                                                                          GGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  <  SRR3722079.518188/100‑1 (MQ=60)
                                                                                          GGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  >  SRR3722079.741646/1‑100 (MQ=60)
                                                                                                  |                                                                                           
GATGAAAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  >  NZ_CP009273/4543005‑4543194

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