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 F2 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,898,828 (T)5→4 coding (130/1878 nt) bglF ← PTS beta‑glucoside transporter subunit IIABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,898,8240T.100.0% 103.7 / NA 25coding (134/1878 nt)bglFPTS beta‑glucoside transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/7);  total (18/7)

CCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898741‑3898907
                                                                                   |                                                                                   
ccAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                               >  1:331617/1‑87 (MQ=255)
             ggCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTgc                                                                 >  1:713965/1‑90 (MQ=255)
                 acaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCgctt                                                             >  1:378783/1‑90 (MQ=255)
                        tACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc                                                      >  1:168895/1‑90 (MQ=255)
                        tACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc                                                      >  1:339084/1‑90 (MQ=255)
                                   accTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatt                                           >  2:636867/1‑90 (MQ=255)
                                    ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                          >  2:191228/1‑90 (MQ=255)
                                         aaCTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAc                                     >  1:457360/1‑90 (MQ=255)
                                             ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa                                 >  1:458407/1‑90 (MQ=255)
                                                   cGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGt                           <  2:508079/90‑1 (MQ=255)
                                                          caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                    <  2:315105/90‑1 (MQ=255)
                                                            ccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg                  >  1:444281/1‑90 (MQ=255)
                                                                 ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCa             >  1:514057/1‑90 (MQ=255)
                                                                  taataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAg            <  1:208672/90‑1 (MQ=255)
                                                                     taataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACt         <  2:130385/90‑1 (MQ=255)
                                                                       ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  1:698375/90‑1 (MQ=255)
                                                                       ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  2:339084/90‑1 (MQ=255)
                                                                         aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  2:142739/1‑90 (MQ=255)
                                                                         aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:217836/1‑90 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:92358/1‑90 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:708451/1‑90 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  <  2:221042/90‑1 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:236792/1‑90 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:641545/1‑90 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:686515/1‑90 (MQ=255)
                                                                            gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:7681/1‑90 (MQ=255)
                                                                                   |                                                                                   
CCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898741‑3898907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  NZ_CP009273/3898731‑3898925
                                                                                                 |                                                                                                 
CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAG                                                                                                >  SRR3722112.335624/1‑100 (MQ=60)
   gtcggcagcgtcagatgtgtataagagacagTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTAC                                                                                             <  SRR3722112.710377/69‑1 (MQ=60)
   ggcagcgtcagatgtgtataagagacagGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTAC                                                                                             <  SRR3722112.312838/72‑1 (MQ=60)
             GGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGC                                                                                   >  SRR3722112.725124/1‑100 (MQ=60)
                 GACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTT                                                                               >  SRR3722112.383702/1‑100 (MQ=60)
                        GCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTC                                                                        >  SRR3722112.170968/1‑100 (MQ=60)
                        GCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTC                                                                        >  SRR3722112.343257/1‑100 (MQ=60)
                                         GACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAC                                                       >  SRR3722112.463797/1‑100 (MQ=60)
                                             ACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAA                                                   >  SRR3722112.464870/1‑100 (MQ=60)
                                                            CCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATG                                    >  SRR3722112.450443/1‑100 (MQ=60)
                                                                 TTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCA                               >  SRR3722112.521585/1‑100 (MQ=60)
                                                                         CCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACA                       >  SRR3722112.220460/1‑100 (MQ=60)
                                                                            TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    >  SRR3722112.93517/1‑100 (MQ=60)
                                                                            TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    >  SRR3722112.719517/1‑100 (MQ=60)
                                                                            TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    >  SRR3722112.239618/1‑100 (MQ=60)
                                                                            TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    <  SRR3722112.211198/100‑1 (MQ=60)
                                                                                 ATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC               <  SRR3722112.709283/100‑1 (MQ=60)
                                                                                               TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  SRR3722112.29771/1‑100 (MQ=60)
                                                                                               TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  SRR3722112.22239/1‑100 (MQ=60)
                                                                                                 |                                                                                                 
CTGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  NZ_CP009273/3898731‑3898925

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