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 2,616,505 C→G A45G (GCC→GGC)  h16_A2404 → predicted esterase of the alpha‑beta hydrolase superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,616,5050CG100.0% 71.8 / NA 21A45G (GCC→GGC) h16_A2404predicted esterase of the alpha‑beta hydrolase superfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (10/11);  total (10/11)

GGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616404‑2616630
                                                                                                     |                                                                                                                             
ggCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCg                                                                                               >  2:64253/1‑134 (MQ=255)
 gccgccCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCgg                                                                                              >  2:881237/1‑134 (MQ=255)
   cgccCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    <  1:188604/110‑1 (MQ=255)
   cgccCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    >  2:188604/1‑110 (MQ=255)
         gctgGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    >  2:659177/1‑104 (MQ=255)
         gctgGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    <  1:659177/104‑1 (MQ=255)
             gCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGgc                                                                                  <  2:101058/134‑1 (MQ=255)
              cGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGgcg                                                                                 <  1:881237/134‑1 (MQ=255)
                     ctCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAg                                                                          <  2:277280/134‑1 (MQ=255)
                        gTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgcg                                                                        >  1:836279/1‑133 (MQ=255)
                        gTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgcg                                                                        >  2:275139/1‑133 (MQ=255)
                            ggCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAg                                                                    <  1:275139/133‑1 (MQ=255)
                                 gcgTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAg                                                                    >  2:354016/1‑128 (MQ=255)
                                 gcgTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAg                                                                    <  1:354016/128‑1 (MQ=255)
                                      ccGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGc                                                          <  2:836279/133‑1 (MQ=255)
                                                                              cGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcg                                                                                                                <  1:872150/39‑1 (MQ=25)
                                                                              cGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcg                                                                                                                >  2:872150/1‑39 (MQ=25)
                                                                               gTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  1:974020/39‑1 (MQ=25)
                                                                               gTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               >  2:974020/1‑39 (MQ=25)
                                                                                      tCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTAc         <  1:64253/134‑1 (MQ=255)
                                                                                              gcgcTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCg  >  2:576091/1‑133 (MQ=255)
                                                                                                     |                                                                                                                             
GGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616404‑2616630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 35 ≤ 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
TCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGGACTGAAT  >  AM260479/2616383‑2616638
                                                                                                                          |                                                                                                                                     
gtgtataagagacagCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCC                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2111:14346:42233/134‑1 (MQ=60)
aagagacagCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCC                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1112:19075:12553/140‑1 (MQ=60)
   cggttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATAT                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2206:7486:59611/69‑1 (MQ=60)
     tcttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATC                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2213:14747:67711/70‑1 (MQ=60)
                              GCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGC                                                                               >  GWNJ‑0478:712:GW2002102894th:2:2205:13734:32352/1‑149 (MQ=60)
                                  GCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGC                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2207:17067:17324/150‑1 (MQ=60)
                                                CGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACC                                                             <  GWNJ‑0478:712:GW2002102894th:2:1203:19535:44626/149‑1 (MQ=60)
                                                CGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACC                                                             <  GWNJ‑0478:712:GW2002102894th:2:1206:19543:64027/149‑1 (MQ=60)
                                                    agacggcatacgagatgtagctccgtctcgtgggctaggagatgtgtataagagacCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTC                                                        <  GWNJ‑0478:712:GW2002102894th:2:1204:20250:63526/94‑1 (MQ=60)
                                                                                                          ATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGGACTGAAT  <  GWNJ‑0478:712:GW2002102894th:2:1104:5343:44907/150‑1 (MQ=60)
                                                                                                                          |                                                                                                                                     
TCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGGACTGAAT  >  AM260479/2616383‑2616638

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