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 AM260479 3,248,064 +A intergenic (+15/‑203) 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,0631.A100.0% 51.4 / NA 15intergenic (+14/‑204)h16_A3001/hemK1transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (5/10);  total (5/10)

TCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGG  >  AM260479/3247925‑3248188
                                                                                                                                           |                                                                                                                             
tCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTc                                                                                                                                <  1:1093797/134‑1 (MQ=255)
                     ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc                                                                                                           <  2:2694/134‑1 (MQ=255)
                     ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcg                                                                                                            >  2:1125456/1‑133 (MQ=255)
                          ggCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCtt                                                                                                      <  1:310863/134‑1 (MQ=255)
                             gTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc                                                                                                           <  2:253968/126‑1 (MQ=255)
                                                                                 ctGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGt                                                   <  2:184044/133‑1 (MQ=255)
                                                                                        gCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAg                                            >  2:990070/1‑134 (MQ=255)
                                                                                                       aCCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGaa                             <  1:1125456/134‑1 (MQ=255)
                                                                                                        ccTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAac                            <  1:990070/134‑1 (MQ=255)
                                                                                                             ttGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTAcc                                                                                                                <  1:1096074/45‑1 (MQ=25)
                                                                                                             ttGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTAcc                                                                                                                >  2:1096074/1‑45 (MQ=25)
                                                                                                                           gACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGc                                                                                  <  1:863610/62‑1 (MQ=255)
                                                                                                                           gACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGc                                                                                  >  2:863610/1‑62 (MQ=255)
                                                                                                                                     gCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACgg  <  1:279801/132‑1 (MQ=255)
                                                                                                                                      ccGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                    >  1:398719/1‑113 (MQ=255)
                                                                                                                                      ccGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                    <  2:398719/113‑1 (MQ=255)
                                                                                                                                           |                                                                                                                             
TCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGG  >  AM260479/3247925‑3248188

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
GACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTAT  >  AM260479/3247918‑3248206
                                                                                                                                                   |                                                                                                                                              
gtataagagacagCGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGG                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2106:13122:77887/136‑1 (MQ=60)
      ATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTAC                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2213:4625:70100/150‑1 (MQ=60)
                gagattagcgagtgtctcgtgggctcggagatgtgtataagagacagctatgatgccgagcgcgaccgcaacctgcgCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGC                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1211:15374:9501/73‑1 (MQ=60)
                     GAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTC                                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1114:10031:30035/1‑150 (MQ=60)
                           ttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCT                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2213:11749:81358/75‑1 (MQ=60)
                                GGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCAT                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1203:17456:52671/150‑1 (MQ=60)
                                                        agacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGC                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2116:11606:67819/92‑1 (MQ=60)
                                                                            CGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAagcctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaa                                                                >  GWNJ‑0478:712:GW2002102894th:2:1201:12998:34342/1‑68 (MQ=60)
                                                                              acggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGC                                                                <  GWNJ‑0478:712:GW2002102894th:2:2207:1497:71774/95‑1 (MQ=60)
                                                                                        cggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaggtcaAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTC                                                      <  GWNJ‑0478:712:GW2002102894th:2:2208:4452:45185/91‑1 (MQ=60)
                                                                                                             CACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAA                                <  GWNJ‑0478:712:GW2002102894th:2:2215:4903:23562/150‑1 (MQ=60)
                                                                                                              ACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAAC                               <  GWNJ‑0478:712:GW2002102894th:2:2207:20365:22187/150‑1 (MQ=60)
                                                                                                                           CAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaca                   >  GWNJ‑0478:712:GW2002102894th:2:2113:14270:70304/1‑68 (MQ=60)
                                                                                                                              ACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACtgtctc                >  GWNJ‑0478:712:GW2002102894th:2:1206:19708:72954/1‑144 (MQ=60)
                                                                                                                                     GCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCct         >  GWNJ‑0478:712:GW2002102894th:2:1113:5869:64869/1‑148 (MQ=60)
                                                                                                                                     GCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTC         >  GWNJ‑0478:712:GW2002102894th:2:1115:9414:65303/1‑150 (MQ=60)
                                                                                                                                         agacaggtcaAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCC     <  GWNJ‑0478:712:GW2002102894th:2:1215:13422:86115/140‑1 (MQ=60)
                                                                                                                                           GGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCC     <  GWNJ‑0478:712:GW2002102894th:2:1116:16323:59510/148‑1 (MQ=60)
                                                                                                                                            GCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaatagaaaaagaagaaacgagggtgg  >  GWNJ‑0478:712:GW2002102894th:2:2211:16264:88522/1‑42 (MQ=52)
                                                                                                                                                   |                                                                                                                                              
GACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCCTAT  >  AM260479/3247918‑3248206

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