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
A4 F8 I1 R1 11 33.4 1306760 96.6% 1262330 122.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 757,731 Δ1 bp coding (43/717 nt) sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913757,7310G.100.0% 67.5 / NA 16coding (43/717 nt)sdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (7/9);  total (7/9)

CCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGC  >  NC_000913/757587‑757832
                                                                                                                                                |                                                                                                     
ccAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTtatgat                                                                                                  >  2:122319/1‑149 (MQ=255)
                           aGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTg                                                                       <  1:127105/149‑1 (MQ=255)
                                          ccGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATg                                                             <  1:103517/144‑1 (MQ=255)
                                          ccGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATg                                                             >  2:103517/1‑144 (MQ=255)
                                             aaaCTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTcgc                                                     <  1:122319/149‑1 (MQ=255)
                                                                 gccgAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgc                                 >  2:178878/1‑149 (MQ=255)
                                                                  ccgAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgcg                                >  7:40223/1‑149 (MQ=255)
                                                                             tacttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgcgc                               >  1:67242/1‑139 (MQ=255)
                                                                             tacttactAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgcgc                               <  2:67242/139‑1 (MQ=255)
                                                                                  actAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaaaa                 <  8:40223/148‑1 (MQ=255)
                                                                                    tAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgcgc                               >  2:188151/1‑132 (MQ=255)
                                                                                    tAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATgcgc                               <  1:188151/132‑1 (MQ=255)
                                                                                                   aaaaTGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCa                        <  5:65032/124‑1 (MQ=255)
                                                                                                   aaaaTGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCa                        >  6:65032/1‑124 (MQ=255)
                                                                                                             tCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGc  >  7:17481/1‑136 (MQ=255)
                                                                                                             tCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGc  <  8:17481/136‑1 (MQ=255)
                                                                                                                                          gATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCtt                             <  8:108185/80‑1 (MQ=255)
                                                                                                                                                |                                                                                                     
CCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGC  >  NC_000913/757587‑757832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCC  >  NC_000913/757613‑757874
                                                                                                                      |                                                                                                                                               
AAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGG                                                                                                                <  NB501016:228:HFK2FAFXY:1:21112:8300:18201/151‑1 (MQ=60)
               ACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGA                                                                                                      <  NB501016:228:HFK2FAFXY:1:21103:17560:18485/146‑1 (MQ=60)
                  GAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCG                                                                                              <  NB501016:228:HFK2FAFXY:1:21110:13384:20159/151‑1 (MQ=60)
                                     CCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACAC                                                                                                                     <  NB501016:228:HFK2FAFXY:4:11407:15918:17149/109‑1 (MQ=60)
                                       GCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGC                                                                         >  NB501016:228:HFK2FAFXY:4:11508:6565:1427/1‑151 (MQ=60)
                                                  GTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCT                                                                        >  NB501016:228:HFK2FAFXY:1:11302:19781:5396/1‑141 (MQ=60)
                                                         CTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCT                                                                        <  NB501016:228:HFK2FAFXY:1:21312:13076:8282/134‑1 (MQ=60)
                                                                        GAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAG                                                                 <  NB501016:228:HFK2FAFXY:3:11603:11194:15299/126‑1 (MQ=60)
                                                                                  CTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCC                                           >  NB501016:228:HFK2FAFXY:4:11409:23002:2730/1‑138 (MQ=60)
                                                                                                               GGATGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTAagatcggaagagcacacgtctgaacccagtcaccgctcattatctcgtatgccgtcttctgcttgaaa  >  NB501016:228:HFK2FAFXY:4:21505:4220:15823/1‑82 (MQ=60)
                                                                                                                  TGTTATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCC                               >  NB501016:228:HFK2FAFXY:4:11608:3370:13681/1‑118 (MQ=60)
                                                                                                                      |                                                                                                                                               
AAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAGTTTTCAATTTATCGCTATAACCCGGATGTTGATGATGCTCCGCGTATGCAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCC  >  NC_000913/757613‑757874

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