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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,003,480 T→C F60S (TTC→TCC)  yedK → SOS response‑associated peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,003,4800TC92.0% 60.5 / ‑0.7 25F60S (TTC→TCC) yedKSOS response‑associated peptidase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (11/12);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

CCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATG  >  NZ_CP009273/2003396‑2003559
                                                                                    |                                                                               
ccATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTggg                                                                            <  1:467677/90‑1 (MQ=255)
 cATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTgggg                                                                           >  1:229635/1‑90 (MQ=255)
           tACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTccc                                                                 >  2:83717/1‑90 (MQ=255)
               aCGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAt                                                             <  1:331048/90‑1 (MQ=255)
                 gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGAtgg                                                           <  1:362445/90‑1 (MQ=255)
                 gTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTTCCGGAtgg                                                           <  1:528143/90‑1 (MQ=255)
                  tCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggt                                                          <  1:411019/90‑1 (MQ=255)
                    gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGAtggtgg                                                        >  1:102196/1‑90 (MQ=255)
                    gcgcCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGAtggtgg                                                        >  1:342490/1‑90 (MQ=255)
                                 aaGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACcgccgc                                           <  1:503871/90‑1 (MQ=255)
                                             tCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGGGGGATAAACCGCCGCTGATTAACGccc                               >  1:401706/1‑90 (MQ=255)
                                             ccAGTGAACGTGATGAACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGccc                               <  2:681469/89‑1 (MQ=255)
                                                gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                            >  2:772367/1‑90 (MQ=255)
                                                gTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCcgcg                            >  2:39489/1‑90 (MQ=255)
                                                     cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                       >  1:652822/1‑90 (MQ=255)
                                                     cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                       <  1:39140/90‑1 (MQ=255)
                                                     cGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGaa                       >  1:350685/1‑90 (MQ=255)
                                                         atgaACCCCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTg                   <  2:187050/90‑1 (MQ=255)
                                                          tgaACCCCTTCATCTGGATCCGGTTTCCGGGGGATATGCCCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                  <  2:692081/90‑1 (MQ=255)
                                                          tgaACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGc                  >  2:605650/1‑90 (MQ=255)
                                                            aaCACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCgg                <  1:667597/90‑1 (MQ=255)
                                                                 cTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAAcc           <  1:605650/90‑1 (MQ=255)
                                                                       cTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGt     <  1:39489/90‑1 (MQ=255)
                                                                          gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg  >  2:11929/1‑90 (MQ=255)
                                                                          gATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATg  >  2:514693/1‑90 (MQ=255)
                                                                                    |                                                                               
CCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATG  >  NZ_CP009273/2003396‑2003559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 27 ≤ 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
ATCCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003387‑2003566
                                                                                             |                                                                                      
ATCCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGG                                                                                  >  SRR3722079.231585/1‑100 (MQ=60)
         CCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCC                                                                         <  SRR3722079.471767/100‑1 (MQ=60)
                   ATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGG                                                               >  SRR3722079.103018/1‑100 (MQ=60)
                   ATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGG                                                               >  SRR3722079.345195/1‑100 (MQ=60)
                        ACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAA                                                          <  SRR3722079.333688/100‑1 (MQ=60)
                          GTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAAC                                                        <  SRR3722079.365279/100‑1 (MQ=60)
                          GTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTTCCGGATGGTGGGATAAAC                                                        <  SRR3722079.533057/100‑1 (MQ=60)
                           TCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACC                                                       <  SRR3722079.414534/100‑1 (MQ=60)
                                          AAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGC                                        <  SRR3722079.508417/100‑1 (MQ=60)
                                            GTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGGGGGATAAACCGCCGCTGATTAACGCCC                                      >  SRR3722079.405095/1‑100 (MQ=60)
                                                    GCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAA                              >  SRR3722079.353446/1‑100 (MQ=60)
                                                    GCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAA                              >  SRR3722079.659581/1‑100 (MQ=60)
                                                              CGTGATGAACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAA                    <  SRR3722079.39433/100‑1 (MQ=60)
                                                                     AACACCTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCG             <  SRR3722079.674562/100‑1 (MQ=60)
                                                                          CTTCATCTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGT        <  SRR3722079.611758/100‑1 (MQ=60)
                                                                                CTGGATCCGGTTTCCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  <  SRR3722079.39783/100‑1 (MQ=60)
                                                                                             |                                                                                      
ATCCCGAACCCATTGGCAGATACAACGTCGCGCCGGGAACCAAAGTCCTACTGCTCAGTGAACGTGATGAACACCTTCATCTGGATCCGGTTTTCTGGGGATATGCTCCCGGATGGTGGGATAAACCGCCGCTGATTAACGCCCGCGTAGAAACTGCGGCAACCAGTCGTATGTTTAAAC  >  NZ_CP009273/2003387‑2003566

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