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 I229 R1 214 17.4 943020 97.5% 919444 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,901,121 A→G D202D (GAT→GAC pstB ← phosphate ABC transporter ATP‑binding protein PstB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,901,1210AG100.0% 52.6 / NA 18D202D (GAT→GACpstBphosphate ABC transporter ATP‑binding protein PstB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/9);  total (9/9)

TCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTC  >  NZ_CP009273/3901041‑3901192
                                                                                |                                                                       
tCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTc                                                                <  1:357255/90‑1 (MQ=255)
  gCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAg                                                              <  1:369886/90‑1 (MQ=255)
       ggTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCgg                                                         >  2:284090/1‑90 (MQ=255)
                     gtgtGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCa                                           >  2:215550/1‑90 (MQ=255)
                     gtgtGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGAAGCCCTGCTTCAGTTCGGGGATCAGCTCTTCa                                           >  1:93763/1‑90 (MQ=255)
                            cGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGcc                                    >  2:37653/1‑90 (MQ=255)
                             ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc                                   <  2:25495/90‑1 (MQ=255)
                              gAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCCCCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCg                                  <  2:305751/90‑1 (MQ=255)
                               aacaacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCgg                                 >  1:363832/1‑90 (MQ=255)
                               aacaacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCgg                                 >  1:51103/1‑90 (MQ=255)
                                 caacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                               <  2:100245/90‑1 (MQ=255)
                                 caacGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                               <  2:327971/90‑1 (MQ=255)
                                   ccgcgCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                               >  2:374775/2‑88 (MQ=255)
                                   ccgcgCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                               <  1:374775/87‑1 (MQ=255)
                                                 gcATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAgcg               <  2:363832/90‑1 (MQ=255)
                                                              gACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTtc  >  1:79763/1‑90 (MQ=255)
                                                                           gtgtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGcc                                    <  1:53623/43‑1 (MQ=255)
                                                                           gtgtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGcc                                    >  2:53623/1‑43 (MQ=255)
                                                                                |                                                                       
TCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTC  >  NZ_CP009273/3901041‑3901192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 33 ≤ 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
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC  >  NZ_CP009273/3901023‑3901215
                                                                                                  |                                                                                              
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGT                                                                                               >  SRR3722113.141443/1‑100 (MQ=60)
  GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                             <  SRR3722113.135941/100‑1 (MQ=60)
  GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                             <  SRR3722113.405690/100‑1 (MQ=60)
       ACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTT                                                                                        <  SRR3722113.267707/100‑1 (MQ=60)
        CTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTC                                                                                       <  SRR3722113.182991/100‑1 (MQ=60)
                  TCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGA                                                                             <  SRR3722113.360910/100‑1 (MQ=60)
                    GCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATC                                                                           <  SRR3722113.373725/100‑1 (MQ=60)
                             CATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGAAGCCCTGCTTCAGTTCGGGGATCAGCTCTTCA                                                                  >  SRR3722113.94607/1‑100 (MQ=60)
                                       GTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGG                                                        >  SRR3722113.367576/1‑100 (MQ=60)
                                       GTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGG                                                        >  SRR3722113.51528/1‑100 (MQ=60)
                                          cgagattcgccttagtctcgtgggctcggagatgtgtataagagacagACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAG                                                     <  SRR3722113.229806/52‑1 (MQ=60)
                                              ttagtctcgtgggctcggagatgtgtataagagacagGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGAT                                                 <  SRR3722113.54077/63‑1 (MQ=60)
                                                   CACCGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGT                                            <  SRR3722113.378663/100‑1 (MQ=60)
                                                                      TGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTC                         >  SRR3722113.80466/1‑100 (MQ=60)
                                                                            GGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAActgtctcttatacacatctccgagcccacgagactaaggcgaatct                   >  SRR3722113.317296/1‑54 (MQ=60)
                                                                                        CCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTC       >  SRR3722113.99427/1‑100 (MQ=60)
                                                                                             GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC  >  SRR3722113.365637/1‑100 (MQ=60)
                                                                                             GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC  >  SRR3722113.24170/1‑100 (MQ=60)
                                                                                             GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC  >  SRR3722113.201882/1‑100 (MQ=60)
                                                                                             GTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC  >  SRR3722113.189740/1‑100 (MQ=60)
                                                                                                  |                                                                                              
TTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGC  >  NZ_CP009273/3901023‑3901215

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