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 I230 R1 226 18.8 1048726 96.7% 1014118 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,002,903 T→C S261P (TCA→CCA)  ygeY → YgeY family selenium metabolism‑linked hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,002,9030TC100.0% 87.2 / NA 27S261P (TCA→CCA) ygeYYgeY family selenium metabolism‑linked hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

TTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAG  >  NZ_CP009273/3002819‑3002984
                                                                                    |                                                                                 
tccgtCGGCAAAGGCACCCTCCCCGTTTCTGAAATCTTGTTCACATCCCAAAGCCGTTGCGCTGTAGCAGACAGCTTCGCAGTTCCAATc                                                                              <  2:293006/86‑1 (MQ=255)
    tCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGAcc                                                                          >  1:485617/1‑90 (MQ=255)
       gCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcgtc                                                                       <  1:50660/90‑1 (MQ=255)
          aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                    <  1:425538/90‑1 (MQ=255)
                    cACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTgggg                                                          >  1:1437/1‑90 (MQ=255)
                     aCCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGc                                                         <  1:399108/90‑1 (MQ=255)
                                     tcttcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcg                                         >  1:307184/1‑90 (MQ=255)
                                      cttcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgc                                        >  1:205784/1‑90 (MQ=255)
                                         cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCACAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt                                     <  2:441606/90‑1 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                       >  1:184873/1‑68 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                       <  2:184873/68‑1 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAgg                                           >  2:344713/1‑80 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAgg                                           <  1:344713/80‑1 (MQ=255)
                                              ccccAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACg                                <  2:205784/90‑1 (MQ=255)
                                                ccAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGaa                              >  2:224295/1‑90 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                    >  2:208228/1‑45 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                    <  1:208228/45‑1 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                       >  1:395948/1‑90 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                       >  1:9879/1‑90 (MQ=255)
                                                            gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                  <  1:135974/90‑1 (MQ=255)
                                                            gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                  <  1:244538/90‑1 (MQ=255)
                                                               gTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcct               >  1:430921/1‑90 (MQ=255)
                                                                 aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc             <  1:161827/90‑1 (MQ=255)
                                                                    aGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGt          <  1:224295/90‑1 (MQ=255)
                                                                         gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc             <  2:225973/82‑1 (MQ=255)
                                                                         gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc             >  1:225973/1‑82 (MQ=255)
                                                                            gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAg  >  2:268852/1‑90 (MQ=255)
                                                                                    |                                                                                 
TTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAG  >  NZ_CP009273/3002819‑3002984

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
TTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTA  >  NZ_CP009273/3002809‑3003002
                                                                                              |                                                                                                   
TTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATC                                                                                                <  SRR3722114.12561/100‑1 (MQ=60)
   TGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGAC                                                                                             <  SRR3722114.151791/100‑1 (MQ=60)
    GACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACC                                                                                            >  SRR3722114.491695/1‑100 (MQ=60)
     ACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCctgtc                                                                                           >  SRR3722114.260717/1‑95 (MQ=60)
                 GCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTG                                                                               <  SRR3722114.51169/100‑1 (MQ=60)
                    ggtctcgtgggctcggagatgtgtataagagacagTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGG                                                                            <  SRR3722114.210242/65‑1 (MQ=60)
                    AAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGG                                                                            >  SRR3722114.1453/1‑100 (MQ=60)
                    AAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGG                                                                            <  SRR3722114.430773/100‑1 (MQ=60)
                               ACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGG                                                                 <  SRR3722114.403958/100‑1 (MQ=60)
                                     TCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCG                                                           >  SRR3722114.310462/1‑100 (MQ=60)
                                      CTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGC                                                          >  SRR3722114.207784/1‑100 (MQ=60)
                                             CTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAActgtctcttata                                                   >  SRR3722114.186694/1‑88 (MQ=60)
                                             CTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGAC                                                   <  SRR3722114.348640/100‑1 (MQ=60)
                                                       TCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCG                                         >  SRR3722114.9977/1‑100 (MQ=60)
                                                       TCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCG                                         >  SRR3722114.400748/1‑100 (MQ=60)
                                                               CCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCT                                 >  SRR3722114.436227/1‑100 (MQ=60)
                                                                      GCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTAC                          <  SRR3722114.246907/100‑1 (MQ=60)
                                                                      GCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTAC                          <  SRR3722114.137377/100‑1 (MQ=60)
                                                                         GTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGA                       >  SRR3722114.228117/1‑100 (MQ=60)
                                                                           AGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAA                     <  SRR3722114.163438/100‑1 (MQ=60)
                                                                              AGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCT                  <  SRR3722114.226424/100‑1 (MQ=60)
                                                                                 CAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAC               >  SRR3722114.222544/1‑100 (MQ=60)
                                                                                              CCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTA  >  SRR3722114.438711/1‑100 (MQ=60)
                                                                                              |                                                                                                   
TTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTA  >  NZ_CP009273/3002809‑3003002

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