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 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,5050CG95.0% 69.1 / ‑3.0 20coding (134/921 nt)h16_A2404predicted esterase of the alpha‑beta hydrolase superfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (7/12);  minor base . (0/1);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CGTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616381‑2616630
                                                                                                                            |                                                                                                                             
cgTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgc                                                                                                                      >  2:561970/1‑134 (MQ=255)
  tCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGgcggcgcgg                                                                                                                    <  1:611430/134‑1 (MQ=255)
        ccTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  1:79907/133‑1 (MQ=255)
        ccTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  2:24945/133‑1 (MQ=255)
         cTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgcgg                                                                                                             <  1:561970/134‑1 (MQ=255)
            ggTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCt                                                                                                           <  2:548133/133‑1 (MQ=255)
                    cgcgGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCata                                                                                                  >  2:795803/1‑134 (MQ=255)
                                gctgGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCCCCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               >  1:593260/1‑109 (MQ=255)
                                gctgGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgc                                                                                                               <  2:593260/109‑1 (MQ=255)
                                            ctCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGa                                                                            <  1:795803/132‑1 (MQ=255)
                                            ctCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         <  2:5979/134‑1 (MQ=255)
                                             tCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         >  1:307433/1‑134 (MQ=255)
                                               gTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         <  2:307433/132‑1 (MQ=255)
                                                      ccGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgcg                                                                                                              >  1:398097/1‑88 (MQ=255)
                                                      ccGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGgcgcggcgcgcg                                                                                                              <  2:398097/88‑1 (MQ=255)
                                                             ccGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGc                                                          <  1:836543/133‑1 (MQ=255)
                                                             ccGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGc                                                          <  1:8072/133‑1 (MQ=255)
                                                                                                          aGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         >  2:680877/1‑73 (MQ=255)
                                                                                                          aGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAgc                                                                         <  1:680877/73‑1 (MQ=255)
                                                                                                                    ggCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCg  >  1:953858/1‑134 (MQ=255)
                                                                                                                            |                                                                                                                             
CGTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCG  >  AM260479/2616381‑2616630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGG  >  AM260479/2616382‑2616631
                                                                                                                           |                                                                                                                              
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGgacgcgggccgaacg                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2102:12779:36365/150‑16 (MQ=60)
      TCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATC                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1107:5827:7127/149‑1 (MQ=60)
       CCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGG                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1216:4254:40147/150‑1 (MQ=60)
          ttttctgagagtggtttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCG                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1113:4247:23039/58‑1 (MQ=60)
                CTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCCCCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCctgtctcttat                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:2103:6778:20482/1‑139 (MQ=60)
                             CTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCG                                                                         >  GWNJ‑0478:712:GW2002102894th:2:1203:20059:49645/1‑150 (MQ=60)
                                      GGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGctgtctcttatacacatctccgagcccacg                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1208:17008:8112/1‑118 (MQ=60)
                                          GCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCA                                                              <  GWNJ‑0478:712:GW2002102894th:2:2112:13609:54056/148‑1 (MQ=60)
                                           gagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGC                                                           <  GWNJ‑0478:712:GW2002102894th:2:2106:16950:48335/104‑1 (MQ=60)
                                                           CCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGG                                            <  GWNJ‑0478:712:GW2002102894th:2:1101:16306:53693/149‑1 (MQ=60)
                                                           CCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGG                                            <  GWNJ‑0478:712:GW2002102894th:2:2114:19673:45764/149‑1 (MQ=60)
                                                                                                GGCCCGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCagcttgttcagctggaagccattcagtc      <  GWNJ‑0478:712:GW2002102894th:2:2107:15463:64751/150‑29 (MQ=60)
                                                                                                    CGTCAAGATCGGCCTGGCGCTGGGCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGG  >  GWNJ‑0478:712:GW2002102894th:2:2209:14527:17971/1‑150 (MQ=60)
                                                                                                                           |                                                                                                                              
GTCATTTCCTCGGTATCTCCGCGGCCGCCCTGCTGGCGGGCTGCTCGTTCGGCCCGCGTCCCGCGCCGCCGCTCGCCACCACCCCGCCCGCGCCGCGGCCCGTCAAGATCGGCCTGGCGCTGGCCGGCGGCGCGGCGCGCGGCTTTGCCCATATCGGCGTGATCAAGGCGCTCGAAGCGCAAGGCATCCATGCCGACCTGGTCACCGGCACCAGCGCCGGCGCGGTCGTGGCGGCGCTTTACGCCAGCGG  >  AM260479/2616382‑2616631

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: