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 154 26.6 2069536 90.2% 1866721 108.9

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% 87.3 / NA 20coding (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 . (9/11);  total (9/11)

CGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGC  >  AM260479/2955063‑2955304
                                                                                                                                    |                                                                                                             
cGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAc                                                                                                               <  1:80804/133‑1 (MQ=255)
  aGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAccag                                                                                                            <  1:894247/134‑3 (MQ=255)
       gAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAccag                                                                                                            <  2:473288/129‑3 (MQ=255)
       gAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAccag                                                                                                            >  1:473288/1‑127 (MQ=255)
                   gCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACgg                                                                                          <  1:148815/134‑1 (MQ=255)
                   gCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACgg                                                                                          <  1:220407/134‑1 (MQ=255)
                           cgccgcCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATgcag                                                                                   >  2:64958/1‑133 (MQ=255)
                             ccgccTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACgg                                                                                          >  1:1005902/1‑124 (MQ=255)
                             ccgccTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACgg                                                                                          <  2:1005902/124‑1 (MQ=255)
                                                 aTCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATcgccg                                                            <  2:322863/134‑1 (MQ=255)
                                                   cGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCcgcg                                                          >  1:246085/1‑134 (MQ=255)
                                                              agaTCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTAc                                               >  1:322413/1‑134 (MQ=255)
                                                                   gCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCg                                           <  2:246085/133‑1 (MQ=255)
                                                                        aCGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGt                                     >  1:447715/1‑134 (MQ=255)
                                                                          gCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTcc                                   <  2:378306/134‑1 (MQ=255)
                                                                            aTGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCcgg                                 <  2:579459/134‑1 (MQ=255)
                                                                               gATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGt                                     >  1:363462/1‑127 (MQ=255)
                                                                               gATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGt                                     <  2:363462/127‑1 (MQ=255)
                                                                                       ccggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGc                      <  1:1001332/134‑1 (MQ=255)
                                                                                        cggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAccag                                                                                                            >  2:321270/1‑46 (MQ=255)
                                                                                        cggcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTAccag                                                                                                            <  1:321270/48‑3 (MQ=255)
                                                                                          gcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                                                  >  1:739222/1‑87 (MQ=255)
                                                                                          gcTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTa                                                                  <  2:739222/87‑1 (MQ=255)
                                                                                                             ggCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAacgcacgc  >  1:204868/1‑132 (MQ=255)
                                                                                                                aCCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCt                                                  >  1:684567/1‑81 (MQ=255)
                                                                                                                aCCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCt                                                  <  2:684567/81‑1 (MQ=255)
                                                                                                                                    |                                                                                                             
CGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGC  >  AM260479/2955063‑2955304

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
AATCGCGCTGTTCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGCCCCTGGCGGTCCGTTT  >  AM260479/2955051‑2955333
                                                                                                                                                  |                                                                                                                                        
AATCGCGCTGTTCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTC                                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1212:8093:55313/150‑1 (MQ=60)
    GCGCTGTTCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCG                                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1215:6879:51028/1‑150 (MQ=60)
           TCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAA                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1107:3169:53861/149‑1 (MQ=60)
             GAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGAC                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2207:18478:96527/150‑1 (MQ=60)
              tttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGAC                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1204:9031:54032/78‑1 (MQ=60)
                          TCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATT                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2215:4993:42750/1‑150 (MQ=60)
                              GGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCA                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1111:16284:77890/150‑1 (MQ=60)
                              GGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCA                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1115:17038:93328/150‑1 (MQ=60)
                                                CATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGC                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1116:19754:40590/1‑150 (MQ=60)
                                                           GCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACG                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1204:7735:61490/1‑150 (MQ=60)
                                                                     CGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTC                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1212:13370:12529/1‑150 (MQ=60)
                                                                            ATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGG                                                           >  GWNJ‑0478:712:GW2002102894th:2:1208:4388:23808/1‑149 (MQ=60)
                                                                                       CATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCctgtctcttatacacatctccgagcccacg                                                  >  GWNJ‑0478:712:GW2002102894th:2:2114:17391:4883/1‑117 (MQ=60)
                                                                                                  GCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACG                                    <  GWNJ‑0478:712:GW2002102894th:2:2215:12844:17929/150‑1 (MQ=60)
                                                                                                          GCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCA                              >  GWNJ‑0478:712:GW2002102894th:2:1115:11632:6937/1‑148 (MQ=60)
                                                                                                             cggagatgtgtataagagtcaGAAGTCGGGCCGTACCATTCGCCGATGGAAGACGGCATGCAGCAGTTCCATGAGTGGTATCGCCGCGCCATGGCTTACTCCGGCTAGTCCTGCGGACTAGTGCTATAACCGCACAACGCACGCAAGAAC                         <  GWNJ‑0478:712:GW2002102894th:2:2114:13674:31169/129‑1 (MQ=60)
                                                                                                             CTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCctgtctcttatacacatctccgagcccacgagacggagc                         >  GWNJ‑0478:712:GW2002102894th:2:2110:12476:91698/1‑111 (MQ=60)
                                                                                                                               AGCGAAGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGctgtctcttatacacatctcc       >  GWNJ‑0478:712:GW2002102894th:2:2210:8579:8290/1‑129 (MQ=60)
                                                                                                                                    AGTCGGGCCGTACCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGCCCCTGGCGGTCCGTTT  >  GWNJ‑0478:712:GW2002102894th:2:1116:4862:90182/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                        
AATCGCGCTGTTCGAGCGCGAGTTCGTCGAGGCCGAGCGCGCCGCCTACATGGAAACCTGCATCGAGGACGACGAGATCGCCGAACGCATGGATGCCGGCCGGCTGGCGCTGCTCAAGCGCGGCACCAGCGAAGTCGGGCCGTACCCAGTCGCCGATGGAAGACGGCATGCAGCATTTCCATGAGTGGTATCGCCGCGCCATGGCTTACGCCGGCTAGTCCGGCGGACTAGTGCGATAACCGCACAACGCACGCAAGAACGGACCGCCCCTGGCGGTCCGTTT  >  AM260479/2955051‑2955333

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: