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 I189 R2 1 26.9 1622034 94.1% 1526333 83.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 37 [0] [0] 37 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)37 (1.370) 27/166 0.1 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
CGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG  >  NZ_CP009273/2399378‑2399558
                                                                                       |                                                                                             
CGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCA                                                                            >  SRR3722065.593350/1‑100 (MQ=60)
    GGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGA                                                                        >  SRR3722065.750280/1‑100 (MQ=60)
    GGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGA                                                                        >  SRR3722065.99833/1‑100 (MQ=60)
         GGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATC                                                                   >  SRR3722065.596705/1‑100 (MQ=60)
            CTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGC                                                                <  SRR3722065.420920/100‑1 (MQ=60)
                      CAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGC                                                      <  SRR3722065.87958/100‑1 (MQ=60)
                                        CGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCT                                    >  SRR3722065.624134/1‑100 (MQ=60)
                                                        CACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT                    <  SRR3722065.2890/100‑1 (MQ=60)
                                                        CACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT                    <  SRR3722065.79817/100‑1 (MQ=60)
                                                            GCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATC                >  SRR3722065.315933/1‑100 (MQ=60)
                                                                     ACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATAC       <  SRR3722065.627138/100‑1 (MQ=60)
                                                                      CGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACA      >  SRR3722065.471041/1‑100 (MQ=60)
                                                                       GGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAA     <  SRR3722065.680299/100‑1 (MQ=60)
                                                                        GCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAG    >  SRR3722065.550375/1‑100 (MQ=60)
                                                                          CGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG  <  SRR3722065.249001/100‑1 (MQ=60)
                                                                                       |                                                                                             
CGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG  >  NZ_CP009273/2399378‑2399558

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