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 F3 I134 R1 6 21.3 1193788 96.3% 1149617 85.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,754,231 G→T R127S (CGC→AGC)  selA ← L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,754,2310GT100.0% 122.2 / NA 32R127S (CGC→AGC) selAL‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/22);  total (10/22)

CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGC  >  NZ_CP009273/3754150‑3754309
                                                                                 |                                                                              
cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        <  2:98772/90‑1 (MQ=255)
cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        <  2:159394/90‑1 (MQ=255)
cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        <  2:388378/90‑1 (MQ=255)
      tATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                  <  2:304859/90‑1 (MQ=255)
      tATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                  >  2:173955/1‑90 (MQ=255)
         tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                               <  2:569234/90‑1 (MQ=255)
           aCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa                                                             <  1:406051/90‑1 (MQ=255)
              aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                          <  1:354133/90‑1 (MQ=255)
              aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                          <  2:465380/90‑1 (MQ=255)
                  aGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                      <  1:75236/90‑1 (MQ=255)
                      aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        <  2:546606/68‑1 (MQ=255)
                      aTCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        >  1:546606/1‑68 (MQ=255)
                         ttCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                               <  2:23595/90‑1 (MQ=255)
                             gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa                                                             >  1:500007/1‑72 (MQ=255)
                             gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa                                                             <  2:500007/72‑1 (MQ=255)
                               ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                         <  1:252812/90‑1 (MQ=255)
                               ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                         <  1:109129/90‑1 (MQ=255)
                                cccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAgg                                                    <  2:415721/78‑1 (MQ=255)
                                cccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAgg                                                    >  1:415721/1‑78 (MQ=255)
                                 ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                                       >  1:123149/1‑90 (MQ=255)
                                    tAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                    <  2:251323/90‑1 (MQ=255)
                                        aCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGc                                >  2:498563/1‑90 (MQ=255)
                                           gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                             >  2:99155/1‑90 (MQ=255)
                                           gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCAGCGTCGTCCAGATCATACTCGAGGGTCACTTGCGAACGCATAGCCtc                             >  1:381266/1‑89 (MQ=255)
                                            cacaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        <  2:193592/46‑1 (MQ=255)
                                            cacaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                        >  1:193592/1‑46 (MQ=255)
                                                gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                        <  2:495152/90‑1 (MQ=255)
                                                   gcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACgg                     >  1:369382/1‑90 (MQ=255)
                                                     tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                   <  2:104141/90‑1 (MQ=255)
                                                             cgcCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACcgc           <  1:323460/90‑1 (MQ=255)
                                                                     cgcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCg   <  1:55460/90‑1 (MQ=255)
                                                                      gcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGc  <  1:503991/90‑1 (MQ=255)
                                                                                 |                                                                              
CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGC  >  NZ_CP009273/3754150‑3754309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ 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
ATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGT  >  NZ_CP009273/3754135‑3754329
                                                                                                |                                                                                                  
ATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCC                                                                                                 <  SRR3722002.236691/100‑1 (MQ=60)
                          ACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGG                                                                       <  SRR3722002.411185/100‑1 (MQ=60)
                           CGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCctgtctcttata                                                                      >  SRR3722002.553879/1‑88 (MQ=60)
                             ATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCA                                                                    <  SRR3722002.358329/100‑1 (MQ=60)
                                 AGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGG                                                                <  SRR3722002.76037/100‑1 (MQ=60)
                                  GGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGctgtctct                                                               >  SRR3722002.506645/1‑92 (MQ=60)
                                     ATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGct                                                            >  SRR3722002.421023/1‑98 (MQ=60)
                                      TCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAC                                                           >  SRR3722002.124563/1‑100 (MQ=60)
                                              CCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCC                                                   <  SRR3722002.110347/100‑1 (MQ=60)
                                              CCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCC                                                   <  SRR3722002.255462/100‑1 (MQ=60)
                                                CCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCAGCGTCGTCCAGATCATACTCGAGGGTCACTTGCGAACGCATAGCCTC                                                 >  SRR3722002.385967/1‑100 (MQ=60)
                                                 CGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCctgtctcttatacacatctccgagcccacgagac                                                >  SRR3722002.195727/1‑66 (MQ=60)
                                                        CGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGG                                         >  SRR3722002.373870/1‑100 (MQ=60)
                                                                            CGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCC                     <  SRR3722002.327220/100‑1 (MQ=60)
                                                                                    CGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGC             <  SRR3722002.56031/100‑1 (MQ=60)
                                                                                     GCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCT            <  SRR3722002.510696/100‑1 (MQ=60)
                                                                                               CTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGT  >  SRR3722002.541617/1‑100 (MQ=60)
                                                                                               CTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGT  >  SRR3722002.104041/1‑100 (MQ=60)
                                                                                                |                                                                                                  
ATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGT  >  NZ_CP009273/3754135‑3754329

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: