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 155 25.8 2285122 88.0% 2010907 104.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,602,688 (G)6→5 intergenic (‑85/+41) h16_A3328 ← / ← h16_A3329 conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,602,6830G.100.0% 100.9 / NA 23intergenic (‑80/+46)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (10/13);  total (10/13)

GCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGC  >  AM260479/3602554‑3602794
                                                                                                                                  |                                                                                                               
gcgGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCggg                                                                                                               <  1:394731/133‑1 (MQ=255)
   gTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggc                                                                                                            >  1:347267/1‑132 (MQ=255)
         aGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCg                                                                                                    <  2:347267/134‑1 (MQ=255)
                            cgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                                                 <  2:735355/134‑1 (MQ=255)
                            cgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                                                 <  2:934721/134‑1 (MQ=255)
                                 aGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGa                                                                            <  2:1141351/134‑1 (MQ=255)
                                                           cGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcg                                                  >  1:799968/1‑134 (MQ=255)
                                                                                  aGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCg                                >  1:418258/1‑129 (MQ=255)
                                                                                  aGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCg                                <  2:418258/129‑1 (MQ=255)
                                                                                       gaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGaaa                                                                          >  2:1140673/3‑82 (MQ=255)
                                                                                       gaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGaaa                                                                          <  1:1140673/80‑1 (MQ=255)
                                                                                       gaaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCg                      <  2:799968/132‑1 (MQ=255)
                                                                                        aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGtgtg                                                                                    >  1:124242/2‑71 (MQ=255)
                                                                                        aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGtgtg                                                                                    <  2:124242/70‑1 (MQ=255)
                                                                                        aaTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCtctc                       >  1:708254/2‑132 (MQ=255)
                                                                                            cAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCtt                 <  1:18753/134‑1 (MQ=255)
                                                                                            cAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCtt                 >  2:18753/1‑134 (MQ=255)
                                                                                                       ggggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACt      >  1:1084071/1‑134 (MQ=255)
                                                                                                        gggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcgcg                                                <  2:1138128/91‑1 (MQ=255)
                                                                                                        gggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcgcg                                                >  1:1138128/1‑91 (MQ=255)
                                                                                                        gggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCtt                 >  2:9511/1‑122 (MQ=255)
                                                                                                        gggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCtt                 <  1:9511/122‑1 (MQ=255)
                                                                                                           aTCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTgcgc  <  2:1084071/134‑1 (MQ=255)
                                                                                                                       ttGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCg                                >  2:877454/1‑92 (MQ=255)
                                                                                                                       ttGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCg                                <  1:877454/92‑1 (MQ=255)
                                                                                                                                  |                                                                                                               
GCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGC  >  AM260479/3602554‑3602794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGCAGCACC  >  AM260479/3602543‑3602828
                                                                                                                                                  |                                                                                                                                            
ggctcggagatgtgtataagagacaGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGAT                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1215:2587:45547/125‑1 (MQ=60)
          CGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGG                                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1206:5376:83520/149‑1 (MQ=60)
                                             gagtgtctcgtgggctcggagatgtgtataagagacagCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCC                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2216:2057:91181/112‑1 (MQ=60)
                                                       GATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGC                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2115:7883:22592/1‑150 (MQ=60)
                                                                              GCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCT                                                            >  GWNJ‑0478:712:GW2002102894th:2:1209:5935:13704/1‑150 (MQ=60)
                                                                                    GGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgt                                                      >  GWNJ‑0478:712:GW2002102894th:2:1108:2863:61625/1‑101 (MQ=60)
                                                                                    GGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCG                                                        >  GWNJ‑0478:712:GW2002102894th:2:2109:3035:32547/1‑148 (MQ=60)
                                                                                        gggctcggagatgtgtataagagacaGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTT                                                   <  GWNJ‑0478:712:GW2002102894th:2:2203:8729:45942/123‑1 (MQ=60)
                                                                                                  AATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTG                                       >  GWNJ‑0478:712:GW2002102894th:2:2214:10891:16651/1‑150 (MQ=60)
                                                                                                    ATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCctgtctcttatacacatctccgagcccac                                      >  GWNJ‑0478:712:GW2002102894th:2:2216:1421:77896/1‑121 (MQ=60)
                                                                                                      ACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGC                                    <  GWNJ‑0478:712:GW2002102894th:2:1101:2878:62954/150‑1 (MQ=60)
                                                                                                      ACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGC                                    <  GWNJ‑0478:712:GW2002102894th:2:1103:12657:23217/150‑1 (MQ=60)
                                                                                                              GTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaacaaaaaagcacgacagtagaaa                              >  GWNJ‑0478:712:GW2002102894th:2:1107:18761:22876/1‑50 (MQ=60)
                                                                                                                                         CAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaataaaa  >  GWNJ‑0478:712:GW2002102894th:2:1107:11679:84448/1‑63 (MQ=60)
                                                                                                                                                  |                                                                                                                                            
ACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGCAGCACC  >  AM260479/3602543‑3602828

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: