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 I196 R1 3 18.2 1013968 97.0% 983548 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,607,750 G→A G204D (GGC→GAC)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,607,7500GA100.0% 79.7 / NA 22G204D (GGC→GAC) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/12);  total (10/12)

TGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTC  >  NZ_CP009273/4607663‑4607833
                                                                                       |                                                                                   
tgAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCggcggacg                                                                                   <  1:325804/90‑1 (MQ=255)
tgAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCggcggacg                                                                                   <  1:287619/90‑1 (MQ=255)
     gcgACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGgcgtg                                                                              <  1:182694/90‑1 (MQ=255)
      cgACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGgcgtgc                                                                             >  1:199776/1‑90 (MQ=255)
                      gcgcATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATgc                                                             <  2:432386/90‑1 (MQ=255)
                                   gAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTAc                                                                         >  1:421573/1‑65 (MQ=255)
                                   gAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTAc                                                                         <  2:421573/65‑1 (MQ=255)
                                   gAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATAtctc                                                >  1:23307/1‑90 (MQ=255)
                                        gATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAt                                           >  1:215275/1‑90 (MQ=255)
                                        gATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAt                                           >  1:265496/1‑90 (MQ=255)
                                         aTCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAtt                                          <  2:320623/90‑1 (MQ=255)
                                         aTCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAtt                                          <  2:236392/90‑1 (MQ=255)
                                             gTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAg                                      >  1:293402/1‑90 (MQ=255)
                                               gATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGAt                                    >  2:343273/1‑90 (MQ=255)
                                                  atGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGaa                                 >  1:13190/1‑90 (MQ=255)
                                                   tGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAAc                                <  1:284032/90‑1 (MQ=255)
                                                   tGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAAc                                <  1:191494/90‑1 (MQ=255)
                                                    gggCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACt                               >  2:168884/1‑90 (MQ=255)
                                                        gTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTc                           <  2:68728/90‑1 (MQ=255)
                                                        gTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTc                           >  1:24971/1‑90 (MQ=255)
                                                        gTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTc                           <  1:223105/90‑1 (MQ=255)
                                                            aaaaaaCCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTg                       <  1:181202/90‑1 (MQ=255)
                                                                       ggTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCa            <  2:316980/90‑1 (MQ=255)
                                                                                 cggcggACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTc  <  2:13190/90‑1 (MQ=255)
                                                                                       |                                                                                   
TGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTC  >  NZ_CP009273/4607663‑4607833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
TAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACT  >  NZ_CP009273/4607652‑4607836
                                                                                                  |                                                                                      
TAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGAC                                                                                       <  SRR3722074.271761/100‑1 (MQ=60)
       GCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGC                                                                                >  SRR3722074.201595/1‑100 (MQ=60)
           TGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTAC                                                                            <  SRR3722074.290498/100‑1 (MQ=60)
           TGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTAC                                                                            <  SRR3722074.329267/100‑1 (MQ=60)
                GCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGG                                                                       <  SRR3722074.184386/100‑1 (MQ=60)
                                    CATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATctgtctcttatacac                                                   >  SRR3722074.426435/1‑85 (MQ=60)
                                    CATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTC                                                   >  SRR3722074.23524/1‑100 (MQ=60)
                                         TGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAT                                              >  SRR3722074.217209/1‑100 (MQ=60)
                                         TGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAT                                              >  SRR3722074.268125/1‑100 (MQ=60)
                                              GAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAG                                         >  SRR3722074.296341/1‑100 (MQ=60)
                                                   GATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAA                                    >  SRR3722074.13303/1‑100 (MQ=60)
                                                         TGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTC                              >  SRR3722074.25200/1‑100 (MQ=60)
                                                              TGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGC                         <  SRR3722074.193257/100‑1 (MQ=60)
                                                              TGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGC                         <  SRR3722074.286879/100‑1 (MQ=60)
                                                                   GTAGAAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACT                    <  SRR3722074.225112/100‑1 (MQ=60)
                                                                       AAAAAACCGTTGGTTTCAAACCGGCGGACGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGC                <  SRR3722074.182888/100‑1 (MQ=60)
                                                                                     tgtataagagacaGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACT  <  SRR3722074.472972/87‑1 (MQ=60)
                                                                                                  |                                                                                      
TAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACT  >  NZ_CP009273/4607652‑4607836

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