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,229,956 A→G K405E (AAA→GAA)  preA → NAD‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,229,9560AG100.0% 65.8 / NA 19K405E (AAA→GAA) preANAD‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/7);  total (12/7)

GAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGA  >  NZ_CP009273/2229869‑2230028
                                                                                       |                                                                        
gagaTATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaa                                                                        >  2:363554/1‑90 (MQ=255)
gagaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaa                                                                        >  1:239671/1‑90 (MQ=255)
gagaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaa                                                                        >  1:25183/1‑90 (MQ=255)
gagaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaa                                                                        >  1:847920/1‑90 (MQ=255)
gagaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGGGAAGTTTAAGAAAGGCGaggaa                                                                        >  2:618670/1‑90 (MQ=255)
      tgtgtgGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAacac                                                                  <  2:569837/90‑1 (MQ=255)
           gggTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGt                                                             >  1:417952/1‑90 (MQ=255)
               tGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACg                                                         <  1:372596/90‑1 (MQ=255)
                 tCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGtt                                                       >  1:639323/1‑90 (MQ=255)
                      tttGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACg                                                         <  2:642412/83‑1 (MQ=255)
                      tttGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACg                                                         >  1:642412/1‑83 (MQ=255)
                              cACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAgaga                                          >  2:851817/1‑90 (MQ=255)
                                      cccGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCg                                  <  1:76162/90‑1 (MQ=255)
                                       ccGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCgg                                 <  1:557756/90‑1 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAacac                                                                  <  1:234787/42‑1 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAacac                                                                  >  2:234787/1‑42 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTc                  >  2:403934/1‑90 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTc                  >  1:53111/1‑90 (MQ=255)
                                                        gCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGc                >  1:765810/1‑90 (MQ=255)
                                                        gCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGc                >  1:129833/1‑90 (MQ=255)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCgg             >  2:337062/1‑90 (MQ=255)
                                                             gggAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGAt           >  2:667236/1‑90 (MQ=255)
                                                                  gtgaagtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGc      >  2:513954/1‑90 (MQ=255)
                                                                      agtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGa  <  1:553273/90‑1 (MQ=255)
                                                                                       |                                                                        
GAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGA  >  NZ_CP009273/2229869‑2230028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
TTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCAT  >  NZ_CP009273/2229859‑2230054
                                                                                                 |                                                                                                  
TTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAA                                                                                                  >  SRR3722072.243688/1‑100 (MQ=60)
TTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAA                                                                                                  >  SRR3722072.25626/1‑100 (MQ=60)
TTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAA                                                                                                  >  SRR3722072.862581/1‑100 (MQ=60)
      TACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACAC                                                                                            <  SRR3722072.536836/100‑1 (MQ=60)
           AGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGT                                                                                       >  SRR3722072.425523/1‑100 (MQ=60)
                ggacgtctcgtgggctcggagatgtgtataagagacagGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGT                                                                                  <  SRR3722072.238705/62‑1 (MQ=60)
                 GTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTT                                                                                 >  SRR3722072.650733/1‑100 (MQ=60)
                      GGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAA                                                                            >  SRR3722072.653866/1‑100 (MQ=60)
                         TGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACC                                                                         <  SRR3722072.379224/100‑1 (MQ=60)
                                                CCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTT                                                  <  SRR3722072.77577/100‑1 (MQ=60)
                                                 CCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTC                                                 <  SRR3722072.567779/100‑1 (MQ=60)
                                                      GGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTC                                            >  SRR3722072.54081/1‑100 (MQ=60)
                                                        GTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGC                                          >  SRR3722072.132175/1‑100 (MQ=60)
                                                        GTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGC                                          >  SRR3722072.779300/1‑100 (MQ=60)
                                                                                AGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATC                  <  SRR3722072.563212/100‑1 (MQ=60)
                                                                                                GGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCAT  >  SRR3722072.322067/1‑100 (MQ=60)
                                                                                                 |                                                                                                  
TTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAAAACCAT  >  NZ_CP009273/2229859‑2230054

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