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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

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% 92.7 / NA 27G71G (GGA→GGGybcODUF1364 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/21);  total (6/21)

AACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568670‑568838
                                                                                  |                                                                                      
aaCCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCaa                                                                                 >  1:325740/1‑90 (MQ=255)
          ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:565688/90‑1 (MQ=255)
          ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:341089/90‑1 (MQ=255)
          ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:183194/90‑1 (MQ=255)
          ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  1:483580/90‑1 (MQ=255)
            tGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgcgc                                                                     <  2:174864/90‑1 (MQ=255)
             gATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCt                                                                    >  2:766803/1‑90 (MQ=255)
                 gCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGaa                                                                <  2:181071/90‑1 (MQ=255)
                     ccaTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTa                                                            <  1:733923/90‑1 (MQ=255)
                       aTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATg                                                          <  2:454407/90‑1 (MQ=255)
                       aTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATg                                                          <  2:325740/90‑1 (MQ=255)
                          gCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  1:523364/90‑1 (MQ=255)
                          gCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  2:447042/90‑1 (MQ=255)
                               ttCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       >  2:330326/1‑85 (MQ=255)
                               ttCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  1:330326/85‑1 (MQ=255)
                                   gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  2:765850/1‑90 (MQ=255)
                                   gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  1:189336/1‑90 (MQ=255)
                                   gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGGAGGCAAAAGAAGGCGCGCTGGAGGGTATGGCGAGAACACAg                                              >  1:399735/1‑90 (MQ=255)
                                                aaaTCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGa                                 <  1:40483/90‑1 (MQ=255)
                                                aaaTCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGa                                 <  2:689334/90‑1 (MQ=255)
                                                     gcccgccgCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAagag                            <  1:649265/88‑1 (MQ=255)
                                                      aCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAgg                           <  1:401857/90‑1 (MQ=255)
                                                               cgcaTTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGttatta                  <  2:28391/90‑1 (MQ=255)
                                                                  aTTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAgg               <  2:66429/90‑1 (MQ=255)
                                                                        ttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGGTGAAAGAGGGGGTTATTAAGGCGTGaa         <  1:588134/90‑1 (MQ=255)
                                                                             gCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCt    <  1:233591/90‑1 (MQ=255)
                                                                               tGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAc  <  2:109235/90‑1 (MQ=255)
                                                                                  |                                                                                      
AACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568670‑568838

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATT  >  NZ_CP009273/568653‑568849
                                                                                                   |                                                                                                 
GTGCGGCCCCGGTACGAAACCGCCAGCCCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGG                                                                                                   <  SRR3722073.61468/100‑1 (MQ=60)
  GCGGCCCCGGTACGAAACCGCCAGCCCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTA                                                                                                 <  SRR3722073.309128/100‑1 (MQ=60)
      cccaggtaaGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCA                                                                                             <  SRR3722073.724501/91‑1 (MQ=60)
      CCCCGGTACAAACCCGCCAGCCCTGATTGCCCCCTTTGCATGTTCTGCCTGCCCCGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCA                                                                                             <  SRR3722073.569527/100‑1 (MQ=60)
       ACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAA                                                                                            >  SRR3722073.332658/1‑100 (MQ=60)
          ggtagcaaccctacagccgtgattgccccctttgccttTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAG                                                                                         <  SRR3722073.9653/62‑1 (MQ=60)
                           CCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGT                                                                        <  SRR3722073.493854/100‑1 (MQ=60)
                                      CCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAAC                                                             <  SRR3722073.748778/100‑1 (MQ=60)
                                          TGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAG                                                         >  SRR3722073.193078/1‑100 (MQ=60)
                                          TGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGGAGGCAAAAGAAGGCGCGCTGGAGGGTATGGCGAGAACACAG                                                         >  SRR3722073.408247/1‑100 (MQ=60)
                                           GCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGG                                                        <  SRR3722073.337341/100‑1 (MQ=60)
                                           GCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGG                                                        <  SRR3722073.534345/100‑1 (MQ=60)
                                                                 AAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGT                                  <  SRR3722073.41466/100‑1 (MQ=60)
                                                                      GCCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTA                             <  SRR3722073.662623/100‑1 (MQ=60)
                                                                       ACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAA                            <  SRR3722073.410421/100‑1 (MQ=60)
                                                                                        GTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGGTATTAAGGCGTGAATACCTACAG           >  SRR3722073.505380/1‑100 (MQ=60)
                                                                                         TTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGGTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGC          <  SRR3722073.600421/100‑1 (MQ=60)
                                                                                              GCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCAC     <  SRR3722073.238480/100‑1 (MQ=60)
                                                                                                 GGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATT  >  SRR3722073.325919/1‑100 (MQ=60)
                                                                                                 GGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATT  >  SRR3722073.87036/1‑100 (MQ=60)
                                                                                                   |                                                                                                 
GTGCGGCACCGGTACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATT  >  NZ_CP009273/568653‑568849

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