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 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,4470TC100.0% 79.2 / NA 26S52P (TCC→CCC) yjjQresponse regulator transcription factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/9);  total (17/9)

AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCAC  >  NZ_CP009273/4593364‑4593528
                                                                                   |                                                                                 
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:215095/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  2:175254/1‑90 (MQ=255)
agGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAt                                                                             >  1:44690/1‑90 (MQ=255)
            cTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTg                                                                 >  2:86580/1‑90 (MQ=255)
              cACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGAt                                                               >  1:76182/1‑90 (MQ=255)
               aCTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGAtt                                                              >  1:22583/1‑90 (MQ=255)
                   cccTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGc                                                          >  2:135740/1‑90 (MQ=255)
                         atatGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGa                                                    <  1:135740/90‑1 (MQ=255)
                            tGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAg                                                 >  2:227260/1‑90 (MQ=255)
                            tGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAg                                                 >  2:145412/1‑90 (MQ=255)
                             gAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGt                                                >  2:197704/1‑90 (MQ=255)
                                 taataTCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAg                                            <  2:88177/90‑1 (MQ=255)
                                 taataTCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAg                                            <  2:101004/90‑1 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  1:3536/1‑90 (MQ=255)
                                          gcgcCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGtt                                   >  1:286574/1‑90 (MQ=255)
                                               tctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGtt                              >  1:259433/1‑90 (MQ=255)
                                                ctGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGttt                             >  1:241993/1‑90 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  1:193233/90‑1 (MQ=255)
                                                      aggaCCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAAc                       <  2:65713/90‑1 (MQ=255)
                                                            ttaccttaTTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACattatt                 <  2:257987/90‑1 (MQ=255)
                                                                 ttattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCt            >  1:284907/1‑90 (MQ=255)
                                                                 ttattaCAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCt            >  1:154703/1‑90 (MQ=255)
                                                                       cAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtc      <  2:3536/90‑1 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc    <  2:76182/90‑1 (MQ=255)
                                                                         aTTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAAtctc    >  1:102976/1‑90 (MQ=255)
                                                                           tACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTcac  <  2:44690/90‑1 (MQ=255)
                                                                                   |                                                                                 
AGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCAC  >  NZ_CP009273/4593364‑4593528

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                                                                                             >  SRR3722113.216789/1‑100 (MQ=60)
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGAT                                                                                             >  SRR3722113.45067/1‑100 (MQ=60)
     AAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTC                                                                                        <  SRR3722113.417505/100‑1 (MQ=60)
              CATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGAT                                                                               >  SRR3722113.76852/1‑100 (MQ=60)
               ATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATT                                                                              >  SRR3722113.22780/1‑100 (MQ=60)
                                   ATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGAT                                                          <  SRR3722113.136877/100‑1 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                                   >  SRR3722113.289190/1‑100 (MQ=60)
                                          ATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTT                                                   >  SRR3722113.3573/1‑100 (MQ=60)
                                               ATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTT                                              >  SRR3722113.261710/1‑100 (MQ=60)
                                                TCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTT                                             >  SRR3722113.244073/1‑100 (MQ=60)
                                                                AGGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTC                             <  SRR3722113.194738/100‑1 (MQ=60)
                                                                 GGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCT                            >  SRR3722113.155985/1‑100 (MQ=60)
                                                                 GGACCTTACCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCT                            >  SRR3722113.287482/1‑100 (MQ=60)
                                                                         CCTTATTACAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTC                    >  SRR3722113.103898/1‑100 (MQ=60)
                                                                                 CAATTACGTCGTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGGGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATC            >  SRR3722113.247296/1‑100 (MQ=60)
                                                                                           GTCCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  SRR3722113.208466/1‑100 (MQ=60)
                                                                                             |                                                                                                 
GGGTTAAAAGAGGTCATGAGGACTCACTTCCCTGAATATGAAATAATATCCAGCGCCTCTGCGGAGGACCTTACCTTATTACAATTACGTCGTTCCGGATTAGTCATTGCTGATTTAGCCGGTGAAAGTGAAGATCCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCA  >  NZ_CP009273/4593354‑4593544

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