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 3,248,092 G→C intergenic (+43/‑175) h16_A3001 → / → hemK1 transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,248,0920GC100.0% 70.5 / NA 19intergenic (+43/‑175)h16_A3001/hemK1transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (9/10);  total (9/10)

GCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAA  >  AM260479/3247966‑3248218
                                                                                                                               |                                                                                                                              
gCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCAt                                                                                                                     >  2:955808/1‑134 (MQ=255)
gCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCAt                                                                                                                     >  1:885104/1‑134 (MQ=255)
                            gcgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAc                                                                                            >  1:958612/5‑134 (MQ=255)
                             cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTtc                                                                                                           <  1:99826/115‑1 (MQ=255)
                             cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTtc                                                                                                           >  2:99826/4‑118 (MQ=255)
                             cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                                                           <  1:67979/131‑1 (MQ=255)
                             cgaccGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                                                           >  2:67979/4‑134 (MQ=255)
                                        ctgcgcGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGt                                                                                 <  2:291762/129‑1 (MQ=255)
                                           cgcGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGtt                                                                                <  2:958612/130‑1 (MQ=255)
                                               gCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGcc                                                                                                      >  2:158569/1‑106 (MQ=255)
                                               gCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGcc                                                                                                      <  1:158569/106‑1 (MQ=255)
                                                      ggccggccACCTGTTTGGCAGCAACGTGACGCACCTGGCTGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGt                                                                   <  1:176019/134‑1 (MQ=255)
                                                                                            gCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACgg                                <  1:955808/132‑1 (MQ=255)
                                                                                              cgtcaagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccg                            <  2:885104/129‑1 (MQ=255)
                                                                                               gtcaagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc                           >  2:752487/5‑134 (MQ=255)
                                                                                                  aagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgccg                         <  1:752487/132‑1 (MQ=255)
                                                                                                     cagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCctc                     <  1:634931/134‑1 (MQ=255)
                                                                                                     cagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCctc                     >  2:1154681/1‑134 (MQ=255)
                                                                                                                        gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCaaa  >  1:1158772/1‑134 (MQ=255)
                                                                                                                               |                                                                                                                              
GCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAA  >  AM260479/3247966‑3248218

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
TGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAACCCAACCATGCATCTG  >  AM260479/3247949‑3248234
                                                                                                                                                |                                                                                                                                              
acgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagaccgcaacctgcgCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2109:3365:10015/88‑1 (MQ=60)
  GGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATC                                                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:2201:4352:100797/1‑150 (MQ=60)
                          GCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCctgtctcttatacacatctccgagcccacg                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2105:6370:48885/1‑120 (MQ=60)
                           CCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACC                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2209:20058:43488/1‑150 (MQ=60)
                              agctatgatgccgagcgcgaccgcaacctgcgCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1106:10096:23413/118‑1 (MQ=60)
                                    gtataagagacagaccgcaacctgcgCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAA                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1111:3516:56722/123‑1 (MQ=60)
                                            GCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGC                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1104:8907:34282/150‑1 (MQ=60)
                                                                      CGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCTGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAG                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1110:2269:51037/150‑1 (MQ=60)
                                                                                                 GTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaacaaaacaaaataaaaaat                                           >  GWNJ‑0478:712:GW2002102894th:2:1102:7496:20375/1‑54 (MQ=60)
                                                                                                 GTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaagaaacacacaaaacc                                           >  GWNJ‑0478:712:GW2002102894th:2:2205:14696:27717/1‑54 (MQ=60)
                                                                                                            GGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCC                                 <  GWNJ‑0478:712:GW2002102894th:2:2209:11462:28439/148‑1 (MQ=60)
                                                                                                                  caAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAA                          <  GWNJ‑0478:712:GW2002102894th:2:2110:15613:35305/147‑1 (MQ=60)
                                                                                                                     GCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAG                      <  GWNJ‑0478:712:GW2002102894th:2:2105:8869:78213/150‑1 (MQ=60)
                                                                                                                          GTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAAC                 >  GWNJ‑0478:712:GW2002102894th:2:2216:7686:57489/1‑150 (MQ=60)
                                                                                                                               TACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGGCGGCCCGCCGCCTCTCCCTATATAAAAAGCAAACCCAAC            >  GWNJ‑0478:712:GW2002102894th:2:2108:9714:60952/1‑150 (MQ=60)
                                                                                                                                         GCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCtgtctcttatacacatctc  >  GWNJ‑0478:712:GW2002102894th:2:1216:9643:69441/1‑131 (MQ=60)
                                                                                                                                                |                                                                                                                                              
TGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTATATAAAAAGCAAACCCAACCATGCATCTG  >  AM260479/3247949‑3248234

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