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 I178 R1 167 28.0 1137820 95.3% 1084342 115.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 609,496 A→C L33F (TTA→TTC ybdZ → enterobactin biosynthesis protein YbdZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273609,4960AC100.0% 91.2 / NA 25L33F (TTA→TTCybdZenterobactin biosynthesis protein YbdZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/13);  total (12/13)

CCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAAT  >  NZ_CP009273/609362‑609614
                                                                                                                                      |                                                                                                                      
ccTCTGGCAACCCCTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCgg                                                                                                                    <  2:262890/139‑1 (MQ=255)
      gCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCt                                                                                                              <  2:56029/139‑1 (MQ=255)
                                      ggCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCg                                                                              >  2:129873/1‑139 (MQ=255)
                                      ggCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCg                                                                              >  2:505973/1‑139 (MQ=255)
                                       gCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGt                                                                             <  2:479331/139‑1 (MQ=255)
                                           tCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTg                                                                         >  1:327112/1‑139 (MQ=255)
                                           tCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTg                                                                         >  1:401911/1‑139 (MQ=255)
                                                 aTCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCcagca                                                                   <  1:129873/139‑1 (MQ=255)
                                                       tCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCt                                                             >  2:369755/1‑139 (MQ=255)
                                                             atCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCgg                                                                                                                    >  2:131896/1‑78 (MQ=255)
                                                             atCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCgg                                                                                                                    <  1:131896/78‑1 (MQ=255)
                                                              tCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGcc                                                      <  2:401911/139‑1 (MQ=255)
                                                               ccGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAg                                                                                       <  1:520715/105‑1 (MQ=255)
                                                               ccGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAg                                                                                       >  2:520715/1‑105 (MQ=255)
                                                                  cAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACt                                                  >  2:124467/1‑139 (MQ=255)
                                                                      gAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCg                                              <  1:415061/139‑1 (MQ=255)
                                                                       aGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCcag                                                                     <  1:234931/115‑1 (MQ=255)
                                                                       aGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCcag                                                                     >  2:234931/1‑115 (MQ=255)
                                                                                      gcgcAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACcgac                               >  2:282233/1‑138 (MQ=255)
                                                                                               gcAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                     <  1:124467/139‑1 (MQ=255)
                                                                                                            aGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCaggagg        <  2:450642/139‑1 (MQ=255)
                                                                                                            aGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCaggagg        <  1:369755/139‑1 (MQ=255)
                                                                                                             gTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                     <  2:386317/125‑1 (MQ=255)
                                                                                                             gTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                     >  1:386317/1‑125 (MQ=255)
                                                                                                                   ggCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAAt  >  1:166920/1‑138 (MQ=255)
                                                                                                                                      |                                                                                                                      
CCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAAT  >  NZ_CP009273/609362‑609614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
ACCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCACAG  >  NZ_CP009273/609361‑609645
                                                                                                                                       |                                                                                                                                                     
atacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagTCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGAC                                                                                                                                        <  SRR3722132.134338/100‑1 (MQ=60)
                               cggagatgtgtataagagacagCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTC                                                                                                           <  SRR3722132.529944/127‑1 (MQ=60)
                                 GAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGC                                                                                                       >  SRR3722132.333219/1‑151 (MQ=60)
                                 GAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGC                                                                                                       >  SRR3722132.409187/1‑151 (MQ=60)
                                                gtataagagacagGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGA                                                                                         <  SRR3722132.239337/137‑1 (MQ=60)
                                                 AATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAG                                                                                       <  SRR3722132.132287/151‑1 (MQ=60)
                                                                      GGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACAC                                                                  <  SRR3722132.422579/151‑1 (MQ=60)
                                                                                               CGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGA                                         <  SRR3722132.126799/151‑1 (MQ=60)
                                                                                                   GGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGActgt                                     >  SRR3722132.393347/1‑147 (MQ=60)
                                                                                                         ATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATG                                >  SRR3722132.170062/1‑150 (MQ=60)
                                                                                                            CAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCC                            <  SRR3722132.376525/151‑1 (MQ=60)
                                                                                                                                      TCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaacaaataagcgtaattataaactctgc  >  SRR3722132.394518/1‑43 (MQ=60)
                                                                                                                                       |                                                                                                                                                     
ACCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCACAG  >  NZ_CP009273/609361‑609645

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