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 I211 R1 223 18.0 1007644 96.3% 970361 85.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,220,044 A→G S236P (TCG→CCG)  yohF ← SDR family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,220,0440AG100.0% 80.8 / NA 25S236P (TCG→CCG) yohFSDR family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/6);  total (19/6)

AATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG  >  NZ_CP009273/2219961‑2220129
                                                                                   |                                                                                     
aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 >  2:284121/1‑90 (MQ=255)
 aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCg                                                                                >  1:93367/1‑90 (MQ=255)
 aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCg                                                                                >  1:401822/1‑90 (MQ=255)
  tGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               <  2:93367/90‑1 (MQ=255)
  tGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               <  1:93919/90‑1 (MQ=255)
                          gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                                       <  2:270715/90‑1 (MQ=255)
                                  gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               >  2:251310/1‑90 (MQ=255)
                                  gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               >  1:165438/1‑90 (MQ=255)
                                  gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGGGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               >  1:113234/1‑90 (MQ=255)
                                       gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:91933/1‑90 (MQ=255)
                                       gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:452865/1‑90 (MQ=255)
                                       gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:329074/1‑90 (MQ=255)
                                       gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:270743/1‑90 (MQ=255)
                                            gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                                     >  1:6510/1‑90 (MQ=255)
                                            gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                                     >  2:462810/1‑90 (MQ=255)
                                                  tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa                               <  1:462810/90‑1 (MQ=255)
                                                  tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa                               <  2:408059/90‑1 (MQ=255)
                                                  tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa                               <  2:452865/90‑1 (MQ=255)
                                                   cGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAt                              >  1:476953/1‑90 (MQ=255)
                                                                  gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc               >  2:364876/1‑90 (MQ=255)
                                                                   ccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCa              >  1:452790/1‑90 (MQ=255)
                                                                            aTCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAAcgccgc     >  2:211658/1‑90 (MQ=255)
                                                                               aaCGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAg  >  1:256123/1‑90 (MQ=255)
                                                                               aaCGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAg  >  1:135003/1‑90 (MQ=255)
                                                                               aaCGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAg  >  2:7928/1‑90 (MQ=255)
                                                                                   |                                                                                     
AATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG  >  NZ_CP009273/2219961‑2220129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
ATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGGAATCG  >  NZ_CP009273/2219952‑2220136
                                                                                            |                                                                                            
ATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCG                                                                                       >  SRR3722094.94722/1‑100 (MQ=60)
ATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCG                                                                                       >  SRR3722094.408765/1‑100 (MQ=60)
           TGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATT                                                                            <  SRR3722094.95284/100‑1 (MQ=60)
                         ccaccgagatctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAC                                                              <  SRR3722094.152387/43‑1 (MQ=60)
                                 GCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAC                                                      >  SRR3722094.167750/1‑100 (MQ=60)
                                 GCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGGGGTGTAATTTGCGCCCTCCGAACAAAGCCAC                                                      >  SRR3722094.114871/1‑100 (MQ=60)
                                      ATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCAC                                                 >  SRR3722094.274779/1‑100 (MQ=60)
                                      ATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCAC                                                 >  SRR3722094.93271/1‑100 (MQ=60)
                                      ATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCAC                                                 >  SRR3722094.460832/1‑100 (MQ=60)
                                           GGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGC                                            >  SRR3722094.6613/1‑100 (MQ=60)
                                                  ACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAT                                     >  SRR3722094.485380/1‑100 (MQ=60)
                                                           TCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGT                            <  SRR3722094.470951/100‑1 (MQ=60)
                                                                  CATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCA                     >  SRR3722094.460755/1‑100 (MQ=60)
                                                                              ATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG         >  SRR3722094.259901/1‑100 (MQ=60)
                                                                              ATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG         >  SRR3722094.136933/1‑100 (MQ=60)
                                                                                   CTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGGAAT    >  SRR3722094.81324/1‑100 (MQ=60)
                                                                                     ATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACActgtctcttatacacatctgacgctgccgacgaccatcttagtgtagatctcggtggtc  >  SRR3722094.390197/1‑41 (MQ=60)
                                                                                     ATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGGAATCG  >  SRR3722094.460063/1‑100 (MQ=60)
                                                                                            |                                                                                            
ATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGGAATCG  >  NZ_CP009273/2219952‑2220136

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: