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 146 27.9 2327082 89.7% 2087392 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 919,006 T→G C159G (TGC→GGC)  h16_B0812 → hydrolase or acyltransferase (alpha/beta hydrolase superfamily)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260480919,0060TG100.0% 81.7 / NA 22C159G (TGC→GGC) h16_B0812hydrolase or acyltransferase (alpha/beta hydrolase superfamily)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/12);  total (10/12)

CGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCG  >  AM260480/918897‑919133
                                                                                                             |                                                                                                                               
cggcctggcCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGcgcc                                                                                                         <  2:623537/134‑1 (MQ=255)
        ccAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGc                                                                                                 >  1:708108/1‑134 (MQ=255)
                      gTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    <  2:1129961/133‑1 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         <  2:1028328/92‑1 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         >  1:1028328/1‑92 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     <  2:759269/112‑1 (MQ=255)
                          aCCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     >  1:759269/1‑112 (MQ=255)
                               ccgGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGct                                                                          >  1:925063/1‑134 (MQ=255)
                               ccgGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGct                                                                          >  2:377560/1‑134 (MQ=255)
                                              gggCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGTCCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGa                                                           <  1:391004/134‑1 (MQ=255)
                                                 cgccgcGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATg                                                        >  1:607232/1‑134 (MQ=255)
                                                  gccgcGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATgg                                                       <  1:817894/134‑1 (MQ=255)
                                                                         ccGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         >  1:586183/1‑45 (MQ=255)
                                                                         ccGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGAc                                                                                                                         <  2:586183/45‑1 (MQ=255)
                                                                          cGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     <  2:1117902/64‑1 (MQ=255)
                                                                          cGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCtt                                                                                                     >  1:1117902/1‑64 (MQ=255)
                                                                              gTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCgtgg                           <  2:708108/134‑1 (MQ=255)
                                                                                              cGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAgg            <  1:377560/133‑1 (MQ=255)
                                                                                               gCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    <  1:712620/60‑1 (MQ=255)
                                                                                               gCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACcg                                                                                    >  2:712620/1‑60 (MQ=255)
                                                                                                ccTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGAtt         >  1:336415/1‑134 (MQ=255)
                                                                                                       cgcCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATcg  <  1:971622/134‑1 (MQ=255)
                                                                                                             |                                                                                                                               
CGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCG  >  AM260480/918897‑919133

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCGCGCCGCGGCGGGCCGGCCTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCGCGACGACATCCATGC  >  AM260480/918873‑919148
                                                                                                                                     |                                                                                                                                              
GCGCGCCGCGGCGGGCCGGCCTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2102:15144:21483/150‑1 (MQ=60)
                 GGCCTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCC                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:2108:17640:64500/1‑150 (MQ=60)
                              GGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACtgtctcttatacacatctccgagcccacgagacactcg                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1116:9675:38115/1‑112 (MQ=60)
                                   GACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCtgtctcttatacacatctccgagccc                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:2211:8000:34728/1‑123 (MQ=60)
                                   GACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCtgtctct                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2110:2472:41415/1‑143 (MQ=60)
                                       GGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCtgtctcttatacacatctccgagcccacgagacactcg                                                                                         >  GWNJ‑0478:712:GW2002102894th:2:1206:2356:91114/1‑112 (MQ=60)
                                        GTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTG                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2204:13076:45057/1‑150 (MQ=60)
                                            gaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCC                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2108:6937:89145/90‑1 (MQ=60)
                                                          GTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGG                                                                      >  GWNJ‑0478:712:GW2002102894th:2:2103:20746:54567/1‑150 (MQ=60)
                                                                     CGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGTCCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTG                                                           <  GWNJ‑0478:712:GW2002102894th:2:1206:7059:28985/150‑1 (MQ=60)
                                                                         CGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAAC                                                       <  GWNJ‑0478:712:GW2002102894th:2:2115:17137:43067/150‑1 (MQ=60)
                                                                                  TATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaa                                               >  GWNJ‑0478:712:GW2002102894th:2:2101:10071:21504/1‑76 (MQ=60)
                                                                                   ATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgt                                             >  GWNJ‑0478:712:GW2002102894th:2:2214:9054:86766/1‑95 (MQ=60)
                                                                                                         TGTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTT                       >  GWNJ‑0478:712:GW2002102894th:2:1205:19241:10355/1‑150 (MQ=60)
                                                                                                          GTGGATTCCATCCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCActgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaacaaaaaaaaaaaaacgaaaata                       >  GWNJ‑0478:712:GW2002102894th:2:2114:9193:96485/1‑46 (MQ=60)
                                                                                                                     CCGCCTGTCGCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCGCGACG            <  GWNJ‑0478:712:GW2002102894th:2:1207:14944:52474/149‑1 (MQ=60)
                                                                                                                              GCGCCAGGGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCGCGACGACATCCATGC  <  GWNJ‑0478:712:GW2002102894th:2:2206:2898:13503/150‑1 (MQ=60)
                                                                                                                                     |                                                                                                                                              
GCGCGCCGCGGCGGGCCGGCCTGGCGGCCTGGCCAGACTGGTGCTGGTCGACCCGCCGGTGTCGGGGCCCGGGCGCCGCGCGTATCCTGCGCAACTGCCGTGGTATGTGGATTCCATCCGCCTGTCGCGCCAGTGCACCGACGTGGAAGGCATGCGCCGCTTCTGCCCGACCTGGACCGACGAGCAGCTGCGCCTGCGCGCCGAATGGCTGCATACCTGCAACGAAGACGCCGTGGTGCGCAGCTTTGAGGATTTCCATCGCGACGACATCCATGC  >  AM260480/918873‑919148

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: