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 582,759 G→T C352F (TGT→TTT)  h16_A0555 → zinc metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479582,7590GT100.0% 73.6 / NA 20C352F (TGT→TTT) h16_A0555zinc metalloprotease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (11/9);  total (11/9)

ATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCC  >  AM260479/582689‑582884
                                                                      |                                                                                                                             
aTTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTgcgcgc                                                                <  2:993875/134‑1 (MQ=255)
  tGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCgcgg                                                              >  2:97566/1‑134 (MQ=255)
              ggACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGcgcacg                                                   >  1:691463/1‑133 (MQ=255)
                  cgcgTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACt                                              >  1:166350/1‑134 (MQ=255)
                    cgTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGc                                            <  1:97566/134‑1 (MQ=255)
                       cAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCa                                                                                      <  2:153362/89‑1 (MQ=255)
                       cAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCa                                                                                      >  1:153362/1‑89 (MQ=255)
                        aGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCgg                                         >  2:339301/1‑133 (MQ=255)
                        aGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCgg                                         >  2:47804/1‑133 (MQ=255)
                        aGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCgg                                         >  2:481119/1‑133 (MQ=255)
                                  gCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCa                              >  1:437414/1‑134 (MQ=255)
                                        gcAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTgcgc                        <  2:437414/134‑1 (MQ=255)
                                           aGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCt                     <  1:47804/134‑1 (MQ=255)
                                                     cacGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTg                                                         <  1:581136/88‑1 (MQ=255)
                                                     cacGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTg                                                         >  2:581136/1‑88 (MQ=255)
                                                     cacGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGc                                            <  1:912232/101‑1 (MQ=255)
                                                     cacGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGc                                            >  2:912232/1‑101 (MQ=255)
                                                         ccTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTa       <  1:481119/134‑1 (MQ=255)
                                                            gAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACgg    >  1:76390/1‑134 (MQ=255)
                                                              aTACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGgcc  <  2:166350/134‑1 (MQ=255)
                                                                      |                                                                                                                             
ATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCC  >  AM260479/582689‑582884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
GCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCACCGAGC  >  AM260479/582610‑582899
                                                                                                                                                     |                                                                                                                                            
gtagctccgtctcgtgggctcggagatgtgtataagagacagGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTT                                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1106:10899:79960/108‑1 (MQ=60)
          CGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATG                                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2215:1621:33426/150‑1 (MQ=60)
                ggctcggagatgtgtataagagacagGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCC                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1113:17230:31040/124‑1 (MQ=60)
                                                                              GATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGC                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2113:19669:51933/1‑149 (MQ=60)
                                                                                  GCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTG                                                            >  GWNJ‑0478:712:GW2002102894th:2:1110:5298:79306/1‑150 (MQ=60)
                                                                                       cgtgggctcggagatgtgtataagagacagGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCG                                                         <  GWNJ‑0478:712:GW2002102894th:2:2104:9244:66819/118‑1 (MQ=60)
                                                                                       CATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCctgtctcttatacacatctccgagcccacg                                                        >  GWNJ‑0478:712:GW2002102894th:2:1111:8579:100082/1‑119 (MQ=60)
                                                                                                  agatgtgtataagagacagGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCAC                                            <  GWNJ‑0478:712:GW2002102894th:2:2211:18219:41312/131‑1 (MQ=60)
                                                                                                  GCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCAC                                            <  GWNJ‑0478:712:GW2002102894th:2:1106:19632:54528/150‑1 (MQ=60)
                                                                                                  GCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCAC                                            >  GWNJ‑0478:712:GW2002102894th:2:1213:17494:65605/1‑150 (MQ=60)
                                                                                                                 ggctcggagatgtgtataagagacagGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGC                             <  GWNJ‑0478:712:GW2002102894th:2:2103:9469:84281/124‑1 (MQ=60)
                                                                                                                         AAGGTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTAC                     <  GWNJ‑0478:712:GW2002102894th:2:1102:6968:25912/150‑1 (MQ=60)
                                                                                                                            GTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAActgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaaagaaaatc                  >  GWNJ‑0478:712:GW2002102894th:2:2108:18484:64465/1‑57 (MQ=60)
                                                                                                                            GTCGCCAACACGCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGC                  >  GWNJ‑0478:712:GW2002102894th:2:1107:5464:20021/1‑150 (MQ=60)
                                                                                                                                       GCCTGAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCAC       <  GWNJ‑0478:712:GW2002102894th:2:1214:11468:8917/150‑1 (MQ=60)
                                                                                                                                           GAATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaa    >  GWNJ‑0478:712:GW2002102894th:2:2111:3093:55868/1‑79 (MQ=60)
                                                                                                                                            AATACGAGTTTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCACCGAGC  >  GWNJ‑0478:712:GW2002102894th:2:1209:19039:20247/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                            
GCGGCTGCAGCGCGCATCGGGGACCGCGGAAAATTCGGTCCCGTCTTATGTGCGGACTCACCCGCTGACCTCCGAGCGGATTGCCGACATGCAGGACCGCGTCAGCCATGCCAGCAAGCGCAAGGTCGCCAACACGCCTGAATACGAGTGTGCGCGGATGCGGGCCCGCGTGATCCAGGAATCGTCGCCCAGCGACCTGCAGAACCTGCGCGCGGCGCTGCGCACGCAACTGCCGGGAGCAACGGCACTGCGCGCTCCGGCGCTGCACTACGGCCTTGCCTTCACCGAGC  >  AM260479/582610‑582899

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