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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,122,677 A→G F152S (TTC→TCC)  wcaE ← colanic acid biosynthesis glycosyltransferase WcaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,122,6770AG100.0% 95.9 / NA 28F152S (TTC→TCC) wcaEcolanic acid biosynthesis glycosyltransferase WcaE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

TATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCC  >  NZ_CP009273/2122594‑2122762
                                                                                   |                                                                                     
tataCATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAata                                                                                 >  2:757802/1‑90 (MQ=255)
  taCATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATc                                                                               <  2:281757/90‑1 (MQ=255)
     aTTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGcc                                                                            <  1:757802/90‑1 (MQ=255)
      ttttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCt                                                                           >  1:135547/1‑90 (MQ=255)
      ttttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCt                                                                           >  2:89529/1‑90 (MQ=255)
      ttttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCt                                                                           >  1:239774/1‑90 (MQ=255)
      ttttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCt                                                                           >  1:340799/1‑90 (MQ=255)
          ggCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCtgatg                                                                       <  1:95695/90‑1 (MQ=255)
                cAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGc                                                                 >  1:138785/1‑90 (MQ=255)
                    gcgTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGc                                                             <  2:784036/90‑1 (MQ=255)
                        aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                                         <  1:214890/90‑1 (MQ=255)
                        aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                                         <  1:689878/90‑1 (MQ=255)
                        aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                                         <  1:281649/90‑1 (MQ=255)
                        aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                                         <  2:135547/90‑1 (MQ=255)
                          tCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCtg                                                       <  2:892708/90‑1 (MQ=255)
                           cGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCtgt                                                      >  1:561091/1‑90 (MQ=255)
                             gaagaaACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGa                                                    <  1:477215/90‑1 (MQ=255)
                             gaagaaACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGa                                                    <  1:882998/90‑1 (MQ=255)
                             gaagaaACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGa                                                    <  2:428704/90‑1 (MQ=255)
                                   aCTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGGAGGCTTGGATAAata                                              >  1:265237/1‑90 (MQ=255)
                                           ttCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGc                                                             <  2:252266/67‑1 (MQ=255)
                                           ttCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGc                                                             >  1:252266/1‑67 (MQ=255)
                                            tCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGccc                                     <  2:576710/90‑1 (MQ=255)
                                              cAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCCGGGTGTGGTAAATATACCAGCCCgg                                   >  1:685713/1‑90 (MQ=255)
                                              cAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCgg                                   >  2:664621/1‑90 (MQ=255)
                                                ggTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATAACAGCCCGGtt                                 >  1:755028/1‑90 (MQ=255)
                                                        cGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTa                         <  2:213869/90‑1 (MQ=255)
                                                                               aGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGcc  <  2:392863/90‑1 (MQ=255)
                                                                                   |                                                                                     
TATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCC  >  NZ_CP009273/2122594‑2122762

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 31 ≤ 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
GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGC  >  NZ_CP009273/2122590‑2122775
                                                                                       |                                                                                                  
GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCT                                                                                        >  SRR3722072.137970/1‑100 (MQ=60)
GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCT                                                                                        >  SRR3722072.243792/1‑100 (MQ=60)
GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCT                                                                                        >  SRR3722072.346805/1‑100 (MQ=60)
         ATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGG                                                                               <  SRR3722072.771168/100‑1 (MQ=60)
          TTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGC                                                                              >  SRR3722072.141255/1‑100 (MQ=60)
              GGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGC                                                                          <  SRR3722072.97466/100‑1 (MQ=60)
                     AGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGT                                                                   >  SRR3722072.571156/1‑100 (MQ=60)
                            AGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAAT                                                            <  SRR3722072.286469/100‑1 (MQ=60)
                            AGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAAT                                                            <  SRR3722072.702216/100‑1 (MQ=60)
                            AGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAAT                                                            <  SRR3722072.218534/100‑1 (MQ=60)
                             GTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGGAGGCTTGGATAAATA                                                           >  SRR3722072.269753/1‑100 (MQ=60)
                                 GAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACC                                                       <  SRR3722072.485895/100‑1 (MQ=60)
                                 GAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACC                                                       <  SRR3722072.898201/100‑1 (MQ=60)
                                     AAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCgggcggtgatctgcctcttttac                                                   >  SRR3722072.256537/1‑77 (MQ=60)
                                        CTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCCGGGTGTGGTAAATATACCAGCCCGG                                                >  SRR3722072.697985/1‑100 (MQ=60)
                                          TTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATAACAGCCCGGTT                                              >  SRR3722072.768355/1‑100 (MQ=60)
                                                                                     GGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAG   <  SRR3722072.676505/100‑1 (MQ=60)
                                                                                      GGAAAATATCGCCTGATGACTGGCGGGCACGCTGTGCTAAATAGACCAGCCCGGGTTGATCCTACGTTTTATTTTATGCCCGTCGCCAAAAACCAGCAGC  >  SRR3722072.458916/1‑100 (MQ=59)
                                                                                       |                                                                                                  
GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGC  >  NZ_CP009273/2122590‑2122775

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