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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

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,0440AG88.5% 60.2 / 2.9 26S236P (TCG→CCG) yohFSDR family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (11/12);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

CAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGT  >  NZ_CP009273/2219960‑2220110
                                                                                    |                                                                  
cAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCATCCACTATCAACGACTGccc                                                              <  2:342943/90‑1 (MQ=255)
 aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              <  1:735416/90‑1 (MQ=255)
 aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGccc                                                              <  2:410327/90‑1 (MQ=255)
  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              <  2:825627/89‑1 (MQ=255)
  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              >  1:825627/1‑89 (MQ=255)
  aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTCTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCg                                                             >  1:96775/1‑90 (MQ=255)
            cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                   >  2:658715/1‑90 (MQ=255)
            cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                   >  2:201853/1‑90 (MQ=255)
            cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                   >  2:278530/1‑90 (MQ=255)
              cGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAttt                                                 >  2:500151/1‑90 (MQ=255)
                 aCCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcg                                              <  2:540699/90‑1 (MQ=255)
                 aCCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcg                                              <  1:22006/90‑1 (MQ=255)
                  ccTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcgc                                             <  1:125013/90‑1 (MQ=255)
                      aCTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTc                                         <  1:419369/90‑1 (MQ=255)
                           gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                    <  2:182510/90‑1 (MQ=255)
                           gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                    <  2:240506/90‑1 (MQ=255)
                                   gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                            <  1:201853/90‑1 (MQ=255)
                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              >  1:275415/1‑51 (MQ=255)
                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                              <  2:275415/51‑1 (MQ=255)
                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                       >  2:780305/1‑90 (MQ=255)
                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                       >  1:265085/1‑90 (MQ=255)
                                        gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                       >  1:391348/1‑90 (MQ=255)
                                             gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                  >  1:701541/1‑90 (MQ=255)
                                                  ttCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa             >  1:182375/1‑90 (MQ=255)
                                                   tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa            <  1:657703/90‑1 (MQ=255)
                                                             aaaGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGt  <  1:249057/90‑1 (MQ=255)
                                                                                    |                                                                  
CAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGT  >  NZ_CP009273/2219960‑2220110

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
ACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGG  >  NZ_CP009273/2219946‑2220131
                                                                                                  |                                                                                       
ACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGC                                                                                        <  SRR3722079.533907/100‑1 (MQ=60)
      ATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCG                                                                                  >  SRR3722079.834464/1‑100 (MQ=60)
      ATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTCTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCG                                                                                  >  SRR3722079.97537/1‑100 (MQ=60)
        ATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGT                                                                                <  SRR3722079.420603/100‑1 (MQ=60)
               AATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAA                                                                         <  SRR3722079.743280/100‑1 (MQ=60)
                               ACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAC                                                         <  SRR3722079.22158/100‑1 (MQ=60)
                                CCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACA                                                        <  SRR3722079.126035/100‑1 (MQ=60)
                                    ACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGC                                                    <  SRR3722079.422980/100‑1 (MQ=60)
                                            ATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAActgtctcttatacacatctccgagcccac                                            >  SRR3722079.277716/1‑71 (MQ=60)
                                            ATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCAC                                            >  SRR3722079.267306/1‑100 (MQ=60)
                                            ATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCAC                                            >  SRR3722079.394604/1‑100 (MQ=60)
                                                 GGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGC                                       >  SRR3722079.708970/1‑100 (MQ=60)
                                                 GGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGC                                       <  SRR3722079.203595/100‑1 (MQ=60)
                                                      GAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCA                                  >  SRR3722079.183970/1‑100 (MQ=60)
                                                                 TCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGT                       <  SRR3722079.664527/100‑1 (MQ=60)
                                                                           AAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAA             <  SRR3722079.251144/100‑1 (MQ=60)
                                                                                      CCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGG  >  SRR3722079.608719/1‑100 (MQ=60)
                                                                                                  |                                                                                       
ACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAGGG  >  NZ_CP009273/2219946‑2220131

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