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 I205 R1 220 17.5 961422 97.1% 933540 86.5

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% 64.9 / NA 22D147D (GAT→GACbglB6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/6);  total (16/6)

CGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896555‑3896707
                                                                  |                                                                                      
cGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCa                                                                 >  2:223447/1‑90 (MQ=255)
   cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             >  2:475866/1‑75 (MQ=255)
   cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  1:475866/75‑1 (MQ=255)
   cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:317744/1‑90 (MQ=255)
         tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATCAGCCCAACCGCCGTa                                                        >  2:310676/1‑90 (MQ=255)
                aCTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTc                                                 >  1:311136/1‑90 (MQ=255)
                aCTTTATGTTGGTAGCGAGTAAAGACCGCGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTc                                                 >  1:196219/1‑90 (MQ=255)
                          ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:156392/1‑90 (MQ=255)
                          ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:195404/1‑90 (MQ=255)
                          ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:276687/1‑90 (MQ=255)
                             aGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATg                                    >  1:104950/1‑90 (MQ=255)
                                  gaaaaGCCGGTGGGGGCATAATGCTCGAAGTGGTCGATCCCCCCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCAtt                               <  2:311136/88‑1 (MQ=255)
                                         ccGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAAt                        >  2:96437/1‑90 (MQ=255)
                                          cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  1:198569/1‑90 (MQ=255)
                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:173640/1‑90 (MQ=255)
                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:255349/1‑90 (MQ=255)
                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:31953/1‑90 (MQ=255)
                                           gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:330140/1‑90 (MQ=255)
                                               gggCGTAGTGCCCGAGGTGGTGAATCCCCGCTCTATTACCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa                  <  1:348311/90‑1 (MQ=255)
                                               gggCGTAATGCTCGAAGTGGTGGATGCCCGCTCGATTAGCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa                  <  1:223080/90‑1 (MQ=255)
                                                ggCGTATTGCCGGAAGGGGTCGATGCCCGCTTGATAACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATaa                 <  2:422438/90‑1 (MQ=255)
                                                               gTGGTCGATGCCCGCTCTATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt  <  2:127919/90‑1 (MQ=255)
                                                                  |                                                                                      
CGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896555‑3896707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ 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
CGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAT  >  NZ_CP009273/3896525‑3896710
                                                                                                |                                                                                         
CGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCG                                                                                        <  SRR3722087.130711/100‑1 (MQ=60)
     GTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAC                                                                                   <  SRR3722087.133787/100‑1 (MQ=60)
     GTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAC                                                                                   <  SRR3722087.176353/100‑1 (MQ=60)
     GTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAC                                                                                   <  SRR3722087.190636/100‑1 (MQ=60)
                  gacagCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCC                                                                      <  SRR3722087.481402/95‑1 (MQ=60)
                       CATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACC                                                                 >  SRR3722087.321143/1‑100 (MQ=60)
                                    CCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTC                                                    >  SRR3722087.314417/1‑100 (MQ=60)
                                    CCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGCGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTC                                                    >  SRR3722087.197944/1‑100 (MQ=60)
                                              ACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCAT                                          >  SRR3722087.157829/1‑100 (MQ=60)
                                              ACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCAT                                          >  SRR3722087.197125/1‑100 (MQ=60)
                                                 TTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATG                                       >  SRR3722087.105957/1‑100 (MQ=60)
                                                              GAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATG                          >  SRR3722087.200318/1‑100 (MQ=60)
                                                               AGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGG                         >  SRR3722087.257824/1‑100 (MQ=60)
                                                                             GGGCGTAATGCTCGAGGGGGTCGATGACCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAG           <  SRR3722087.231350/100‑1 (MQ=60)
                                                                             GGGCGTAATGCTCGAAGTGGTGGATGCCCGCTCGATTAGCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAG           <  SRR3722087.225085/100‑1 (MQ=60)
                                                                                     TGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGA   >  SRR3722087.259348/1‑100 (MQ=60)
                                                                                      GCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAT  >  SRR3722087.53668/1‑100 (MQ=60)
                                                                                                |                                                                                         
CGCGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAT  >  NZ_CP009273/3896525‑3896710

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