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 1,191,771 Δ1 bp coding (5/390 nt) h16_A1095 ← Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,191,7690C.100.0% 68.0 / NA 16coding (7/390 nt)h16_A1095Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/7);  total (9/7)

TCGCGTCGCGCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGC  >  AM260479/1191642‑1191874
                                                                                                                               |                                                                                                         
tcgcgtcgcgCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACg                                                                                                    >  2:1011582/1‑134 (MQ=255)
    gtcgcgCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCg                                                                                                <  2:647741/134‑1 (MQ=255)
    gtcgcgCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCg                                                                                                <  1:78496/134‑1 (MQ=255)
                                    gCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGa                                                                >  1:519248/1‑134 (MQ=255)
                                                   gCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCg                                                   >  2:973306/1‑132 (MQ=255)
                                                   gCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCg                                                   <  1:973306/132‑1 (MQ=255)
                                                       tGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGc                                             >  2:357680/1‑134 (MQ=255)
                                                        ggCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGActct                                                                              >  2:1001419/1‑100 (MQ=255)
                                                        ggCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGActct                                                                              <  1:1001419/100‑1 (MQ=255)
                                                                      gcgcggcgAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCg                                                                                                <  1:294635/68‑1 (MQ=255)
                                                                      gcgcggcgAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCg                                                                                                >  2:294635/1‑68 (MQ=255)
                                                                       cgcggcgAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGccgca                                                                                                       <  1:387323/61‑4 (MQ=255)
                                                                       cgcggcgAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGccgca                                                                                                       >  2:387323/1‑58 (MQ=255)
                                                                                       aCCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCg                                                                                                <  1:4407/51‑1 (MQ=255)
                                                                                       aCCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCg                                                                                                >  2:4407/1‑51 (MQ=255)
                                                                                               gccgGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCcagcgcag                                                                    <  2:271937/71‑1 (MQ=255)
                                                                                               gccgGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCcagcgcag                                                                    >  1:271937/1‑71 (MQ=255)
                                                                                                  ggCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGc  >  2:1139783/1‑134 (MQ=255)
                                                                                                                               |                                                                                                         
TCGCGTCGCGCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGC  >  AM260479/1191642‑1191874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ 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
GCATCGCGTCGCGCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGCGGTCGGCCGCGTTCACGACGGCGTCTCCTGCGAGCCCGGCACGC  >  AM260479/1191639‑1191918
                                                                                                                                    |                                                                                                                                                   
aagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCacaacag                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1206:9328:43649/91‑8 (MQ=60)
      catacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagatcgccagcgggcCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGAC                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1201:10192:66796/85‑1 (MQ=60)
      CGTCGCGCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGAC                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1104:18113:98221/150‑1 (MQ=60)
       ttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCacaacaggaggag                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1101:16216:30154/81‑14 (MQ=60)
        aggaagctcgctgatttttgtgttattttattttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTC                                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2101:20001:48884/40‑1 (MQ=60)
                        gagatgtgtataagagacagTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGAC                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2207:11757:85694/130‑1 (MQ=60)
                        GCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGAC                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1215:15091:41040/1‑150 (MQ=60)
                                                   ACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGC                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2206:14654:23929/150‑1 (MQ=60)
                                                                                   GGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgt                                                >  GWNJ‑0478:712:GW2002102894th:2:1116:20724:38324/1‑102 (MQ=60)
                                                                                                    gctcggagatgtgtataagagacaggttgtgCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGCGGTCGGCCGCGTTC                                <  GWNJ‑0478:712:GW2002102894th:2:1103:17621:42225/119‑1 (MQ=60)
                                                                                                                       CTCCTGTTGTGCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCctgtctcttatacacatctccgagcccacgagacactcgctaatc            >  GWNJ‑0478:712:GW2002102894th:2:2112:6682:76278/1‑105 (MQ=60)
                                                                                                                                   CCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaacgt    >  GWNJ‑0478:712:GW2002102894th:2:1112:10518:49887/1‑59 (MQ=60)
                                                                                                                                   CCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGCGGTCGGCCGCGTTCACGctgtctcttatacacatctccgagccc  >  GWNJ‑0478:712:GW2002102894th:2:2215:16071:49611/1‑122 (MQ=60)
                                                                                                                                    |                                                                                                                                                   
GCATCGCGTCGCGCGGCCTGCGCCGCCTGCAGGCGCGCGGCTTCTTCCTGCACCGCCATGGCCAGCGTTGGCGGCGCGGCGAAGGGCAGGACCGCCGCGCCGGCCCGCTGGCGATCCTCCTCCTGTTGTGCCCGCACGTGCGGCCGCGTGCACCGACTCTGCCAGCGCAGCGGACGGCGCGTAAGCGGGCAGCTCCCGCAGCAGCGCCGAAAGCCGGTCCTTGCCGTCGATAAAGCGGTCGGCCGCGTTCACGACGGCGTCTCCTGCGAGCCCGGCACGC  >  AM260479/1191639‑1191918

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: