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 1,941,604 A→C G10G (GGT→GGG ruvC ← crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,941,6040AC100.0% 62.9 / NA 21G10G (GGT→GGGruvCcrossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/10);  total (11/10)

TTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941524‑1941686
                                                                                |                                                                                  
ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 >  2:457163/1‑84 (MQ=255)
ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 <  1:457163/84‑1 (MQ=255)
ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCa                                                                           >  2:222306/1‑90 (MQ=255)
ttGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCa                                                                           >  1:427277/1‑90 (MQ=255)
   gTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATg                                                                        >  2:123394/1‑90 (MQ=255)
      cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  1:130496/90‑1 (MQ=255)
      cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  2:466024/90‑1 (MQ=255)
      cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg                                                                     <  1:527451/90‑1 (MQ=255)
       ggATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCga                                                                    <  2:534335/90‑1 (MQ=255)
        gATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCgag                                                                   >  2:577300/1‑90 (MQ=255)
        gATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCgag                                                                   <  2:427277/90‑1 (MQ=255)
                    aCCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCAt                                                       >  1:171419/1‑90 (MQ=255)
                              gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCgag                                                                   >  2:39296/1‑68 (MQ=255)
                              gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCgag                                                                   <  1:39296/68‑1 (MQ=255)
                              gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtctc                                             <  1:135567/90‑1 (MQ=255)
                              gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGtctc                                             >  1:527917/1‑90 (MQ=255)
                                  gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt                                         >  1:19387/1‑90 (MQ=255)
                                       ctacctacCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 <  1:536858/45‑1 (MQ=255)
                                       ctacctacCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAcccc                                                                                 >  2:536858/1‑45 (MQ=255)
                                       ctacctacCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCt                                    >  2:496486/1‑90 (MQ=255)
                                               ctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCa                            >  2:243547/1‑90 (MQ=255)
                                                tGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  2:150127/90‑1 (MQ=255)
                                                tGGAGGTTGACGCCGTCACCGGTCGCTCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAg                           <  1:126298/90‑1 (MQ=255)
                                                                    ggTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTTGCTg       >  1:58583/1‑90 (MQ=255)
                                                                         cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc  <  2:469404/90‑1 (MQ=255)
                                                                                |                                                                                  
TTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC  >  NZ_CP009273/1941524‑1941686

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 33 ≤ 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
ATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACC  >  NZ_CP009273/1941514‑1941703
                                                                                          |                                                                                                   
ATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCA                                                                                            >  SRR3722111.433634/1‑100 (MQ=60)
    attcgtcggcagcgtcagatgtgtataagagacagGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGC                                                                                        <  SRR3722111.545293/65‑1 (MQ=60)
    TCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGC                                                                                        <  SRR3722111.464104/100‑1 (MQ=60)
                CGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAG                                                                            <  SRR3722111.132131/100‑1 (MQ=60)
                CGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAG                                                                            <  SRR3722111.535728/100‑1 (MQ=60)
                  tataagagacagACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCC                                                                          <  SRR3722111.39769/88‑1 (MQ=60)
                    TGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCAT                                                                        >  SRR3722111.173560/1‑100 (MQ=60)
                              ACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTC                                                              >  SRR3722111.536204/1‑100 (MQ=60)
                                  AGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTT                                                          >  SRR3722111.19622/1‑100 (MQ=60)
                                        GACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTG                                                    <  SRR3722111.137273/100‑1 (MQ=60)
                                                          TGGAGGTTGACGCCGTCACCGGTCGCTCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAG                                  <  SRR3722111.127881/100‑1 (MQ=60)
                                                                    CGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTTGCTG                        >  SRR3722111.59284/1‑100 (MQ=60)
                                                                                   CGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGA         >  SRR3722111.437842/1‑100 (MQ=60)
                                                                                          CCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACC  >  SRR3722111.544212/1‑100 (MQ=60)
                                                                                          |                                                                                                   
ATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACC  >  NZ_CP009273/1941514‑1941703

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