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 I195 R1 222 27.1 1551108 96.0% 1489063 86.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,1030AG100.0% 94.1 / NA 27E129G (GAA→GGA) uxuBfructuronate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/16);  total (11/16)

CGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543025‑4543187
                                                                              |                                                                                    
cGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACt                                                                           >  1:229331/1‑90 (MQ=255)
cGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACt                                                                           >  2:205319/1‑90 (MQ=255)
      tCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTg                                                                     >  2:71190/1‑90 (MQ=255)
           ggCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCg                                                                >  1:673040/1‑90 (MQ=255)
            gCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCgg                                                               >  2:585116/1‑90 (MQ=255)
              tGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCa                                                             <  2:472527/90‑1 (MQ=255)
                cGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAg                                                           <  2:402513/90‑1 (MQ=255)
                    ggTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGt                                                       <  2:262711/90‑1 (MQ=255)
                         tctcAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACt                                                  <  1:614366/90‑1 (MQ=255)
                         tctcAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACt                                                  <  2:257762/90‑1 (MQ=255)
                             aaCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAt                                              <  1:377528/90‑1 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCTAGCGGTCAACTGGAtc                                             <  1:402798/90‑1 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:114841/1‑90 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:26483/1‑90 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             <  1:660403/90‑1 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  1:343074/1‑90 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  1:248206/1‑90 (MQ=255)
                              aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             <  1:237379/90‑1 (MQ=255)
                                          gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  2:196070/90‑1 (MQ=255)
                                           tCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCt                                <  2:163022/90‑1 (MQ=255)
                                           tCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCt                                <  1:657214/90‑1 (MQ=255)
                                               cAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATc                            <  1:626874/90‑1 (MQ=255)
                                                aaaCGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCa                           <  2:179039/90‑1 (MQ=255)
                                                   cGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGc                        <  2:712883/90‑1 (MQ=255)
                                                     gCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCAc                      >  2:248883/1‑90 (MQ=255)
                                                                      ggTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGAc     <  1:249172/90‑1 (MQ=255)
                                                                         cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:201128/1‑90 (MQ=255)
                                                                              |                                                                                    
CGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543025‑4543187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  >  NZ_CP009273/4543014‑4543194
                                                                                         |                                                                                           
AGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAC                                                                                   <  SRR3722073.241516/100‑1 (MQ=60)
AGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAC                                                                                   <  SRR3722073.463106/100‑1 (MQ=60)
 GCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACT                                                                                  >  SRR3722073.234114/1‑100 (MQ=60)
            GGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCG                                                                       >  SRR3722073.686799/1‑100 (MQ=60)
                               GGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATC                                                    >  SRR3722073.253411/1‑100 (MQ=60)
                               GGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATC                                                    >  SRR3722073.350343/1‑100 (MQ=60)
                                    TCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAAT                                               <  SRR3722073.627147/100‑1 (MQ=60)
                                        AACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACC                                           <  SRR3722073.385569/100‑1 (MQ=60)
                                         ACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCTAGCGGTCAACTGGATCTCAATAACCC                                          <  SRR3722073.411386/100‑1 (MQ=60)
                                         ACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCC                                          <  SRR3722073.242354/100‑1 (MQ=60)
                                         ACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCC                                          <  SRR3722073.673941/100‑1 (MQ=60)
                                                      TCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCAC                             <  SRR3722073.670691/100‑1 (MQ=60)
                                                          CAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATC                         <  SRR3722073.639876/100‑1 (MQ=60)
                                                                          CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722073.205204/1‑100 (MQ=60)
                                                                                 GGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  <  SRR3722073.254403/100‑1 (MQ=60)
                                                                                         |                                                                                           
AGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  >  NZ_CP009273/4543014‑4543194

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