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 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,602,688 (G)6→5 intergenic (‑85/+41) h16_A3328 ← / ← h16_A3329 conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,602,6830G.100.0% 68.5 / NA 16intergenic (‑80/+46)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/7);  total (9/7)

GGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCC  >  AM260479/3602559‑3602798
                                                                                                                             |                                                                                                                   
gggCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggcgat                                                                                                             >  1:337231/1‑130 (MQ=255)
    aGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCg                                                                                                        >  1:905766/1‑134 (MQ=255)
     ggTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGt                                                                                                       >  1:253689/1‑134 (MQ=255)
         gCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCgg                                                                                                   <  2:337231/134‑1 (MQ=255)
                   aGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgtgt                                                                                         >  1:240048/1‑134 (MQ=255)
                      gcgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTgg                                                                                       >  1:592073/1‑133 (MQ=255)
                         aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGc                                                                                   <  2:592073/134‑1 (MQ=255)
                           cAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCggg                                                                                                                   <  2:252946/101‑1 (MQ=255)
                           cAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCggg                                                                                                                   >  1:252946/1‑101 (MQ=255)
                                                            cAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcgcgagcg                                                <  2:240048/134‑1 (MQ=255)
                                                                                  gaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                          <  2:253689/132‑1 (MQ=255)
                                                                                   aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                                                     >  1:555103/2‑74 (MQ=255)
                                                                                   aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                                                     <  2:555103/73‑1 (MQ=255)
                                                                                   aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCtctc                           >  1:419925/2‑132 (MQ=255)
                                                                                   aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCtctc                           >  2:588226/2‑132 (MQ=255)
                                                                                        aGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCtt                     >  1:520249/1‑133 (MQ=255)
                                                                                                          tGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCgccgcc  <  2:520249/134‑1 (MQ=255)
                                                                                                             gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                          >  2:923793/1‑107 (MQ=255)
                                                                                                             gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                          <  1:923793/107‑1 (MQ=255)
                                                                                                                             |                                                                                                                   
GGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCC  >  AM260479/3602559‑3602798

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGC  >  AM260479/3602544‑3602822
                                                                                                                                                 |                                                                                                                                      
CACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATG                                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:1207:1410:59356/1‑150 (MQ=60)
 atgttagtttttttttttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaggaATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGG                                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2207:12243:56076/52‑1 (MQ=60)
 ACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGG                                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1207:17750:52914/150‑1 (MQ=60)
    ACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGT                                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:2208:17342:70589/1‑150 (MQ=60)
     CCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTG                                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1116:8902:89012/1‑150 (MQ=60)
         gatcgcgctgcacttcagcgacctgcccaccgcgatggtggtgtCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGT                                                                                                                          >  GWNJ‑0478:712:GW2002102894th:2:2106:1837:95842/45‑150 (MQ=60)
                   agcgtgctgctGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTG                                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1116:14881:6051/12‑150 (MQ=60)
                      TCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGC                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:2107:13692:17667/1‑149 (MQ=60)
                           GAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGctgtctcttatacacatct                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1116:4359:84194/1‑131 (MQ=60)
                                                                                   GGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCctgtctcttatacacatctccgagcccacgagacggagctacatct                                                >  GWNJ‑0478:712:GW2002102894th:2:2102:20561:61992/1‑104 (MQ=60)
                                                                                   GGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCG                                                  >  GWNJ‑0478:712:GW2002102894th:2:1210:13834:55545/1‑148 (MQ=60)
                                                                                        ATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTC                                            >  GWNJ‑0478:712:GW2002102894th:2:2101:18626:11579/1‑149 (MQ=60)
                                                                                                gtataagagacAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAAC                                   <  GWNJ‑0478:712:GW2002102894th:2:2211:17806:78442/139‑1 (MQ=60)
                                                                                                                                ATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGG   >  GWNJ‑0478:712:GW2002102894th:2:1109:13461:47472/1‑150 (MQ=60)
                                                                                                                                  TGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGC  >  GWNJ‑0478:712:GW2002102894th:2:1105:9517:87525/1‑149 (MQ=60)
                                                                                                                                                 |                                                                                                                                      
CACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGC  >  AM260479/3602544‑3602822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: