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 2,526,050 A→G T95A (ACC→GCC)  cysK → cysteine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,526,0500AG100.0% 89.9 / NA 28T95A (ACC→GCC) cysKcysteine synthase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/11);  total (17/11)

GACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGC  >  NZ_CP009273/2525971‑2526135
                                                                               |                                                                                     
gACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTAt                                                                             >  2:433954/1‑90 (MQ=255)
 aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  >  1:401335/1‑84 (MQ=255)
 aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:401335/84‑1 (MQ=255)
            tACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCaa                                                                 <  2:26837/90‑1 (MQ=255)
                gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  >  1:229665/1‑69 (MQ=255)
                gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:229665/69‑1 (MQ=255)
                gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctg                                                             >  1:159974/1‑90 (MQ=255)
                 ggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAGCGCCGCAAgctgc                                                            >  1:509672/1‑90 (MQ=255)
                  gATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctgct                                                           >  1:570256/1‑90 (MQ=255)
                              cTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:531624/55‑1 (MQ=255)
                              cTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  >  1:531624/1‑55 (MQ=255)
                                     gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTGTTGAACGCCGCGAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  2:237191/1‑90 (MQ=255)
                                     gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  1:306078/1‑90 (MQ=255)
                                     gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  1:119124/1‑90 (MQ=255)
                                       tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGaa                                                                         >  1:503322/1‑55 (MQ=255)
                                       tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGaa                                                                         <  2:503322/55‑1 (MQ=255)
                                        gCCGCTCGCGGTTACAAATTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTg                                     >  1:252187/1‑90 (MQ=255)
                                          cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                   >  1:429402/1‑90 (MQ=255)
                                          cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                   <  2:306078/90‑1 (MQ=255)
                                                    tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                         <  2:252187/90‑1 (MQ=255)
                                                     aCAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAg                        <  2:570256/90‑1 (MQ=255)
                                                            cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 >  1:450580/1‑90 (MQ=255)
                                                            cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 <  1:376945/90‑1 (MQ=255)
                                                            cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 <  2:54450/90‑1 (MQ=255)
                                                            cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 >  1:314768/1‑90 (MQ=255)
                                                                  gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt           >  1:547497/1‑90 (MQ=255)
                                                                       tGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAg      >  1:531951/1‑90 (MQ=255)
                                                                            gAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGgcgc  <  2:520228/89‑1 (MQ=255)
                                                                               |                                                                                     
GACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGC  >  NZ_CP009273/2525971‑2526135

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
GCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCA  >  NZ_CP009273/2525952‑2526136
                                                                                                  |                                                                                      
GCGTTGAACTGGTTGACCCGCCCAGCGGTAATCCCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGC                                                                                       <  SRR3722111.361203/100‑1 (MQ=60)
      AACTGGTTGAACCGACCAGCGGTAATCCCGGGATTGCACGGGCCTATGTAGCTGCCGCGCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAG                                                                                 <  SRR3722111.544241/100‑1 (MQ=60)
          GGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATT                                                                             >  SRR3722111.407185/1‑100 (MQ=60)
                         CGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACctgtctcttat                                                              >  SRR3722111.232429/1‑89 (MQ=60)
                         CGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTG                                                              >  SRR3722111.161995/1‑100 (MQ=60)
                          GGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAGCGCCGCAAGCTGC                                                             >  SRR3722111.517563/1‑100 (MQ=60)
                           GTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCT                                                            >  SRR3722111.579219/1‑100 (MQ=60)
                                       TTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACctgtctcttatacacatctgacgct                                                >  SRR3722111.539958/1‑75 (MQ=60)
                                              GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAC                                         >  SRR3722111.310053/1‑100 (MQ=60)
                                              GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAC                                         >  SRR3722111.120627/1‑100 (MQ=60)
                                                CCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCctgtctcttatacacatctgacgct                                       >  SRR3722111.511090/1‑75 (MQ=60)
                                                 CTATGTAGCTGCCGCTCGCGGTTACAAATTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTG                                      >  SRR3722111.255153/1‑100 (MQ=60)
                                                   ATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGT                                    >  SRR3722111.435800/1‑100 (MQ=60)
                                                                     GTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAA                  >  SRR3722111.318907/1‑100 (MQ=60)
                                                                     GTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAA                  >  SRR3722111.457416/1‑100 (MQ=60)
                                                                           AACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAT            >  SRR3722111.556105/1‑100 (MQ=60)
                                                                               CACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA        <  SRR3722111.382330/100‑1 (MQ=60)
                                                                                ACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAG       >  SRR3722111.540292/1‑100 (MQ=60)
                                                                                     GACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCA  >  SRR3722111.528358/1‑100 (MQ=60)
                                                                                                  |                                                                                      
GCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCA  >  NZ_CP009273/2525952‑2526136

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