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 F3 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,084,653 A→G G358G (GGT→GGC frvB ← PTS fructose‑like transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,084,6530AG100.0% 118.6 / NA 38G358G (GGT→GGCfrvBPTS fructose‑like transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/14);  total (24/14)

GACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAA  >  NZ_CP009273/4084569‑4084739
                                                                                    |                                                                                      
gACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCccgc                                                                                   >  2:153735/1‑90 (MQ=255)
          ccGACCACAATGGCACTACTGGCAGCCCCTTTTTCCTGGCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggc                                                                         <  1:978130/90‑1 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                       >  1:184339/1‑90 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                       >  1:287293/1‑90 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                       >  2:575661/1‑90 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                       >  1:852596/1‑90 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                       >  1:800395/1‑90 (MQ=255)
            gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                       >  1:83748/1‑90 (MQ=255)
                acaATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATc                                                                   <  2:847111/90‑1 (MQ=255)
                      gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                             >  1:657338/1‑90 (MQ=255)
                      gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                             >  2:565905/1‑90 (MQ=255)
                        actactGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGAc                                                           <  1:8509/90‑1 (MQ=255)
                           actGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAAc                                                        >  2:169196/1‑90 (MQ=255)
                           actGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAAc                                                        >  1:731846/1‑90 (MQ=255)
                              ggCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGcc                                                     >  2:257558/1‑90 (MQ=255)
                                cAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcga                                                   >  2:151625/1‑90 (MQ=255)
                                 ggCCCCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgac                                                  <  1:432076/89‑1 (MQ=255)
                                 aGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgac                                                  <  2:825879/90‑1 (MQ=255)
                                              tGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAg                                     >  1:522666/1‑90 (MQ=255)
                                               gTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGt                                    >  2:536433/1‑90 (MQ=255)
                                                tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGAc                                                                             <  2:987259/48‑1 (MQ=255)
                                                tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGAc                                                                             >  1:987259/1‑48 (MQ=255)
                                                tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTc                                   >  1:368124/1‑90 (MQ=255)
                                                    gTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATaa                               <  1:15011/90‑1 (MQ=255)
                                                             ttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtttt                      >  2:726017/1‑90 (MQ=255)
                                                               gggCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                               >  1:748534/1‑63 (MQ=255)
                                                               gggCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                               <  2:748534/63‑1 (MQ=255)
                                                                     aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag              >  1:47404/1‑90 (MQ=255)
                                                                     aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag              >  2:829296/1‑90 (MQ=255)
                                                                       tGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTg            <  1:153735/90‑1 (MQ=255)
                                                                        gAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTga           <  2:657338/90‑1 (MQ=255)
                                                                           cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa        <  2:368124/90‑1 (MQ=255)
                                                                           cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa        <  1:709127/90‑1 (MQ=255)
                                                                           cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa        <  2:800395/90‑1 (MQ=255)
                                                                              cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCa     >  1:994328/1‑90 (MQ=255)
                                                                              cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCa     >  2:271911/1‑90 (MQ=255)
                                                                              cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCa     >  1:320096/1‑90 (MQ=255)
                                                                                 aaGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCCGCAGTGAGAAGCagaa  <  1:754816/90‑1 (MQ=255)
                                                                                    |                                                                                      
GACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAA  >  NZ_CP009273/4084569‑4084739

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ 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
CCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTT  >  NZ_CP009273/4084556‑4084752
                                                                                                 |                                                                                                   
agagacaGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCC                                                                                                   <  SRR3722236.467573/93‑1 (MQ=60)
      ttcgtcggcagcgtcagatgtgtataagagacagACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGC                                                                                             <  SRR3722236.7417/66‑1 (MQ=60)
               CTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGG                                                                                    >  SRR3722236.867521/1‑100 (MQ=60)
               CTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGG                                                                                    >  SRR3722236.187066/1‑100 (MQ=60)
               CTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGG                                                                                    >  SRR3722236.291453/1‑100 (MQ=60)
               CTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGG                                                                                    >  SRR3722236.84910/1‑100 (MQ=60)
               CTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGG                                                                                    >  SRR3722236.814296/1‑100 (MQ=60)
                       gcgtcagatgtgtataagagacagGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCG                                                                            >  SRR3722236.864215/25‑100 (MQ=60)
                       CCGACCACAATGGCACTACTGGCAGCCCCTTTTTCCTGGCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCG                                                                            <  SRR3722236.995373/100‑1 (MQ=60)
                         GACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCG                                                                          >  SRR3722236.668224/1‑100 (MQ=60)
                              CAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAAC                                                                     >  SRR3722236.744313/1‑100 (MQ=60)
                                     ACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACG                                                              <  SRR3722236.8625/100‑1 (MQ=60)
                                              ggcccCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTAC                                                     <  SRR3722236.438130/95‑1 (MQ=60)
                                                 CTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAG                                                  >  SRR3722236.530656/1‑100 (MQ=60)
                                                   CTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCctgtctcttatacacatctgacgctgccgacg                                                >  SRR3722236.1004655/1‑68 (MQ=60)
                                                   CTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTC                                                >  SRR3722236.373320/1‑100 (MQ=60)
                                                        cccggcCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAA                                           <  SRR3722236.968764/94‑1 (MQ=60)
                                                                 GTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTT                                  <  SRR3722236.15215/100‑1 (MQ=60)
                                                                  TAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGctgtctcttatacacat                                 >  SRR3722236.761357/1‑83 (MQ=60)
                                                                        GTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAG                           >  SRR3722236.48086/1‑100 (MQ=60)
                                                                                 CAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCA                  >  SRR3722236.1011834/1‑100 (MQ=60)
                                                                                 CAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCA                  >  SRR3722236.324677/1‑100 (MQ=60)
                                                                                    TGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAA               <  SRR3722236.155982/100‑1 (MQ=60)
                                                                                        CGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAAC           <  SRR3722236.721138/100‑1 (MQ=60)
                                                                                         GTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACC          >  SRR3722236.905567/1‑100 (MQ=60)
                                                                                              AAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCCGCAGTGAGAAGCAGAAGAACCAGGCG     <  SRR3722236.767764/100‑1 (MQ=60)
                                                                                                GGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGT   >  SRR3722236.635259/1‑100 (MQ=60)
                                                                                                 GCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTT  >  SRR3722236.39558/1‑100 (MQ=60)
                                                                                                 |                                                                                                   
CCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTT  >  NZ_CP009273/4084556‑4084752

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