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 AM260479 160,041 G→C G102A (GGG→GCG)  h16_A0146 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479160,0410GC100.0% 62.1 / NA 18G102A (GGG→GCG) h16_A0146conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (10/8);  total (10/8)

GCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGGGGTACGACCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCCGCG  >  AM260479/159937‑160166
                                                                                                        |                                                                                                                              
gcgcGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATg                                                                                                   >  1:405978/1‑134 (MQ=255)
                    ggaaGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGcgacga                                                                                                     <  2:117045/109‑1 (MQ=255)
                    ggCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGcgacga                                                                                                     >  1:117045/1‑112 (MQ=255)
                    ggCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCaa                                                                               >  2:108627/1‑134 (MQ=255)
                                    gCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAAcc                                                                 <  1:228014/132‑1 (MQ=255)
                                    gCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAAcc                                                                 >  2:228014/1‑132 (MQ=255)
                                       gcgcGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTc                                                                                 >  2:730734/1‑113 (MQ=255)
                                       gcgcGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTc                                                                                 <  1:730734/113‑1 (MQ=255)
                                                         ccgcTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAAc                                                                  <  2:918703/110‑1 (MQ=255)
                                                         ccgcTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAAc                                                                  >  1:918703/1‑110 (MQ=255)
                                                                           gcggGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTAc                        >  2:358254/1‑134 (MQ=255)
                                                                                     gggCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCg                                                                          >  1:395365/1‑74 (MQ=255)
                                                                                     gggCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCg                                                                          <  2:395365/74‑1 (MQ=255)
                                                                                     gggCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATc                                                   <  1:75931/97‑1 (MQ=255)
                                                                                     gggCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATc                                                   >  2:75931/1‑97 (MQ=255)
                                                                                       gcggcgTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTc            <  1:269868/134‑1 (MQ=255)
                                                                                               gATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCcg    >  2:1012217/1‑134 (MQ=255)
                                                                                                 tGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCcgcg  <  1:154591/134‑1 (MQ=255)
                                                                                                        |                                                                                                                              
GCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGGGGTACGACCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCCGCG  >  AM260479/159937‑160166

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
CCACCAGCCCGCAGGGCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGGGGTACGACCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCCGCGTGCTGTGGTTCCGCT  >  AM260479/159922‑160181
                                                                                                                       |                                                                                                                                             
CCACCAGCCCGCAGGGCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGC                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1211:12837:75634/1‑150 (MQ=60)
                    ACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTctgtctct                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1108:1879:85630/1‑142 (MQ=60)
                                gagacagGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAAC                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2115:18666:59223/143‑1 (MQ=60)
                                      tcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGA                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1205:18261:100133/83‑1 (MQ=60)
                                               gtcgtagcgcacgccatcgacGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCctgtctcttatacacatctccgagcccacgagacggagc                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2113:8523:93500/22‑111 (MQ=60)
                                                gcGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCC                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1114:15014:37592/148‑1 (MQ=60)
                                                         CGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCtgtctctt                                                         >  GWNJ‑0478:712:GW2002102894th:2:2211:3342:49309/1‑141 (MQ=60)
                                                               cggagatgtgtataagagacagGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGAC                                                   <  GWNJ‑0478:712:GW2002102894th:2:1107:5342:22557/127‑1 (MQ=60)
                                                                                     GTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGctgtctcttatacacatctccgagcccacgagacggagctacatct                            >  GWNJ‑0478:712:GW2002102894th:2:1211:9885:16044/1‑104 (MQ=60)
                                                                                                     GGCGGCGTCGATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCCGCGTGCTG            <  GWNJ‑0478:712:GW2002102894th:2:1203:13697:3890/150‑1 (MQ=60)
                                                                                                               ATGGCGTGCGCTACGACCCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCCGCGTGCTGTGGTTCCGCT  <  GWNJ‑0478:712:GW2002102894th:2:1110:12116:11743/150‑1 (MQ=60)
                                                                                                                       |                                                                                                                                             
CCACCAGCCCGCAGGGCGCGACCGAGATGGCGTCGGGCAGTGGCGGCGAGGGCAGCGCGGCCAGGCCTGCCGCCGCTGCCGCGCCGTCCGGCGGGCCGCAGGGCGGCGTCGATGGCGTGGGGTACGACCGCCGCCGTCGGCGACGATGCACTACGCCAGCTTCGTCAACGGCGTGCAGAACCCCGACGGCCTGATCCGCTGGGAGCAGGACGGCAGCCGCTACCGGCTGGCGGTCGAAACCCGCGTGCTGTGGTTCCGCT  >  AM260479/159922‑160181

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