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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC100.0% 96.3 / NA 29G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (14/15);  total (14/15)

CCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785605‑3785772
                                                                                  |                                                                                     
ccATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAc                                                                                <  2:63658/90‑1 (MQ=255)
ccATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAc                                                                                <  2:308152/90‑1 (MQ=255)
 cATCCCCAGGAAGCCCGCCAGTTTTTGTTCAAACTCTTTATGCCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAcc                                                                               <  1:3771/90‑1 (MQ=255)
       cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  1:558133/1‑90 (MQ=255)
       cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  2:306735/1‑90 (MQ=255)
                 gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                               >  2:314675/1‑90 (MQ=255)
                 gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                               >  2:21114/1‑90 (MQ=255)
                           gTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgccg                                                     >  1:507860/1‑90 (MQ=255)
                              cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAACCCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  1:240470/90‑1 (MQ=255)
                               aaGCTCTTTATGGCTGCCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCcgccgc                                                 <  1:392720/90‑1 (MQ=255)
                                        aTGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:203389/90‑1 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      >  1:743968/1‑54 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      <  2:743968/54‑1 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:771764/1‑90 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:511178/1‑90 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:49864/1‑90 (MQ=255)
                                              gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:494054/1‑90 (MQ=255)
                                                   gAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCg                             >  2:439014/1‑90 (MQ=255)
                                                        ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:526072/1‑90 (MQ=255)
                                                          gCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  1:436297/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:114867/1‑90 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:435545/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:518329/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:578896/90‑1 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:187806/1‑90 (MQ=255)
                                                          gCAAATAAAACGCACCGAAGCCACCCCGAAACCGTGAGAATCCATCCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  1:334222/90‑1 (MQ=255)
                                                                            aGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGcaca    <  2:85583/90‑1 (MQ=255)
                                                                             gCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  1:685614/90‑1 (MQ=255)
                                                                              ccATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaaa  <  1:657519/90‑1 (MQ=255)
                                                                                  |                                                                                     
CCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785605‑3785772

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGAC  >  NZ_CP009273/3785588‑3785782
                                                                                                   |                                                                                               
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                 >  SRR3722073.590989/1‑100 (MQ=60)
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                 >  SRR3722073.438342/1‑100 (MQ=60)
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                 >  SRR3722073.529233/1‑100 (MQ=60)
              CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                   >  SRR3722073.569779/1‑100 (MQ=60)
                  CATCCCCAGGAAGCCCGCCAGTTTTTGTTCAAACTCTTTATGCCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGACAATCC                                                                               <  SRR3722073.3877/100‑1 (MQ=60)
                                  GCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCG                                                               >  SRR3722073.518571/1‑100 (MQ=60)
                                               CAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAACCCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                  <  SRR3722073.245515/100‑1 (MQ=60)
                                                AAGCTCTTTATGGCTGCCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCA                                                 <  SRR3722073.401084/100‑1 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCctgtctcttatacacatctgacgctg                                            >  SRR3722073.759043/1‑74 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                            >  SRR3722073.787331/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                            >  SRR3722073.521952/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                            >  SRR3722073.51037/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                            >  SRR3722073.504514/1‑100 (MQ=60)
                                                         ATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                        <  SRR3722073.207515/100‑1 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                  >  SRR3722073.537105/1‑100 (MQ=60)
                                                                 CCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCC                                >  SRR3722073.191501/1‑100 (MQ=60)
                                                                 CCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCC                                >  SRR3722073.117332/1‑100 (MQ=60)
                                                                           GCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                      <  SRR3722073.445628/100‑1 (MQ=60)
                                                                           GCAAATAAAACGCACCGAAGCCACCCCGAAACCGTGAGAATCCATCCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                      <  SRR3722073.341314/100‑1 (MQ=60)
                                                                                       CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCGCAAAAG          >  SRR3722073.637903/1‑100 (MQ=60)
                                                                                         CCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTT        >  SRR3722073.743278/1‑100 (MQ=60)
                                                                                              GCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA   <  SRR3722073.699585/100‑1 (MQ=60)
                                                                                               CCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGAC  <  SRR3722073.670999/100‑1 (MQ=60)
                                                                                                   |                                                                                               
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGAC  >  NZ_CP009273/3785588‑3785782

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