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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 900,176 T→G V43G (GTT→GGT)  ybjQ → heavy metal‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273900,1760TG100.0% 81.3 / NA 23V43G (GTT→GGT) ybjQheavy metal‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (7/16);  total (7/16)

TCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCG  >  NZ_CP009273/900089‑900261
                                                                                       |                                                                                     
tCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTgg                                                                                    >  1:251023/1‑91 (MQ=255)
  gTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCg                                                                                  <  2:1116078/91‑1 (MQ=255)
    tGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGa                                                                                <  2:798134/91‑1 (MQ=255)
    tGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGa                                                                                <  2:251023/91‑1 (MQ=255)
     gAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGAc                                                                               >  2:137168/1‑91 (MQ=255)
            gCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCgg                                                                        >  1:1040439/1‑91 (MQ=255)
            gCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCgg                                                                        <  1:427671/91‑1 (MQ=255)
              ggTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTg                                                                      <  1:968278/91‑1 (MQ=255)
                   ggtgACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTAt                                                                 <  2:901127/91‑1 (MQ=255)
                   ggtgACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTAt                                                                 <  2:501390/91‑1 (MQ=255)
                   ggtgACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTAt                                                                 <  2:546503/91‑1 (MQ=255)
                   ggtgACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTAt                                                                 <  2:654005/91‑1 (MQ=255)
                    gtgAGCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATg                                                                >  1:905567/1‑91 (MQ=255)
                           gCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGa                                                         >  1:893078/1‑91 (MQ=255)
                             gAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAAc                                                       <  1:15423/91‑1 (MQ=255)
                                    ttttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTaa                                                <  1:7437/91‑1 (MQ=255)
                                    ttttAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTaa                                                <  2:188742/91‑1 (MQ=255)
                                          gTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTCTGAAAAAGAACTGCGTAAAGCACg                                          <  2:1078342/91‑1 (MQ=255)
                                          gTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACg                                          <  1:222064/91‑1 (MQ=255)
                                                    tttCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGcc                                >  2:158397/1‑91 (MQ=255)
                                                         gTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGa                           <  2:921160/91‑1 (MQ=255)
                                                                             cgcgATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGgcgc       <  2:1040221/91‑1 (MQ=255)
                                                                                  tatCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGgcg  >  1:1092300/1‑91 (MQ=255)
                                                                                       |                                                                                     
TCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCG  >  NZ_CP009273/900089‑900261

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ 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
GGCCAGACCATCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCC  >  NZ_CP009273/900079‑900276
                                                                                                 |                                                                                                    
GGCCAGACCATCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGG                                                                                                   >  SRR3722125.256079/1‑101 (MQ=60)
   CAGACCATCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGG                                                                                                <  SRR3722125.1050917/101‑1 (MQ=60)
     GACCATCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGAC                                                                                              <  SRR3722125.343067/101‑1 (MQ=60)
            GTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGG                                                                                       >  SRR3722125.1063767/1‑101 (MQ=60)
                    TTGCGGTGTGGTGAGCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATG                                                                               >  SRR3722125.925636/1‑101 (MQ=60)
                      GCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAA                                                                             <  SRR3722125.436686/101‑1 (MQ=60)
                        GGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAA                                                                           <  SRR3722125.989875/101‑1 (MQ=60)
                           GTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGA                                                                        >  SRR3722125.912844/1‑101 (MQ=60)
                                       GAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGC                                                            <  SRR3722125.15708/101‑1 (MQ=60)
                                              TTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAG                                                     <  SRR3722125.7578/101‑1 (MQ=60)
                                                    GTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCC                                               <  SRR3722125.226458/101‑1 (MQ=60)
                                                                                  GTATCCGCGATATCGGTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCG                 >  SRR3722125.1116878/1‑101 (MQ=60)
                                                                                                 GTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCtaattc  >  SRR3722125.291513/1‑95 (MQ=60)
                                                                                                 |                                                                                                    
GGCCAGACCATCGTTGAATATTGCGGTGTGGTGACCGGCGAAGCGATTTTAGGTGCCAATATTTTCCGTGATTTCTTTGCCGGTATCCGCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCC  >  NZ_CP009273/900079‑900276

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