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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

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% 142.3 / NA 40intergenic (‑24/+65)blc/ampClipocalin Blc/BlaEC family class C beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/27);  total (13/27)

GAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATC  >  NZ_CP009273/4367476‑4367649
                                                                                       |                                                                                      
gagaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCt                                                                                      <  1:66650/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  2:765727/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  1:331852/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  2:293571/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  1:374359/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  1:400144/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  2:148204/90‑1 (MQ=255)
 agaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTa                                                                                     <  1:722349/90‑1 (MQ=255)
    aCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgt                                                                                  >  1:420786/1‑90 (MQ=255)
      tGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAgtgt                                                                                <  2:324779/90‑1 (MQ=255)
       gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTa                                                                               >  1:507203/1‑90 (MQ=255)
         aGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGa                                                                             >  2:32782/1‑90 (MQ=255)
           gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATg                                                                           <  2:654766/90‑1 (MQ=255)
             aaCGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGAc                                                                         <  2:273642/90‑1 (MQ=255)
                  ccAGAAATGCAGCTGTCGCTGCGGTAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCaa                                                                    <  1:171862/90‑1 (MQ=255)
                     gAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGa                                                                 >  2:748032/1‑90 (MQ=255)
                          gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGGAGATGACAGCAAGGAAAAGc                                                            >  1:373292/1‑90 (MQ=255)
                               tGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGgaga                                                       >  1:628269/1‑90 (MQ=255)
                               tGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGgaga                                                       >  1:407089/1‑90 (MQ=255)
                                              gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAtt                                        >  1:765005/1‑90 (MQ=255)
                                                  gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGa                                    >  2:613155/1‑90 (MQ=255)
                                                   ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                   <  2:361216/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  2:69494/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  2:272928/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  2:720473/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  2:310888/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  1:337538/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  1:263285/90‑1 (MQ=255)
                                                        aGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                              <  2:507203/90‑1 (MQ=255)
                                                           cgcATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGa                           <  1:588028/90‑1 (MQ=255)
                                                             cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                         <  2:199815/90‑1 (MQ=255)
                                                             cATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                         <  1:617332/90‑1 (MQ=255)
                                                              aTAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt                        >  2:653566/1‑90 (MQ=255)
                                                                       tCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag               <  2:7754/90‑1 (MQ=255)
                                                                       tCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag               >  2:520936/1‑90 (MQ=255)
                                                                          ttACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc            <  1:748032/90‑1 (MQ=255)
                                                                          ttACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc            <  1:396361/90‑1 (MQ=255)
                                                                           tACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCg           >  1:600543/1‑90 (MQ=255)
                                                                             cTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGtt         >  2:54359/1‑90 (MQ=255)
                                                                              tGGTTTATTCCTAAGTGTAGAGGACAGGAAGGAAAAGAGGAGAAAAGGTTCGAAAATTCGGACCCGATGCAATTTTCCTGTCGAGCGTTa        <  1:295801/90‑1 (MQ=255)
                                                                                    tttcCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATc  <  1:244149/90‑1 (MQ=255)
                                                                                       |                                                                                      
GAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATC  >  NZ_CP009273/4367476‑4367649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG  >  NZ_CP009273/4367470‑4367659
                                                                                             |                                                                                                
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGT                                                                                            >  SRR3722073.429795/1‑100 (MQ=60)
   TAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTA                                                                                         >  SRR3722073.517906/1‑100 (MQ=60)
      GAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGAT                                                                                      <  SRR3722073.68212/100‑1 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722073.736945/100‑1 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722073.338904/100‑1 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722073.382306/100‑1 (MQ=60)
       AGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATG                                                                                     <  SRR3722073.408667/100‑1 (MQ=60)
                      GACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGGAGATGACAGCAAGGAAAAGC                                                                      >  SRR3722073.381216/1‑100 (MQ=60)
                        CCAGAAATGCAGCTGTCGCTGCGGTAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGG                                                                    <  SRR3722073.175228/100‑1 (MQ=60)
                           GAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGA                                                                 >  SRR3722073.415780/1‑100 (MQ=60)
                           GAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGA                                                                 >  SRR3722073.641293/1‑100 (MQ=60)
                                          CTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATT                                                  >  SRR3722073.780451/1‑100 (MQ=60)
                                                              AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC                              <  SRR3722073.268841/100‑1 (MQ=60)
                                                              AGGCGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAC                              <  SRR3722073.344707/100‑1 (MQ=60)
                                                                 CGCATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGT                           <  SRR3722073.600312/100‑1 (MQ=60)
                                                                 CGCATAAATGTTGCCGTACTGGGTTTTTGCTAAGGGGAGTGGACAGCAAGGAAAAGCGGAGAAAAGGCCCGAAGAGTCGGACCCGgtggggttctaatgt                           >  SRR3722073.428151/1‑85 (MQ=17)
                                                                   CATAAATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAG                         <  SRR3722073.630164/100‑1 (MQ=60)
                                                                       AATGTTTCCTTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCG                     >  SRR3722073.613067/1‑100 (MQ=60)
                                                                                TTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATC            <  SRR3722073.404804/100‑1 (MQ=60)
                                                                                TTACTGGTTTTTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATC            <  SRR3722073.763198/100‑1 (MQ=60)
                                                                                    TGGTTTATTCCTAAGTGTAGAGGACAGGAAGGAAAAGAGGAGAAAAGGTTCGAAAATTCGGACCCGATGCAATTTTCCTGTCGAGCGTTAAGAATCTGCC        <  SRR3722073.302060/100‑1 (MQ=60)
                                                                                          TTTCCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG  <  SRR3722073.249269/100‑1 (MQ=60)
                                                                                             |                                                                                                
GCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGG  >  NZ_CP009273/4367470‑4367659

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