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 152 24.9 2083062 91.6% 1908084 109.2

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% 68.7 / NA 19A45G (GCC→GGC) h16_A2404predicted esterase of the alpha‑beta hydrolase superfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (9/10);  total (9/10)

CGTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616381‑2616630
                                                                                                                            |                                                                                                                             
cgTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgc                                                                                                                      >  2:554644/1‑134 (MQ=255)
  tCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    <  1:554644/134‑1 (MQ=255)
  tCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGACAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    <  1:26445/134‑1 (MQ=255)
        ccTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  1:769771/133‑1 (MQ=255)
                            cccTGCTGGCGGGCTGCTAGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGt                                                                                           >  2:223057/1‑133 (MQ=255)
                                                                   cgccgcTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTg                                                   <  1:1001141/134‑1 (MQ=255)
                                                                        cTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCAc                                              >  2:236891/1‑134 (MQ=255)
                                                                               ccaccCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCAcc                                       <  1:223057/134‑1 (MQ=255)
                                                                                    ccgcccgcGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCa                                         >  1:415285/1‑127 (MQ=255)
                                                                                    ccgcccgcGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCa                                         <  2:415285/127‑1 (MQ=255)
                                                                                      gcccgcGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATc                                                               >  1:858528/1‑103 (MQ=255)
                                                                                      gcccgcGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATc                                                               <  2:858528/103‑1 (MQ=255)
                                                                                       cccgcGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  1:832635/54‑1 (MQ=255)
                                                                                       cccgcGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               >  2:832635/1‑54 (MQ=255)
                                                                                           cgccgcgGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGgcg                                                                                 <  1:990534/80‑1 (MQ=255)
                                                                                           cgccgcgGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGgcg                                                                                 >  2:990534/1‑80 (MQ=255)
                                                                                                                gCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAg                                                                    >  1:724030/1‑72 (MQ=255)
                                                                                                                gCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAg                                                                    <  2:724030/72‑1 (MQ=255)
                                                                                                                     gcgcTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCg  >  1:58458/1‑133 (MQ=255)
                                                                                                                            |                                                                                                                             
CGTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616381‑2616630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGG  >  AM260479/2616382‑2616631
                                                                                                                           |                                                                                                                              
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2102:15521:26384/150‑1 (MQ=60)
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGACAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1103:1778:87314/150‑1 (MQ=60)
     caagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATC                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2204:2293:80069/85‑1 (MQ=60)
      TCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATC                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2116:15886:50940/149‑1 (MQ=60)
                                   agctccgtctcgtgggctcggagatgtgtataagagacaGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGC                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2214:2206:29545/111‑1 (MQ=60)
                                                                 CCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGC                                     <  GWNJ‑0478:712:GW2002102894th:2:2214:13381:86430/150‑1 (MQ=60)
                                                                    CCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCC                                  >  GWNJ‑0478:712:GW2002102894th:2:1210:3570:14182/1‑150 (MQ=60)
                                                                      GCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCctgtctcttatacacat                                >  GWNJ‑0478:712:GW2002102894th:2:2209:3517:48361/1‑133 (MQ=60)
                                                                             ACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTC                         <  GWNJ‑0478:712:GW2002102894th:2:1115:5190:98917/150‑1 (MQ=60)
                                                                                                GGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCctgtctcttatacacatctccgagcccacgagacggagctacatctcg      >  GWNJ‑0478:712:GW2002102894th:2:2112:19697:96175/1‑102 (MQ=60)
                                                                                                     GTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGG  >  GWNJ‑0478:712:GW2002102894th:2:1104:14134:7083/1‑149 (MQ=60)
                                                                                                                           |                                                                                                                              
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGG  >  AM260479/2616382‑2616631

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