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,642 +A intergenic (‑39/+87) 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,6411.A100.0% 70.8 / NA 20intergenic (‑38/+88)h16_A3328/h16_A3329conserved hypothetical membrane protein/unassigned peptidase, M61 subfamily (glycyl aminopeptidase)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (8/12);  total (8/12)

CATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCG  >  AM260479/3602514‑3602764
                                                                                                                                |                                                                                                                           
cATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                       <  2:919914/134‑6 (MQ=255)
cATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaataca                                                                                                                       <  2:251004/134‑6 (MQ=255)
      gCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaatac                                                                                                                        <  1:420725/127‑5 (MQ=255)
      gCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaatac                                                                                                                        >  2:420725/1‑123 (MQ=255)
                  cgcCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGt                                                                                                      <  2:57116/134‑1 (MQ=255)
                    cccccATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTgc                                                                                                    >  1:104376/1‑134 (MQ=255)
                          tCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                              <  2:104376/134‑1 (MQ=255)
                                        gcgGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCggg                                                                                 <  1:394731/133‑1 (MQ=255)
                                           gTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCgggggc                                                                              >  1:347267/1‑132 (MQ=255)
                                                 aGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCg                                                                      <  2:347267/134‑1 (MQ=255)
                                                                   gcgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaata                                                                                                                         <  2:751568/65‑4 (MQ=255)
                                                                   gcgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGaata                                                                                                                         >  1:751568/1‑62 (MQ=255)
                                                                    cgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                   <  2:934721/134‑1 (MQ=255)
                                                                    cgATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCa                                                   <  2:735355/134‑1 (MQ=255)
                                                                      aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGaa                                                                                                                    >  2:85880/1‑68 (MQ=255)
                                                                      aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGaa                                                                                                                    <  1:85880/68‑1 (MQ=255)
                                                                      aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATtgca                                                                                          >  2:522709/1‑94 (MQ=255)
                                                                      aTCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATtgca                                                                                          <  1:522709/94‑1 (MQ=255)
                                                                         aGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGa                                              <  2:1141351/134‑1 (MQ=255)
                                                                                                   cGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                              >  2:549807/1‑61 (MQ=255)
                                                                                                   cGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                              <  1:549807/61‑1 (MQ=255)
                                                                                                   cGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                              <  2:86157/61‑1 (MQ=255)
                                                                                                   cGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTAt                                                                                              >  1:86157/1‑61 (MQ=255)
                                                                                                   cGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcg                    >  1:799968/1‑134 (MQ=255)
                                                                                                                          aGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCg  <  2:418258/129‑1 (MQ=255)
                                                                                                                          aGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCg  >  1:418258/1‑129 (MQ=255)
                                                                                                                                |                                                                                                                           
CATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCG  >  AM260479/3602514‑3602764

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
CCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGC  >  AM260479/3602513‑3602792
                                                                                                                                  |                                                                                                                                                      
CCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGA                                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1209:17701:26947/149‑1 (MQ=60)
 CATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAatacaggctgtctcttatac                                                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1112:8520:38293/1‑129 (MQ=60)
    catacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTG                                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1107:6195:37683/99‑1 (MQ=60)
      GGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCG                                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1106:11601:46062/1‑150 (MQ=60)
                          aagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGG                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1214:17393:73003/91‑1 (MQ=60)
                             GACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGcg                                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:1204:14670:6672/1‑147 (MQ=60)
                              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)
                                                     TCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccg                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2110:15291:67048/1‑95 (MQ=60)
                                                                           gagtgtctcgtgggctcggagatgtgtataagagacagCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCC                                                         <  GWNJ‑0478:712:GW2002102894th:2:2216:2057:91181/112‑1 (MQ=60)
                                                                                     GATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtctt                                                >  GWNJ‑0478:712:GW2002102894th:2:1107:15356:39132/1‑91 (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)
                                                                                                                                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)
                                                                                                                                  |                                                                                                                                                      
CCATCAGGCTGGTGCACAGCGCCCCCATCGACACCACCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGC  >  AM260479/3602513‑3602792

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: