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
A4 F25 I1 R1 60 73.8 3348830 98.7% 3305295 145.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,586,137 2 bp→TC intergenic (+113/+101) PP_4061 → / ← PP_4063 hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,586,1370CT100.0% 88.0 / NA 26intergenic (+113/+102)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/20);  total (6/20)
*NC_002947_CJ‑RC4,586,1380TC100.0% 88.1 / NA 26intergenic (+114/+101)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/20);  total (6/20)

GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAC  >  NC_002947_CJ‑RC/4585994‑4586210
                                                                                                                                               ||                                                                          
gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGGCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgg                                                                       <  7:74005/149‑1 (MQ=255)
                   gTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAg                                                     <  5:318879/149‑1 (MQ=255)
                       aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAc                                                <  3:47372/149‑1 (MQ=255)
                       aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAc                                                <  8:246342/149‑1 (MQ=255)
                           ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAggg                                                   >  3:293770/1‑143 (MQ=255)
                           ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAggg                                                   <  4:293770/143‑1 (MQ=255)
                           ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAccgtg                                            >  7:87410/1‑145 (MQ=39)
                              aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGcc                                           <  8:39044/148‑1 (MQ=39)
                                ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGc                                        <  1:86597/149‑1 (MQ=39)
                                    gACCACGGCACCGGCGTGAAGAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                     <  5:384786/148‑1 (MQ=38)
                                       catggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                     <  5:147707/142‑1 (MQ=38)
                                       catggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                     >  6:147707/4‑145 (MQ=38)
                                          ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                              <  8:65068/149‑1 (MQ=39)
                                                 gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGc                       <  8:87410/149‑1 (MQ=39)
                                                        cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                <  3:311301/149‑1 (MQ=39)
                                                           cTGATCTGCGCTGCTCCATTGAATGGTTGAGGTTGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCACCTGGTAg              <  5:405115/148‑1 (MQ=38)
                                                           cTGATCTGCGCTGCTCAATTGAAAGGTTGAGGATGACAGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGTCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCATGAAAGCGGCGCAGCCCCTGGTAg              >  6:405115/1‑148 (MQ=37)
                                                                   cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGt      <  8:84622/148‑1 (MQ=38)
                                                                   cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCggg       <  2:20565/147‑1 (MQ=38)
                                                                    gctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGttt    <  6:358852/149‑1 (MQ=38)
                                                                     ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTa   <  5:255635/149‑1 (MQ=38)
                                                                        ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc  <  7:341388/147‑1 (MQ=38)
                                                                                aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgg                                                                       <  7:113106/69‑1 (MQ=255)
                                                                                aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgg                                                                       >  8:113106/1‑69 (MQ=255)
                                                                                 aaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc          >  4:81661/1‑130 (MQ=38)
                                                                                 aaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc          <  3:81661/130‑1 (MQ=38)
                                                                                                                                               ||                                                                          
GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAC  >  NC_002947_CJ‑RC/4585994‑4586210

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACC  >  NC_002947_CJ‑RC/4585987‑4586211
                                                                                                                                                      |                                                                            
GCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTT                                                                              <  NB501016:282:H2V7GAFX2:1:11104:22685:1465/151‑1 (MQ=60)
                         GGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGG                                                     <  NB501016:282:H2V7GAFX2:3:21507:19237:15328/151‑1 (MQ=60)
                             AAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC                                                <  NB501016:282:H2V7GAFX2:2:11109:14622:6458/151‑1 (MQ=60)
                              AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG                                                                        <  NB501016:282:H2V7GAFX2:4:21503:2560:5674/126‑1 (MQ=60)
                                 CTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGC                                             >  NB501016:282:H2V7GAFX2:4:11503:15943:4441/1‑151 (MQ=60)
                                 CTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGA                                                   >  NB501016:282:H2V7GAFX2:2:21205:23150:16108/1‑145 (MQ=60)
                                      GGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCC                                        <  NB501016:282:H2V7GAFX2:1:11202:20434:19048/151‑1 (MQ=60)
                                          AGACCACGGCACCGGCGTGAAGAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTG                                     <  NB501016:282:H2V7GAFX2:3:21606:4299:14106/150‑1 (MQ=60)
                                             CCATGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTG                                     <  NB501016:282:H2V7GAFX2:3:11602:15586:2460/147‑1 (MQ=60)
                                                              ACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTA                <  NB501016:282:H2V7GAFX2:2:21208:11277:16092/151‑1 (MQ=60)
                                                                 GCTGATCTGCGCTGCTCCATTGAATGGTTGAGGTTGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCACCTGGTAGG              <  NB501016:282:H2V7GAFX2:3:21610:7808:2389/150‑1 (MQ=60)
                                                                           GCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAC   <  NB501016:282:H2V7GAFX2:3:21408:22698:11155/151‑1 (MQ=60)
                                                                              GCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACC  <  NB501016:282:H2V7GAFX2:4:21509:17513:9783/149‑1 (MQ=60)
                                                                                       AAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCG          <  NB501016:282:H2V7GAFX2:2:11203:8970:5253/132‑1 (MQ=60)
                                                                                                                                                      |                                                                            
GCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACC  >  NC_002947_CJ‑RC/4585987‑4586211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: