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 I203 R1 4 36.0 2027912 97.0% 1967074 86.4

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 32 [0] [0] 32 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)32 (0.890) 23/170 0.3 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
CCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG  >  NZ_CP009273/2399366‑2399558
                                                                                                   |                                                                                             
gagacagGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGA                                                                                               <  SRR3722085.643398/93‑1 (MQ=60)
                GGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGA                                                                        >  SRR3722085.597800/1‑100 (MQ=60)
                          CATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTT                                                              >  SRR3722085.892374/1‑100 (MQ=60)
                               TCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTA                                                         <  SRR3722085.1027123/100‑1 (MQ=60)
                               TCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTA                                                         <  SRR3722085.407622/100‑1 (MQ=60)
                               TCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAAAGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTA                                                         <  SRR3722085.1003840/100‑1 (MQ=60)
                                        GCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCC                                                >  SRR3722085.26013/1‑100 (MQ=60)
                                          CTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAG                                              >  SRR3722085.1025745/1‑100 (MQ=60)
                                             gatgtgtataagagacaGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAC                                           <  SRR3722085.143768/83‑1 (MQ=60)
                                             GCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCAC                                           <  SRR3722085.436238/100‑1 (MQ=60)
                                                  CACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATAT                                      <  SRR3722085.710155/100‑1 (MQ=60)
                                                    CGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGGCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCG                                    >  SRR3722085.983038/1‑100 (MQ=60)
                                                                    CACTGCCGCACGAACGGCCGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGAT                    <  SRR3722085.165710/100‑1 (MQ=60)
                                                                                      CGGAAGCGTCGAG‑‑‑‑‑‑‑AAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG  <  SRR3722085.571859/100‑1 (MQ=60)
                                                                                                   |                                                                                             
CCACTCAGCACACGCAGGTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCAATACAAGAG  >  NZ_CP009273/2399366‑2399558

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