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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,176,805 T→G D109A (GAC→GCC)  yegX ← glycoside hydrolase family 25 protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,176,8050TG100.0% 72.2 / NA 22D109A (GAC→GCC) yegXglycoside hydrolase family 25 protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/10);  total (12/10)

GGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTG  >  NZ_CP009273/2176720‑2176887
                                                                                     |                                                                                  
ggATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGGTGCCAGTTACGCGAAAAATAGGGGGccac                                                                                >  1:796349/1‑90 (MQ=255)
  aTACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGccacca                                                                              >  2:1126083/1‑90 (MQ=255)
     ccGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCt                                                                           >  2:205337/1‑90 (MQ=255)
           gggAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGc                                                                     >  2:1131148/1‑90 (MQ=255)
            ggAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGcc                                                                    <  1:347543/90‑1 (MQ=255)
               gaAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTc                                                                 >  2:418355/1‑90 (MQ=255)
               gaAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCACCTTc                                                                 >  1:1232067/1‑90 (MQ=255)
                       gATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCt                                                         <  1:1131148/90‑1 (MQ=255)
                          aCGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTtaa                                                      >  1:742920/1‑90 (MQ=255)
                           cGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTtaat                                                     <  2:1227506/90‑1 (MQ=255)
                             ccccACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTtaataa                                                   <  2:270346/90‑1 (MQ=255)
                             ccccACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTtaataa                                                   >  2:290752/1‑90 (MQ=255)
                             ccccACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTtaataa                                                   <  2:742920/90‑1 (MQ=255)
                                 aCGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGc                                               <  1:1011835/90‑1 (MQ=255)
                                      gcagGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACt                                          >  1:371810/1‑90 (MQ=255)
                                         ggCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTa                                       <  2:708438/90‑1 (MQ=255)
                                              ttttCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCgg                                  <  1:447449/90‑1 (MQ=255)
                                              ttttCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCgg                                  >  2:1234369/1‑90 (MQ=255)
                                                     ggCTTAGTTGCCAGTTACGGGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCAt                           <  2:59609/90‑1 (MQ=255)
                                                                cAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCAtt                >  1:678037/1‑90 (MQ=255)
                                                                  gTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCAtttt              >  2:444338/1‑90 (MQ=255)
                                                                              aTAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTg  <  1:244192/90‑1 (MQ=255)
                                                                                     |                                                                                  
GGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTG  >  NZ_CP009273/2176720‑2176887

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATC  >  NZ_CP009273/2176707‑2176897
                                                                                                  |                                                                                            
CCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGC                                                                                             <  SRR3722089.636239/100‑1 (MQ=60)
   GAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGGTGCCAGTTACGCGAAAAATAGGGGGCCAC                                                                                          >  SRR3722089.807329/1‑100 (MQ=60)
       TGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGC                                                                                      <  SRR3722089.287761/100‑1 (MQ=60)
                  CCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCACCTTC                                                                           >  SRR3722089.1250798/1‑100 (MQ=60)
                         GGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCC                                                                    <  SRR3722089.351863/100‑1 (MQ=60)
                             AAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAA                                                                >  SRR3722089.752687/1‑100 (MQ=60)
                                    GATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGC                                                         <  SRR3722089.1148350/100‑1 (MQ=60)
                                         GCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACT                                                    >  SRR3722089.376372/1‑100 (MQ=60)
                                              ACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAG                                               <  SRR3722089.1026906/100‑1 (MQ=60)
                                                           TTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGT                                  <  SRR3722089.452814/100‑1 (MQ=60)
                                                                   GCTTAGTTGCCAGTTACGCGAAAAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATT                          >  SRR3722089.686833/1‑100 (MQ=60)
                                                                                         AAATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGA    >  SRR3722089.953488/1‑100 (MQ=60)
                                                                                           ATAGGGGGCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATC  <  SRR3722089.247302/100‑1 (MQ=60)
                                                                                                  |                                                                                            
CCTGAACTGAAGCGGATACCGACGGGGAGAAATAATGATACGCCCCACGCAGCAGGCCATTTTCGCGGCTTAGTTGCCAGTTACGCGAAAAATAGGGGTCCACCAGCTTTTCGCCTTCCGTCGCCTTAATAAAAGCAAACTGTAAGCGGATACCATTGTCGCGCATTTTTGCCACCCGCTGCCAGTCGATC  >  NZ_CP009273/2176707‑2176897

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