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 F3 I95 R1 196 14.2 808436 96.0% 776098 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,896,621 A→G D147D (GAT→GAC bglB ← 6‑phospho‑beta‑glucosidase BglB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,896,6210AG100.0% 52.5 / NA 17D147D (GAT→GACbglB6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/3);  total (14/3)

TGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAC  >  NZ_CP009273/3896539‑3896693
                                                                                   |                                                                        
tGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGa                                                                   >  1:211018/1‑90 (MQ=255)
         cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                          >  1:164652/1‑90 (MQ=255)
         cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                          <  1:167279/90‑1 (MQ=255)
         cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                          >  1:131443/1‑90 (MQ=255)
           ttAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                               >  2:348509/1‑83 (MQ=255)
           ttAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                               <  1:348509/83‑1 (MQ=255)
                               cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                                    >  2:45720/1‑90 (MQ=255)
                               cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                                    >  2:385341/1‑90 (MQ=255)
                                   ttATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCAcc                                >  1:69511/1‑90 (MQ=255)
                                          ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                         >  1:41830/1‑90 (MQ=255)
                                          ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                         >  2:402748/1‑90 (MQ=255)
                                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg        >  2:2586/1‑90 (MQ=255)
                                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg        >  2:302241/1‑90 (MQ=255)
                                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg        >  1:259495/1‑90 (MQ=255)
                                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg        >  1:175417/1‑90 (MQ=255)
                                                               gggCATAATGCCGAAGTGGGTCGGGGCCCGCTCAATAACCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa    <  1:402748/90‑1 (MQ=255)
                                                                 gCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAc  >  1:25913/1‑90 (MQ=255)
                                                                                   |                                                                        
TGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAC  >  NZ_CP009273/3896539‑3896693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ 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
TGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGC  >  NZ_CP009273/3896529‑3896720
                                                                                            |                                                                                                   
TGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGA                                                                                              >  SRR3722235.213627/1‑100 (MQ=60)
 gcctcttgtctcgtgggctcggagatgtgtataagagacagACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAC                                                                                             <  SRR3722235.305514/59‑1 (MQ=60)
         ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAT                                                                                     >  SRR3722235.132969/1‑100 (MQ=60)
         ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAT                                                                                     >  SRR3722235.166550/1‑100 (MQ=60)
              GATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCC                                                                                <  SRR3722235.353530/100‑1 (MQ=60)
                   CATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACC                                                                           <  SRR3722235.169201/100‑1 (MQ=60)
                                   TAACCCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACC                                                           >  SRR3722235.70300/1‑100 (MQ=60)
                                          ACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCAT                                                    >  SRR3722235.42292/1‑100 (MQ=60)
                                             tttttTTGGTAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAGGGGGTCGAGGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATG                                                 <  SRR3722235.317579/95‑1 (MQ=44)
                                             tttaaaagagaaagagaaaagcccggggGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATG                                                 <  SRR3722235.213061/72‑1 (MQ=60)
                                                 gatgtttaaaagagaccgccggggGGGCGTAATGCTAGAAGGGGTCGATGCCCGCCCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCAT                                             <  SRR3722235.375779/76‑1 (MQ=56)
                                                           AGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGG                                   >  SRR3722235.262918/1‑100 (MQ=60)
                                                           AGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGG                                   >  SRR3722235.177431/1‑100 (MQ=60)
                                                            GTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCtgtctcttatacacatctccgagcccacgagacaagaggcaatatcgtatgccgtctt                                  >  SRR3722235.272165/1‑42 (MQ=60)
                                                                 GACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAC                             >  SRR3722235.26212/1‑100 (MQ=60)
                                                                                            GTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGC  >  SRR3722235.26156/1‑100 (MQ=60)
                                                                                            GTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGC  >  SRR3722235.108902/1‑100 (MQ=60)
                                                                                            |                                                                                                   
TGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGC  >  NZ_CP009273/3896529‑3896720

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