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 I230 R1 226 18.8 1048726 96.7% 1014118 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,318,468 A→C A609A (GCT→GCG ftsH ← ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,318,4680AC100.0% 75.1 / NA 25A609A (GCT→GCGftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (14/11);  total (14/11)

GACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCAT  >  NZ_CP009273/3318382‑3318553
                                                                                      |                                                                                     
gACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAcgcg                                                                                    >  2:260546/1‑90 (MQ=255)
   aTGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 <  1:251061/90‑1 (MQ=255)
   aTGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 <  2:132365/90‑1 (MQ=255)
     ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  1:305208/1‑90 (MQ=255)
        gtTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttc                                                                            >  1:282203/1‑90 (MQ=255)
          tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  1:287405/1‑90 (MQ=255)
              cGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCa                                                                      >  1:228458/1‑90 (MQ=255)
               gggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCAGTGTCGCCAGAATTGTTAGACGCGCCTGGGTCTTCCCAg                                                                     >  2:179519/1‑90 (MQ=255)
                ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  1:132253/1‑90 (MQ=255)
                ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  1:341836/1‑90 (MQ=255)
                ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:36250/1‑90 (MQ=255)
                ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:370273/1‑90 (MQ=255)
                    cGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGc                                                                >  2:239268/1‑90 (MQ=255)
                         tACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCg                                                           >  2:372902/1‑90 (MQ=255)
                             cgGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACg                                                       <  2:132253/90‑1 (MQ=255)
                               gTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTa                                                     >  1:125724/1‑90 (MQ=255)
                                         cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  1:464525/90‑1 (MQ=255)
                                         cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  2:154046/90‑1 (MQ=255)
                                          ggACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGt                                          <  1:170460/90‑1 (MQ=255)
                                                     cTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtc                               >  2:62509/1‑90 (MQ=255)
                                                                   tGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTg                 <  2:160065/90‑1 (MQ=255)
                                                                          aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  1:99975/90‑1 (MQ=255)
                                                                          aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  1:505342/90‑1 (MQ=255)
                                                                          aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  2:452477/90‑1 (MQ=255)
                                                                                gTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGtctc    >  1:403706/1‑90 (MQ=255)
                                                                                  tAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCat  <  1:239268/90‑1 (MQ=255)
                                                                                      |                                                                                     
GACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCAT  >  NZ_CP009273/3318382‑3318553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ 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
GCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATATTTCATGAG  >  NZ_CP009273/3318377‑3318563
                                                                                           |                                                                                               
GCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGG                                                                                         >  SRR3722114.308462/1‑100 (MQ=60)
  TCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTT                                                                                       <  SRR3722114.24745/100‑1 (MQ=60)
   CTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTC                                                                                      >  SRR3722114.285165/1‑100 (MQ=60)
   CTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTC                                                                                      <  SRR3722114.323044/100‑1 (MQ=60)
     GACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTT                                                                                    >  SRR3722114.290445/1‑100 (MQ=60)
        ATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCC                                                                                 <  SRR3722114.253534/100‑1 (MQ=60)
         TGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCA                                                                                >  SRR3722114.230622/1‑100 (MQ=60)
           GTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGC                                                                              >  SRR3722114.133626/1‑100 (MQ=60)
           GTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGC                                                                              >  SRR3722114.345716/1‑100 (MQ=60)
                          GGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTA                                                               >  SRR3722114.127041/1‑100 (MQ=60)
                                              CGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGG                                           <  SRR3722114.470327/100‑1 (MQ=60)
                                               GGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGT                                          <  SRR3722114.172155/100‑1 (MQ=60)
                                                                           CGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTC              >  SRR3722114.408624/1‑100 (MQ=60)
                                                                               tgtataagagacAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATAT          <  SRR3722114.12269/88‑1 (MQ=60)
                                                                               AGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATAT          <  SRR3722114.101027/100‑1 (MQ=60)
                                                                               AGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATAT          <  SRR3722114.511674/100‑1 (MQ=60)
                                                                                       TAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATATTTCATGAG  <  SRR3722114.241551/100‑1 (MQ=60)
                                                                                           |                                                                                               
GCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATATTTCATGAG  >  NZ_CP009273/3318377‑3318563

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