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 I194 R1 197 32.2 1845014 95.2% 1756453 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,1030AG94.4% 112.2 / ‑0.6 36E129G (GAA→GGA) uxuBfructuronate reductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (16/18);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543019‑4543187
                                                                                    |                                                                                    
tGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAg                                                                                 <  2:896875/90‑1 (MQ=255)
     ccGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            <  2:597433/90‑1 (MQ=255)
                  gCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCgg                                                               >  1:717019/1‑90 (MQ=255)
                     gCGAGGGTATTCTCAACGCGATGGCGGGTCCGCAAACGGCGATTGTCCCGCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCaa                                                            <  2:46694/90‑1 (MQ=255)
                     gCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            >  1:760905/1‑74 (MQ=255)
                     gCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            <  2:760905/74‑1 (MQ=255)
                         aGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCgg                                                        <  1:305455/90‑1 (MQ=255)
                          ggTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGt                                                       <  2:433778/90‑1 (MQ=255)
                           gTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTc                                                      <  2:621287/90‑1 (MQ=255)
                            tATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCa                                                     >  2:237484/1‑90 (MQ=255)
                               tctcAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACt                                                  <  2:664980/90‑1 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            >  2:494324/1‑59 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAc                                                                            <  1:494324/59‑1 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:683114/1‑90 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  1:823967/1‑90 (MQ=255)
                                    aCGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtc                                             >  2:76366/1‑90 (MQ=255)
                                     cgcgATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAtct                                            >  1:96236/1‑90 (MQ=255)
                                          tGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATa                                       >  2:917204/1‑90 (MQ=255)
                                          tGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATa                                       <  1:237484/90‑1 (MQ=255)
                                          tGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATa                                       <  1:140787/90‑1 (MQ=255)
                                                gTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGc                                 <  2:854001/90‑1 (MQ=255)
                                                     cAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATc                            <  1:87156/90‑1 (MQ=255)
                                                         cGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGc                        <  1:426066/90‑1 (MQ=255)
                                                          ggCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCa                       <  1:917204/90‑1 (MQ=255)
                                                                 gTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTg                >  2:293255/1‑90 (MQ=255)
                                                                     ctctAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGaaa            <  2:96236/90‑1 (MQ=255)
                                                                           cGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCGGa      >  2:191774/1‑90 (MQ=255)
                                                                            ggTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGAc     <  1:333845/90‑1 (MQ=255)
                                                                            ggTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGAc     <  1:244480/90‑1 (MQ=255)
                                                                             gTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGAcc    >  1:868671/1‑89 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:500799/1‑90 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:238086/1‑90 (MQ=255)
                                                                               cACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:190653/1‑90 (MQ=255)
                                                                               cACGGGAAAAGGCTACCGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTgc  >  1:658976/1‑90 (MQ=255)
                                                                               cACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCt                                <  2:411490/60‑1 (MQ=255)
                                                                               cACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCt                                >  1:411490/1‑60 (MQ=255)
                                                                                    |                                                                                    
TGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC  >  NZ_CP009273/4543019‑4543187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ 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                                                                                            <  SRR3722072.220114/100‑1 (MQ=60)
     AAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGG                                                                                       <  SRR3722072.884478/100‑1 (MQ=60)
      aagcagaagacggcatacgagataatcggacgtctcgtgggctcggagatgtgtataagagacagCGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGC                                                                                      <  SRR3722072.475669/35‑1 (MQ=59)
      AGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACAGGAAAAGGC                                                                                      <  SRR3722072.749723/100‑1 (MQ=60)
         AGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTAC                                                                                   <  SRR3722072.230833/100‑1 (MQ=60)
                   ggagatgtgtataagagacagGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATG                                                                         <  SRR3722072.503258/79‑1 (MQ=60)
                      GAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGG                                                                      >  SRR3722072.729778/1‑100 (MQ=60)
                         ATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCtgtct                                                                   >  SRR3722072.774317/1‑95 (MQ=60)
                                       AGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGAT                                                     <  SRR3722072.310748/100‑1 (MQ=60)
                                        GGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATC                                                    >  SRR3722072.838293/1‑100 (MQ=60)
                                         GTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCT                                                   >  SRR3722072.98020/1‑100 (MQ=60)
                                                        TGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGAT                                    <  SRR3722072.241454/100‑1 (MQ=60)
                                                        TGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGAT                                    <  SRR3722072.143307/100‑1 (MQ=60)
                                                                   CAAACGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATC                         <  SRR3722072.88780/100‑1 (MQ=60)
                                                                       CGGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGA                     <  SRR3722072.433817/100‑1 (MQ=60)
                                                                        GGCGATTGTCTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAA                    <  SRR3722072.932961/100‑1 (MQ=60)
                                                                                 CTCTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACC           >  SRR3722072.883637/1‑100 (MQ=60)
                                                                                   CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722072.509844/1‑100 (MQ=60)
                                                                                   CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722072.242066/1‑100 (MQ=60)
                                                                                   CTCTAACGGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722072.193961/1‑100 (MQ=60)
                                                                                   CTCTAACGGTCACGGGAAAAGGCTACCGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGC         >  SRR3722072.670734/1‑100 (MQ=60)
                                                                                   CTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACctgtctcttatacacatctc         >  SRR3722072.418941/1‑80 (MQ=60)
                                                                                          GGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  <  SRR3722072.339713/100‑1 (MQ=60)
                                                                                          GGTCACGGGAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  <  SRR3722072.248584/100‑1 (MQ=60)
                                                                                                  |                                                                                           
GATGAAAGAAGCGCTGCATCCGGAAATCGATGGCTGCGAAGGTATTCTCAACGCGATGGCGCGTCCGCAAACGGCGATTGTCTCTCTAACGGTCACGGAAAAAGGCTACTGCGCTGATGCGGCAAGCGGTCAACTGGATCTCAATAACCCGCTGATCAAGCACGATCTGGAAAACCCGACTGCGCCGAAG  >  NZ_CP009273/4543005‑4543194

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