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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,907,305 2 bp→AT coding (181‑182/1221 nt) h16_A2693 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,907,3050GA100.0% 107.5 / NA 29D61I (GAC→ATC) h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/16);  total (13/16)
*AM2604792,907,3060AT100.0% 104.6 / NA 29D61I (GAC→ATC) h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/16);  total (13/16)

CGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGG  >  AM260479/2907183‑2907431
                                                                                                                          ||                                                                                                                             
cgccgcTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGaacagc                                                                                                                     <  1:134310/134‑5 (MQ=255)
cgccgcTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGaacagc                                                                                                                     <  1:890879/134‑5 (MQ=255)
                           cTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTAc                                                                                       <  1:898926/134‑1 (MQ=255)
                             gccgccGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCt                                                                                     >  2:924090/1‑134 (MQ=255)
                             gccgccGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCt                                                                                     >  1:757676/1‑134 (MQ=255)
                              ccgccgGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTg                                                                                    <  1:191909/134‑1 (MQ=255)
                                              cACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACg                                                                    <  1:169776/134‑1 (MQ=255)
                                                          gctgctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTg                                                                                    <  1:925306/106‑1 (MQ=255)
                                                          gctgctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTg                                                                                    >  2:925306/1‑106 (MQ=255)
                                                                  ccGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGcc                                                <  2:405601/134‑1 (MQ=255)
                                                                       atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGc                                                                                                            <  2:105105/71‑1 (MQ=255)
                                                                       atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGc                                                                                                            >  1:105105/1‑71 (MQ=255)
                                                                                    acacCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCgg                              >  1:945560/1‑134 (MQ=255)
                                                                                            gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcggc                      >  2:796408/1‑134 (MQ=255)
                                                                                              gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               >  1:178533/1‑91 (MQ=255)
                                                                                              gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               <  2:178533/91‑1 (MQ=255)
                                                                                                           agctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTcgc                                                          >  2:11432/1‑83 (MQ=255)
                                                                                                           agctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTcgc                                                          <  1:11432/83‑1 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                  >  1:985803/1‑54 (MQ=19)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                  <  2:985803/54‑1 (MQ=19)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               >  2:25169/1‑73 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               <  1:25169/73‑1 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCt                            <  1:913945/108‑1 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCt                            >  2:913945/1‑108 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc               >  2:257430/1‑121 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc               <  1:257430/121‑1 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt              <  2:310689/122‑1 (MQ=255)
                                                                                                                cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt              >  1:310689/1‑122 (MQ=255)
                                                                                                                   ggAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAgg  <  1:924090/131‑1 (MQ=255)
                                                                                                                          ||                                                                                                                             
CGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGG  >  AM260479/2907183‑2907431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
GCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC  >  AM260479/2907182‑2907446
                                                                                                                           |                                                                                                                                             
GCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGcgac                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1109:17597:86071/150‑5 (MQ=60)
GCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGcgac                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2208:7946:11351/150‑5 (MQ=60)
               ggttccgtgGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTG                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2116:18998:45646/10‑150 (MQ=60)
                   agatgtagctccgtctcgtgggctcggagatgtgtataagagacagCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGC                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1108:7394:7974/104‑1 (MQ=60)
                           GCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCC                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2208:16312:60580/150‑1 (MQ=60)
                              GCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAAC                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1112:13636:32593/150‑1 (MQ=60)
                               gtataagagacagGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAAC                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2210:9941:18655/136‑1 (MQ=60)
                                              CCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGC                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1110:14613:99641/150‑1 (MQ=60)
                                                   gagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCC                                                               <  GWNJ‑0478:712:GW2002102894th:2:1103:18786:74729/103‑1 (MQ=60)
                                                         ccgtctcgtgggctcggagatgtgtataagagacAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGC                                                          <  GWNJ‑0478:712:GW2002102894th:2:1101:12063:82268/115‑1 (MQ=60)
                                                         CTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGcgactacgaggcctgtctcttatacacatctccgagcccacgagacggagctacatctcgt                                                           >  GWNJ‑0478:712:GW2002102894th:2:1108:17017:3946/1‑89 (MQ=60)
                                                                      GCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGC                                            >  GWNJ‑0478:712:GW2002102894th:2:2212:18185:38482/1‑150 (MQ=60)
                                                                                CAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCctgtctcttatacacatctccgagcccac                                  >  GWNJ‑0478:712:GW2002102894th:2:1113:6072:52496/1‑121 (MQ=60)
                                                                                       ataagagacagGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTC                            <  GWNJ‑0478:712:GW2002102894th:2:2211:19492:51654/138‑1 (MQ=60)
                                                                                                  GTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttg                >  GWNJ‑0478:712:GW2002102894th:2:2214:8455:65945/1‑84 (MQ=60)
                                                                                                  GTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGC                >  GWNJ‑0478:712:GW2002102894th:2:1205:20663:95239/1‑150 (MQ=60)
                                                                                                    ACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCC               <  GWNJ‑0478:712:GW2002102894th:2:1201:9304:35418/149‑1 (MQ=60)
                                                                                                                   GGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC  <  GWNJ‑0478:712:GW2002102894th:2:2210:4280:11554/147‑1 (MQ=60)
                                                                                                                           |                                                                                                                                             
GCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC  >  AM260479/2907182‑2907446

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