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 3,728,024 A→C S315S (TCA→TCC xylH → xylose ABC transporter permease XylH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,728,0240AC100.0% 83.4 / NA 24S315S (TCA→TCCxylHxylose ABC transporter permease XylH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/8);  total (16/8)

ACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCG  >  NZ_CP009273/3727941‑3728109
                                                                                   |                                                                                     
aCCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCt                                                                                 >  2:685613/1‑90 (MQ=255)
    aaCTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCgg                                                                             <  1:583774/90‑1 (MQ=255)
         gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAat                                                                        >  1:336426/1‑90 (MQ=255)
         gCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAat                                                                        >  2:713869/1‑90 (MQ=255)
                        aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGAc                                                         >  2:558141/1‑90 (MQ=255)
                          cGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGc                                                       <  2:672047/90‑1 (MQ=255)
                             aTTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAAt                                                    >  2:644741/1‑90 (MQ=255)
                                  tGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTgcag                                               <  2:428626/90‑1 (MQ=255)
                                           tCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGa                                      >  2:404037/1‑90 (MQ=255)
                                           tCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGa                                      >  2:247668/1‑90 (MQ=255)
                                            cGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGAt                                     <  2:273787/90‑1 (MQ=255)
                                                       tCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCa                          <  1:240706/90‑1 (MQ=255)
                                                           tAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCct                      >  2:658530/1‑90 (MQ=255)
                                                           tAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCct                      >  1:501208/1‑90 (MQ=255)
                                                               tctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggct                  <  1:634828/90‑1 (MQ=255)
                                                                ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggctg                 >  1:6709/1‑90 (MQ=255)
                                                                ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggctg                 >  2:510910/1‑90 (MQ=255)
                                                                ctcGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCctggctg                 >  2:757572/1‑90 (MQ=255)
                                                                  cGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCtggctggc               >  2:44166/1‑90 (MQ=255)
                                                                       tGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGgtgt          >  1:662630/1‑90 (MQ=255)
                                                                         gcgcTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGGGTgg        >  1:431425/1‑90 (MQ=255)
                                                                              ggTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAc                            <  2:505490/65‑1 (MQ=255)
                                                                              ggTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCAc                            >  1:505490/1‑65 (MQ=255)
                                                                               gTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCg  <  2:39396/90‑1 (MQ=255)
                                                                                   |                                                                                     
ACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCG  >  NZ_CP009273/3727941‑3728109

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
TAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGC  >  NZ_CP009273/3727928‑3728113
                                                                                                |                                                                                         
TAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCT                                                                                        <  SRR3722073.624043/100‑1 (MQ=60)
            CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAAT                                                                            >  SRR3722073.343565/1‑100 (MQ=60)
                 AACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGC                                                                       <  SRR3722073.595957/100‑1 (MQ=60)
                                                              GATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCT                          >  SRR3722073.511810/1‑100 (MQ=60)
                                                                   ATCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTG                     >  SRR3722073.6881/1‑100 (MQ=60)
                                                                    TCCTTAGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGG                    <  SRR3722073.245758/100‑1 (MQ=60)
                                                                          GTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGT              >  SRR3722073.676211/1‑100 (MQ=60)
                                                                            TCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGG            <  SRR3722073.647946/100‑1 (MQ=60)
                                                                            TCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGGGTGG            >  SRR3722073.440654/1‑100 (MQ=60)
                                                                                 ACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTctgtctcttatacac       >  SRR3722073.516164/1‑85 (MQ=60)
                                                                                      GCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGGGTGGGAAGCGTTGC  >  SRR3722073.187931/1‑100 (MQ=60)
                                                                                                |                                                                                         
TAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGC  >  NZ_CP009273/3727928‑3728113

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: