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 2,303,559 A→C S217R (AGT→AGG h16_A2116 ← Putative DNA polymerase bacteriophage‑type,uracil‑DNA glycosylase superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,303,5590AC100.0% 82.5 / NA 22S217R (AGT→AGGh16_A2116Putative DNA polymerase bacteriophage‑type,uracil‑DNA glycosylase superfamily
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/10);  total (12/10)

TAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACACTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCCGC  >  AM260479/2303443‑2303686
                                                                                                                    |                                                                                                                               
tAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAgcgc                                                                                                                >  2:266919/1‑134 (MQ=255)
        caccacGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTggc                                                                                                        <  2:776250/134‑1 (MQ=255)
                         cATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGAt                                                                                                             <  2:872672/112‑1 (MQ=255)
                         cATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGAt                                                                                                             >  1:872672/1‑112 (MQ=255)
                                            gcgcAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGgcgc                                                                                                   <  2:844631/103‑1 (MQ=255)
                                            gcgcAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGgcgc                                                                                                   >  1:844631/1‑103 (MQ=255)
                                            gcgcAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCgg                                                                    >  2:527300/1‑134 (MQ=255)
                                            gcgcAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCgg                                                                    >  1:833406/1‑134 (MQ=255)
                                                 gCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGgcgc                                                                                                   <  1:178442/98‑1 (MQ=255)
                                                 gCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGgcgc                                                                                                   >  2:178442/1‑98 (MQ=255)
                                                                    gctgcGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcgggcg                                            >  1:125486/1‑134 (MQ=255)
                                                                     ctgcGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcgggcgg                                           <  1:527300/134‑1 (MQ=255)
                                                                                 gACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAg                               >  1:778547/1‑134 (MQ=255)
                                                                                 gACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAg                               >  2:966466/1‑134 (MQ=255)
                                                                                      cgccgcAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGgcgc                                                                                                   <  2:386057/61‑1 (MQ=255)
                                                                                      cgccgcAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGgcgc                                                                                                   >  1:386057/1‑61 (MQ=255)
                                                                                      cgccgcAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGt                          <  1:793956/134‑1 (MQ=255)
                                                                                           cAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcg                                                >  2:524725/1‑107 (MQ=255)
                                                                                           cAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACcg                                                <  1:524725/107‑1 (MQ=255)
                                                                                            aaaGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGa                    <  2:778547/134‑1 (MQ=255)
                                                                                                       gCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTg         >  2:60093/1‑134 (MQ=255)
                                                                                                              gATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCcgc  <  1:60093/134‑1 (MQ=255)
                                                                                                                    |                                                                                                                               
TAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACACTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCCGC  >  AM260479/2303443‑2303686

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
AGGGTGCGCAGCAGGTAAGCCGGATGGTAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACACTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCCGCGCGCCAGCCCGAGCG  >  AM260479/2303416‑2303701
                                                                                                                                               |                                                                                                                                              
acgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1106:13504:27618/101‑1 (MQ=60)
      acagGCAGGTAAGCCGGATGGTAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTACCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACC                                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1111:7713:45587/146‑1 (MQ=60)
                                     CCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACtgtctct                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2206:5994:62751/1‑143 (MQ=60)
                                       cggagatgtgtataagagacagGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCAC                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1113:18523:52571/128‑1 (MQ=60)
                                                        GCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACctgtctcttatacacat                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2205:5613:86142/1‑133 (MQ=60)
                                                        GCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGG                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2205:14862:13300/1‑150 (MQ=60)
                                                                                GCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGG                                                          >  GWNJ‑0478:712:GW2002102894th:2:1109:2936:39231/1‑150 (MQ=60)
                                                                                          gtataagagacagGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTG                                                <  GWNJ‑0478:712:GW2002102894th:2:2101:17730:44041/137‑1 (MQ=60)
                                                                                             GTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGC                                             >  GWNJ‑0478:712:GW2002102894th:2:2201:20233:38958/1‑150 (MQ=60)
                                                                                               GCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCAC                                           <  GWNJ‑0478:712:GW2002102894th:2:2101:9476:62322/150‑1 (MQ=60)
                                                                                                  GCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACctttcgcttatacacatatccgaacctacgagaccggggtccacctcgtatgccgtctt                                        >  GWNJ‑0478:712:GW2002102894th:2:1209:15542:61833/1‑91 (MQ=60)
                                                                                                                GCGCCGCAAAGCGGCCCAGCACCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACAT                          <  GWNJ‑0478:712:GW2002102894th:2:2203:19057:46117/150‑1 (MQ=60)
                                                                                                                                     CCACGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCCGCGCGCCAGCCCG      >  GWNJ‑0478:712:GW2002102894th:2:1110:4024:72203/1‑149 (MQ=60)
                                                                                                                                        CGATCACCCTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCCGCGCGCCAGCCCGAGCG  <  GWNJ‑0478:712:GW2002102894th:2:1104:16003:16211/150‑1 (MQ=60)
                                                                                                                                               |                                                                                                                                              
AGGGTGCGCAGCAGGTAAGCCGGATGGTAGGTCACCACCACGGGGATACCTTCATAGCTATGTACGGTGCCGCGCAGCTTGCCGATGGGCGTGGTGCTGCGCAGCAGCGACTGCGCCGCAAAGCGGCCCAGCACCACGATCACACTGGGCTTCACCAGCGCGATCTGGCGGCGCAGGAACGGCTCGCACTGCGCCACCTCCTCGGGCTCCGGGTTGCGGTTACCGGGCGGGCGGCACTTGAGCACGTTGGCGATAAACACATTGCGCCCGCGCGCCAGCCCGAGCG  >  AM260479/2303416‑2303701

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: