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 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 858,772 T→C I66V (ATT→GTT)  ybiY ← glycyl‑radical enzyme activating protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273858,7720TC100.0% 116.2 / NA 31I66V (ATT→GTT) ybiYglycyl‑radical enzyme activating protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/19);  total (12/19)

GTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGC  >  NZ_CP009273/858685‑858855
                                                                                       |                                                                                   
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                  >  1:196604/1‑91 (MQ=255)
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                  >  2:716535/1‑91 (MQ=255)
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                  >  2:198065/1‑91 (MQ=255)
gtggCGGACAGCGGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTCTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                  >  1:813518/1‑91 (MQ=255)
        cagcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcg                                                                          <  2:442593/91‑1 (MQ=255)
        cagcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcg                                                                          <  2:229886/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                         <  2:512232/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                         <  2:374816/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                         <  1:96752/91‑1 (MQ=255)
          gcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGcgcg                                                                        >  1:213454/1‑91 (MQ=255)
                  ttAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCCATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAgc                                                                >  2:561235/1‑91 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGc                                                     <  2:706677/89‑1 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGc                                                     >  1:706677/1‑89 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                   <  2:182491/91‑1 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                   <  2:229728/91‑1 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                   <  1:314476/91‑1 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                   <  1:716535/91‑1 (MQ=255)
                                    ccGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTc                                              >  2:169939/1‑91 (MQ=255)
                                                tttCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCg                                  <  2:176564/91‑1 (MQ=255)
                                                tttCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCg                                  <  2:815767/91‑1 (MQ=255)
                                                   cccGATGAATAAGCAAACCATTCGGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGc                               <  2:342238/91‑1 (MQ=255)
                                                   cccGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGc                               <  1:10781/91‑1 (MQ=255)
                                                   cccGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGc                               <  1:33124/91‑1 (MQ=255)
                                                                   aCCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATcc               <  2:202796/91‑1 (MQ=255)
                                                                     cATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCtg             >  2:661758/1‑91 (MQ=255)
                                                                           gcgcgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                   >  1:787340/1‑47 (MQ=255)
                                                                           gcgcgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                   <  2:787340/47‑1 (MQ=255)
                                                                             gcgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGc                               >  1:562213/1‑65 (MQ=255)
                                                                             gcgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGc                               <  2:562213/65‑1 (MQ=255)
                                                                              cgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATa                         <  2:426211/70‑1 (MQ=255)
                                                                              cgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATa                         >  1:426211/1‑70 (MQ=255)
                                                                                cgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGgcgc  <  1:537557/91‑1 (MQ=255)
                                                                                cgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGgcgc  >  1:322618/1‑91 (MQ=255)
                                                                                       |                                                                                   
GTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGC  >  NZ_CP009273/858685‑858855

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ 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
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGG  >  NZ_CP009273/858673‑858865
                                                                                                   |                                                                                             
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACC                                                                                              >  SRR3722118.206848/1‑101 (MQ=60)
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACC                                                                                              >  SRR3722118.42430/1‑101 (MQ=60)
  GTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAC                                                                                            >  SRR3722118.200502/1‑101 (MQ=60)
  GTTAATGCCTGTGGCGGACAGCGGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTCTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAC                                                                                            >  SRR3722118.831413/1‑101 (MQ=60)
            GTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCG                                                                                  >  SRR3722118.217745/1‑101 (MQ=60)
                     AGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCG                                                                         <  SRR3722118.98520/101‑1 (MQ=60)
                                 ACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAG                                                             >  SRR3722118.722116/1‑101 (MQ=60)
                                           ATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGAC                                                   <  SRR3722118.320952/101‑1 (MQ=60)
                                           ATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCACAC                                                   <  SRR3722118.732252/101‑1 (MQ=60)
                                                               CCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGC                               <  SRR3722118.33718/101‑1 (MQ=60)
                                                               CCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGC                               <  SRR3722118.10996/101‑1 (MQ=60)
                                                                             AAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACctgtctcttatacacatctccgagcccacgagac                 >  SRR3722118.804624/1‑67 (MQ=60)
                                                                               ACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCctgtctcttatacaca               >  SRR3722118.574133/1‑85 (MQ=60)
                                                                                CCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCctgtctcttat              >  SRR3722118.435206/1‑90 (MQ=60)
                                                                                  ATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGC            >  SRR3722118.329261/1‑101 (MQ=60)
                                                                                       GCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTT       >  SRR3722118.393219/1‑101 (MQ=60)
                                                                                            CGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGG  <  SRR3722118.548963/101‑1 (MQ=60)
                                                                                            CGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGG  >  SRR3722118.107059/1‑101 (MQ=60)
                                                                                                   |                                                                                             
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGG  >  NZ_CP009273/858673‑858865

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