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 I231 R1 218 12.8 703680 97.2% 683976 86.9

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% 62.7 / NA 19Y292S (TAC→TCC) waaLO‑antigen ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/4);  total (15/4)

ATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791098‑3791263
                                                                                    |                                                                                 
aTCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGaaa                                                                              >  1:101985/1‑90 (MQ=255)
aTCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGaaa                                                                              >  1:148179/1‑90 (MQ=255)
aTCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGaaa                                                                              >  1:164480/1‑90 (MQ=255)
aTCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGaaa                                                                              >  1:51979/1‑90 (MQ=255)
                        aCTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGt                                                      >  2:158510/1‑90 (MQ=255)
                         cTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTc                                                     >  2:269454/1‑90 (MQ=255)
                            aaaCAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCAcc                                                  <  1:246366/90‑1 (MQ=255)
                                   tataCCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCa                                           >  1:32158/1‑90 (MQ=255)
                                         aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                                     >  1:163212/1‑90 (MQ=255)
                                           tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATc                                   >  1:13927/1‑90 (MQ=255)
                                           tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATc                                   >  2:188306/1‑90 (MQ=255)
                                                     gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTc                                                     <  1:15565/62‑1 (MQ=255)
                                                     gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTc                                                     >  2:15565/1‑62 (MQ=255)
                                                     gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:131117/1‑90 (MQ=255)
                                                     gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  1:9545/1‑90 (MQ=255)
                                                             ccTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGa                 <  1:150601/90‑1 (MQ=255)
                                                             ccTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGa                 >  2:314120/1‑90 (MQ=255)
                                                              cTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGaa                <  1:149836/90‑1 (MQ=255)
                                                                            ggCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGtt  >  1:255178/1‑90 (MQ=255)
                                                                                    |                                                                                 
ATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791098‑3791263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ 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
AAACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791085‑3791263
                                                                                                 |                                                                                 
AAACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCG                                                                                 >  SRR3722116.226239/1‑100 (MQ=60)
   CCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAA                                                                              >  SRR3722116.103237/1‑100 (MQ=60)
   CCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAA                                                                              >  SRR3722116.149884/1‑100 (MQ=60)
   CCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAA                                                                              >  SRR3722116.166370/1‑100 (MQ=60)
   CCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAA                                                                              >  SRR3722116.52620/1‑100 (MQ=60)
                                      gatgtgtataagagacagTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCA                                           <  SRR3722116.15775/82‑1 (MQ=60)
                                      CTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCA                                           >  SRR3722116.32536/1‑100 (MQ=60)
                                         AAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTT                                        <  SRR3722116.249806/100‑1 (MQ=60)
                                            CAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGA                                     >  SRR3722116.165082/1‑100 (MQ=60)
                                              GTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATC                                   >  SRR3722116.14113/1‑100 (MQ=60)
                                                        TGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCA                         >  SRR3722116.132657/1‑100 (MQ=60)
                                                        TGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCA                         >  SRR3722116.9663/1‑100 (MQ=60)
                                                                          CCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAAT       <  SRR3722116.152323/100‑1 (MQ=60)
                                                                           CTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATT      <  SRR3722116.151553/100‑1 (MQ=60)
                                                                               GTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  SRR3722116.258802/1‑100 (MQ=60)
                                                                                                 |                                                                                 
AAACCAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791085‑3791263

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