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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 635,449 T→C S81P (TCC→CCC)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273635,4490TC100.0% 83.4 / NA 24S81P (TCC→CCC) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/9);  total (15/9)

AGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAG  >  NZ_CP009273/635370‑635526
                                                                               |                                                                             
agaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGc                                                                    >  1:151691/1‑91 (MQ=255)
          ttGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCa                                                          >  1:672188/1‑91 (MQ=255)
          ttGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCa                                                          >  1:763070/1‑91 (MQ=255)
           tGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCAc                                                         >  2:199116/1‑91 (MQ=255)
            gCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCAcc                                                        >  2:901354/1‑91 (MQ=255)
               ggTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCg                                                     <  2:235237/91‑1 (MQ=255)
                          cGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGc                                          >  2:502121/1‑91 (MQ=255)
                            tCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGt                                        <  1:887589/91‑1 (MQ=255)
                            tCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGt                                        <  2:1084193/91‑1 (MQ=255)
                            tCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGt                                        >  2:342813/1‑91 (MQ=255)
                                  cTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCtgt                                  <  2:763070/91‑1 (MQ=255)
                                       cACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGAcc                             <  2:428709/91‑1 (MQ=255)
                                              ccAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTc                      >  2:405337/1‑91 (MQ=255)
                                                aGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCAt                    >  1:879058/1‑91 (MQ=255)
                                                aGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCAt                    >  2:343708/1‑91 (MQ=255)
                                                 ggTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATc                   <  1:839369/91‑1 (MQ=255)
                                                 ggTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATc                   <  2:829348/91‑1 (MQ=255)
                                                      cAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCg              >  2:146136/1‑91 (MQ=255)
                                                      cAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCg              >  1:725285/1‑91 (MQ=255)
                                                       aaCCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGa             >  1:482858/1‑91 (MQ=255)
                                                           aGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGAcc                             >  2:120957/1‑71 (MQ=255)
                                                           aGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGAcc                             <  1:120957/71‑1 (MQ=255)
                                                            ggggCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAaagaa        >  1:941222/1‑91 (MQ=255)
                                                                  aCGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAg  <  2:756236/91‑1 (MQ=255)
                                                                               |                                                                             
AGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAG  >  NZ_CP009273/635370‑635526

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCC  >  NZ_CP009273/635349‑635548
                                                                                                    |                                                                                                   
ctcggagatgtgtataagagacagTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCC                                                                                                     <  SRR3722125.440195/77‑1 (MQ=60)
         CACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGG                                                                                            <  SRR3722125.582438/101‑1 (MQ=60)
           CCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGC                                                                                          >  SRR3722125.154583/1‑101 (MQ=60)
                     AGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCA                                                                                >  SRR3722125.686660/1‑101 (MQ=60)
                     AGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCA                                                                                >  SRR3722125.779604/1‑101 (MQ=60)
                                                 TCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGAC                                                    <  SRR3722125.907232/101‑1 (MQ=60)
                                                           TCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCAT                                          >  SRR3722125.898502/1‑101 (MQ=60)
                                                            taagagacagGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATC                                         <  SRR3722125.123227/91‑1 (MQ=60)
                                                                 ACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCG                                    >  SRR3722125.740917/1‑101 (MQ=60)
                                                                  CCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGA                                   >  SRR3722125.493048/1‑101 (MQ=60)
                                                                      GGTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGA                               <  SRR3722125.857854/101‑1 (MQ=60)
                                                                       GTTCCAACCAGGGGCCACGTTTTGCAGGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAA                              >  SRR3722125.962138/1‑101 (MQ=60)
                                                                                                 GGCCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGCGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCG    >  SRR3722125.742726/1‑101 (MQ=60)
                                                                                                   CCCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCtgtctcttatacacatctccgagcccacgagactaggc  >  SRR3722125.263774/1‑63 (MQ=60)
                                                                                                    |                                                                                                   
AGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACGAGTTCACCTCGCTGGTACTGGCGTTGCTGTGGACCGGTGGTCATCCGTCGAAAGAAGCGCAGTCTCTGCTGGAGCAGATTCGCC  >  NZ_CP009273/635349‑635548

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