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 834,220 A→G D148D (GAT→GAC ybiX ← PKHD‑type hydroxylase YbiX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273834,2200AG100.0% 112.4 / NA 31D148D (GAT→GACybiXPKHD‑type hydroxylase YbiX
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/14);  total (17/14)

ACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAA  >  NZ_CP009273/834135‑834307
                                                                                     |                                                                                       
aCTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGcc                                                                                    <  2:187750/91‑1 (MQ=255)
  tGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGCTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTg                                                                                  >  1:257886/1‑91 (MQ=255)
    gATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGcc                                                                                >  1:119645/1‑91 (MQ=255)
    gATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGcc                                                                                >  1:214722/1‑91 (MQ=255)
          caTAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGt                                                                          >  1:507170/1‑91 (MQ=255)
          caTAAATGATGCCACTCGTACGCCGCGGGTTAAGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTCACGGGAGt                                                                          >  1:97647/1‑91 (MQ=255)
                  atgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCg                                                                  >  2:86713/1‑91 (MQ=255)
                   tgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                                 >  1:581341/1‑91 (MQ=255)
                   tgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                                 >  1:104762/1‑91 (MQ=255)
                   tgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                                 >  1:116499/1‑91 (MQ=255)
                      cACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgt                                                              <  1:53065/91‑1 (MQ=255)
                       aCTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgtt                                                             >  1:622973/1‑91 (MQ=255)
                            tACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCg                                                        >  1:228532/1‑91 (MQ=255)
                                  gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGTAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGGCCGAAGgtg                                                  <  2:49188/91‑1 (MQ=255)
                                  gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGgtg                                                  >  1:423210/1‑91 (MQ=255)
                                                        tGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgccgc                            <  2:257886/91‑1 (MQ=255)
                                                              ctgctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCcgtcgt                      >  1:320701/1‑91 (MQ=255)
                                                              ctgctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCcgtcgt                      <  1:490672/91‑1 (MQ=255)
                                                                gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                    <  1:267784/91‑1 (MQ=255)
                                                                gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                    <  1:546867/91‑1 (MQ=255)
                                                                gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                    <  1:718750/91‑1 (MQ=255)
                                                                gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                    <  2:440981/91‑1 (MQ=255)
                                                                 ctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGc                   >  1:338443/1‑91 (MQ=255)
                                                                        gATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTgg            >  1:202837/1‑91 (MQ=255)
                                                                        gATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTgg            >  1:370110/1‑91 (MQ=255)
                                                                           acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc         <  2:19818/91‑1 (MQ=255)
                                                                           acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc         <  2:284169/91‑1 (MQ=255)
                                                                           acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc         <  2:423210/91‑1 (MQ=255)
                                                                           acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc         <  1:668995/91‑1 (MQ=255)
                                                                           acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTCCGCCGCCGTCGTAGCTTTGTGGATc         <  2:214722/91‑1 (MQ=255)
                                                                           acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGGTCGCCGCCGTCGTAGCTTTGTGGATc         <  2:663562/91‑1 (MQ=255)
                                                                              acacGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCt      >  1:166534/1‑91 (MQ=255)
                                                                                 cGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTaa   >  2:588146/1‑91 (MQ=255)
                                                                                  gaggTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTaaa  >  2:194270/1‑91 (MQ=255)
                                                                                  gaggTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTaaa  >  2:659513/1‑91 (MQ=255)
                                                                                     |                                                                                       
ACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAA  >  NZ_CP009273/834135‑834307

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ 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
GATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGC  >  NZ_CP009273/834127‑834318
                                                                                             |                                                                                                  
GATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGCTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTG                                                                                             >  SRR3722118.263125/1‑101 (MQ=60)
  TCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCC                                                                                           >  SRR3722118.121881/1‑101 (MQ=60)
  TCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCC                                                                                           >  SRR3722118.219032/1‑101 (MQ=60)
        ACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGT                                                                                     >  SRR3722118.517920/1‑101 (MQ=60)
        ACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTAAGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTCACGGGAGT                                                                                     >  SRR3722118.99437/1‑101 (MQ=60)
                 ACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGA                                                                            >  SRR3722118.118662/1‑101 (MQ=60)
                 ACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGA                                                                            >  SRR3722118.593664/1‑101 (MQ=60)
                 ACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGA                                                                            >  SRR3722118.106689/1‑101 (MQ=60)
                     AAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTT                                                                        >  SRR3722118.636357/1‑101 (MQ=60)
                          ATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCG                                                                   >  SRR3722118.233149/1‑101 (MQ=60)
                              CACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGG                                                               <  SRR3722118.54038/101‑1 (MQ=60)
                                CTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTG                                                             >  SRR3722118.432156/1‑101 (MQ=60)
                                                            GCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGT                                 >  SRR3722118.327289/1‑101 (MQ=60)
                                                               ATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGC                              >  SRR3722118.345414/1‑101 (MQ=60)
                                                                      CTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGG                       >  SRR3722118.377832/1‑101 (MQ=60)
                                                                      CTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGG                       <  SRR3722118.501074/101‑1 (MQ=60)
                                                                      CTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGG                       >  SRR3722118.206891/1‑101 (MQ=60)
                                                                        GCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                     <  SRR3722118.273229/101‑1 (MQ=60)
                                                                        GCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                     <  SRR3722118.558465/101‑1 (MQ=60)
                                                                        GCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                     <  SRR3722118.734522/101‑1 (MQ=60)
                                                                            GAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCT                 >  SRR3722118.169793/1‑101 (MQ=60)
                                                                                   ACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAC          <  SRR3722118.683492/101‑1 (MQ=60)
                                                                                           AGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGC  <  SRR3722118.414240/101‑1 (MQ=60)
                                                                                             |                                                                                                  
GATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGC  >  NZ_CP009273/834127‑834318

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