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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 581,821 +C coding (117/1734 nt) h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479581,8201.C100.0% 53.6 / NA 15S39T (AGC→ACC) h16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (6/9);  total (6/9)

CCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGT  >  AM260479/581702‑581939
                                                                                                                       |                                                                                                                       
ccATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTggcgg                                                                                                           <  1:846068/134‑1 (MQ=255)
         cgcTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCcg                                                                                                  <  2:118790/134‑1 (MQ=255)
                 gcgcATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAgg                                                                                          <  2:66833/134‑1 (MQ=255)
                                                                              ggCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTg                              >  2:676957/1‑133 (MQ=255)
                                                                              ggCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGa                             <  2:124596/134‑1 (MQ=255)
                                                                                           ggTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGaa                <  1:673831/134‑1 (MQ=255)
                                                                                              aCGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGc                                                                                <  1:514879/67‑1 (MQ=255)
                                                                                              aCGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGc                                                                                >  2:514879/1‑67 (MQ=255)
                                                                                                 ccgcGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCggg          <  1:676957/134‑1 (MQ=255)
                                                                                                  cgcgCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCgg                                                                             >  1:178233/1‑66 (MQ=255)
                                                                                                  cgcgCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCgg                                                                             <  2:178233/66‑1 (MQ=255)
                                                                                                         ggCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGcc                                                                     >  1:351323/1‑67 (MQ=255)
                                                                                                         ggCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGcc                                                                     <  2:351323/67‑1 (MQ=255)
                                                                                                         ggCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGt  >  1:484096/1‑134 (MQ=255)
                                                                                                         ggCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGt  >  2:970205/1‑134 (MQ=255)
                                                                                                                       |                                                                                                                       
CCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGT  >  AM260479/581702‑581939

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 35 ≤ 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
GCCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGAGTTCGGCCTGCGCACCAACAACG  >  AM260479/581701‑581966
                                                                                                                         |                                                                                                                                                 
GCCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGggggcctggc                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2209:16475:31991/150‑11 (MQ=60)
                           gcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCT                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2101:12766:8684/98‑1 (MQ=60)
                                                                         GGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCA                                                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2103:18898:98080/1‑49 (MQ=60)
                                                                                    GGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgcc                                   >  GWNJ‑0478:712:GW2002102894th:2:1113:14636:50515/1‑96 (MQ=60)
                                                                                           GGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgc                            >  GWNJ‑0478:712:GW2002102894th:2:1206:17061:56280/1‑97 (MQ=60)
                                                                                           GGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTC                            >  GWNJ‑0478:712:GW2002102894th:2:1214:6610:25646/1‑150 (MQ=60)
                                                                                           GGGTACGCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTC                            <  GWNJ‑0478:712:GW2002102894th:2:2110:2237:48911/150‑1 (MQ=60)
                                                                                                 GCCGCGCACGGCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGAGTT                      <  GWNJ‑0478:712:GW2002102894th:2:2110:13087:67430/150‑1 (MQ=60)
                                                                                                           GCCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGAGTTCGGCCTGCGC            >  GWNJ‑0478:712:GW2002102894th:2:2103:20467:51575/1‑150 (MQ=60)
                                                                                                            CCCCGGCGCCAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGAGTTCGGCCTGCGC            >  GWNJ‑0478:712:GW2002102894th:2:1113:8532:17050/1‑149 (MQ=60)
                                                                                                                     CAGCCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGAGTTCGGCCTGCGCACCAACAACG  <  GWNJ‑0478:712:GW2002102894th:2:2213:8547:79548/150‑1 (MQ=60)
                                                                                                                         |                                                                                                                                                 
GCCATTGCCGCGCTGGCGGCGCATGGACGCCATGGCGGCGGCCGCTGGCGGCGCTGCTGGCGCTGACGTTCGCGGCGCCGGCCTGGCCCCAGGGTACGCCGCGCACGGCCCCGGCGCCAGCACGCCGGTGGCGGCTGCGCCCGCGGCCAGGCCCCCGCAGCCGGTGCCCGCCATCTCGGGCGAGGCCAGCGACCAGGTCTATGACAACCTGAACCGCAGCGTGAAGGCGGGCCAGAAGTCCGAGTTCGGCCTGCGCACCAACAACG  >  AM260479/581701‑581966

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: