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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

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% 86.5 / NA 27K330N (AAA→AACglcAglycolate permease GlcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/19);  total (8/19)

ACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGA  >  NZ_CP009273/3113573‑3113730
                                                                            |                                                                                 
acacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCa                                                                      <  2:559738/90‑1 (MQ=255)
 cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                     <  1:232489/90‑1 (MQ=255)
 cacTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                     <  2:42062/90‑1 (MQ=255)
  acTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATg                                                                    <  2:401322/90‑1 (MQ=255)
            gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                          <  1:286451/90‑1 (MQ=255)
            gCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATgg                                                          <  2:528929/90‑1 (MQ=255)
                    gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                  <  2:174622/90‑1 (MQ=255)
                    gggATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                  <  2:385517/90‑1 (MQ=255)
                           ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                           <  2:82752/90‑1 (MQ=255)
                           ggAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTa                                           >  2:579592/1‑90 (MQ=255)
                                          gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                     <  1:419319/49‑1 (MQ=255)
                                          gTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCAt                                                                     >  2:419319/1‑49 (MQ=255)
                                             aaTAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGAcc                         <  1:471377/90‑1 (MQ=255)
                                                    cgcgCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCg                  >  1:135525/1‑90 (MQ=255)
                                                      cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                  <  2:9531/56‑1 (MQ=255)
                                                      cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGc                                                  >  1:9531/1‑56 (MQ=255)
                                                      cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                >  1:3416/1‑90 (MQ=255)
                                                      cgccgcCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGaa                >  1:338850/1‑90 (MQ=255)
                                                           cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt           >  1:458343/1‑90 (MQ=255)
                                                           cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt           <  1:38292/90‑1 (MQ=255)
                                                           cccGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGAtt           <  1:142088/90‑1 (MQ=255)
                                                                aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc      <  2:3416/90‑1 (MQ=255)
                                                                aGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGccc      <  2:47887/90‑1 (MQ=255)
                                                                  cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa    <  1:579592/90‑1 (MQ=255)
                                                                  cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa    <  2:582691/90‑1 (MQ=255)
                                                                  cAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCa    <  1:309978/90‑1 (MQ=255)
                                                                    aaTAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGa  >  2:242321/1‑90 (MQ=255)
                                                                            |                                                                                 
ACACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGA  >  NZ_CP009273/3113573‑3113730

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTC  >  NZ_CP009273/3113574‑3113738
                                                                           |                                                                                         
gtcggcagcgtcagatgtgtataagagacagTTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATG                                                                   <  SRR3722111.425511/69‑1 (MQ=60)
CACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATG                                                                   <  SRR3722111.235278/100‑1 (MQ=60)
           GCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCAC                                                        <  SRR3722111.290090/100‑1 (MQ=60)
                                         GTGAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCG                          >  SRR3722111.137231/1‑100 (MQ=60)
                                           GAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGActgtctcttatacacatctgacgc                        >  SRR3722111.9652/1‑76 (MQ=60)
                                           GAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAA                        >  SRR3722111.343386/1‑100 (MQ=60)
                                           GAATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAA                        >  SRR3722111.3459/1‑100 (MQ=60)
                                            AATAAAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAAT                       <  SRR3722111.478579/100‑1 (MQ=60)
                                                AAACGCGCCGCCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATT                   >  SRR3722111.465308/1‑100 (MQ=60)
                                                          CCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC         <  SRR3722111.143881/100‑1 (MQ=60)
                                                          CCCGGAGCAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGAC         <  SRR3722111.38744/100‑1 (MQ=60)
                                                                 CAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTC  <  SRR3722111.314038/100‑1 (MQ=60)
                                                                 CAAATAACGCGTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTC  <  SRR3722111.588682/100‑1 (MQ=60)
                                                                           |                                                                                         
CACTTGTTGATGCAAATGAGGGATCTGGAAATTAATCACCAGTGAATAAAACGCGCCGCCCGGAGCAAATAACGCTTTAAACGGCTTCATGGTCCAGATGGTGACCAGCACCGTTAAGATTAAAAACGGTGACCACGCTCGAATGATTTGCCCCAGACTATATTC  >  NZ_CP009273/3113574‑3113738

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