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 F2 I222 R1 6 25.9 1502420 97.3% 1461854 86.7

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 28 [0] [0] 28 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)28 (1.040) 15/172 0.4 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
CATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGA  >  NZ_CP009273/2399395‑2399561
                                                                      |                                                                                                
CATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTAT                                                              >  SRR3722106.543853/1‑100 (MQ=60)
                GCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAC                                              <  SRR3722106.337485/100‑1 (MQ=60)
                GCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAC                                              <  SRR3722106.342856/100‑1 (MQ=60)
                     CACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATAT                                         <  SRR3722106.456562/100‑1 (MQ=60)
                       CGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCG                                       >  SRR3722106.646529/1‑100 (MQ=60)
                                                       GCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAG       <  SRR3722106.677093/100‑1 (MQ=60)
                                                         CGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG     <  SRR3722106.308884/100‑1 (MQ=60)
                                                         CGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG     <  SRR3722106.64366/100‑1 (MQ=60)
                                                         CGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG     <  SRR3722106.729665/100‑1 (MQ=60)
                                                            AAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGA  <  SRR3722106.515342/100‑1 (MQ=60)
                                                                      |                                                                                                
CATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAGGGA  >  NZ_CP009273/2399395‑2399561

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: ‑