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,791,182 A→C Y292S (TAC→TCC)  waaL → O‑antigen ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,791,1820AC100.0% 66.6 / NA 21Y292S (TAC→TCC) waaLO‑antigen ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/6);  total (15/6)

GAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAA  >  NZ_CP009273/3791096‑3791249
                                                                                      |                                                                   
gAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGa                                                                  >  1:199/1‑90 (MQ=255)
gAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGa                                                                  >  1:325096/1‑90 (MQ=255)
gAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGa                                                                  >  1:487620/1‑90 (MQ=255)
    cGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAAtt                                                              <  1:401780/90‑1 (MQ=255)
     gCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTg                                                             >  1:251954/1‑90 (MQ=255)
          aaTGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTa                                                        >  2:482815/1‑90 (MQ=255)
               aGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAatat                                                   <  1:525061/90‑1 (MQ=255)
                gCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAatata                                                  >  2:123169/1‑90 (MQ=255)
                    tAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCa                                              <  2:487620/90‑1 (MQ=255)
                                   gTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCttttt                               <  1:110691/90‑1 (MQ=255)
                                     tataCCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCa                             >  2:274238/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                       >  1:141073/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                       >  1:283907/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                       >  2:283247/1‑90 (MQ=255)
                                             tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAAtt                                                              >  2:490905/1‑49 (MQ=255)
                                             tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAAtt                                                              <  1:490905/49‑1 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa           >  2:512248/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa           >  2:535128/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa           >  1:188140/1‑90 (MQ=255)
                                                                 tAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGaa  <  2:498318/89‑1 (MQ=255)
                                                                 tAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGaa  >  1:498318/1‑89 (MQ=255)
                                                                                      |                                                                   
GAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAA  >  NZ_CP009273/3791096‑3791249

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  >  NZ_CP009273/3791086‑3791279
                                                                                                |                                                                                                 
AACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGA                                                                                                >  SRR3722111.495108/1‑100 (MQ=60)
AACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGA                                                                                                >  SRR3722111.329434/1‑100 (MQ=60)
AACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGA                                                                                                >  SRR3722111.201/1‑100 (MQ=60)
   CAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAAT                                                                                             <  SRR3722111.82144/100‑1 (MQ=60)
     ATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTG                                                                                           >  SRR3722111.254916/1‑100 (MQ=60)
     ATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTG                                                                                           >  SRR3722111.326238/1‑100 (MQ=60)
              gtcggcagcgtcagatgtgtataagagacagGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATA                                                                                  <  SRR3722111.498449/69‑1 (MQ=60)
              atctacacattgacattcgtcggcagcgtcagatgtgtataagagacagGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATA                                                                                  <  SRR3722111.51779/51‑1 (MQ=60)
              CGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATA                                                                                  <  SRR3722111.407638/100‑1 (MQ=60)
                  cgaccaccgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATT                                                                              <  SRR3722111.357977/42‑1 (MQ=60)
                         AGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAG                                                                       <  SRR3722111.533280/100‑1 (MQ=60)
                                           CAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGA                                                     >  SRR3722111.287494/1‑100 (MQ=60)
                                           CAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGA                                                     >  SRR3722111.142853/1‑100 (MQ=60)
                                             GTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATC                                                   <  SRR3722111.112082/100‑1 (MQ=60)
                                                       TGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCA                                         >  SRR3722111.190488/1‑100 (MQ=60)
                                                                 GTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAA                               >  SRR3722111.505987/1‑100 (MQ=60)
                                                                                              GTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  >  SRR3722111.266792/1‑100 (MQ=60)
                                                                                                |                                                                                                 
AACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  >  NZ_CP009273/3791086‑3791279

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