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 2,174,381 T→G C270G (TGC→GGC)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,174,3810TG100.0% 117.2 / NA 35C270G (TGC→GGC) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (16/19);  total (16/19)

AGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTT  >  NZ_CP009273/2174295‑2174467
                                                                                      |                                                                                      
aGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCgcggcg                                                                                     >  2:411913/1‑90 (MQ=255)
  gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc                                                                                   >  2:7553/1‑90 (MQ=255)
       tGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCt                                                                              <  1:1067975/90‑1 (MQ=255)
                   cAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCt                                                                  >  1:901316/1‑90 (MQ=255)
                    aGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCtt                                                                 >  2:1153262/1‑90 (MQ=255)
                        cTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcg                                                             >  1:313625/1‑90 (MQ=255)
                        cTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcg                                                             >  2:517210/1‑90 (MQ=255)
                         tATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcgg                                                            >  2:320718/1‑90 (MQ=255)
                          aTTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcggc                                                           >  1:8522/1‑90 (MQ=255)
                          aTTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcggc                                                           >  1:79656/1‑90 (MQ=255)
                              aGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGaca                                                       <  2:313625/90‑1 (MQ=255)
                               gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                      >  2:1158851/1‑90 (MQ=255)
                                  cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                                   >  1:39790/1‑90 (MQ=255)
                                   gttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGac                                                  <  1:277027/90‑1 (MQ=255)
                                    ttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGaca                                                 <  2:975378/90‑1 (MQ=255)
                                       cgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCa                                              >  2:715646/1‑90 (MQ=255)
                                         tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                            <  2:804553/90‑1 (MQ=255)
                                         tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                            <  2:293975/90‑1 (MQ=255)
                                         tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                            <  1:517210/90‑1 (MQ=255)
                                         tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                            <  2:693486/90‑1 (MQ=255)
                                         tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                            >  2:416540/1‑90 (MQ=255)
                                                      tGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGa                               <  1:1006179/90‑1 (MQ=255)
                                                       ggTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGAc                              >  1:613646/1‑90 (MQ=255)
                                                         ttGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACgg                            <  1:1206767/90‑1 (MQ=255)
                                                            aaCTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCaa                         >  2:401262/1‑90 (MQ=255)
                                                                 gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTTCTATGGCGACGGCAATTTgc                    >  1:1118184/1‑90 (MQ=255)
                                                                   acaAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgcgg                  <  1:653186/90‑1 (MQ=255)
                                                                          gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTg           >  2:716380/1‑90 (MQ=255)
                                                                             ccGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt        <  1:320718/90‑1 (MQ=255)
                                                                             ccGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt        <  2:1160814/90‑1 (MQ=255)
                                                                             ccGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt        <  2:8522/90‑1 (MQ=255)
                                                                              cGACTCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg       <  2:950621/90‑1 (MQ=255)
                                                                              cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg       <  2:39790/90‑1 (MQ=255)
                                                                              cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg       <  1:416540/90‑1 (MQ=255)
                                                                              cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg       <  1:1071525/90‑1 (MQ=255)
                                                                              cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg       <  1:187470/90‑1 (MQ=255)
                                                                                   cgcgGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGtt  <  1:996648/90‑1 (MQ=255)
                                                                                      |                                                                                      
AGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTT  >  NZ_CP009273/2174295‑2174467

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
TACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCG  >  NZ_CP009273/2174291‑2174477
                                                                                          |                                                                                                
tcggcagcgtcagatgtgtataagagacagATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCC                                                                                         <  SRR3722089.659841/70‑1 (MQ=60)
           TGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAAC                                                                              <  SRR3722089.1084010/100‑1 (MQ=60)
             GCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCT                                                                            >  SRR3722089.914390/1‑100 (MQ=60)
                  GGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCG                                                                       >  SRR3722089.317601/1‑100 (MQ=60)
                    TACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGC                                                                     >  SRR3722089.8620/1‑100 (MQ=60)
                    TACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGC                                                                     >  SRR3722089.80629/1‑100 (MQ=60)
                            CTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGA                                                             >  SRR3722089.40286/1‑100 (MQ=60)
                                       GTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTG                                                  <  SRR3722089.280568/100‑1 (MQ=60)
                                             TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                            <  SRR3722089.523283/100‑1 (MQ=60)
                                                 CCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGAC                                        >  SRR3722089.621307/1‑100 (MQ=60)
                                                          TGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTG                               <  SRR3722089.1021159/100‑1 (MQ=60)
                                                           GGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGctgtctcttatacacatctgacgctgccgacgaccatcttagtgtagatctcggtgg                              >  SRR3722089.1089562/1‑43 (MQ=60)
                                                           GGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTTCTATGGCGACGGCAATTTGC                              >  SRR3722089.1135138/1‑100 (MQ=60)
                                                             TTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGG                            <  SRR3722089.1225114/100‑1 (MQ=60)
                                                                       ACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAT                  <  SRR3722089.661578/100‑1 (MQ=60)
                                                                              GACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTA           >  SRR3722089.356737/1‑100 (MQ=60)
                                                                                 CCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACG        <  SRR3722089.324763/100‑1 (MQ=60)
                                                                                  CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC       <  SRR3722089.421577/100‑1 (MQ=60)
                                                                                  CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC       <  SRR3722089.189813/100‑1 (MQ=60)
                                                                                  CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC       <  SRR3722089.1087617/100‑1 (MQ=60)
                                                                                       CGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCG  <  SRR3722089.1011453/100‑1 (MQ=60)
                                                                                          |                                                                                                
TACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGCTATCG  >  NZ_CP009273/2174291‑2174477

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