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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,113,649 T→G K330N (AAA→AAC glcA ← glycolate permease GlcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,113,6490TG100.0% 82.4 / NA 26K330N (AAA→AACglcAglycolate permease GlcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/16);  total (10/16)

CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC  >  NZ_CP009273/3113564‑3113731
                                                                                     |                                                                                  
ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                >  1:50794/1‑90 (MQ=255)
ccGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTaa                                                                                <  1:537606/90‑1 (MQ=255)
  gCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAAc                                                                              <  2:330728/90‑1 (MQ=255)
     ttCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGc                                                                           >  2:102785/1‑90 (MQ=255)
                ttgATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCa                                                                <  2:408033/90‑1 (MQ=255)
                     gCAAATGAGGGATCTGGAAATTGATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                           <  1:106563/90‑1 (MQ=255)
                     gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                           <  1:309688/90‑1 (MQ=255)
                      cAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGt                                                          >  2:307692/1‑90 (MQ=255)
                          tGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGAcc                                                      <  1:364734/90‑1 (MQ=255)
                          tGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGAcc                                                      <  2:363500/90‑1 (MQ=255)
                             gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                   <  2:540169/90‑1 (MQ=255)
                                aTCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAcc                                                >  2:408575/1‑90 (MQ=255)
                                 tCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCg                                               >  1:135696/1‑90 (MQ=255)
                                       aaTTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGa                                         >  1:445915/1‑90 (MQ=255)
                                           aaTCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTaa                                     >  2:221676/1‑90 (MQ=255)
                                                   gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                             >  1:418772/1‑58 (MQ=255)
                                                   gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGAt                                                             <  2:418772/58‑1 (MQ=255)
                                                          aaaCGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGc                      >  2:406540/1‑90 (MQ=255)
                                                               cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                 <  1:49560/90‑1 (MQ=255)
                                                                    cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt            <  1:63278/90‑1 (MQ=255)
                                                                       ggAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGc         <  2:434484/90‑1 (MQ=255)
                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc       <  1:552482/90‑1 (MQ=255)
                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc       <  2:208025/90‑1 (MQ=255)
                                                                         aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc       <  1:406540/90‑1 (MQ=255)
                                                                             aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa   >  1:44142/1‑90 (MQ=255)
                                                                              aTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAc  <  1:307692/90‑1 (MQ=255)
                                                                                     |                                                                                  
CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC  >  NZ_CP009273/3113564‑3113731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
ACAATGGGTGCCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTCTGAAGGCA  >  NZ_CP009273/3113554‑3113746
                                                                                               |                                                                                                 
ACAATGGGTGCCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAA                                                                                               >  SRR3722109.51395/1‑100 (MQ=60)
          CCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAT                                                                                     <  SRR3722109.545783/100‑1 (MQ=60)
                    ccaccgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacagGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATG                                                                           <  SRR3722109.57589/43‑1 (MQ=60)
                               GCAAATGAGGGATCTGGAAATTGATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAC                                                                <  SRR3722109.107898/100‑1 (MQ=60)
                               GCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAC                                                                <  SRR3722109.313752/100‑1 (MQ=60)
                                 AAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCG                                                              >  SRR3722109.137372/1‑100 (MQ=60)
                                    TGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTA                                                           <  SRR3722109.369818/100‑1 (MQ=60)
                                       GGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGA                                                        >  SRR3722109.452546/1‑100 (MQ=60)
                                                   TTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCctgtctcttatacacatctgac                                            >  SRR3722109.424911/1‑78 (MQ=60)
                                                                         CGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCC                      <  SRR3722109.50143/100‑1 (MQ=60)
                                                                             GCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGA                  >  SRR3722109.44666/1‑100 (MQ=60)
                                                                              CCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC                 <  SRR3722109.64026/100‑1 (MQ=60)
                                                                                   AGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATAT            <  SRR3722109.412453/100‑1 (MQ=60)
                                                                                   AGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATAT            <  SRR3722109.560885/100‑1 (MQ=60)
                                                                                        ATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTCTGA       <  SRR3722109.311709/100‑1 (MQ=60)
                                                                                             GCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATctgtctctta  >  SRR3722109.492709/1‑90 (MQ=60)
                                                                                               |                                                                                                 
ACAATGGGTGCCGCTTTCAACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTCTGAAGGCA  >  NZ_CP009273/3113554‑3113746

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