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 I187 R2 6 23.1 1435966 95.1% 1365603 84.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NZ_CP009273 2,399,465 Δ7 bp coding (650‑656/939 nt) lrhA ← transcriptional regulator LrhA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NZ_CP009273 2399465 2399471 7 22 [0] [0] 22 lrhA transcriptional regulator LrhA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 23994640 (0.000)22 (0.900) 19/166 0.2 100% coding (657/939 nt) lrhA transcriptional regulator LrhA
?NZ_CP009273 2399472 = 0 (0.000)coding (649/939 nt) lrhA transcriptional regulator LrhA

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGG  >  NZ_CP009273/2399386‑2399560
                                                                               |                                                                                               
CGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATAT                                                                      <  SRR3722063.110932/100‑1 (MQ=60)
      CATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTT                                                                >  SRR3722063.443728/1‑100 (MQ=60)
           TCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTA                                                           <  SRR3722063.687909/100‑1 (MQ=60)
               AACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCG                                                       <  SRR3722063.412776/100‑1 (MQ=60)
                   GGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGC                                                   >  SRR3722063.161420/1‑100 (MQ=60)
                         agagacagGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAC                                             <  SRR3722063.88238/92‑1 (MQ=60)
                                          GGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGC                            >  SRR3722063.591970/1‑100 (MQ=60)
                                               TCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGA                       <  SRR3722063.488704/100‑1 (MQ=60)
                                                CACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT                      <  SRR3722063.171870/100‑1 (MQ=60)
                                                             ACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATAC         <  SRR3722063.335470/100‑1 (MQ=60)
                                                                    GAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCGGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGG  <  SRR3722063.15056/100‑1 (MQ=60)
                                                                               |                                                                                               
CGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGG  >  NZ_CP009273/2399386‑2399560

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