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
A0 F0 I1 R1 149 26.7 2337422 87.5% 2045244 105.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 582,589 T→G H295Q (CAT→CAG h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479582,5890TG100.0% 103.5 / NA 27H295Q (CAT→CAGh16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/14);  total (13/14)

CGGCGCTGGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTC  >  AM260479/582465‑582712
                                                                                                                            |                                                                                                                           
cggcgCTGGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGccc                                                                                                                    >  1:661347/1‑134 (MQ=255)
                      gggCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAgc                                                                                              >  1:968194/1‑134 (MQ=255)
                                 tCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCggg                                                                                    >  1:568412/1‑133 (MQ=255)
                                     cAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACcg                                                                               <  2:726871/134‑1 (MQ=255)
                                            cTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGaaaa                                                                        <  2:240194/134‑1 (MQ=255)
                                               tCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                                      <  2:568412/133‑1 (MQ=255)
                                                      cgcgcgGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGaaaa                                                                        <  2:619685/124‑1 (MQ=255)
                                                      cgcgcgGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGaaaa                                                                        >  1:619685/1‑124 (MQ=255)
                                                       gcgcgGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCg                                                              <  2:968194/133‑1 (MQ=255)
                                                                cGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCa                                                                                                <  1:850076/90‑1 (MQ=255)
                                                                cGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCa                                                                                                >  2:850076/1‑90 (MQ=255)
                                                                           ccGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATttctt                                                                                                            <  1:297452/67‑1 (MQ=255)
                                                                           ccGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATttctt                                                                                                            >  2:297452/1‑67 (MQ=255)
                                                                                            cAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATttctt                                                                                                            <  1:1137677/50‑1 (MQ=255)
                                                                                            cAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATttctt                                                                                                            >  2:1137677/1‑50 (MQ=255)
                                                                                             aGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt                       <  1:29266/134‑1 (MQ=255)
                                                                                             aGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGAtt                       <  2:887735/134‑1 (MQ=255)
                                                                                                     aCCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGaaaa                                                                        >  1:354804/1‑77 (MQ=255)
                                                                                                     aCCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGaaaa                                                                        <  2:354804/77‑1 (MQ=255)
                                                                                                      ccGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATg              >  1:453360/1‑134 (MQ=255)
                                                                                                       cgcgGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGc             <  1:729168/134‑1 (MQ=255)
                                                                                                         cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                                      <  2:864816/75‑1 (MQ=255)
                                                                                                         cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAAtt                                                                      >  1:864816/1‑75 (MQ=255)
                                                                                                         cggcggGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAg           >  2:454682/1‑134 (MQ=255)
                                                                                                              gggTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACcg      >  1:116627/1‑134 (MQ=255)
                                                                                                               ggTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACcgc     >  1:185439/1‑134 (MQ=255)
                                                                                                                   tGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTc  <  1:454682/133‑1 (MQ=255)
                                                                                                                            |                                                                                                                           
CGGCGCTGGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTC  >  AM260479/582465‑582712

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
GCCCGGACGCCGCGGCGGCGCTGGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAG  >  AM260479/582450‑582727
                                                                                                                                           |                                                                                                                                          
GCCCGGACGCCGCGGCGGCGCTGGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCG                                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2107:13267:34255/1‑150 (MQ=60)
                      gcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCG                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1201:8003:84930/98‑1 (MQ=60)
                      GGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCG                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2209:12882:94504/1‑150 (MQ=60)
                       tttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCG                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2215:19133:65782/80‑1 (MQ=60)
                                 GGCTGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGG                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1216:11767:69543/1‑149 (MQ=60)
                                  cgtgggctcggagatgtgtataagagacagGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGAC                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2116:2989:10504/120‑1 (MQ=60)
                                                      ACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGT                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2102:17410:87782/1‑150 (MQ=60)
                                                                                                     GGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTActgtctcttatacacatctccgagcccacgagacactcgcta                              >  GWNJ‑0478:712:GW2002102894th:2:1204:6039:49265/1‑107 (MQ=60)
                                                                                                      GCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGC                            >  GWNJ‑0478:712:GW2002102894th:2:1210:14810:11200/1‑150 (MQ=60)
                                                                                                        cggagatgtgtataagagacaGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGC                            <  GWNJ‑0478:712:GW2002102894th:2:2204:17663:15184/127‑1 (MQ=60)
                                                                                                         TCCAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGctgtctcttatacacatctccgagcccacgagacactcgctaatc                         >  GWNJ‑0478:712:GW2002102894th:2:2116:8941:77203/1‑105 (MQ=60)
                                                                                                           CAGATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGAC                       <  GWNJ‑0478:712:GW2002102894th:2:1103:18387:88127/150‑1 (MQ=60)
                                                                                                              ATCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGC                    >  GWNJ‑0478:712:GW2002102894th:2:1108:10745:24486/1‑150 (MQ=60)
                                                                                                               TCCTGACCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCG                   >  GWNJ‑0478:712:GW2002102894th:2:1113:10073:3880/1‑150 (MQ=60)
                                                                                                                     CCGCGGCGGGGTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCC             <  GWNJ‑0478:712:GW2002102894th:2:2111:18696:15567/150‑1 (MQ=60)
                                                                                                                               GTTTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCtgtctcttatacacatctc   >  GWNJ‑0478:712:GW2002102894th:2:1112:18270:4676/1‑131 (MQ=60)
                                                                                                                                 TTGACCCGCAGGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAG  <  GWNJ‑0478:712:GW2002102894th:2:1210:19966:18131/149‑1 (MQ=60)
                                                                                                                                           |                                                                                                                                          
GCCCGGACGCCGCGGCGGCGCTGGCCATGGGCGGGCAGGGCGCGGCCATCGCCAACCAGCTGTCGTTCTCGCGCGGGGCCGAGCGCGAGGCCGACCGGGTCGGCTTCCAGATCCTGACCGCGGCGGGGTTTGACCCGCATGGCATGCCCGATTTCTTCATGCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAG  >  AM260479/582450‑582727

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