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 I96 R1 227 23.3 882352 98.0% 864704 126.9

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% 58.7 / NA 17I66V (ATT→GTT) ybiYglycyl‑radical enzyme activating protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/8);  total (9/8)

CTCAACGCTTTTCACTTCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAG  >  NZ_CP009273/858647‑858904
                                                                                                                              |                                                                                                                                    
ctcAACGCTTTTCACTTCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcg                                                                                                                           <  1:191071/137‑1 (MQ=255)
                    cacacaCGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCa                                                                                                     <  1:342171/139‑1 (MQ=255)
                    cacacaCGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCa                                                                                                     <  2:405728/139‑1 (MQ=255)
                           ggTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTc                                                                                               >  1:36481/1‑139 (MQ=255)
                                aTGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGaca                                                                                         >  2:3764/1‑139 (MQ=255)
                                       tgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGt                                                                                  >  2:83917/1‑139 (MQ=255)
                                       tgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGt                                                                                  >  2:412385/1‑139 (MQ=255)
                                           ggACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCg                                                                               <  2:191417/138‑1 (MQ=255)
                                           ggACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCg                                                                               <  1:3764/138‑1 (MQ=255)
                                                 cagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTGAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATa                                                                          >  2:249665/1‑137 (MQ=255)
                                                  agTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATAcag                                                                       <  1:47275/139‑1 (MQ=255)
                                                  agTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATAcag                                                                       <  1:210353/139‑1 (MQ=255)
                                                                 tCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCtgcgtgcg                                                        >  2:288176/1‑139 (MQ=255)
                                                                                                              ttCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCgg                                         >  2:190948/1‑110 (MQ=255)
                                                                                                              ttCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCgg                                         <  1:190948/110‑1 (MQ=255)
                                                                                                                cAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAg          >  2:331213/1‑139 (MQ=255)
                                                                                                                         cTCAACCACTTCCGGCGCAGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAg  >  2:293327/1‑138 (MQ=255)
                                                                                                                              |                                                                                                                                    
CTCAACGCTTTTCACTTCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACCAGCGGCAGCCCAGCGAACAGCCTTTAAG  >  NZ_CP009273/858647‑858904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ 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
CCTCAACGCTTTTCACTTCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACC  >  NZ_CP009273/858646‑858876
                                                                                                                               |                                                                                                        
CCTCAACGCTTTTCACTTCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCG                                                                                    <  SRR3722119.193118/149‑1 (MQ=60)
                 TCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCC                                                                  >  SRR3722119.36801/1‑151 (MQ=60)
                    CCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGAC                                                              <  SRR3722119.345658/151‑1 (MQ=60)
                                           CGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAG                                        <  SRR3722119.3796/150‑1 (MQ=60)
                                                  CAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCG                                <  SRR3722119.212613/151‑1 (MQ=60)
                                                  CAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCG                                <  SRR3722119.47691/151‑1 (MQ=60)
                                                                                agatgtgtataagagacagTAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACC  <  SRR3722119.192994/132‑1 (MQ=60)
                                                                                                                               |                                                                                                        
CCTCAACGCTTTTCACTTCACCACACACGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTGACACC  >  NZ_CP009273/858646‑858876

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