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 I210 R1 226 19.6 1130256 95.4% 1078264 84.4

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% 81.8 / NA 26G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (15/11);  total (15/11)

CCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785611‑3785771
                                                                            |                                                                                    
ccAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                         >  1:97067/1‑90 (MQ=255)
ccAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                         <  2:303760/90‑1 (MQ=255)
  aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGaa                                                                       >  2:476658/1‑90 (MQ=255)
      aGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCa                                                                   >  2:512630/1‑90 (MQ=255)
                   ttGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgc                                                      <  2:501354/90‑1 (MQ=255)
                        cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                 <  1:21578/90‑1 (MQ=255)
                             tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:256696/90‑1 (MQ=255)
                             tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:79243/90‑1 (MQ=255)
                             tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:329423/90‑1 (MQ=255)
                                  aTGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                       <  1:383156/90‑1 (MQ=255)
                                  aTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                       <  1:383097/90‑1 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:1943/1‑90 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:464236/1‑90 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:443428/1‑90 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:375503/1‑90 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:210312/1‑90 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:555035/1‑90 (MQ=255)
                                        gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:24770/1‑90 (MQ=255)
                                             gAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGGCTTTGCCGCCGCAATCAGATCAGGATGATTCg                            >  2:479945/1‑90 (MQ=255)
                                                  ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       >  2:450065/1‑90 (MQ=255)
                                                  ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       <  1:375503/90‑1 (MQ=255)
                                                  ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       >  2:545411/1‑90 (MQ=255)
                                                  ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       >  2:6026/1‑90 (MQ=255)
                                                    gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     <  1:450065/90‑1 (MQ=255)
                                                                cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt         >  1:452922/1‑90 (MQ=255)
                                                                       gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  2:18437/90‑1 (MQ=255)
                                                                            |                                                                                    
CCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785611‑3785771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ 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
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA  >  NZ_CP009273/3785595‑3785781
                                                                                            |                                                                                              
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                         <  SRR3722092.182873/100‑1 (MQ=60)
      TCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAG                                                                                   >  SRR3722092.98729/1‑100 (MQ=60)
                                        CAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                 <  SRR3722092.21939/100‑1 (MQ=60)
                                             TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                            <  SRR3722092.260883/100‑1 (MQ=60)
                                             TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                            <  SRR3722092.80604/100‑1 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722092.213774/1‑100 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722092.25194/1‑100 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                           >  SRR3722092.566168/1‑100 (MQ=60)
                                                  ATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                       <  SRR3722092.390439/100‑1 (MQ=60)
                                                  ATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                       <  SRR3722092.390377/100‑1 (MQ=60)
                                                                  CCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGT                       <  SRR3722092.382601/100‑1 (MQ=60)
                                                                    GCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                     <  SRR3722092.458931/100‑1 (MQ=60)
                                                                      AAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTT                   >  SRR3722092.461855/1‑100 (MQ=60)
                                                                                       GCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA  >  SRR3722092.349412/1‑100 (MQ=60)
                                                                                            |                                                                                              
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA  >  NZ_CP009273/3785595‑3785781

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