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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,849,816 T→C E117G (GAA→GGA)  ygbA ← nitrous oxide‑stimulated promoter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,849,8160TC100.0% 133.9 / NA 39E117G (GAA→GGA) ygbAnitrous oxide‑stimulated promoter family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/16);  total (23/16)

CACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAG  >  NZ_CP009273/2849732‑2849896
                                                                                    |                                                                                
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  2:981022/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  2:806760/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  1:195208/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  1:366073/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  2:1056257/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  1:426181/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  1:460123/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                             >  1:941587/1‑90 (MQ=255)
cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCAGCCTGTCCGGCATACTATCCATGa                                                                             >  2:207072/1‑90 (MQ=255)
     aTGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGttt                                                                        >  2:350459/1‑90 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGAc                                                                <  2:85055/90‑1 (MQ=255)
             gCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCACACTATCCATGAGGTTTTTTCGGAc                                                                <  2:290568/90‑1 (MQ=255)
                           aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                  <  1:1039517/90‑1 (MQ=255)
                            aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                 <  2:941587/90‑1 (MQ=255)
                                          aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                 >  1:787946/1‑76 (MQ=255)
                                          aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                 <  2:787946/76‑1 (MQ=255)
                                              gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCt                               >  2:752190/1‑90 (MQ=255)
                                               cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                              >  1:1003823/1‑90 (MQ=255)
                                               cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                              >  1:633681/1‑90 (MQ=255)
                                               cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                              >  1:923370/1‑90 (MQ=255)
                                                aaaTTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTg                             >  1:1049280/1‑90 (MQ=255)
                                                   ttGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCa                          <  1:908913/90‑1 (MQ=255)
                                                        ggACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGAACCGGACGCTTGTCATCgat                     <  2:408166/90‑1 (MQ=255)
                                                         gACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                    <  1:406037/90‑1 (MQ=255)
                                                          aCAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatga                   <  2:88951/90‑1 (MQ=255)
                                                            aGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgaga                 <  2:1195597/90‑1 (MQ=255)
                                                             gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGACATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat                >  1:366414/1‑90 (MQ=255)
                                                               tCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagatga              >  2:1063797/1‑90 (MQ=255)
                                                                  gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc           >  1:1106934/1‑90 (MQ=255)
                                                                   ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa          <  2:93413/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc         >  2:530121/1‑90 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc         <  1:806760/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc         <  1:543549/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc         <  2:1001329/90‑1 (MQ=255)
                                                                    cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc         <  1:424724/90‑1 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt      >  1:40571/1‑90 (MQ=255)
                                                                       tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt      >  2:1165221/1‑90 (MQ=255)
                                                                         cGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTa    <  2:426181/90‑1 (MQ=255)
                                                                           gCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAg  >  2:934927/1‑90 (MQ=255)
                                                                                    |                                                                                
CACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAG  >  NZ_CP009273/2849732‑2849896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  >  NZ_CP009273/2849720‑2849914
                                                                                                |                                                                                                  
CGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCAT                                                                                                 >  SRR3722089.413205/1‑100 (MQ=60)
 GACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATG                                                                                                <  SRR3722089.594381/100‑1 (MQ=60)
 GACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATG                                                                                                <  SRR3722089.582236/100‑1 (MQ=60)
  ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722089.955380/1‑100 (MQ=60)
  ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722089.465619/1‑100 (MQ=60)
  ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722089.431324/1‑100 (MQ=60)
  ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722089.197667/1‑100 (MQ=60)
  ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722089.370582/1‑100 (MQ=60)
  ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722089.281821/1‑100 (MQ=60)
   CCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAG                                                                                              <  SRR3722089.115207/100‑1 (MQ=60)
                                       AATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCA                                                          <  SRR3722089.1055076/100‑1 (MQ=60)
                                            GCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCtgt                                                     >  SRR3722089.798756/1‑97 (MQ=60)
                                                 GATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTT                                                >  SRR3722089.1018762/1‑100 (MQ=60)
                                                 GATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTT                                                >  SRR3722089.936840/1‑100 (MQ=60)
                                                 GATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTT                                                >  SRR3722089.641745/1‑100 (MQ=60)
                                                  ATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTG                                               >  SRR3722089.1065010/1‑100 (MQ=60)
                                                               TTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGAT                                  <  SRR3722089.922131/100‑1 (MQ=60)
                                                               TTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGACATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGAT                                  >  SRR3722089.370925/1‑100 (MQ=60)
                                                                    GGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGC                             >  SRR3722089.1123672/1‑100 (MQ=60)
                                                                     gtgtataagagacagCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCA                            <  SRR3722089.372777/85‑1 (MQ=60)
                                                                     GACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCA                            <  SRR3722089.410958/100‑1 (MQ=60)
                                                                         GATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGT                        >  SRR3722089.41081/1‑100 (MQ=60)
                                                                         GATCTGCCTGTCCGGCATACTATCCATGAGGTATTTTCGGACGATATTTTTCCtgtctcttatacacatctgacgctgccgacgaccatcttagtgtaga                        >  SRR3722089.343885/1‑53 (MQ=60)
                                                                                CTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATC                 <  SRR3722089.429851/100‑1 (MQ=60)
                                                                                CTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATC                 <  SRR3722089.817988/100‑1 (MQ=60)
                                                                                CTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATC                 <  SRR3722089.549877/100‑1 (MQ=60)
                                                                                    CCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAT             >  SRR3722089.35554/1‑100 (MQ=60)
                                                                                            CTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCAC     >  SRR3722089.149712/1‑100 (MQ=60)
                                                                                              ATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAG   >  SRR3722089.246359/1‑100 (MQ=60)
                                                                                               TCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  <  SRR3722089.504539/100‑1 (MQ=60)
                                                                                               TCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  >  SRR3722089.1225900/1‑100 (MQ=60)
                                                                                               TCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  >  SRR3722089.143376/1‑100 (MQ=60)
                                                                                                |                                                                                                  
CGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  >  NZ_CP009273/2849720‑2849914

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