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
A4 F30 I1 R1 17 32.7 1267382 95.3% 1207815 125.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC100.0% 53.4 / NA 16I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/7);  total (9/7)

GTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132654‑2132926
                                                                                                                                     |                                                                                                                                           
gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                               <  4:103220/148‑1 (MQ=255)
   gagGGTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGATGTACCCGGAGAATTTTTTCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  3:20368/149‑1 (MQ=255)
                              ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                                >  2:157414/1‑149 (MQ=255)
                                                          cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtg                                                                    >  2:17461/1‑149 (MQ=255)
                                                             tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCCCCACCATCCGCAg                                                                                      >  8:101859/1‑128 (MQ=255)
                                                                        gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:157414/149‑1 (MQ=255)
                                                                           cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               >  5:4579/1‑89 (MQ=255)
                                                                           cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  6:4579/89‑1 (MQ=255)
                                                                            ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                   >  3:48354/1‑148 (MQ=255)
                                                                                   aTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                            <  4:48354/148‑1 (MQ=255)
                                                                                               aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTg                               <  1:17461/149‑1 (MQ=255)
                                                                                                         cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                     >  2:46027/1‑149 (MQ=255)
                                                                                                             aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  >  5:128843/1‑100 (MQ=255)
                                                                                                             aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  <  6:128843/100‑1 (MQ=255)
                                                                                                                  aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      >  2:70208/1‑107 (MQ=255)
                                                                                                                            atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGtt  >  3:35369/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                           
GTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132654‑2132926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTG  >  NC_000913/2132648‑2132927
                                                                                                                                           |                                                                                                                                            
gtgtgtatggagttctgacgtgtgctcttccgatgtgtCGCGGTCGAATTTGTCTTTCTGATTGCGGTAATTGTGTAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCC                                                                                                                                   <  NB501016:228:HFK2FAFXY:4:21507:26657:13143/113‑1 (MQ=60)
        AGAGGGTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGATGTACCCGGAGAATTTTTTCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAC                                                                                                                           <  NB501016:228:HFK2FAFXY:2:11111:24507:8870/151‑1 (MQ=60)
                                             cgacgtgggcgcttccgagctgtaacAGTGGAAGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGGCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGG                                                                                      <  NB501016:228:HFK2FAFXY:4:21503:2170:18407/125‑1 (MQ=60)
                                                                             AGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTC                                                      <  NB501016:228:HFK2FAFXY:1:21302:16766:12206/151‑1 (MQ=60)
                                                                              atcgctgtgactggagttcagacgtggtgctcttccgatcTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTC                                                      <  NB501016:228:HFK2FAFXY:1:11304:20186:4365/110‑1 (MQ=60)
                                                                                ACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTC                                                                                                               >  NB501016:228:HFK2FAFXY:3:11402:2227:18523/1‑91 (MQ=60)
                                                                                 CCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGT                                                   >  NB501016:228:HFK2FAFXY:2:11301:17602:13260/1‑150 (MQ=60)
                                                                                                    GACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGA                               <  NB501016:228:HFK2FAFXY:1:11107:6927:11934/151‑1 (MQ=60)
                                                                                                                  CATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGC                                                                  >  NB501016:228:HFK2FAFXY:3:21508:2491:15190/1‑102 (MQ=60)
                                                                                                                                 GATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTG  >  NB501016:228:HFK2FAFXY:2:11207:17698:4802/1‑151 (MQ=60)
                                                                                                                                           |                                                                                                                                            
GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTG  >  NC_000913/2132648‑2132927

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑