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 1,081,436 +C coding (437/447 nt) h16_A0989 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,081,4331.C100.0% 44.4 / NA 13G145A (GGC→GCC) h16_A0989Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (6/7);  total (6/7)

AGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCC  >  AM260479/1081304‑1081517
                                                                                                                                  |                                                                                    
aGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGcccc                                                                                  <  2:227010/134‑2 (MQ=255)
           aTTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGa                                                                         <  2:545326/133‑1 (MQ=255)
               cgcTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCATGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTAc                                                                    <  2:1109391/134‑1 (MQ=255)
                                  cgtcGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCgggg                                                  >  1:1057023/1‑133 (MQ=255)
                                  cgtcGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCgggg                                                  >  1:932751/1‑133 (MQ=255)
                                   gtcGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGcccc                                                                                  >  1:313020/1‑98 (MQ=255)
                                   gtcGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGcccc                                                                                  <  2:313020/99‑2 (MQ=255)
                                      gACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGa                                             <  2:1051891/134‑1 (MQ=255)
                                        cGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACg                                           >  2:859687/1‑134 (MQ=255)
                                                                                 gcaATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGcc  <  1:859687/134‑1 (MQ=255)
                                                                                         cgcCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAg                                               <  1:91397/81‑1 (MQ=255)
                                                                                         cgcCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAg                                               >  2:91397/1‑81 (MQ=255)
                                                                                                cccGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATg                                                                <  1:553382/57‑1 (MQ=39)
                                                                                                cccGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATg                                                                >  2:553382/1‑57 (MQ=39)
                                                                                                cccGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGcc            <  1:580893/109‑1 (MQ=255)
                                                                                                cccGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGcc            >  2:580893/1‑109 (MQ=255)
                                                                                                                                  |                                                                                    
AGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCC  >  AM260479/1081304‑1081517

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ 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
AACGGCAGCCGGCAGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCTTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGACGAGGGCCAGATC  >  AM260479/1081291‑1081580
                                                                                                                                                  |                                                                                                                                                
agcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGG                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2112:14846:98111/86‑1 (MQ=60)
       GCCGGCAGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGA                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1111:7929:95148/150‑1 (MQ=60)
       GCCGGCAGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGA                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2114:7913:9173/150‑1 (MQ=60)
                            CGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCctggacctgtctcttatacacatctccgagcc                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1101:9465:23459/1‑118 (MQ=60)
                               gcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGC                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1214:7870:99866/87‑1 (MQ=60)
                                GGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGC                                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2210:14318:82712/1‑149 (MQ=60)
                                GGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGC                                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2204:16674:92759/1‑149 (MQ=60)
                                 GCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACctgtctcttatacacatctcc                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1203:17406:83520/1‑129 (MQ=60)
                                                cgagtgtctcgtgggctcggagatgtgtataagagacaGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACC                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1107:19184:73795/112‑1 (MQ=60)
                                                                gccgaggcttcgacgcgggcatGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACctgtctcttatacacatctccgagccc                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2208:19826:61313/23‑122 (MQ=60)
                                                                                   ataagagacagGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGC                                                            <  GWNJ‑0478:712:GW2002102894th:2:2101:7175:41548/139‑1 (MQ=60)
                                                                                             AGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTG                                                  <  GWNJ‑0478:712:GW2002102894th:2:2116:15205:53309/150‑1 (MQ=60)
                                                                                                gagattagcgagtgtctcgtgggctcggagatgtgtataagagacagcCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTGACG                                               <  GWNJ‑0478:712:GW2002102894th:2:2101:21071:84593/102‑1 (MQ=60)
                                                                                                                                  GCGGCACGCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGGCG             >  GWNJ‑0478:712:GW2002102894th:2:1216:20560:98335/1‑150 (MQ=60)
                                                                                                                                         GCAAGGCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGACGAGGGCC       >  GWNJ‑0478:712:GW2002102894th:2:2208:12090:89388/1‑149 (MQ=60)
                                                                                                                                             GCCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaa  >  GWNJ‑0478:712:GW2002102894th:2:1206:2135:62725/1‑77 (MQ=60)
                                                                                                                                              CCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGACGAGGGCC       >  GWNJ‑0478:712:GW2002102894th:2:2211:10953:83368/1‑143 (MQ=60)
                                                                                                                                              CCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGACGAGGGCC       >  GWNJ‑0478:712:GW2002102894th:2:1114:17832:99045/1‑143 (MQ=60)
                                                                                                                                              CCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCCTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGACGAGGGCC       >  GWNJ‑0478:712:GW2002102894th:2:1107:17676:96318/1‑143 (MQ=60)
                                                                                                                                                  |                                                                                                                                                
AACGGCAGCCGGCAGGAGCTGGTGATTGCGCTGGCGCAGGCCTCGTCCGTCGACGTGCTGCGCTTCCGCGGCCTGTCCGATGCCGGGGCGCGCAGCAATGCGCGCCATGCCCGCGTCGAAGCCTCGGCCGGCGGCACGCAAGGCCCTGGACCGTGATCTACGATGCGGACCTGCCGGGGCAGGACGGCTACGATACCCGCGTGCCGCTGGCCCAGCCCGTCAAGGCCGGGGCCCTGCGCTTGACGCTTACCGGCAACCGCTGGCCGGGCGACCAGGACGAGGGCCAGATC  >  AM260479/1081291‑1081580

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