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,329 Δ1 bp coding (205/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,3270G.100.0% 109.6 / NA 26coding (203/1221 nt)h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (13/13);  total (13/13)

CTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCC  >  AM260479/2907210‑2907442
                                                                                                                     |                                                                                                                   
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                                                                                               >  2:925306/1‑106 (MQ=255)
                               gctgctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTg                                                                                               <  1:925306/106‑1 (MQ=255)
                                       ccGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGcc                                                           <  2:405601/134‑1 (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                                                                     <  1:11432/83‑1 (MQ=255)
                                                                                agctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTcgc                                                                     >  2:11432/1‑83 (MQ=255)
                                                                                     cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                             <  2:985803/54‑1 (MQ=19)
                                                                                     cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                             >  1:985803/1‑54 (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)
                                                                                                 tcaagaaCAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGccc  >  2:655661/6‑134 (MQ=255)
                                                                                                                     |                                                                                                                   
CTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCC  >  AM260479/2907210‑2907442

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
CAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC  >  AM260479/2907197‑2907446
                                                                                                                                    |                                                                                                                     
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)
                                                       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)
                                                                                                                                    |                                                                                                                     
CAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC  >  AM260479/2907197‑2907446

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