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 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% 59.2 / 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)

GAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791096‑3791263
                                                                                      |                                                                                 
gAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGa                                                                                >  2:458328/1‑90 (MQ=255)
   tCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAAt                                                                             <  2:655159/90‑1 (MQ=255)
     gCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTg                                                                           >  1:77256/1‑90 (MQ=255)
               aGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAatat                                                                 <  1:400090/90‑1 (MQ=255)
                          aCTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGt                                                      >  1:91069/1‑90 (MQ=255)
                              aaaCAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCAcc                                                  <  1:458202/90‑1 (MQ=255)
                                    ttATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTc                                            >  2:67681/1‑90 (MQ=255)
                                     tataCCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGGCTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCa                                           >  1:297414/1‑90 (MQ=255)
                                     tataCCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCa                                           >  2:531614/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTTATAAAGTCACCTTTTTCATTTAGa                                     >  1:332965/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                                     >  1:512114/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                                     >  2:652336/1‑90 (MQ=255)
                                           aaTGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGa                                     <  1:447518/90‑1 (MQ=255)
                                             tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATc                                   >  2:553575/1‑90 (MQ=255)
                                             tGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATc                                   >  2:653610/1‑90 (MQ=255)
                                                 aaTAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATcagca                               >  2:230232/1‑90 (MQ=255)
                                                 aaTAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATcagca                               >  1:411704/1‑90 (MQ=255)
                                                       gtTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCa                         >  2:218921/1‑90 (MQ=255)
                                                                              ggCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGtt  >  2:653443/1‑90 (MQ=255)
                                                                                      |                                                                                 
GAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTT  >  NZ_CP009273/3791096‑3791263

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
CAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  >  NZ_CP009273/3791089‑3791279
                                                                                             |                                                                                                 
CAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAAT                                                                                             <  SRR3722088.400262/100‑1 (MQ=60)
  ATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTG                                                                                           >  SRR3722088.78039/1‑100 (MQ=60)
                      AGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAG                                                                       <  SRR3722088.404438/100‑1 (MQ=60)
                       GCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGT                                                                      >  SRR3722088.91992/1‑100 (MQ=60)
                                  CTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGGCTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCA                                                           >  SRR3722088.300307/1‑100 (MQ=60)
                                     AAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTT                                                        <  SRR3722088.463560/100‑1 (MQ=60)
                                        CAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTTATAAAGTCACCTTTTTCATTTAGA                                                     >  SRR3722088.336427/1‑100 (MQ=60)
                                        CAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGA                                                     >  SRR3722088.518294/1‑100 (MQ=60)
                                              TACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCA                                               >  SRR3722088.416252/1‑100 (MQ=60)
                                                  AATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGT                                           <  SRR3722088.452707/100‑1 (MQ=60)
                                                     GCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTCCGAAATctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgc                                        >  SRR3722088.501464/1‑47 (MQ=60)
                                                                         AGGTGCAAGACTGGCAATGTCCGAAATTGGTTTAAATATATTCATActgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgcc                    >  SRR3722088.620404/1‑46 (MQ=60)
                                                                                 GACTGGCAATGTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGC            >  SRR3722088.598390/1‑100 (MQ=60)
                                                                                           GTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  >  SRR3722088.163785/1‑100 (MQ=60)
                                                                                           GTCCGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  <  SRR3722088.661541/100‑1 (MQ=60)
                                                                                             |                                                                                                 
CAATACAGAATCGCTATAATGAAGCATTAAATGACTTAAACAGTTATACCAATGCTAATAGTGTTACTTCCCTAGGTGCAAGACTGGCAATGTACGAAATTGGTTTAAATATATTCATAAAGTCACCTTTTTCATTTAGATCAGCAGAGTCACGCGCTGAAAGTATGAATTTGTTAGTTGCAGAACACAAT  >  NZ_CP009273/3791089‑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: