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
A2 F5 I184 R1 278 29.6 1166432 97.3% 1134938 124.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 857,539 T→C E175G (GAA→GGA)  ybiW ← glycyl radical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273857,5390TC100.0% 54.0 / NA 16E175G (GAA→GGA) ybiWglycyl radical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/7);  total (9/7)

TTCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTTCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTC  >  NZ_CP009273/857409‑857673
                                                                                                                                  |                                                                                                                                      
ttCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCcagcag                                                                                                                                 <  2:228280/138‑1 (MQ=255)
 tCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCcagcag                                                                                                                                 >  1:344847/1‑137 (MQ=255)
  cACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCcagcagca                                                                                                                               <  2:428322/138‑1 (MQ=255)
                 gCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGc                                                                                                               <  2:430475/139‑1 (MQ=255)
                  cACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGcc                                                                                                              <  1:443596/139‑1 (MQ=255)
                                   tCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCg                                                                                               >  2:32569/1‑137 (MQ=255)
                                          ttGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCata                                                                                      >  2:565405/1‑139 (MQ=255)
                                                      tAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGc                                                                           <  2:571913/138‑1 (MQ=255)
                                                              ccAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGAt                                                                  >  2:42593/1‑139 (MQ=255)
                                                                       tCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGc                                                         >  1:43320/1‑139 (MQ=255)
                                                                               aTGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGAc                                                 >  1:265778/1‑139 (MQ=255)
                                                                               aTGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGAc                                                 >  2:35890/1‑139 (MQ=255)
                                                                                                     cTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACa                            <  2:136926/138‑1 (MQ=255)
                                                                                                          tCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCg                       >  1:489497/1‑138 (MQ=255)
                                                                                                             cgcAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCg                                                              <  1:569408/96‑1 (MQ=255)
                                                                                                             cgcAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCg                                                              >  2:569408/1‑96 (MQ=255)
                                                                                                                              ttttcCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTc  >  1:375154/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TTCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTTCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTC  >  NZ_CP009273/857409‑857673

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTCAATGTGTTCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTTCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTCTGACCGCGCCACCAC  >  NZ_CP009273/857399‑857688
                                                                                                                                            |                                                                                                                                                     
GTTCAATGTGTTCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGC                                                                                                                                               >  SRR3722139.349031/1‑149 (MQ=60)
                           GCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATC                                                                                                                  <  SRR3722139.448802/151‑1 (MQ=60)
                                                                      TTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCC                                                                       >  SRR3722139.43824/1‑151 (MQ=60)
                                                                              cgtgggctcggagatgtgtataagagacaGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAG                                                                  <  SRR3722139.575927/119‑1 (MQ=60)
                                                                              CCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACC                                                               >  SRR3722139.269114/1‑151 (MQ=60)
                                                                                                         CCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGT                                     >  SRR3722139.495178/1‑150 (MQ=60)
                                                                                                                             CCATCAAGCCCTTTTCCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTCT                >  SRR3722139.379643/1‑151 (MQ=60)
                                                                                                                                            CCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTCTGACCGCGCCACCAC  <  SRR3722139.458479/150‑1 (MQ=60)
                                                                                                                                            |                                                                                                                                                     
GTTCAATGTGTTCACTGACTGCCACCAGCACGATATCAATCGCTTTCAGGAATTGCTCACCGTGTAAATCTTCCAGCACCGTCAGGTTGATGCGCGAGCGACGTTCCGCTACTTCCTCGCGCAGACCATCAAGCCCTTTTTCCAGCAGCAGCGGGAAATTCACCGCCAGGTGCGCATCGCCGGAGGTCATATTGCCTTCCGCTTTAATGATTCCGGTCGCCAGCAGACCTTTTTGCTCATCGGTAAACATGCCGTAGCAGCGATCCTGTACGGTCTGACCGCGCCACCAC  >  NZ_CP009273/857399‑857688

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