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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,927,010 T→C K2E (AAA→GAA)  kdgA ← bifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,927,0100TC100.0% 116.1 / NA 34K2E (AAA→GAA) kdgAbifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/19);  total (15/19)

ACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACC  >  NZ_CP009273/1926925‑1927085
                                                                                     |                                                                           
aCCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatc                                                                         <  2:36731/90‑1 (MQ=255)
   gcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCaga                                                                      >  2:728814/1‑90 (MQ=255)
       gtTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGttt                                                                  >  2:505284/1‑90 (MQ=255)
        ttCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtttt                                                                 >  1:409245/1‑90 (MQ=255)
             gTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTtctctc                                                            >  1:754346/1‑90 (MQ=255)
                tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                         <  2:564110/90‑1 (MQ=255)
                tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCt                                                        <  2:334108/90‑1 (MQ=255)
                                   tACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcgtc                                      >  2:121866/1‑90 (MQ=255)
                                         cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                   >  2:432496/1‑55 (MQ=255)
                                         cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                   <  1:432496/55‑1 (MQ=255)
                                         cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                <  2:863983/90‑1 (MQ=255)
                                         cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                <  2:730224/90‑1 (MQ=255)
                                         cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                <  2:186357/90‑1 (MQ=255)
                                          gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                               >  2:1149984/1‑90 (MQ=255)
                                           ggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                              <  1:644383/90‑1 (MQ=255)
                                            gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGt                             >  1:919655/1‑90 (MQ=255)
                                                gtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATa                         <  2:343684/90‑1 (MQ=255)
                                                gtggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATa                         <  2:356254/90‑1 (MQ=255)
                                                  ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                       <  1:543489/90‑1 (MQ=255)
                                                       ggATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTg                  <  1:1138166/90‑1 (MQ=255)
                                                          ttgattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTgcgc               <  2:1070783/90‑1 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:551014/90‑1 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:180758/90‑1 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:754346/90‑1 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  2:828861/90‑1 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  1:943933/90‑1 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             >  1:481266/1‑90 (MQ=255)
                                                            gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc             <  1:27300/90‑1 (MQ=255)
                                                                 tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTTGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  1:259287/1‑90 (MQ=255)
                                                                 tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  1:1028933/1‑90 (MQ=255)
                                                                 tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  1:406233/1‑90 (MQ=255)
                                                                 tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        <  1:270821/90‑1 (MQ=255)
                                                                 tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc        >  2:97304/1‑90 (MQ=255)
                                                                  gCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCg       >  2:51665/1‑90 (MQ=255)
                                                                       tGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAcc  >  2:569620/1‑90 (MQ=255)
                                                                                     |                                                                           
ACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACC  >  NZ_CP009273/1926925‑1927085

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGTTTTT  >  NZ_CP009273/1926923‑1927096
                                                                                       |                                                                                      
GCACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTT                                                                            >  SRR3722089.414200/1‑100 (MQ=60)
    tctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCT                                                                        <  SRR3722089.496856/52‑1 (MQ=60)
     GCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTC                                                                       >  SRR3722089.764383/1‑100 (MQ=60)
        agcgtcagatgtgtataagagacagCCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCC                                                                    <  SRR3722089.437695/75‑1 (MQ=60)
                               AACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTctgtctcttatacacatctgacgctgccgacgac                                             >  SRR3722089.529448/1‑66 (MQ=60)
                                    GTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGT                                        >  SRR3722089.933073/1‑100 (MQ=60)
                                             GGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGT                               <  SRR3722089.652623/100‑1 (MQ=60)
                                                    GGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCC                        >  SRR3722089.486983/1‑100 (MQ=60)
                                                    GGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCC                        <  SRR3722089.549816/100‑1 (MQ=60)
                                                         GGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTTGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTC                   >  SRR3722089.262614/1‑100 (MQ=60)
                                                         GGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTC                   <  SRR3722089.1155482/100‑1 (MQ=60)
                                                         GGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTC                   >  SRR3722089.1044295/1‑100 (MQ=60)
                                                         GGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTC                   >  SRR3722089.411155/1‑100 (MQ=60)
                                                              GATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAC              <  SRR3722089.957767/100‑1 (MQ=60)
                                                              GATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAC              <  SRR3722089.27630/100‑1 (MQ=60)
                                                                   TGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAC         <  SRR3722089.274282/100‑1 (MQ=60)
                                                                          GTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGTTTTT  >  SRR3722089.451517/1‑100 (MQ=60)
                                                                                       |                                                                                      
GCACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGTTTTT  >  NZ_CP009273/1926923‑1927096

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