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 149 26.7 2337422 87.5% 2045244 105.0

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% 86.1 / NA 20intergenic (‑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/10);  total (10/10)

GCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGC  >  AM260479/3602561‑3602794
                                                                                                                           |                                                                                                               
gCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATgg                                                                                                      >  2:283650/1‑134 (MQ=255)
   ggTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGt                                                                                                   >  1:119855/1‑134 (MQ=255)
     tCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                                                 <  1:167866/134‑1 (MQ=255)
     tCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGc                                                                                                 <  2:414618/134‑1 (MQ=255)
        cTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGt                                                                                              >  1:1093281/1‑134 (MQ=255)
                                 cgcCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAg                                                                     <  1:517923/134‑1 (MQ=255)
                                                      aGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCagcgcg                                                <  2:1093281/134‑1 (MQ=255)
                                                                          tAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGtgtg                                                                                    >  1:327503/1‑78 (MQ=255)
                                                                          tAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGtgtg                                                                                    <  2:327503/78‑1 (MQ=255)
                                                                           aGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgtg                                                                                      <  2:1014493/75‑1 (MQ=255)
                                                                           aGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgtg                                                                                      >  1:1014493/1‑75 (MQ=255)
                                                                                  aTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATgg                                                                                                      <  1:410798/52‑1 (MQ=255)
                                                                                  aTACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATgg                                                                                                      >  2:410798/1‑52 (MQ=255)
                                                                                                 gggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTg                                   >  1:471012/1‑104 (MQ=255)
                                                                                                 gggATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTg                                   <  2:471012/104‑1 (MQ=255)
                                                                                                    aTCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTgcgc  >  2:326835/1‑134 (MQ=255)
                                                                                                           gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgtg                                                                                      >  2:408873/1‑43 (MQ=38)
                                                                                                           gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGgtg                                                                                      <  1:408873/43‑1 (MQ=38)
                                                                                                           gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGaaaa                                                                         >  1:799203/1‑56 (MQ=255)
                                                                                                           gcgTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGaaaa                                                                         <  2:799203/56‑1 (MQ=255)
                                                                                                                           |                                                                                                               
GCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGC  >  AM260479/3602561‑3602794

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
CCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGCAG  >  AM260479/3602549‑3602824
                                                                                                                                            |                                                                                                                                        
CCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTG                                                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1108:1956:42489/1‑150 (MQ=60)
     GCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTG                                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2213:11731:63163/1‑150 (MQ=60)
           gcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCTGGGGCGATGGCGTGCGGTGGCGGGT                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1201:3663:15984/87‑1 (MQ=60)
            tcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGT                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1208:15898:78684/83‑1 (MQ=60)
                GTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGC                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1110:11534:5955/150‑1 (MQ=60)
                            ttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACT                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1208:7284:67682/75‑1 (MQ=60)
                                            GCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCC                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1215:7187:35474/150‑1 (MQ=60)
                                               gtttttgttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagcagcGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCC                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2104:17317:72121/59‑1 (MQ=60)
                                                                       CGCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGctgtctcttatacacatctccgagcccacgagacactcgct                                                          >  GWNJ‑0478:712:GW2002102894th:2:1202:18184:79189/1‑108 (MQ=60)
                                                                        GCCCCCGGTAGATGTAGTGCGAATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTctgtctcttatacacatctccgagcccacgagacactcgctaatc                                                        >  GWNJ‑0478:712:GW2002102894th:2:2208:6152:58420/1‑105 (MQ=60)
                                                                                              ATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTgtctcttatacac                                   >  GWNJ‑0478:712:GW2002102894th:2:1213:8587:4629/1‑136 (MQ=60)
                                                                                                        GTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgct                        >  GWNJ‑0478:712:GW2002102894th:2:2112:6032:70909/1‑86 (MQ=60)
                                                                                                                      CGCGTATTGCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACctgtctcttatacacatctccgagccc          >  GWNJ‑0478:712:GW2002102894th:2:1115:16337:61470/1‑123 (MQ=60)
                                                                                                                              GCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCctgtctcttatacacatctccgagcccacgagacactcgctaatctc  >  GWNJ‑0478:712:GW2002102894th:2:1102:17244:39520/1‑103 (MQ=60)
                                                                                                                              GCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGC    >  GWNJ‑0478:712:GW2002102894th:2:1211:7660:44290/1‑148 (MQ=60)
                                                                                                                              GCATGCAGCGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGC    >  GWNJ‑0478:712:GW2002102894th:2:1211:7871:99644/1‑148 (MQ=60)
                                                                                                                                            |                                                                                                                                        
CCATCGCGGTGGGCAGGTCGCTGAAGTGCAGCGCGATCAGCGTGGCGCCGATCACCCCGGTGGCCGAGCGCAGCCCCCGGTAGATGTAGTGCGATACAGGAAGGTGCGGGGATCGTGCGCGTATTGCATGCAGCGGGGGGCGATGGCGTGCGGTGGCGGGTGTGGCAGCGGAAAACTAGTTCACTGGCCCAGCCAGCGCGAGCGCGCCGCGTGCCGTCCGCTCTCGAGCTTCACCTTGAACTGCGCCGCCGGCTCGCGTGCCAGCGTCACCGGCAG  >  AM260479/3602549‑3602824

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