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 I210 R1 226 19.6 1130256 95.4% 1078264 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 444,519 G→A R527C (CGT→TGT)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273444,5190GA100.0% 82.6 / NA 23R527C (CGT→TGT) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/11);  total (12/11)

TATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTG  >  NZ_CP009273/444447‑444593
                                                                        |                                                                          
tataGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgcggtcgc                                                           >  2:264650/1‑90 (MQ=255)
        ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                   >  2:9057/1‑90 (MQ=255)
        ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                   >  2:169597/1‑90 (MQ=255)
        ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                   >  1:261488/1‑90 (MQ=255)
        ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                   >  1:308449/1‑90 (MQ=255)
            gCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGt                                               <  1:135683/90‑1 (MQ=255)
              ggggAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTAc                                             >  1:176801/1‑90 (MQ=255)
              ggggAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTAc                                             <  1:199577/90‑1 (MQ=255)
                             cACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAgg                              <  1:9057/90‑1 (MQ=255)
                               cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                            <  1:223753/90‑1 (MQ=255)
                               cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                            <  2:176801/90‑1 (MQ=255)
                                         cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATc                                          >  1:39421/1‑66 (MQ=255)
                                         cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATc                                          <  2:39421/66‑1 (MQ=255)
                                         cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                  >  1:38506/1‑90 (MQ=255)
                                             cACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcg              <  1:169597/90‑1 (MQ=255)
                                             cACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcg              <  2:485942/90‑1 (MQ=255)
                                                cccACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCaa           <  2:407161/90‑1 (MQ=255)
                                                  cACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc         <  1:264650/90‑1 (MQ=255)
                                                     gggTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc         >  1:518220/1‑87 (MQ=255)
                                                     gggTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc         <  2:518220/87‑1 (MQ=255)
                                                     gggTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGt      >  2:143848/1‑90 (MQ=255)
                                                        tCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACt   >  1:422481/1‑90 (MQ=255)
                                                         cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTg  >  1:394149/1‑90 (MQ=255)
                                                                        |                                                                          
TATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTG  >  NZ_CP009273/444447‑444593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ 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
TGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAC  >  NZ_CP009273/444422‑444596
                                                                                                 |                                                                             
caccgagatctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagCTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGT                                                                             <  SRR3722092.245639/44‑1 (MQ=60)
  AACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTC                                                                           <  SRR3722092.321052/100‑1 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                      >  SRR3722092.265797/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                      >  SRR3722092.313888/1‑100 (MQ=60)
                             GAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTAC                                                >  SRR3722092.179752/1‑100 (MQ=60)
                                     GCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAAT                                        <  SRR3722092.137963/100‑1 (MQ=60)
                                       GGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAA                                      <  SRR3722092.202874/100‑1 (MQ=60)
                                                      CACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAC                       <  SRR3722092.9221/100‑1 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGActgtctcttataca                     >  SRR3722092.40106/1‑86 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                     <  SRR3722092.227409/100‑1 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                     >  SRR3722092.39180/1‑100 (MQ=60)
                                                                    GCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGT         >  SRR3722092.528610/1‑100 (MQ=60)
                                                                      CACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAC       <  SRR3722092.172442/100‑1 (MQ=60)
                                                                       ACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACT      >  SRR3722092.430664/1‑100 (MQ=60)
                                                                        CCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTG     >  SRR3722092.401696/1‑100 (MQ=60)
                                                                           CACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAC  <  SRR3722092.269053/100‑1 (MQ=60)
                                                                                                 |                                                                             
TGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAC  >  NZ_CP009273/444422‑444596

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: