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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,367,563 T→C intergenic (‑24/+65) blc ← / ← ampC lipocalin Blc/BlaEC family class C beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,367,5630TC100.0% 108.2 / NA 33intergenic (‑24/+65)blc/ampClipocalin Blc/BlaEC family class C beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/24);  total (9/24)

AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA  >  NZ_CP009273/4367477‑4367647
                                                                                      |                                                                                    
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:66492/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:578553/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:285245/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:130063/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  2:114921/90‑1 (MQ=255)
agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                   <  1:483101/90‑1 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                                >  2:175029/1‑90 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                                >  1:367435/1‑90 (MQ=255)
   aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                                >  1:588723/1‑90 (MQ=255)
      gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTa                                                                             >  1:182400/1‑90 (MQ=255)
          gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                         <  1:138316/90‑1 (MQ=255)
          gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                         <  2:86295/90‑1 (MQ=255)
          gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                         <  1:603246/90‑1 (MQ=255)
          gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                         <  1:550119/90‑1 (MQ=255)
                 ccAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCaa                                                                  <  1:622010/90‑1 (MQ=255)
                         gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGc                                                          >  2:362936/1‑90 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                            <  1:535214/90‑1 (MQ=255)
                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                            >  2:135682/1‑90 (MQ=255)
                                                        ggCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATg                           >  1:137141/1‑90 (MQ=255)
                                                            cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                       <  2:238601/90‑1 (MQ=255)
                                                            cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                       <  2:355002/90‑1 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                     <  1:380627/90‑1 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                     <  2:182400/90‑1 (MQ=255)
                                                              tAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                     <  2:381564/90‑1 (MQ=255)
                                                                 aTGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                  <  1:175029/90‑1 (MQ=255)
                                                                     ttCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa              >  1:408193/1‑90 (MQ=255)
                                                                     ttCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa              <  1:601631/90‑1 (MQ=255)
                                                                     ttCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAACTCGGACCCGATGGAATTTTACTGTa              >  1:179676/1‑90 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga            <  1:270978/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga            <  2:367435/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga            <  2:647131/90‑1 (MQ=255)
                                                                       ccTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga            <  1:140479/90‑1 (MQ=255)
                                                                             tGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa      <  2:170588/90‑1 (MQ=255)
                                                                                 tttttcCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGaa  <  1:16332/90‑1 (MQ=255)
                                                                                      |                                                                                    
AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA  >  NZ_CP009273/4367477‑4367647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ 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
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGC  >  NZ_CP009273/4367470‑4367657
                                                                                             |                                                                                              
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGT                                                                                          >  SRR3722088.595949/1‑100 (MQ=60)
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGT                                                                                          >  SRR3722088.371343/1‑100 (MQ=60)
   TAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTA                                                                                       >  SRR3722088.184225/1‑100 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                   <  SRR3722088.488853/100‑1 (MQ=60)
                 GCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGA                                                                         <  SRR3722088.139736/100‑1 (MQ=60)
                 GCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGA                                                                         <  SRR3722088.610672/100‑1 (MQ=60)
                 GCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGA                                                                         <  SRR3722088.556842/100‑1 (MQ=60)
                        CCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGG                                                                  <  SRR3722088.629688/100‑1 (MQ=60)
                                                     AGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATG                                     >  SRR3722088.138555/1‑100 (MQ=60)
                                                              AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC                            <  SRR3722088.541714/100‑1 (MQ=60)
                                                                  GCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTA                        >  SRR3722088.412686/1‑100 (MQ=60)
                                                                  GCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAACTCGGACCCGATGGAATTTTACTGTA                        >  SRR3722088.181477/1‑100 (MQ=60)
                                                                     TAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAG                     <  SRR3722088.384668/100‑1 (MQ=60)
                                                                        ATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGT                  <  SRR3722088.176783/100‑1 (MQ=60)
                                                                            TTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAG              <  SRR3722088.609039/100‑1 (MQ=60)
                                                                              CCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA            <  SRR3722088.273596/100‑1 (MQ=60)
                                                                              CCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAA            <  SRR3722088.141923/100‑1 (MQ=60)
                                                                                        TTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGC  <  SRR3722088.16510/100‑1 (MQ=60)
                                                                                             |                                                                                              
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGC  >  NZ_CP009273/4367470‑4367657

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: