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 2,955,197 Δ1 bp coding (1034/1071 nt) h16_A2735 → phenylpropionate dioxygenase or related ring‑hydroxylating dioxygenase, large terminal subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,955,1950C.100.0% 132.7 / NA 30coding (1032/1071 nt)h16_A2735phenylpropionate dioxygenase or related ring‑hydroxylating dioxygenase, large terminal subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/17);  total (13/17)

GAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGCCCCTG  >  AM260479/2955064‑2955322
                                                                                                                                   |                                                                                                                               
gAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAcc                                                                                                                                <  2:308022/133‑1 (MQ=255)
         ttCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGAt                                                                                                                     <  1:974292/134‑1 (MQ=255)
           cgtcgAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATg                                                                                                                    >  2:1000645/1‑133 (MQ=255)
                  gCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACgg                                                                                                            <  1:1000645/134‑1 (MQ=255)
                     gAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCAt                                                                                                         <  1:80716/134‑1 (MQ=255)
                       gcgcgcCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgc                                                                                                       <  1:313224/134‑1 (MQ=255)
                       gcgcgcCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgc                                                                                                       <  2:85764/134‑1 (MQ=255)
                                              gCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgc                                                                                >  1:957584/1‑134 (MQ=255)
                                                aTCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccg                                                                              >  1:804563/1‑134 (MQ=255)
                                                aTCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccg                                                                              <  1:254329/134‑1 (MQ=255)
                                                 tCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccgc                                                                             >  2:362591/1‑134 (MQ=255)
                                                 tCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccgc                                                                             <  2:957584/134‑1 (MQ=255)
                                                  cGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                                                                    <  1:304287/126‑1 (MQ=255)
                                                  cGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                                                                    >  2:304287/1‑126 (MQ=255)
                                                               aTCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCt                                                                    >  2:521647/1‑129 (MQ=255)
                                                               aTCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCt                                                                    <  1:521647/129‑1 (MQ=255)
                                                                              gATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcggcgg                                                <  2:804563/134‑1 (MQ=255)
                                                                                     gccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgcag                                                                                                     >  1:184868/1‑74 (MQ=255)
                                                                                     gccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgcag                                                                                                     <  2:329744/74‑1 (MQ=255)
                                                                                     gccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgcag                                                                                                     <  2:184868/74‑1 (MQ=255)
                                                                                     gccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgcag                                                                                                     >  1:329744/1‑74 (MQ=255)
                                                                                      ccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGc                                        <  1:944180/134‑1 (MQ=255)
                                                                                                               aCCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAg                                                        <  2:593766/93‑1 (MQ=255)
                                                                                                               aCCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAg                                                        >  1:593766/1‑93 (MQ=255)
                                                                                                                 cAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcggcg                                                 >  2:365953/1‑98 (MQ=255)
                                                                                                                 cAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcggcg                                                 <  1:365953/98‑1 (MQ=255)
                                                                                                                 cAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAacg                         <  1:835812/122‑1 (MQ=255)
                                                                                                                 cAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAacg                         >  2:835812/1‑122 (MQ=255)
                                                                                                                 cAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACg             >  1:987518/1‑134 (MQ=255)
                                                                                                                 cAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACg             >  2:330720/1‑134 (MQ=255)
                                                                                                                            ggCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGCCCCTg  <  2:221283/134‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
GAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGCCCCTG  >  AM260479/2955064‑2955322

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 36 ≤ 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
TGTTCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGAC  >  AM260479/2955059‑2955314
                                                                                                                                          |                                                                                                                     
ccgtctcgtgggctcggagatgtgtataagagacaGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGAT                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1111:19273:47169/113‑1 (MQ=60)
             GTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCA                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2214:3641:2860/150‑1 (MQ=60)
                      GGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCA                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2215:19150:49317/150‑1 (MQ=60)
                         CGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTCCATGAG                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1107:12960:47263/150‑1 (MQ=60)
                           AGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGT                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1204:18726:12152/150‑1 (MQ=60)
                                    CCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCC                                                                       >  GWNJ‑0478:712:GW2002102894th:2:2213:9332:8645/1‑150 (MQ=60)
                                      TACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGC                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2202:16969:60724/1‑150 (MQ=60)
                                              AACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGC                                                             <  GWNJ‑0478:712:GW2002102894th:2:1205:20293:54798/150‑1 (MQ=60)
                                                    CATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGC                                                       <  GWNJ‑0478:712:GW2002102894th:2:1201:16634:13911/150‑1 (MQ=60)
                                                              GACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCC                                             <  GWNJ‑0478:712:GW2002102894th:2:2101:16945:56386/150‑1 (MQ=60)
                                                                           AACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGctgtctcttatacacatctccgagcccacgagacggagctacatct                                >  GWNJ‑0478:712:GW2002102894th:2:1207:2277:21071/1‑104 (MQ=60)
                                                                           AACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGctgtctcttatacacatctccgagcccacgagacggagctacatct                                >  GWNJ‑0478:712:GW2002102894th:2:1113:8761:90412/1‑104 (MQ=60)
                                                                                 cggagatgtgtataagagacagGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAAC                          <  GWNJ‑0478:712:GW2002102894th:2:1208:15903:48014/128‑1 (MQ=60)
                                                                                  catacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaggtaCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAAC                          <  GWNJ‑0478:712:GW2002102894th:2:1203:1980:93390/95‑1 (MQ=60)
                                                                                          GCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACG                 <  GWNJ‑0478:712:GW2002102894th:2:2212:2885:30982/150‑1 (MQ=60)
                                                                                                     CTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCtgtctcttatacacatctccgagccc      >  GWNJ‑0478:712:GW2002102894th:2:2107:12999:48592/1‑124 (MQ=60)
                                                                                                       GCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGG    >  GWNJ‑0478:712:GW2002102894th:2:2214:12447:76064/1‑150 (MQ=60)
                                                                                                          CAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGAC  <  GWNJ‑0478:712:GW2002102894th:2:2205:16925:65707/149‑1 (MQ=60)
                                                                                                                                          |                                                                                                                     
TGTTCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGAC  >  AM260479/2955059‑2955314

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