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 I227 R1 223 21.3 1167702 97.0% 1132670 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,593,447 T→C S52P (TCC→CCC)  yjjQ → response regulator transcription factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,593,4470TC88.5% 60.8 / 2.9 26S52P (TCC→CCC) yjjQresponse regulator transcription factor
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (16/7);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAA  >  NZ_CP009273/4593364‑4593530
                                                                                   |                                                                                   
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                               >  1:450612/1‑90 (MQ=255)
 ggTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAtt                                                                              >  1:317977/1‑90 (MQ=255)
  gTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTa                                                                             >  2:398051/1‑90 (MQ=255)
                       gAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGt                                                        <  2:450090/90‑1 (MQ=255)
                        aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTg                                                       >  2:66306/1‑90 (MQ=255)
                        aaTATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGACTAGTCATTGCTGATTTAGCCGGTg                                                       >  1:295866/1‑90 (MQ=255)
                          tatGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGaa                                                     <  2:226605/90‑1 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                     >  2:320824/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                     >  1:230802/1‑90 (MQ=255)
                                           cgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTc                                    <  1:169543/90‑1 (MQ=255)
                                           cgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTc                                    <  1:141193/90‑1 (MQ=255)
                                            gcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   <  1:289140/90‑1 (MQ=255)
                                            gcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   >  2:297860/1‑90 (MQ=255)
                                               tctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGtt                                >  2:359104/1‑90 (MQ=255)
                                                       ggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACa                        >  2:107388/1‑90 (MQ=255)
                                                          ccttaccttaTTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   >  2:224056/1‑76 (MQ=255)
                                                          ccttaccttaTTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCt                                   <  1:224056/76‑1 (MQ=255)
                                                                 ttattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCt              >  1:295496/1‑90 (MQ=255)
                                                                 ttattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCt              >  2:514749/1‑90 (MQ=255)
                                                                  tattaCAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCtt             >  1:233673/1‑90 (MQ=255)
                                                                   attaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCttt            >  1:49085/1‑90 (MQ=255)
                                                                        aaTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtct       >  1:200724/1‑90 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc      <  1:276918/90‑1 (MQ=255)
                                                                             cgtcgtCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACaa  >  2:137427/1‑90 (MQ=255)
                                                                             cgtcgtCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACaa  >  2:204917/1‑90 (MQ=255)
                                                                             cgtcgtCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACaa  <  1:221464/90‑1 (MQ=255)
                                                                                   |                                                                                   
AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAA  >  NZ_CP009273/4593364‑4593530

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  NZ_CP009273/4593354‑4593544
                                                                                             |                                                                                                 
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722111.457448/1‑100 (MQ=60)
 GGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATT                                                                                            <  SRR3722111.450591/100‑1 (MQ=60)
 GGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATT                                                                                            >  SRR3722111.322179/1‑100 (MQ=60)
     AAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTC                                                                                        <  SRR3722111.335645/100‑1 (MQ=60)
                        CACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGACTAGTCATTGCTGATTTAGCCGGTG                                                                     >  SRR3722111.299664/1‑100 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                                   >  SRR3722111.233580/1‑100 (MQ=60)
                                                     CGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAA                                        <  SRR3722111.171664/100‑1 (MQ=60)
                                                     CGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAA                                        <  SRR3722111.142975/100‑1 (MQ=60)
                                                      acagCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAC                                       <  SRR3722111.226770/96‑1 (MQ=60)
                                                      GCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAC                                       <  SRR3722111.292827/100‑1 (MQ=60)
                                                                 GGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCT                            >  SRR3722111.299285/1‑100 (MQ=60)
                                                                  GACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTT                           >  SRR3722111.236468/1‑100 (MQ=60)
                                                                   ACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTT                          >  SRR3722111.49679/1‑100 (MQ=60)
                                                                        ACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCT                     >  SRR3722111.203220/1‑100 (MQ=60)
                                                                                   ATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGG          <  SRR3722111.280376/100‑1 (MQ=60)
                                                                                       CGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAA      <  SRR3722111.224146/100‑1 (MQ=60)
                                                                                           GTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  SRR3722111.42024/1‑100 (MQ=60)
                                                                                             |                                                                                                 
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  NZ_CP009273/4593354‑4593544

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