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 2,554,121 A→C N169T (AAC→ACC)  BW25113_RS12740 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,554,1210AC100.0% 131.7 / NA 37N169T (AAC→ACC) BW25113_RS12740hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (14/23);  total (14/23)

TCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCA  >  NZ_CP009273/2554037‑2554203
                                                                                    |                                                                                  
tCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACgg                                                                               >  2:739546/1‑90 (MQ=255)
tCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACgg                                                                               >  1:497811/1‑90 (MQ=255)
  gCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGtt                                                                             >  2:357008/1‑90 (MQ=255)
   cTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTa                                                                            <  1:275905/90‑1 (MQ=255)
               aaaaCTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCg                                                                >  1:404809/1‑90 (MQ=255)
                aaaCTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGa                                                               <  2:393990/90‑1 (MQ=255)
                           tGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCa                                                    >  1:874939/1‑90 (MQ=255)
                              gATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGCCTGGTCAAcc                                                 <  1:634888/90‑1 (MQ=255)
                                 ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt                                              >  1:198109/1‑90 (MQ=255)
                                 ggTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACt                                              >  2:400331/1‑90 (MQ=255)
                                     cAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCaa                                          <  2:270682/90‑1 (MQ=255)
                                      aGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCaaa                                         <  1:659744/90‑1 (MQ=255)
                                       gggTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAAc                                        <  1:831790/90‑1 (MQ=255)
                                            ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  1:474095/1‑90 (MQ=255)
                                            ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  1:182250/1‑90 (MQ=255)
                                            ggTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgtt                                   >  1:3200/1‑90 (MQ=255)
                                                 gggAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCg                              >  1:47274/1‑90 (MQ=255)
                                                   gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgttg                                  >  1:448929/1‑84 (MQ=255)
                                                   gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACttgttg                                  <  2:448929/84‑1 (MQ=255)
                                                   gAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTa                            <  2:523917/90‑1 (MQ=255)
                                                    aaCAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTAc                           <  1:577814/90‑1 (MQ=255)
                                                        aTCGACAAACGAGCCAGGCTATCAAAAACCCCGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTCCTGaa                       <  2:3200/90‑1 (MQ=255)
                                                        aTCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGaa                       <  2:49478/90‑1 (MQ=255)
                                                         tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg                      <  2:391950/90‑1 (MQ=255)
                                                         tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg                      <  1:739830/90‑1 (MQ=255)
                                                         tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg                      <  1:400331/90‑1 (MQ=255)
                                                         tCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAg                      <  1:468091/90‑1 (MQ=255)
                                                             cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta                  <  1:382698/90‑1 (MQ=255)
                                                             cAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCta                  <  1:61571/90‑1 (MQ=255)
                                                              aaaCGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCtat                 <  2:198109/90‑1 (MQ=255)
                                                              aaaCGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCtat                 <  2:474095/90‑1 (MQ=255)
                                                              aaaCGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCtat                 <  2:483722/90‑1 (MQ=255)
                                                                  gAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATAc             <  1:739546/90‑1 (MQ=255)
                                                                          cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg     <  2:182250/90‑1 (MQ=255)
                                                                          cTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCg     <  2:661690/90‑1 (MQ=255)
                                                                            aTCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGc   >  1:824134/1‑90 (MQ=255)
                                                                             tCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCa  >  2:196544/1‑90 (MQ=255)
                                                                                    |                                                                                  
TCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCA  >  NZ_CP009273/2554037‑2554203

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
CTTTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAG  >  NZ_CP009273/2554027‑2554212
                                                                                              |                                                                                           
CTTTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGG                                                                                        >  SRR3722072.506803/1‑100 (MQ=60)
      TACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAA                                                                                  <  SRR3722072.550479/100‑1 (MQ=60)
             CTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCT                                                                           <  SRR3722072.280625/100‑1 (MQ=60)
               GGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCG                                                                         >  SRR3722072.412114/1‑100 (MQ=60)
                           AACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCA                                                             >  SRR3722072.890008/1‑100 (MQ=60)
                                 CAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACT                                                       >  SRR3722072.201525/1‑100 (MQ=60)
                                        GATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGCCTGGTCAACCACTTCAAACT                                                <  SRR3722072.646217/100‑1 (MQ=60)
                                            GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT                                            >  SRR3722072.185405/1‑100 (MQ=60)
                                            GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT                                            >  SRR3722072.482740/1‑100 (MQ=60)
                                            GTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTT                                            >  SRR3722072.3278/1‑100 (MQ=60)
                                                AGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGC                                        <  SRR3722072.671516/100‑1 (MQ=60)
                                                 GGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCG                                       <  SRR3722072.846216/100‑1 (MQ=60)
                                                 GGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCG                                       >  SRR3722072.48136/1‑100 (MQ=60)
                                                   GTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTA                                     >  SRR3722072.457149/1‑100 (MQ=60)
                                                              AACAATCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTAT                          <  SRR3722072.588134/100‑1 (MQ=60)
                                                                   TCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACC                     <  SRR3722072.407550/100‑1 (MQ=60)
                                                                   TCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACC                     <  SRR3722072.752934/100‑1 (MQ=60)
                                                                   TCGACAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACC                     <  SRR3722072.476632/100‑1 (MQ=60)
                                                                       CAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGA                 <  SRR3722072.62697/100‑1 (MQ=60)
                                                                       CAAACGAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGA                 <  SRR3722072.389539/100‑1 (MQ=60)
                                                                            GAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGC            <  SRR3722072.752646/100‑1 (MQ=60)
                                                                            GAGCCAGGCTATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGC            >  SRR3722072.838464/1‑100 (MQ=60)
                                                                                      ATCAAAAACCACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCtgtctcttataca  >  SRR3722072.353149/1‑87 (MQ=60)
                                                                                              |                                                                                           
CTTTGTTACTTCGCTGGTACTGGCAAAAACTCACAGTTGGGATGGTGCAGGGTTGGTAGGGGAACAATCGACAAACGAGCCAGGCTATCAAAAAACACGGTTAAAAAATGCCTCGATCGTCTGGTCAACCACTTCAAACTTGTTGAGCGTACTGAAGGCTATATACCAGGATCGGCAGAAAGGGAG  >  NZ_CP009273/2554027‑2554212

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