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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,011,010 A→G E531G (GAA→GGA)  ygfK → putative selenate reductase subunit YgfK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,011,0100AG100.0% 88.3 / NA 28E531G (GAA→GGA) ygfKputative selenate reductase subunit YgfK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/13);  total (15/13)

CACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGT  >  NZ_CP009273/3010931‑3011095
                                                                               |                                                                                     
cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                             >  2:271664/1‑90 (MQ=255)
cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                             >  2:198793/1‑90 (MQ=255)
  ccAGTGCAAATACAGCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                           <  1:83119/90‑1 (MQ=255)
              cAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                               <  1:38795/90‑1 (MQ=255)
              cAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                               <  1:461394/90‑1 (MQ=255)
               aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                          <  1:459630/78‑1 (MQ=255)
               aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                          >  2:459630/1‑78 (MQ=255)
                  tGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc                                                           <  2:285661/90‑1 (MQ=255)
                   gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt                                                          >  2:297087/1‑90 (MQ=255)
                        cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                     >  2:521222/1‑90 (MQ=255)
                                 tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                            <  1:397633/90‑1 (MQ=255)
                                  aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                           >  2:627298/1‑90 (MQ=255)
                                  aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                           >  1:285974/1‑90 (MQ=255)
                                       aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                      >  1:616952/1‑90 (MQ=255)
                                       aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                      >  1:438885/1‑90 (MQ=255)
                                        gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGtt                                                                               >  1:616542/1‑48 (MQ=255)
                                        gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGtt                                                                               <  2:616542/48‑1 (MQ=255)
                                        gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                     >  2:182679/1‑90 (MQ=255)
                                        gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                     >  1:564983/1‑90 (MQ=255)
                                                aaTATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGc                             <  2:564983/90‑1 (MQ=255)
                                                    tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                         >  2:598208/1‑90 (MQ=255)
                                                       gcgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCgg                      >  1:488910/1‑90 (MQ=255)
                                                        cgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGt                     <  1:297087/90‑1 (MQ=255)
                                                             tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg                <  2:285974/90‑1 (MQ=255)
                                                              gAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGt               >  1:493928/1‑90 (MQ=255)
                                                                  aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt           <  1:478501/90‑1 (MQ=255)
                                                                       cgcgCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGc      <  1:474185/90‑1 (MQ=255)
                                                                           cTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGt  <  2:28330/90‑1 (MQ=255)
                                                                               |                                                                                     
CACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGT  >  NZ_CP009273/3010931‑3011095

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ 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
GTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTCCGGCG  >  NZ_CP009273/3010917‑3011101
                                                                                             |                                                                                           
gtaattttgGCGATCACCAGTGCCAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGG                                                                                       <  SRR3722088.325274/91‑1 (MQ=60)
GTCATTTTTGCGATCCCCAGTGCCAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGG                                                                                       <  SRR3722088.96204/100‑1 (MQ=60)
GTCATTTTGGCGATCACCAGTCCAAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGG                                                                                       <  SRR3722088.285499/100‑1 (MQ=60)
                CCAGTGCAAATACAGCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAG                                                                       <  SRR3722088.83952/100‑1 (MQ=60)
                 gggcccaaaaacacctgtcCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGC                                                                      <  SRR3722088.465006/81‑1 (MQ=60)
                            CAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCAC                                                           <  SRR3722088.39202/100‑1 (MQ=60)
                            CAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCAC                                                           <  SRR3722088.466792/100‑1 (MQ=60)
                                      CGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCG                                                 >  SRR3722088.288730/1‑100 (MQ=60)
                                           GGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCT                                            >  SRR3722088.443933/1‑100 (MQ=60)
                                           GGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCT                                            >  SRR3722088.624559/1‑100 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATActgtctcttatacacatctccgagcccacgag                                           >  SRR3722088.624142/1‑68 (MQ=60)
                                            GATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTG                                           >  SRR3722088.571911/1‑100 (MQ=60)
                                               TACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTT                                        <  SRR3722088.401935/100‑1 (MQ=60)
                                                           CTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATctgtctcttatacacatctccgagcccacgagaccgtactagatct                            >  SRR3722088.14748/1‑54 (MQ=60)
                                                           CTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGG                            >  SRR3722088.494735/1‑100 (MQ=60)
                                                                  TCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGT                     >  SRR3722088.499830/1‑100 (MQ=60)
                                                                      CGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT                 <  SRR3722088.299973/100‑1 (MQ=60)
                                                                                AAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTC       <  SRR3722088.484191/100‑1 (MQ=60)
                                                                                     CGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTCCGGCG  <  SRR3722088.479810/100‑1 (MQ=60)
                                                                                             |                                                                                           
GTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGTCCGGCG  >  NZ_CP009273/3010917‑3011101

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