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 AM260480 2,431,093 Δ1 bp coding (2380/2469 nt) h16_B2148 → membrane carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604802,431,0920T.100.0% 38.8 / NA 10coding (2379/2469 nt)h16_B2148membrane carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (4/6);  total (4/6)

CCCGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATTCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGCC  >  AM260480/2430963‑2431168
                                                                                                                                 |                                                                            
cccGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGAtcggg                                                                          <  2:1025819/134‑5 (MQ=255)
cccGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGAtcggg                                                                          <  2:339024/134‑5 (MQ=255)
  cGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCggggg                                                                       >  1:598920/1‑134 (MQ=255)
   ggCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCggggg                                                                       >  2:622386/1‑133 (MQ=255)
       ccgcgccCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGAc                                                                  >  2:376519/1‑134 (MQ=255)
             ccGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGAt                                                            <  1:376519/134‑1 (MQ=255)
                     gcgTCGGCCGGCGCTTCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGc                                                    <  2:987796/134‑1 (MQ=255)
                     gcgTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGc                                                    <  1:784039/134‑1 (MQ=255)
                          ggccggcGCTGCCCACGCCGGAGGCCGAGGTCCCGCCCCCGGCCGGCGCGGCATCGCCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCGGGCGCTGGATCGGGGGGGGACGGTGATTGCGGGGCCGCAg                                               >  1:987796/1‑134 (MQ=255)
                                                             ccgccgGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCAcc            <  2:88127/134‑1 (MQ=255)
                                                                    ccggcGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCATCAg     <  1:622386/134‑1 (MQ=255)
                                                                       gcgcGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGcc  <  2:598920/134‑1 (MQ=255)
                                                                                                                                 |                                                                            
CCCGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATTCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGCC  >  AM260480/2430963‑2431168

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ 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
AGGCCAGCGAGACCCCGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATTCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGCCGAACCAGCGCCTAACGCCGCCCCTGCCTCCCCGGCTTCCGCAAATCCGCCCGCACCGCCTGCGCCCC  >  AM260480/2430950‑2431235
                                                                                                                                               |                                                                                                                                              
AGGCCAGCGAGACCCCGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGG                                                                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:2107:7367:61324/1‑150 (MQ=60)
                        GCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCCCGCCCCCGGCCGGCGCGGCATCGCCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCGGGCGCTGGATCGGGGGGGGACGGTGATTGCGGGGCCGCAGG                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2214:11859:43406/1‑150 (MQ=60)
                         CCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGA                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1209:5347:3217/150‑1 (MQ=60)
                            tcaggctcagtatgtctagtttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACA                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1210:9352:31702/53‑1 (MQ=60)
                                 CGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCG                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2201:6498:78277/150‑1 (MQ=60)
                                                       tttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGC                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1215:7331:8951/80‑1 (MQ=60)
                                                                                GCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCATCAGGCCGAACCAGCGCCT                                                         <  GWNJ‑0478:712:GW2002102894th:2:2108:11139:15620/150‑1 (MQ=60)
                                                                                                 gggctcggagatgtgtataagagacagGTCCAGCCTGGCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGCCGAACCAGCGCCTAACGCCGCCCCTGCCTC                                        <  GWNJ‑0478:712:GW2002102894th:2:2204:2023:83927/123‑1 (MQ=60)
                                                                                                                           AGTCCAGCCTGGCGCTGGATCGGTGGGAGACGGTGATTGCGGTGCCGCAGGACActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaagtcggaacacaatcatt              >  GWNJ‑0478:712:GW2002102894th:2:2215:14692:28283/1‑54 (MQ=60)
                                                                                                                                      GCGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaataaaaaaggacagaccagctccagaccat   >  GWNJ‑0478:712:GW2002102894th:2:1103:14762:28170/1‑43 (MQ=60)
                                                                                                                                       CGCTGGATCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGCCGAACCAGCGCCTgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtc  >  GWNJ‑0478:712:GW2002102894th:2:2210:5845:12998/1‑95 (MQ=60)
                                                                                                                                               |                                                                                                                                              
AGGCCAGCGAGACCCCGGCGCCGCGCCCGGCACCGCGTCGGCCGGCGCTGCCCACGCCGGAGGCCGAGGTCGCGCCGCCGGCCGGCGCGGCATCGTCTGCCGCGGTGGATGTCGGCCCGGATGAGTCCAGCCTGGCGCTGGATTCGGGGGGAGACGGTGATTGCGGTGCCGCAGGACATGAGCGCGGGCAGCGCTCCGGTTGCGCCACCCCCTCAGGCCGAACCAGCGCCTAACGCCGCCCCTGCCTCCCCGGCTTCCGCAAATCCGCCCGCACCGCCTGCGCCCC  >  AM260480/2430950‑2431235

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