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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 568,752 A→G G71G (GGA→GGG ybcO → DUF1364 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273568,7520AG100.0% 79.4 / NA 25G71G (GGA→GGGybcODUF1364 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/18);  total (7/18)

TACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568665‑568838
                                                                                       |                                                                                      
tccGAACCCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGgta                                                                                      <  1:10584/88‑1 (MQ=255)
          ccAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGaa                                                                            <  1:168138/90‑1 (MQ=255)
                 tGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgcgc                                                                     >  2:528372/1‑90 (MQ=255)
                     tGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGa                                                                 <  2:352286/90‑1 (MQ=255)
                       ccaccaTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAg                                                               <  2:30247/90‑1 (MQ=255)
                               gCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  1:379107/90‑1 (MQ=255)
                                    ttCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAac                                                  <  2:388438/90‑1 (MQ=255)
                                    ttCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAac                                                  >  2:3063/1‑90 (MQ=255)
                                        gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  1:110440/1‑90 (MQ=255)
                                        gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              <  2:442490/90‑1 (MQ=255)
                                           tgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAg                                                               <  2:493995/70‑1 (MQ=255)
                                           tgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAg                                                               >  1:493995/1‑70 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  1:25180/90‑1 (MQ=255)
                                                          gACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAagag                            <  2:218614/90‑1 (MQ=255)
                                                           aCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAgg                           <  1:528372/90‑1 (MQ=255)
                                                            ccgccgCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAggg                          <  2:274545/90‑1 (MQ=255)
                                                                         tttGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCga                                                      <  2:149780/49‑1 (MQ=255)
                                                                         tttGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCga                                                      >  1:50739/1‑49 (MQ=255)
                                                                         tttGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCga                                                      >  1:149780/1‑49 (MQ=255)
                                                                         tttGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCga                                                      <  2:50739/49‑1 (MQ=255)
                                                                          ttgttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGt            <  1:189804/90‑1 (MQ=255)
                                                                              tgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAAt        <  1:173582/90‑1 (MQ=255)
                                                                               gACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATa       <  2:389212/90‑1 (MQ=255)
                                                                                aCGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATAc      >  2:239616/1‑90 (MQ=255)
                                                                                    tGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAc  <  2:49110/90‑1 (MQ=255)
                                                                                    tGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAc  <  1:107123/90‑1 (MQ=255)
                                                                                       |                                                                                      
TACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568665‑568838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ 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
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTAC  >  NZ_CP009273/568653‑568851
                                                                                                   |                                                                                                   
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGG                                                                                                     >  SRR3722109.30615/1‑100 (MQ=60)
    GGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATG                                                                                                 <  SRR3722109.554672/100‑1 (MQ=60)
          GGTACAAAACCGCCAGCCCTGATTGCCACCTTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAG                                                                                           <  SRR3722109.335950/100‑1 (MQ=60)
            TCCGAACCCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAA                                                                                         <  SRR3722109.10728/100‑1 (MQ=60)
                      CCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGG                                                                               <  SRR3722109.170215/100‑1 (MQ=60)
                                          TGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAG                                                           >  SRR3722109.111818/1‑100 (MQ=60)
                                           GCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGG                                                          <  SRR3722109.384479/100‑1 (MQ=60)
                                             ATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGctgtctctta                                                        >  SRR3722109.501455/1‑90 (MQ=60)
                                                          CACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAG                                           <  SRR3722109.25496/100‑1 (MQ=60)
                                                                       ACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAA                              <  SRR3722109.536414/100‑1 (MQ=60)
                                                                           CCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTctgtctcttatacacatctgacgctgccgac                          >  SRR3722109.151631/1‑69 (MQ=60)
                                                                           CCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTctgtctcttatacacatctgacgctgccgac                          >  SRR3722109.51337/1‑69 (MQ=60)
                                                                                  CATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGctgtctcttatacacatctgacgctgccgacgaatgtcaatgtgtagat                   >  SRR3722109.313513/1‑51 (MQ=60)
                                                                                      TTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC               <  SRR3722109.192137/100‑1 (MQ=60)
                                                                                          TGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCA           <  SRR3722109.175724/100‑1 (MQ=60)
                                                                                                TGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACAT     <  SRR3722109.108461/100‑1 (MQ=60)
                                                                                                   GTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTAC  <  SRR3722109.517285/100‑1 (MQ=60)
                                                                                                   |                                                                                                   
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTAC  >  NZ_CP009273/568653‑568851

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