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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,613,222 C→G R238P (CGC→CCC)  lplA ← lipoate‑‑protein ligase LplA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,613,2220CG100.0% 50.1 / NA 15R238P (CGC→CCC) lplAlipoate‑‑protein ligase LplA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (7/8);  total (7/8)

CCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT  >  NZ_CP009273/4613139‑4613303
                                                                                   |                                                                                 
ccACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCaa                                                                             >  1:681718/1‑90 (MQ=255)
                 aaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGttt                                                            >  2:911671/1‑90 (MQ=255)
                     gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                        >  1:17660/1‑90 (MQ=255)
                     gTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCa                                                        <  2:861128/90‑1 (MQ=255)
                       tCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAg                                                      >  1:484491/1‑90 (MQ=255)
                       tCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAg                                                      >  1:497058/1‑90 (MQ=255)
                              gcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCg                                               >  1:469386/1‑90 (MQ=255)
                                 gATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGttt                                            <  1:77777/90‑1 (MQ=255)
                                   tGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTg                                          >  1:730839/1‑90 (MQ=255)
                                      gagaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTc                                       <  2:469386/90‑1 (MQ=255)
                                                               tCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCAcgc              <  1:182274/90‑1 (MQ=255)
                                                                 cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc            <  1:614785/90‑1 (MQ=255)
                                                                 cATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACgcgc            <  1:882986/90‑1 (MQ=255)
                                                                          cTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATaa   <  2:17660/90‑1 (MQ=255)
                                                                           tACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAt  <  1:470072/90‑1 (MQ=255)
                                                                                   |                                                                                 
CCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT  >  NZ_CP009273/4613139‑4613303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG  >  NZ_CP009273/4613124‑4613314
                                                                                                  |                                                                                            
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGG                                                                                             <  SRR3722072.213784/100‑1 (MQ=60)
     AAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAA                                                                                        >  SRR3722072.693913/1‑100 (MQ=60)
                          CAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCA                                                                   >  SRR3722072.17989/1‑100 (MQ=60)
                            GGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAG                                                                 >  SRR3722072.493261/1‑100 (MQ=60)
                            GGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAG                                                                 >  SRR3722072.506032/1‑100 (MQ=60)
                                   CGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCG                                                          >  SRR3722072.477940/1‑100 (MQ=60)
                                        ATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTG                                                     >  SRR3722072.743803/1‑100 (MQ=60)
                                                GATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGA                                             <  SRR3722072.79220/100‑1 (MQ=60)
                                                                              TCCATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATA               <  SRR3722072.185429/100‑1 (MQ=60)
                                                                                ataagagacagACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722072.502273/89‑1 (MQ=60)
                                                                                agacagCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722072.897977/94‑1 (MQ=60)
                                                                                CATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722072.625749/100‑1 (MQ=60)
                                                                                CATTCCCAGCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAAT             <  SRR3722072.898189/100‑1 (MQ=60)
                                                                                        GCTACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAA     >  SRR3722072.655271/1‑100 (MQ=60)
                                                                                          TACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAA   <  SRR3722072.478637/100‑1 (MQ=60)
                                                                                           ACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATctgtctcttat  >  SRR3722072.433743/1‑89 (MQ=60)
                                                                                           ACTCTGGGGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG  >  SRR3722072.366668/1‑100 (MQ=60)
                                                                                                  |                                                                                            
CGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAG  >  NZ_CP009273/4613124‑4613314

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