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 I227 R1 223 21.3 1167702 97.0% 1132670 86.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% 68.0 / NA 22K405E (AAA→GAA) preANAD‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/7);  total (15/7)

GAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGG  >  NZ_CP009273/2229869‑2230041
                                                                                       |                                                                                     
gagaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaa                                                                                     >  2:271031/1‑90 (MQ=255)
  gaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaaga                                                                                   >  2:70339/1‑90 (MQ=255)
  gaAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGaggaaga                                                                                   >  2:60791/1‑90 (MQ=255)
      tgtgtgGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAacac                                                                               <  1:562882/90‑1 (MQ=255)
       gtgtgGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACAcc                                                                              >  1:163002/1‑90 (MQ=255)
               tGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACg                                                                      >  1:226370/1‑90 (MQ=255)
               tGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACg                                                                      >  2:288631/1‑90 (MQ=255)
                          tgGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCa                                                           >  2:402070/1‑90 (MQ=255)
                          tgGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCa                                                           >  2:341320/1‑90 (MQ=255)
                                   cTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCAt                                                  >  2:28463/1‑90 (MQ=255)
                                    tGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATc                                                 >  2:576638/1‑90 (MQ=255)
                                          gTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCAt                                           >  2:524940/1‑90 (MQ=255)
                                          gTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCAt                                           >  2:342579/1‑90 (MQ=255)
                                                      gAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTc                               >  1:207320/1‑90 (MQ=255)
                                                        gCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGc                             >  2:71865/1‑90 (MQ=255)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAac                                                                                 <  2:434770/35‑1 (MQ=38)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAac                                                                                 >  1:434770/1‑35 (MQ=38)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCgg                          >  2:359503/1‑90 (MQ=255)
                                                           cGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCgg                          >  2:94809/1‑90 (MQ=255)
                                                                      agtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGa               <  1:60791/90‑1 (MQ=255)
                                                                      agtTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGa               <  1:524940/90‑1 (MQ=255)
                                                                                aGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATc     <  2:406789/90‑1 (MQ=255)
                                                                                aGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATc     <  1:285532/90‑1 (MQ=255)
                                                                                 ggCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATcc    >  1:9185/1‑90 (MQ=255)
                                                                                   cGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCgg  <  2:226370/90‑1 (MQ=255)
                                                                                       |                                                                                     
GAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGG  >  NZ_CP009273/2229869‑2230041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ 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
ATTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAA  >  NZ_CP009273/2229858‑2230048
                                                                                                  |                                                                                            
accgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacagGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGA                                                                                             <  SRR3722111.567035/45‑1 (MQ=60)
       TACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACAC                                                                                      <  SRR3722111.330099/100‑1 (MQ=60)
        ACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACC                                                                                     >  SRR3722111.165057/1‑100 (MQ=60)
                ATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACG                                                                             >  SRR3722111.229102/1‑100 (MQ=60)
                 TGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGT                                                                            <  SRR3722111.571732/100‑1 (MQ=60)
                                                       GGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTC                                      >  SRR3722111.209869/1‑100 (MQ=60)
                                                            GTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACctgtctcttatacacatctgacgctgccgacgaatgtcaatgtgt                                 >  SRR3722111.441271/1‑55 (MQ=60)
                                                                                 AGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATC            <  SRR3722111.533158/100‑1 (MQ=60)
                                                                                 AGTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATC            <  SRR3722111.61524/100‑1 (MQ=60)
                                                                                  GTTTAAGAAAGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCC           >  SRR3722111.9301/1‑100 (MQ=60)
                                                                                           AGGCGAGGAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAA  <  SRR3722111.289152/100‑1 (MQ=60)
                                                                                                  |                                                                                            
ATTGTAATACCGAGAAATGTGTGGGTTGTCTGCTTTGTGGTCACGTCTGCCCGGTGGGTTGTATTGAGCTCGGGGAAGTGAAGTTTAAGAAAGGCGAGAAAGAACACCCGGTAACGTTGTAAAACCAGAGACGCATCCGGCATTTGGTTCCTGTCGCCGGATGCGGCTTGAACGCCTTATCCGGCCTACAA  >  NZ_CP009273/2229858‑2230048

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