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 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,4680AC96.0% 72.7 / ‑3.4 25A609A (GCT→GCGftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (14/10);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCAT  >  NZ_CP009273/3318382‑3318553
                                                                                      |                                                                                     
gACATGGTGTTACCCGGGTTCGGCGTGCTCGGTTCATCAACCGGACGAGGAGCCTTTGGACTTCCATTGTCGCCAGAATTGTTAGAcgcg                                                                                    >  1:250518/1‑90 (MQ=255)
   aTGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 <  1:192292/90‑1 (MQ=255)
   aTGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCt                                                                                 <  1:573669/90‑1 (MQ=255)
     ggTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTgg                                                                               >  2:311304/1‑90 (MQ=255)
          tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  2:373415/1‑90 (MQ=255)
          tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGttctt                                                                          >  2:190361/1‑90 (MQ=255)
          tACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGttctt                                                                          >  1:352523/1‑90 (MQ=255)
                ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:99491/1‑90 (MQ=255)
                ggTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGc                                                                    >  2:350540/1‑90 (MQ=255)
                    cGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGc                                                                >  2:18710/1‑90 (MQ=255)
                                ttCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTAc                                                    >  2:166141/1‑90 (MQ=255)
                                  catcaACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACAt                                                  >  2:295816/1‑90 (MQ=255)
                                     caACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATcgc                                               >  2:185150/1‑90 (MQ=255)
                                         cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTTTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTGCATCGCGACg                                           >  1:66390/1‑90 (MQ=255)
                                         cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           <  1:373415/90‑1 (MQ=255)
                                         cGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACg                                           >  2:393699/1‑90 (MQ=255)
                                              gaggagCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCa                                      <  1:250885/90‑1 (MQ=255)
                                                ggagCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATc                                    <  2:346779/90‑1 (MQ=255)
                                                ggagCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATc                                    <  1:311304/90‑1 (MQ=255)
                                                     cTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGtc                               >  2:566635/1‑90 (MQ=255)
                                                                    gTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGc                >  1:360238/1‑90 (MQ=255)
                                                                      cGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGt              >  2:536330/1‑90 (MQ=255)
                                                                          aGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGAt          <  1:315204/90‑1 (MQ=255)
                                                                                gTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGtctc    <  1:349756/90‑1 (MQ=255)
                                                                                  tAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCat  <  2:505097/90‑1 (MQ=255)
                                                                                  tAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCat  <  2:505098/90‑1 (MQ=255)
                                                                                      |                                                                                     
GACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCAT  >  NZ_CP009273/3318382‑3318553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ 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
TAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATATTTCATGAGAGC  >  NZ_CP009273/3318372‑3318566
                                                                                                |                                                                                                  
TAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTGCTCGGTTCATCAACCGGACGAGGAGCCTTTGGACTTCCATTGTCGCCAGAATTGTTAGACGCG                                                                                                 >  SRR3722111.253464/1‑100 (MQ=60)
          GACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTT                                                                                       >  SRR3722111.357338/1‑100 (MQ=60)
             ATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCC                                                                                    <  SRR3722111.194693/100‑1 (MQ=60)
             ATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCC                                                                                    <  SRR3722111.582672/100‑1 (MQ=60)
                                         GTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTTTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTGCATCGCGACG                                                        >  SRR3722111.67174/1‑100 (MQ=60)
                                                   CGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGG                                              <  SRR3722111.378715/100‑1 (MQ=60)
                                                        GAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATC                                         <  SRR3722111.253835/100‑1 (MQ=60)
                                                          GGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAA                                       <  SRR3722111.315390/100‑1 (MQ=60)
                                                                    GACTACCATTGTCGCCAGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGC                             >  SRR3722111.365226/1‑100 (MQ=60)
                                                                                    AGAATTGTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATAT             <  SRR3722111.319348/100‑1 (MQ=60)
                                                                                          GTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCctgtctcttatacacatctgacgctgccgacgaatgtcaatgtgtagatctcg       >  SRR3722111.532725/1‑47 (MQ=60)
                                                                                          GTTAGACGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATATTTCATG       <  SRR3722111.354497/100‑1 (MQ=60)
                                                                                               ACGCGCCTGGTTCTTCCCGGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCAACAATCTGCGGTGCGTCGATAGTCTCATATTTCATGAGAGC  <  SRR3722111.533693/100‑1 (MQ=60)
                                                                                                |                                                                                                  
TAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAGAATTGTTAGAAGCGCCTGGTTCTTCCCAGCCCGCTGGCGGACGTACATCGCGACGTGCCATCAGGTCATCAATCTGCGGTGCGTCGATAGTCTCATATTTCATGAGAGC  >  NZ_CP009273/3318372‑3318566

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