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 3,529,607 +C intergenic (‑121/‑183) secA ← / → h16_A3265 preprotein translocase subunit SecA/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,529,6061.C100.0% 107.1 / NA 29intergenic (‑120/‑184)secA/h16_A3265preprotein translocase subunit SecA/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (15/14);  total (15/14)

CGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGC  >  AM260479/3529477‑3529730
                                                                                                                                   |                                                                                                                            
cGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCAcggccg                                                                                                                           <  1:492959/134‑3 (MQ=255)
                   tCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGt                                                                                                          <  1:314128/133‑1 (MQ=255)
                              gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCAcggccg                                                                                                                           >  1:210932/1‑102 (MQ=255)
                              gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCAcggccg                                                                                                                           <  2:210932/104‑3 (MQ=255)
                              gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGcc                                                                                                              >  2:741699/1‑118 (MQ=255)
                              gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGcc                                                                                                              <  1:741699/118‑1 (MQ=255)
                                  cGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGggcaggca                                                                                          <  1:629565/134‑1 (MQ=255)
                                        gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTg                                                                                    >  2:73654/1‑134 (MQ=255)
                                                        ggCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCAcg                                                                    <  1:648253/134‑1 (MQ=255)
                                                          ctgctCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGc                                                                                   >  2:863475/1‑117 (MQ=255)
                                                          ctgctCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGc                                                                                   <  1:863475/117‑1 (MQ=255)
                                                          ctgctCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCAcgcg                                                                  >  1:591198/1‑134 (MQ=255)
                                                                 ggcggcCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCaa                                                                                                   <  2:992743/94‑1 (MQ=255)
                                                                 ggcggcCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCaa                                                                                                   >  1:992743/1‑94 (MQ=255)
                                                                                              cgcAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGcc                                                                                                              >  1:888901/1‑54 (MQ=255)
                                                                                              cgcAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGcc                                                                                                              <  2:888901/54‑1 (MQ=255)
                                                                                              cgcAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTgcg                              <  1:148378/134‑1 (MQ=255)
                                                                                                  aaaCGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCa                          >  2:870350/1‑134 (MQ=255)
                                                                                                       ccGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCg                                                                          <  2:326621/81‑1 (MQ=255)
                                                                                                       ccGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCg                                                                          >  1:326621/1‑81 (MQ=255)
                                                                                                           gTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCca                   >  2:896401/1‑132 (MQ=255)
                                                                                                                        cGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgcc    >  1:960007/1‑134 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGAtt                                                   <  2:38439/86‑1 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGAtt                                                   >  1:38439/1‑86 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCAt                         >  2:533301/1‑112 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCAt                         <  1:533301/112‑1 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTg                       >  1:599217/1‑114 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTg                       <  2:599217/114‑1 (MQ=255)
                                                                                                                         ggggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgccg   <  1:870350/134‑1 (MQ=255)
                                                                                                                          gtgTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGgccgc  >  1:80257/3‑134 (MQ=255)
                                                                                                                          gggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGAtt                                                   <  2:549212/85‑1 (MQ=255)
                                                                                                                          gggTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGAtt                                                   >  1:549212/1‑85 (MQ=255)
                                                                                                                                   |                                                                                                                            
CGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGC  >  AM260479/3529477‑3529730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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
CCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGCCGCAGCTTCAATCC  >  AM260479/3529476‑3529744
                                                                                                                                     |                                                                                                                                         
CCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCC                                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1214:18090:73891/150‑1 (MQ=60)
              agGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGC                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2114:21273:49489/148‑1 (MQ=60)
                GCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCtgtctcttatacac                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:1115:17586:43397/1‑135 (MQ=60)
                   TTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGC                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1204:14380:18175/149‑1 (MQ=60)
                                  CCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCcttgcgg                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2108:8238:80596/150‑8 (MQ=60)
                                         cagccccagGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACG                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2206:14594:38123/141‑1 (MQ=60)
                                            CCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGT                                                                               >  GWNJ‑0478:712:GW2002102894th:2:2107:8159:32578/1‑150 (MQ=60)
                                                   GATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCttgcctgctgtctcttatacacatctccgagcc                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2215:2001:6637/1‑117 (MQ=60)
                                                        GGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGA                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2109:20968:97414/150‑1 (MQ=60)
                                                                                GAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttg                                           >  GWNJ‑0478:712:GW2002102894th:2:2207:21221:97259/1‑84 (MQ=60)
                                                                                         GATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCtgtctcttatacacatctccgagcccacgagacggagc                                  >  GWNJ‑0478:712:GW2002102894th:2:1204:2899:99487/1‑112 (MQ=60)
                                                                                              GCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGT                             <  GWNJ‑0478:712:GW2002102894th:2:1111:7192:71493/150‑1 (MQ=60)
                                                                                                   agagacaGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGG                        <  GWNJ‑0478:712:GW2002102894th:2:2103:6879:35602/143‑1 (MQ=60)
                                                                                                          GAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCG                 >  GWNJ‑0478:712:GW2002102894th:2:2213:11412:21893/1‑150 (MQ=60)
                                                                                                           AGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCtgtctcttatacacatctccgagcccacgagac                >  GWNJ‑0478:712:GW2002102894th:2:1105:10349:63054/1‑117 (MQ=60)
                                                                                                           AGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCctgtct                >  GWNJ‑0478:712:GW2002102894th:2:2107:7467:82823/1‑144 (MQ=60)
                                                                                                            GTTTATCATGTGCCGGTGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGCC               >  GWNJ‑0478:712:GW2002102894th:2:1107:9001:44351/1‑150 (MQ=60)
                                                                                                            GTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCtgtctcttatacacatctccgagcccacgagacg               >  GWNJ‑0478:712:GW2002102894th:2:2102:3050:49267/1‑116 (MQ=60)
                                                                                                                         CGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGCCGCAGCTTCAATCC  <  GWNJ‑0478:712:GW2002102894th:2:2206:14317:82281/150‑1 (MQ=60)
                                                                                                                                     |                                                                                                                                         
CCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGTGGTCCGCAGATTGCCCGCGTGATTGCCCGTGCGCCCATTGGCCACAGTTTGGGGCGGCCGCCGCAGCTTCAATCC  >  AM260479/3529476‑3529744

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