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,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% 66.8 / NA 21D202D (GAT→GACpstBphosphate ABC transporter ATP‑binding protein PstB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/12);  total (9/12)

ATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAG  >  NZ_CP009273/3901035‑3901204
                                                                                      |                                                                                   
atcaatTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTcc                                                                                  >  2:381231/1‑90 (MQ=255)
     ttCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCtt                                                                             >  1:29835/1‑90 (MQ=255)
      tCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTc                                                                            <  2:268087/90‑1 (MQ=255)
           cAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTc                                                                       <  2:497813/90‑1 (MQ=255)
                     aaCGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGc                                                             <  2:492766/90‑1 (MQ=255)
                     aaCGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGc                                                             >  1:258384/1‑90 (MQ=255)
                           gtgtGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCa                                                       >  1:357235/1‑90 (MQ=255)
                                   ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCCCCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc                                               <  1:150434/90‑1 (MQ=255)
                                   ggAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc                                               <  2:223965/90‑1 (MQ=255)
                                                     ctgcATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAg                             <  1:44853/90‑1 (MQ=255)
                                                            ttgtGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc                                               <  1:471688/65‑1 (MQ=255)
                                                            ttgtGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGccc                                               >  2:471688/1‑65 (MQ=255)
                                                              gtgGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAaca                    >  2:464776/1‑90 (MQ=255)
                                                                gggTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGa                              >  1:432500/1‑78 (MQ=255)
                                                                gggTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGa                              <  2:432500/78‑1 (MQ=255)
                                                                    gACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                                           <  2:457867/61‑1 (MQ=255)
                                                                    gACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                                           >  1:457867/1‑61 (MQ=255)
                                                                            ccacGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagca      <  2:357235/90‑1 (MQ=255)
                                                                                ggTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGagcagcag  <  2:29222/90‑1 (MQ=255)
                                                                                 gtgtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                                           >  2:371291/1‑48 (MQ=255)
                                                                                 gtgtAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTa                                           <  1:371291/48‑1 (MQ=255)
                                                                                      |                                                                                   
ATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAG  >  NZ_CP009273/3901035‑3901204

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  NZ_CP009273/3901025‑3901220
                                                                                                |                                                                                                   
tacactaagatggtcgtcggcagcgtcagatgtgtataagagacaGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722092.230941/55‑1 (MQ=60)
gagatctaccctaagatggtcgtcggcagcgtcagatgtgtataagagacagCGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722092.288054/48‑1 (MQ=60)
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722092.382304/100‑1 (MQ=60)
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722092.199809/100‑1 (MQ=60)
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGAAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCC                                                                                                  <  SRR3722092.29184/100‑1 (MQ=60)
     ACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTT                                                                                             >  SRR3722092.30352/1‑100 (MQ=60)
                     CAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGC                                                                             >  SRR3722092.262611/1‑100 (MQ=60)
                           CATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCA                                                                       >  SRR3722092.363792/1‑100 (MQ=60)
                                             gtgtataagagacagCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA                                                     <  SRR3722092.481048/85‑1 (MQ=60)
                                             GGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCCCCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATA                                                     <  SRR3722092.152965/100‑1 (MQ=60)
                                                 cgtcggcagcgtcagatgtgtataagagacagGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGT                                                 <  SRR3722092.378285/68‑1 (MQ=60)
                                                               CTGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACAC                                   <  SRR3722092.45617/100‑1 (MQ=60)
                                                                TGCATGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACT                                  >  SRR3722092.440932/1‑100 (MQ=60)
                                                                    TGTTGTGGGTGACGATCACCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCtgtctcttatacacatct                              >  SRR3722092.466911/1‑82 (MQ=60)
                                                                                      CCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACctgtctcttatacacatctgacg            >  SRR3722092.247964/1‑77 (MQ=60)
                                                                                      CCACGGTGTAGTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTC            >  SRR3722092.502236/1‑100 (MQ=60)
                                                                                                GTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  SRR3722092.160877/1‑100 (MQ=60)
                                                                                                GTCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  SRR3722092.1406/1‑100 (MQ=60)
                                                                                                |                                                                                                   
GCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAGCCTGCTGCATGTTGTGGGTGACGATCACCACGGTGTAATCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATG  >  NZ_CP009273/3901025‑3901220

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