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 I194 R1 197 32.2 1845014 95.2% 1756453 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% 110.1 / NA 33D147D (GAT→GACbglB6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/8);  total (25/8)

ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896538‑3896707
                                                                                   |                                                                                      
aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  >  2:5895/1‑90 (MQ=255)
aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                                  >  1:563141/1‑90 (MQ=255)
 tGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGa                                                                                 >  1:789864/1‑90 (MQ=255)
   ttGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAcc                                                                               >  2:406989/1‑90 (MQ=255)
        ttCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCg                                                                          >  1:774756/1‑90 (MQ=255)
         tCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGa                                                                         <  2:789864/90‑1 (MQ=255)
         tCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGa                                                                         <  1:758219/90‑1 (MQ=255)
          cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                                        >  1:128097/1‑90 (MQ=255)
          cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                                        >  1:847869/1‑90 (MQ=255)
                                   tttATGTTGGAAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCAc                                               <  1:702141/90‑1 (MQ=255)
                                     tttGTTGGTAGCGAAAAAAGCCCGGGCGGGCGTAATGCTCGAGGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCa                                             <  1:41336/88‑1 (MQ=255)
                                     tATGTTGGTAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCa                                             <  1:303966/90‑1 (MQ=255)
                                      ttGTTGGTGGCGAGAAAAGCCCGGGGGGGCGTAATGCTGGAGGGGGTCGATGCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAg                                            <  2:40196/89‑1 (MQ=255)
                                           ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:562288/1‑90 (MQ=255)
                                           ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:37039/1‑90 (MQ=255)
                                            gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata                                      >  2:278154/1‑90 (MQ=255)
                                             tAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCatat                                     >  2:671375/1‑90 (MQ=255)
                                                 gAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCt                                                                            >  1:575675/1‑47 (MQ=255)
                                                 gAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCt                                                                            <  2:575675/47‑1 (MQ=255)
                                                     aaaGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTc                             >  1:449749/1‑90 (MQ=255)
                                                           cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  1:248978/1‑90 (MQ=255)
                                                           cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:317103/1‑90 (MQ=255)
                                                           cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:234377/1‑90 (MQ=255)
                                                           cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  1:542837/1‑90 (MQ=255)
                                                           cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  1:414873/1‑90 (MQ=255)
                                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:122093/1‑90 (MQ=255)
                                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:714217/1‑90 (MQ=255)
                                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:338887/1‑90 (MQ=255)
                                                             tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:278778/1‑90 (MQ=255)
                                                             tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:208438/1‑90 (MQ=255)
                                                             tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:116052/1‑90 (MQ=255)
                                                             tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  2:640852/1‑90 (MQ=255)
                                                                                gTGGTCGACGCCCGCTCGATCACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt  <  2:550397/90‑1 (MQ=255)
                                                                                   |                                                                                      
ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896538‑3896707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ 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
CGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACG  >  NZ_CP009273/3896527‑3896694
                                                                                              |                                                                         
CGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGAT                                                                      <  SRR3722072.536217/100‑1 (MQ=60)
 GTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATG                                                                     >  SRR3722072.573239/1‑100 (MQ=60)
  TGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGA                                                                    >  SRR3722072.803693/1‑100 (MQ=60)
   GTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGAC                                                                   <  SRR3722072.261295/100‑1 (MQ=60)
         ACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCG                                                             >  SRR3722072.788363/1‑100 (MQ=60)
           ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAT                                                           >  SRR3722072.862530/1‑100 (MQ=60)
           ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAT                                                           >  SRR3722072.130415/1‑100 (MQ=60)
                    TCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAC                                                  <  SRR3722072.771595/100‑1 (MQ=60)
                                    aaaaccccttTTTATTTTGGTAGGGAGAAAAGCCCGTGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCAC                                  <  SRR3722072.791678/90‑1 (MQ=57)
                                            ACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCAT                          >  SRR3722072.572374/1‑100 (MQ=60)
                                              TTTATGTTGGAAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATAT                        <  SRR3722072.714676/100‑1 (MQ=60)
                                                TTTGTTGGTAGCGAAAAAAGCCCGGGCGGGCGTAATGCTCGAGGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGG                      <  SRR3722072.42092/100‑1 (MQ=60)
                                                TATGTTGGTAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGG                      <  SRR3722072.309232/100‑1 (MQ=60)
                                                  TGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCctgtctcttatacacatctccgagcccacgaga                    >  SRR3722072.585961/1‑67 (MQ=60)
                                                      GGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTC                >  SRR3722072.457983/1‑100 (MQ=60)
                                                            GAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATG          >  SRR3722072.552611/1‑100 (MQ=60)
                                                            GAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATG          >  SRR3722072.422381/1‑100 (MQ=60)
                                                            GAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATG          >  SRR3722072.253183/1‑100 (MQ=60)
                                                             AGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGG         >  SRR3722072.344847/1‑100 (MQ=60)
                                                                    cccggGCGGGCGTAATGCTCGAAGTGGTCAATCCCCGCTGGATTAGCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACG  <  SRR3722072.656708/95‑1 (MQ=60)
                                                                                              |                                                                         
CGTGTAACGACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACG  >  NZ_CP009273/3896527‑3896694

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