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 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,634,151 Δ1 bp coding (164/534 nt) rrrQ ← lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,634,1500G.100.0% 129.8 / NA 31coding (165/534 nt)rrrQlysozyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (11/20);  total (11/20)

CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634062‑1634235
                                                                                        |                                                                                     
ccTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgc                                                                                      <  2:952038/90‑2 (MQ=255)
ccTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgc                                                                                      <  1:366518/90‑2 (MQ=255)
  ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    >  1:156041/1‑87 (MQ=255)
  ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    <  2:2335/90‑4 (MQ=255)
         gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACa                                                                            >  2:698269/1‑90 (MQ=255)
           tCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGa                                                                          <  1:941171/90‑1 (MQ=255)
            cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGAt                                                                         <  1:343323/90‑1 (MQ=255)
              aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATgg                                                                       <  1:547261/90‑1 (MQ=255)
              aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATgg                                                                       <  2:399435/90‑1 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  <  2:673535/90‑1 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  <  1:317564/90‑1 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  <  2:205421/90‑1 (MQ=255)
                      gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATg                                                               >  2:497951/1‑90 (MQ=255)
                           tgTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:995055/2‑90 (MQ=255)
                           ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGa                                                          <  2:590144/90‑1 (MQ=255)
                           ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:609504/1‑90 (MQ=255)
                           ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:753755/1‑90 (MQ=255)
                             tCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAAc                                                        <  1:833856/90‑1 (MQ=255)
                                   ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  1:118621/90‑1 (MQ=255)
                                   ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  2:411689/90‑1 (MQ=255)
                                     ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  1:828740/90‑1 (MQ=255)
                                             aCAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCAt                                        <  1:448748/90‑1 (MQ=255)
                                             aCAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCAt                                        <  1:547471/90‑1 (MQ=255)
                                               aGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtg                                      <  1:109342/90‑1 (MQ=255)
                                                       atTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTa                              <  1:455590/90‑1 (MQ=255)
                                                          gggAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  1:417920/1‑90 (MQ=255)
                                                          gggAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  2:759916/1‑90 (MQ=255)
                                                          gggAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  1:927112/1‑90 (MQ=255)
                                                               aaaCGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCttt                      >  1:671691/1‑90 (MQ=255)
                                                                 aCGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCttttt                    <  2:675431/90‑1 (MQ=255)
                                                                            atccaCCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         <  1:698269/90‑1 (MQ=255)
                                                                              ccaccaCTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTg       >  2:335281/1‑90 (MQ=255)
                                                                                 ccacTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTc    >  2:514336/1‑90 (MQ=255)
                                                                                   acTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCaa  <  1:524919/90‑1 (MQ=255)
                                                                                   acTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCaa  <  2:885881/90‑1 (MQ=255)
                                                                                        |                                                                                     
CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634062‑1634235

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGG  >  NZ_CP009273/1634052‑1634245
                                                                                                   |                                                                                              
CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGC                                                                                                <  SRR3722236.170969/100‑1 (MQ=60)
  TGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCAC                                                                                             <  SRR3722236.47472/100‑1 (MQ=60)
  TGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCAC                                                                                             >  SRR3722236.158317/1‑100 (MQ=60)
        CGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCAAC                                                                                       <  SRR3722236.502725/100‑1 (MQ=60)
          CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAG                                                                                     <  SRR3722236.371696/100‑1 (MQ=60)
                     TCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATAT                                                                          <  SRR3722236.957781/100‑1 (MQ=60)
                      CAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATG                                                                         <  SRR3722236.348208/100‑1 (MQ=60)
                        ATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCC                                                                       <  SRR3722236.555696/100‑1 (MQ=60)
                           GCGTTGACCTTGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGA                                                                    >  SRR3722236.1012574/1‑100 (MQ=60)
                           GCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGA                                                                    >  SRR3722236.766683/1‑100 (MQ=60)
                           GCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGA                                                                    >  SRR3722236.619247/1‑100 (MQ=60)
                             GTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAAC                                                                  <  SRR3722236.322119/100‑1 (MQ=60)
                                       TCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTA                                                        <  SRR3722236.848419/100‑1 (MQ=60)
                                             TTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCAT                                                  <  SRR3722236.120278/100‑1 (MQ=60)
                                               TTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTG                                                <  SRR3722236.843197/100‑1 (MQ=60)
                                                  ccttcgtcggcagcgtcagatgtgtataagagacagATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGT                                             <  SRR3722236.417575/64‑1 (MQ=60)
                                                       ACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTA                                        <  SRR3722236.455217/100‑1 (MQ=60)
                                                       ACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTA                                        <  SRR3722236.555910/100‑1 (MQ=60)
                                                         AGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACC                                      <  SRR3722236.110853/100‑1 (MQ=60)
                                                          GTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCT                                     >  SRR3722236.423764/1‑100 (MQ=60)
                                                          GTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCT                                     >  SRR3722236.943483/1‑100 (MQ=60)
                                                               ATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTT                                >  SRR3722236.682876/1‑100 (MQ=60)
                                                                 ATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTT                              <  SRR3722236.462208/100‑1 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         >  SRR3722236.40339/1‑100 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         <  SRR3722236.52817/100‑1 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         <  SRR3722236.710041/100‑1 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         >  SRR3722236.350179/1‑100 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         >  SRR3722236.848062/1‑100 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         >  SRR3722236.890075/1‑100 (MQ=60)
                                                                                      ATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGACATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTA         >  SRR3722236.694132/1‑100 (MQ=60)
                                                                                        CCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACctgtctcttata       >  SRR3722236.207255/1‑88 (MQ=60)
                                                                                        CCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATC       >  SRR3722236.134208/1‑100 (MQ=60)
                                                                                             ACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGG  <  SRR3722236.532945/100‑1 (MQ=60)
                                                                                                   |                                                                                              
CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGG  >  NZ_CP009273/1634052‑1634245

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