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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 354,025 T→C E126G (GAA→GGA)  cynR ← transcriptional regulator CynR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273354,0250TC100.0% 118.1 / NA 35E126G (GAA→GGA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/24);  total (11/24)

CTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAA  >  NZ_CP009273/353940‑354101
                                                                                     |                                                                            
cTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTg                                                                          <  1:727049/90‑1 (MQ=255)
 tCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGt                                                                         <  1:1171903/90‑1 (MQ=255)
  ccGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                           >  2:152316/1‑87 (MQ=255)
  ccGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCt                                                                           <  1:152316/87‑1 (MQ=255)
  ccGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTa                                                                        <  1:30726/90‑1 (MQ=255)
    ggCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGc                                                                      >  2:949797/1‑90 (MQ=255)
       gAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTgg                                                                   <  1:484658/90‑1 (MQ=255)
             acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                             <  1:240662/90‑1 (MQ=255)
             acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                             <  2:650607/90‑1 (MQ=255)
             acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                             <  2:1054248/90‑1 (MQ=255)
             acaGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTg                                                             <  2:546326/90‑1 (MQ=255)
                 gcgcGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGc                                                         <  1:224430/90‑1 (MQ=255)
                  cgcgAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCt                                                        <  1:351168/90‑1 (MQ=255)
                   gcgAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTg                                                       >  1:1238687/1‑90 (MQ=255)
                   gcgAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTg                                                       >  2:940766/1‑90 (MQ=255)
                   gcgAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTg                                                       >  2:636143/1‑90 (MQ=255)
                          cAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAgc                                                >  1:112804/1‑90 (MQ=255)
                           aaTCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAgcg                                               <  2:261566/90‑1 (MQ=255)
                            aTCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAgcgc                                              >  1:341096/1‑90 (MQ=255)
                              cccAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAgcgcgc                                            <  2:800510/90‑1 (MQ=255)
                                  aCGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGAGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAg                                        >  2:111824/1‑90 (MQ=255)
                                        aaCTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATcc                                  <  1:605247/90‑1 (MQ=255)
                                        aaCTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATcc                                  <  2:1230053/90‑1 (MQ=255)
                                        aaCTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATcc                                  <  2:802957/90‑1 (MQ=255)
                                           tcgtcgCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCGGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGcc                               <  1:952424/90‑1 (MQ=255)
                                               cgcgGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTa                           <  2:14678/90‑1 (MQ=255)
                                               cgcgGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTa                           <  1:234724/90‑1 (MQ=255)
                                                    cAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAgggg                      <  1:571900/90‑1 (MQ=255)
                                                      aaGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGcc                    <  2:527393/90‑1 (MQ=255)
                                                        gCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGa                  <  2:688836/90‑1 (MQ=255)
                                                                  gATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCt        >  2:1084801/1‑90 (MQ=255)
                                                                      ttCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTa    <  1:159503/90‑1 (MQ=255)
                                                                       tctcCTGCCACATTCCCTGTAGCGGGCGCGTGATGCTGGGATAGCGCGCCTAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaa   <  1:717531/90‑1 (MQ=255)
                                                                        ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaaa  >  1:1189850/1‑90 (MQ=255)
                                                                        ctcCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTaaa  >  1:1223117/1‑90 (MQ=255)
                                                                                     |                                                                            
CTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAA  >  NZ_CP009273/353940‑354101

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
GGAATTGCCTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGCGATACGC  >  NZ_CP009273/353926‑354123
                                                                                                   |                                                                                                  
GGAATTGCCTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTC                                                                                                    >  SRR3722089.1070070/1‑100 (MQ=60)
        CTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGC                                                                                            <  SRR3722089.1224874/100‑1 (MQ=60)
        CTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGC                                                                                            <  SRR3722089.340191/100‑1 (MQ=60)
             GCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAG                                                                                       <  SRR3722089.154209/100‑1 (MQ=60)
              CTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGC                                                                                      <  SRR3722089.736507/100‑1 (MQ=60)
               TCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCG                                                                                     <  SRR3722089.1189741/100‑1 (MQ=60)
                CCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGT                                                                                    <  SRR3722089.31097/100‑1 (MQ=60)
                     GAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGC                                                                               <  SRR3722089.490411/100‑1 (MQ=60)
                       ATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTG                                                                             >  SRR3722089.1257537/1‑100 (MQ=60)
                           ACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGAT                                                                         <  SRR3722089.243721/100‑1 (MQ=60)
                              GGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGC                                                                      >  SRR3722089.114186/1‑100 (MQ=60)
                               GCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCG                                                                     <  SRR3722089.227278/100‑1 (MQ=60)
                                CGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGC                                                                    <  SRR3722089.355523/100‑1 (MQ=60)
                                CGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGC                                                                    >  SRR3722089.345347/1‑100 (MQ=60)
                                                      AACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGG                                              <  SRR3722089.612732/100‑1 (MQ=60)
                                                         TCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCGGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGC                                           <  SRR3722089.966414/100‑1 (MQ=60)
                                                             CGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGAT                                       <  SRR3722089.237706/100‑1 (MQ=60)
                                                                  CAAAGCATATCCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGT                                  <  SRR3722089.578672/100‑1 (MQ=60)
                                                                            CCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAA                        >  SRR3722089.1241705/1‑100 (MQ=60)
                                                                            CCTCGATTTTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAA                        >  SRR3722089.1207945/1‑100 (MQ=60)
                                                                                    TTCTCCTGCGACATTCCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGG                <  SRR3722089.161482/100‑1 (MQ=60)
                                                                                     TCTCCTGCCACATTCCCTGTAGCGGGCGCGTGATGCTGGGATAGCGCGCCTAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGT               <  SRR3722089.726860/100‑1 (MQ=60)
                                                                                                  TCCCTGTAGCTGGAGGGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGCGATACGC  <  SRR3722089.1038921/100‑1 (MQ=60)
                                                                                                   |                                                                                                  
GGAATTGCCTCCAGCTCCGGCGAATGCACAGGCGCGAAGGCAATCCCAACGTCCAACTCGTCGCGGCAAAGCATATCCTCGATTTTCTCCTGCGACATTTCCTGTAGCTGGAGCGTGATGCTGGGATAGCGCGCATAGAAATCCGCCATTAAGGGGCCGATAAAGTAGCTCGTAAAGGTGGGGGTGACGGCGATACGC  >  NZ_CP009273/353926‑354123

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